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');

此外,您可以使用近似 pdfhistc

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 回答