瑞利分布理论和近似概率密度函数

%瑞利分布概率密度函数
clear
clc
N=100000;
x=randn(1,N);
y=randn(1,N);
r=sqrt(0.5*(x.^2+y.^2));
sigma=0.5; %方差取0.5
step=0.1;
range=0:step:3;
h=hist(r,range);%h为落入区间的次数
%bar(h);
pr_approx=h/(step*sum(h));
pr=(range/0.5).*exp(-range.^2/(2*sigma)); %(sigma^2=0.5)
plot(range,pr_approx,'ko',range,pr,'k');
grid
xlabel('r')
ylabel('p(r)')
legend('近似概率密度函数','理论概率密度函数','Northeast')