%%?模擬產(chǎn)生勻速運動目標??
%?先做多普勒模糊補償，后做CZT補償距離徙動????
clear?all;clc;?close?all;

?%%?雷達參數(shù)
?fs=15e6;???????????????????????????%?基帶采樣率
?B=10e6;????????????????????????????%?基帶帶寬
?ts=1/fs;
?tp=200e-6;?????????????????????????%?脈寬
?kr=B/tp;???????????????????????????%?調(diào)頻斜率
?Tr=1e-3;???????????????????????????%?脈沖重復周期
?PRF=1/Tr;
?c=3e8;
?fc=3e8;????????????????????????????%?載頻
?lamda=c/fc;
?N=200;?????????????????????????????%?脈沖數(shù)
?%%?目標參數(shù)
?N01=0:N-1;
?R0=5e4;????????????????????????????%?初始距離?m
?v=499;?????????????????????????????%?徑向運動速度m/s
?R=R0-v*N01*Tr;?????????????????????%?目標與雷達徑向距離?m
?fd=2*v/lamda;??????????????????????%?多普勒頻率??Hz
?
?%%回波信號
?t0=0:ts:tp-ts;?????????????????????????????%發(fā)射序列
?Tt=exp（j*pi*kr*（t0-tp/2）.^2）;??????????????%發(fā)射波形
?tt=tp:ts:（Tr-ts）;??????????????????????????%采樣序列
?n_fft=2^nextpow2（length（tt）-1）;????????????%頻域脈壓FFT點數(shù)
?st=zeros（length（tt）N）;????????????????????%回波波形
for?m=1:N
?????tm（m）=2*R（m）/c-floor（2*R（m）/c/Tr）*Tr;
?????st（ceil（tm（m）*fs）:ceil（tm（m）*fs）+tp*fs-1m）=exp（j*pi*kr*（t0-tm（m）+ceil（tm（m）*fs）*ts-tp/2）.^2?-j*2*pi*2*R（m）/lamda）;???%+j*2*pi*fd*t0
end
%?SNR=-18;?%輸入信噪比
%?st=noisegen（stSNR）;
%figure;plot（real（st（:1）））;
?
%%?脈壓
fft_ht=fft（Ttn_fft）;???????????????????????%匹配函數(shù)
fft_st=fft（stn_fft）;
Array_ones=ones（N1）;
fft_ht_N=Array_ones*fft_ht;
xt=fft_st.*conj（fft_ht_N.‘）;????????????????%共軛相乘
%?yt=ifft（xt）;????????????????????????????????%返回到時域，脈壓
%?figure;contour（abs（yt））;?
%?%?xlabel（‘距離單元‘‘fontsize‘mesh_TextFontSize）;ylabel（‘脈沖數(shù)‘‘fontsize‘mesh_TextFontSize）
%?%?set（gca‘fontname‘‘Helvetica‘‘fontsize‘mesh_NumFontSize）;
%?%?set（gcf‘windowstyle‘‘normal‘）;
%?%?set（gcf‘Position‘[500300mesh_widthmesh_height]）;
%?title（‘脈壓結果‘）

%%?MTD
%?inte_xt=fft（yt.‘）;??????????????????????????%脈間FFT
%?figure;mesh（abs（inte_xt））;
%?%?xlabel（‘距離單元‘‘fontsize‘mesh_TextFontSize）;ylabel（‘多普勒單元‘‘fontsize‘mesh_TextFontSize）;
%?%?zlabel（‘幅度‘‘fontsize‘mesh_TextFontSize）
%?%?set（gca‘fontname‘‘Helvetica‘‘fontsize‘mesh_NumFontSize）;
%?%?grid?on;
%?%?set（gcf‘windowstyle‘‘normal‘）;
%?%?set（gcf‘Position‘[500300mesh_widthmesh_height]）;
%?title（‘未補償‘）

%%?keystone距離徙動補償
N01=0:N-1;
k1=floor（fd/PRF）;?????%?多普勒模糊數(shù)，速度模糊部分的補償
%%%%?頻譜正半部分補償
for?l=0:min（floor（n_fft*B*1.15/fs/2）n_fft/2）-1??????%?1:min（floor（n_fft*B*1.1/fs/2）n_fft/2）-1??選擇頻譜主要能量?1.1*n_fft*B/fs/2
????phi=（fc+l*fs/n_fft）/fc;?????????????????????????%?（fc+l*fs/n_fft）/fc?為對應快時間頻率（l）處的系數(shù)，大小=1左右
????f_amb_xt（l+1:）=xt（l+1:）.*exp（-j*2*pi*k1*phi*N01）;???%先做多普勒模糊補償????????
????w=exp（-j*2*pi*phi/N）;
????Czt_KT（l+1:）?=?czt（f_amb_xt（l+1:）Nw）;???????%?慢時間維CZT實現(xiàn)距離徙動補償?
end
%%%%?頻譜付半部分補償
for?l=max（-floor（n_fft*B*1.15/fs/2）-n_fft/2）:-1?????%?max（-floor（n_fft*B*1.1/fs/2）-n_fft/2）:-1??選擇頻譜主要能量?1.1*n_fft*B/fs/2
????phi=（fc+l*fs/n_fft）/fc;???????????????
????f_amb_xt（l+1+n_fft:）=xt（l+1+n_fft:）.*exp（-j*2*pi*k1*phi*N01）;?????%先做多普勒模糊補償????
????w=exp（-j*2*pi*phi/N）;
????Czt_KT（l+1+n_fft:）?=?czt（f_a

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????102504??2017-04-15?05:39??CZT仿真.docx
?????文件?????????817??2017-04-15?05:39??tryCZT2.m
?????文件????????3450??2017-04-15?05:39??keystone1.m