0

から一様分布行列の確率密度関数をプロットする必要があります。

U = rand (1,1000)

しかし、ksdensity 関数は使用できません。私はこれを試しました:

term = 1000;
U = rand (1,term);
x=0:0.001:1;
for j = 2:term;
    u_height(j) = u_height(j-1)+((abs(x(j)-U(j))<0.01/2)/0.01)/term;
    n_height(j) = n_height(j-1)+((abs(x(j)-N(j))<0.01/2)/0.01)/term;
end

しかし、正しくプロットされていません

4

1 に答える 1

2

関数は正常ksdensityに動作するはずです。PDF の範囲が有限であることを指定し、ボックス カーネルを使用する必要があります。

u = rand(10000,1);

ksdensity(u, 'Support', [0 1], 'kernel', 'box');

さらに、次を使用してpdfを概算できますhistc

u = rand(10000,1);
bins = 0:.05:1;
counts = hist(u, bins);
p = counts ./ trapz(bins, counts);

plot(bins, p);
于 2013-02-12T14:28:25.143 に答える