%==產生LFM信號-分析其頻譜====%?
clc;clear?all;close?all;warning?off;??????????????
snr=0;???????????????????????????%信噪比
fc=0e10;?????????????????????????%?載波頻率;
c=3e8;???????????????????????????%?真空光速;
Tr=750e-6;???????????????????????%?脈沖重復周期s;
fd=1.6e6;??????????????????????????%?DOPPLOR?頻率
td=2e-6;?????????????????????????%?時間延遲s
Tp=1e-5;?????????????????????????%?LFM調制脈沖寬度s；
B=8e6;???????????????????????????%?LFM調制脈沖帶寬hz；
fs=2*B;??????????????????????????%?LFM采樣頻率；
t=-Tp/2:1/fs:Tp/2;N=length（t）;
lfm=10^（snr/10）*exp（j*2*pi*（fc*t+0.5*B/Tp*t.^2））;
fft_lfm=fftshift（fft（lfmN））;
LFM_Back=10^（snr/10）*exp（j*2*pi*（fc*（t-td）+0.5*B/Tp*（t-td）.^2））.*exp（-j*2*pi*fd*t）;;
figure（1）
subplot（211）plot（-Tp/2+Tp/N:Tp/N:Tp/2real（lfm））;xlabel（‘LFM脈沖寬度（s）‘）;?ylabel（‘幅度‘）;title（‘實LFM信號波形‘）;axis（[-TpTp-22]）
subplot（212）plot（-Tp/2+Tp/N:Tp/N:Tp/2real（LFM_Back））;xlabel（‘LFM脈沖寬度（s）‘）;?ylabel（‘幅度‘）;title（‘實LFM信號回波波形‘）;axis（[-TpTp-22]）

figure（2）
?????subplot（211）plot（-0.5+1/N:1/N:0.5abs（fft_lfm（1:N）））;xlabel（‘歸一化頻率‘）;?ylabel（‘幅度‘）;title（‘LFM信號頻域幅度‘）;
?????subplot（212）plot（-0.5+1/N:1/N:0.5angle（fft_lfm（1:N）））;xlabel（‘歸一化頻率‘）;?ylabel（‘相角‘）;title（‘LFM信號頻域相角‘）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%?頻域匹配濾波?%?放大可查看半功率寬度=1/B
out_mat=fftshift（ifft（fft（（LFM_Back））.*conj（fft（（lfm）））））?;
figure（3）
subplot（211）plot（-Tp/2:1/fs:Tp/220*log10（abs（out_mat）/max（abs（out_mat））））;title（‘匹配輸出LFM‘）;axis（[-Tp/2Tp/2-2000]）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%?頻域匹配濾波?%?放大可查看半功率寬度=1/B???加窗
out_mat=fftshift（ifft（fft（（LFM_Back））.*conj（fft（（lfm.*hamming（N）‘）））））?;
subplot（212）plot（-Tp/2:1/fs:Tp/220*log10（abs（out_mat）/max（abs（out_mat））））;title（‘匹配輸出LFM‘）;axis（[-Tp/2Tp/2-2000]）;

%%%%%%%%%%%%%%%%%%%%%%%%%%%?時域匹配濾波?
matcher=conj（fliplr（lfm））;????????????????%?時域匹配濾波?
weight=conj（fliplr（lfm.*hamming（N）‘））;????%?時域匹配濾波?加窗
out_matcher=conv（LFM_Backmatcher）;
out_weight=conv（LFM_Backweight）;
figure（4）
subplot（211）plot（-Tp:1/fs:Tp20*log10（abs（out_matcher）/max（abs（out_matcher））））;title（‘時間匹配輸出LFM‘）;axis（[-Tp/2Tp/2-2000]）;
subplot（212）plot（-Tp:1/fs:Tp20*log10（abs（out_weight）/max（abs（out_weight））））;title（‘時間加窗匹配輸出LFM‘）;axis（[-Tp/2Tp/2-2000]）;



%%%%%%%%%%%%%%%%%%%%%%%%%%%?頻率時間聯合估計
fdd=（B/Tp）*Tp;
ii=0;
NN=20;%?頻域采樣數目
out_mat=zeros（12*N-1）;
for?fdx=-fdd:fdd/NN:fdd;
????ii=ii+1;
????matcher=conj（fliplr（lfm））.*exp（-j*2*pi*fdx*t）;?
????out_mat（ii:）=conv（LFM_Backmatcher）;
end
figure（5）;[FT]=meshgrid（-Tp:1/fs:Tp-fdd:fdd/NN:fdd）;mesh（FTabs（out_mat））;title（‘頻率時間聯合估計‘）;axis（[-TpTp-fddfddmin（min（abs（out_mat）））max（max（abs（out_mat）））]）;grid?on;
figure（6）;contour（FTa