陷阱频率密度被trapz（x，f）低估

Question

我想从下面的数据集中绘制频率密度图，其中bin宽度= 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是一个积分计算，通过first-order hold重建信号，而sum(10*f)使用zero-order hold。这就是您的比率不是1的原因。有关详细信息，您还可以查看Orders of approximation。