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;??