matlab实现正态分布到瑞利分布随机变量.doc

clear all
a=1; %方差
mean1=0; %均值
dt=0.2; %直方图的横坐标间隔
%由正态随机变量得出瑞利随机变量
x=mean1+sqrt(a)*randn(1,10000);
y=mean1+sqrt(a)*randn(1,10000);
X=x.^2;
Y=y.^2;
Z=sqrt(X+Y);
%瑞利随机变量直方图
t=0:dt:10;
n=hist(Z,t);
%直方图到概率密度曲线
N=n./(10000*dt);
hold on
plot(t,N,'r.-')
%理想瑞利分布概率密度曲线
k=0:0.1:10;
f=k/a.*exp(-k.^2/(2*a^2));
plot(k,f,'b:')
legend('正态随机变量得出的瑞利分布','瑞利分布概率密度曲线')