莱斯衰落概率密度

clear
clc
N=100000;%表示获取r的个数
x=randn（1N）;
y=randn（1N）;
r=sqrt（0.5*（x.^2+y.^2））;
A=1;
sigma=0.5; %方差取0.5，代表多径信号分量的功率
step=0.1;
range=0:step:10;
h=hist（rrange）;%h为r落入区间的次数
%bar（h）;
pr_approx=h/（step*sum（h））;%仿真的
pr=@（sigmarangeA）（range./sigma）.*exp（-（A.^2