clc;
close?all;
clear?all;

T=10.e-6;??????????????????????????????????????????%脈寬10us
B=100.e6;??????????????????????????????????????????%信號帶寬100MHz
fs=120.e6;?????????????????????????????????????????%采樣頻率
K=B/T;?????????????????????????????????????????????%調頻斜率
r1=9850;??????????????????????????????????????????%目標1距離
r2=10150;??????????????????????????????????????????%目標2距離
c=3.e8;????????????????????????????????????????????%光速
tao1=2*r1/c;
tao2=2*r2/c;
N=2*T*fs;??????????????????????????????????????????%采樣點數
t=linspace（-TTN）;
f=（-N/2:N/2-1）.‘*fs/（N-1）;
tran_signal=exp（pi*j*K*t.^2）.*（abs（t）<=T/2）;???????%產生線性調頻信號
tran_signal1=（abs（t）<=T/2）;?
t0=linspace（-T+tao1T+tao1N）;
rec_signal1=exp（pi*j*K*（t0-tao1）.^2）.*（abs（t0-tao1）<=T/2）;
rec_signal2=exp（pi*j*K*（t0-tao2）.^2）.*（abs（t0-tao2）<=T/2）;
rec_signal=rec_signal1+rec_signal2;????????????????%產生回波信號??

?figure;
?plot（f/1e6fftshift（abs（fft（rec_signal））））;grid?on;xlabel（‘MHZ‘）;title（‘頻譜‘）;

pulse_compression1=fftshift（ifft（fft（rec_signal1）.*conj（fft（tran_signal））））;
pulse_compression2=fftshift（ifft（fft（rec_signal2）.*conj（fft（tran_signal））））;
pulse_compression=fftshift（ifft（fft（rec_signal）.*conj（fft（tran_signal））））;
pulse_compression3=ifft（fftshift（fft（rec_signal））.*（exp（pi*j*f.^2/K）.*（abs（f）
figure
plot（t0*c/2rec_signal）;
figure;
plot（t0*c/2abs（pulse_compression））;grid?on;xlabel（‘r/m‘）;title（‘脈壓后‘）;