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