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