clutter 模拟杂波的生成

模拟杂波的生成，包括高斯杂波，韦布尔杂波和K分布杂波，对数正态分布杂波

azi_num=2000;
fr=1000;
lamda0=0.05;
sigmav=1.0;
sigmaf=2*sigmav/lamda0;
rand（‘state‘sum（100*clock））;
d1=rand（1azi_num）;
rand（‘state‘7*sum（100*clock）+3）;
d2=rand（1azi_num）;
xi=1*（sqrt（-2*log（d1））.*cos（2*pi*d2））;
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
%Gaussion clutter generation
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）;
vmuc=2;
xisigmac=std（xxi）;
ximuc=mean（xxi）;
xxi=（xxi-ximuc）/xisigmac;
xqsigmac=std（xxq）;
xqmuc=mean（xxq）;
xxq=（xxq-xqmuc）/xqsigmac;
xdata=xxi+1j*xxq;
tmpdat=randn（1azi_num）;
[ba]=butter（50.01）;
sk_dat=filter（batmpdat）;
sk_dat=sk_dat/std（sk_dat）;

%solving nonlinear equations next
max_z=6;
step=0.005;
table_z=0:step:max_z;
table_s=nonline_eq_sirp（table_zvmuc）;
for n=1:azi_num
    index=floor（abs（sk_dat（n））/max_z*length（table_z）+1）;
    sk_dat（n）=table_s（index）;
end
ydata=xdata.*sk_dat;

figuresubplot（211）plot（real（ydata））;xlabel（‘时间‘）;ylabel（‘幅度‘）;title（‘k分布杂波时域图，实部‘）;
subplot（212）plot（imag（ydata））;xlabel（‘时间‘）;ylabel（‘幅度‘）;title（‘k分布杂波时域图，虚部‘）;

%PDF
num=100;
maxdat=max（abs（ydata））;
mindat=min（abs（ydata））;
NN=hist（abs（xdata）num）;
xpdf1=num*NN/（（sum（NN））*（maxdat-mindat））;
%xpdf1=getpdf（abs（ydata）nummaxdatmindat）;
xaxis1=mindat:（maxdat-mindat）/num:maxdat-（maxdat-mindat）/num;
alpha=sqrt（std（xdata）.^2./（2*vmuc））;
%th_val=2*（（xaxis1/（2*alpha））.^vmuc）.*ESSELK（（vmuc-1）xaxis1/alpha）./（alpha*gamma（vmuc））;

figuresubplot（211）plot（xaxis1xpdf1）;
%holdplot（xaxis1th_val‘:r‘）;
title（‘杂波幅度分布‘）;xlabel（‘幅度‘）;ylabel（‘概率密度‘）;

signal=ydata;
signal=signal-mean（ydata）;

M=256;
psd_dat=real（pburg（real（signal）16Mfr））;
psd_dat=psd_dat/（max（psd_dat））;
freqx=0:0.5*M;
freqx=freqx*fr/M;
subplot（212）plot（freqxpsd_dat）;
title（‘杂波频谱‘）;xlabel（‘频率Hz‘）;ylabel（‘功率谱密度‘）;
%理想高斯谱线
powerf=exp（-freqx.^2/（2*sigmaf.^2））;
hold;
plot（freqxpowerf‘:r‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2232  2014-05-29 15:37  clutter_generate.m

     文件       1771  2014-05-29 15:37  Gaussion.m

     文件        295  2014-05-19 13:38  glp.m

     文件       1501  2014-05-29 15:37  LN.m

     文件        757  2014-05-17 00:40  nonline_eq_sirp.m

     文件       1683  2014-05-29 15:38  weibull.m

----------- ---------  ---------- -----  ----

                 8239                    6