%Rayleigh?distribution

%采樣點數N=32768；
%標準差sigma_f=20Hz；
%lambda=2；
%采樣頻率fs=1000；

clc
clear?all
close?all

Num=32768;
sigma_f=15;
lambda=1;

N=20;
fs=1000;
f=fs/Num*[-Num/2:Num/2-1];

n1=randn（1Num）;
n2=randn（1Num）;

L=2*N+1;%濾波器長度
h=zeros（1L）;%濾波器時域系數
n=0:N;
C=2*sigma_f*1/fs*sqrt（pi）*exp（-4*sigma_f^2*pi^2*（1/fs）^2*n.^2）;
for?i=1:2*N+1
????if?i<=N+1
????????h（i）=1/2*C（N+2-i）;
????else
????????h（i）=1/2*C（i-N）;
????end
end

Hw=fftshift（fft（hNum））;
abs_Hw=abs（Hw）;
x1=ifft（fft（n1Num）.*Hw）;
x2=ifft（fft（n2Num）.*Hw）;
%*****************************%未標準正態化處理
x1_mu=mean（x1）;????%平均值
x2_mu=mean（x2）;
x1_sigma=std（x1）;??%標準差
x2_sigma=std（x2）;
x1=（x1-x1_mu）/x1_sigma;
x2=（x2-x2_mu）/x2_s