%%demo?of?chirp?signal
clear?all
tao=30e-6;?????????????????%信號持續時間

B=5e6;?????????????????????%chirp?frequency?modulation?bandwidth?30MHz
f0=30e6;
K=B/tao;??????????????????????%調頻斜率
fs=40e6;??????????????????????%采樣頻率
Ts=1/fs;??????????????????????%采樣周期
N=tao*fs;?????????????????????%采樣點數?
t=-tao/2:Ts:tao/2;

%線性調頻信號的時頻譜圖分析
st=cos（pi*K*t.^2+2*pi*f0*t）;??%產生線性調頻信號
figure（1）
subplot（211）
plot（tst）;
xlabel（‘Time‘）;
title（‘線性調頻信號的實部‘）;
grid?on;
subplot（212）
freq=linspace（0fsN）;
plot（freqabs（fft（stN）））;
xlabel（‘Frequency?in?MHz‘）;
title（‘線性調頻信號頻譜‘）;
grid?on;

%數字下變頻?輸入I路和Q路正交信號
y1=cos（2*pi*f0*t）;
y2=sin（2*pi*f0*t）;
xi=st.*y1;
xq=st.*y2;
fft_xi=fftshift（abs（fft（xiN）））;
fft_xq=fftshift（abs（fft（xqN）））;
figure（2）
subplot（211）
plot（txi）;
title（‘輸入I信號時域