雷达杂波瑞利分布模型MATLAB仿真

%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