matlab - 绘制被 trapz(x,f) 低估的频率密度

Question

我想从下面的数据集中绘制频率密度图，bin width=10 个单位：

[f x]=hist(p,1:10:200); bar(x,f/trapz(x,f));

然而sum(10*f)/trapz(x,f)=1.0461。面积不应该是trapz=1吗？

`p=[23.7648
40.5013
34.7466
24.8569
68.2174
40.8811
78.1818
18.6804
22.3158
80.4086
15.0927
44.1342
163.6714
6.4450
22.9075
17.5615
8.8536
26.3300
32.6789
42.4225
22.3832
3.0519
16.6581
32.9496
10.3606]`

Answer 1

trapz是一种积分计算，它通过一阶保持重构信号，而sum(10*f)使用零阶保持。这就是为什么您的比率不是 1。有关更多信息，您还可以查看Orders of approximation。