matlab - trapz（x、f）によって過小評価された周波数密度のプロット

Question

以下のデータセットから、ビン幅=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は、ゼロ次ホールドを使用するのに対し、1次ホールドを介して信号を再構築する積分計算です。そのため、比率は1ではありません。詳細については、近似の順序を確認することもできます。sum(10*f)