%%雷達信號基本處理流程?脈壓
%%平臺：R2016a

clear?all;close?all;clc;

%%?基本參數設置
f0?=?10e9;???????%載頻
Tp?=?10e-6;??????%脈沖寬度
B??=?10e6;???????%信號帶寬
fs?=?100e6;??????%采樣率
R0?=?3000;???????%目標初始距離
c??=?3e8;????????%光速
tr?=?2*R0/c;?????%目標點延時
k??=?B/Tp;???????%調頻斜率
N??=?4096;???????%采樣點數
t??=?（0:N-1）/fs;?%采樣時間步進

%%?發射信號
S0?=?rectpuls（t-Tp/2Tp）?.*?exp（1i*pi*k*（t-Tp/2）.^2）;
figure（1）;
subplot（211）;plot（t*c/2real（S0））;title（‘發射信號實部‘）;xlabel（‘距離/m‘）;
subplot（212）;plot（t*c/2imag（S0））;title（‘發射信號虛部‘）;xlabel（‘距離/m‘）;
%subplot（313）;plot（abs（fft（S0）））;

%%?發射信號頻譜
S0_fft?=?fft（S0N）;

%%?回波構造
S1?=?rectpuls（t-tr-Tp/2Tp）?.*?exp（1i*pi*k*（t-tr-Tp/2）.^2）?.*?exp（-1i*2*pi*f0*tr）;
figure（2）;
subplot（211）;plot（t*c/2real（S1））;title（‘回波函數實部‘）;xlabel（‘距離/m‘）;
subplot（212）;plot（t*c/2imag（S1））;title（‘回波函數虛部‘）;xlabel（‘距離/m‘）;
%subplot（313）;plot（abs（fft（S1）））;

%%?回波頻譜
f?=?fs/N?*?（-N/2:N/2?-?1）;
S1_fft?=?fft（S1N）;
figure（3）;
subplot（211）;plot（abs（S1_fft））;title（‘回波頻譜‘）;xlabel（‘頻率‘）;
subplot（212）;plot（fftshift（abs（S1_fft）））;title（‘回波頻譜‘）;xlabel（‘頻率‘）;

%%?脈沖壓縮
S_F?=?ifft（?S1_fft?.*?conj（S0_fft））;
figure（4）;
subplot（211）;plot（t*c/2abs（S_F））;title（‘脈壓‘）;xlabel（‘距離/m‘）;
subplot（212）;plot（t*c/2db（abs（S_F）/max（S_F）））;title（‘歸一化‘）;xlabel（‘距離/m‘）;

%%?頻域加窗
win?=?hamming（410）;
window?=?[zeros（11843）