雷达杂波瑞利模型仿真

clear all；close all;
azi_num=2000;   %取2000个点
fr=1000;        %雷达重复频率
lamda0=0.05;   %杂波波长
sigmav=1.0;     %杂波方差
sigmaf=2*sigmav/lamda0;  
rand（‘state‘sum（100*clock））; %产生服从U（01）分布的随机序列
d1=rand（1azi_num）;            
rand（‘state‘7*sum（100*clock）+3）;
d2=rand（1azi_num）;
xi=2*sqrt（-2*log（d1））.*cos（2*pi*d2）;  %正交且独立的高斯序列N（0，1）
xq=2*sqrt（-2*log（d1））.*sin（2*pi*d2）;
%形成滤波器频率响应
coe_num=12;           %求滤波器系数，用傅里叶级数展开法
for n=0:coe_num
    coeff（n+1）=2*sigmaf*sqrt（pi）*exp（-4*sigmaf^2*pi^2*n^2/fr^2）/fr;  
end
for n=1:2*coe_num+1
    if n<=coe_num+1
        b（n）=1/2*coeff（coe_num+2-n）;
    else
        b（n）=1/2*coeff（n-coe_num）;
    end
end
%生成高斯谱杂波
xxi=conv（bxi）;   
xxq=conv（bxq）;   
xxi=xxi（coe_num*2+1:azi_num+coe_num*2）;%目的是去掉暂态响应
xxq=xxq（coe_num*2+1:azi_num+coe_num*2）;
xisigmac=std（xxi）;     
ximuc=mean（xxi）;       
yyi=（xxi-ximuc）/xisigmac;    
xqsigmac=std（xxq）;     
xqmuc=mean（xxq）;       
yyq=（xxq-xqmuc）/xqsigmac;    %归一化
sigmac=1.2 ;           %杂波的标准差
yyi=sigmac*yyi;        %使瑞利分布杂波具有指定的标准差
yyq=sigmac*yyq;