% 雷達回波線性調頻信號及其頻譜圖
close all
clear all
clc;
tao=1*10^(-4);    
f0=0;
B=1*10^6;
u=B/tao;
t=-tao/2:tao/1023:tao/2;
R=5000;                                           %雷達與物體距離 
c=3e8;
t0=2*R/c;
s=exp(j*2*pi*(f0*(t-t0)+0.5*u*(t-t0).^2));                                           %信號表達式
subplot(2,1,1)
plot(t,s,'r'),axis tight; 
xlabel('t(單位：秒)'),ylabel('y(單位：伏)'),title('回波信號時域');
grid on
S0=fft(s,1024);
S1=abs(S0);
S2=fftshift(S1);                                                     %信號幅頻譜
f1=(1023/tao)*(0:511)/1024;
f=[-1023/(tao*2),-fliplr(f1),f1(2:length(f1))];                               %頻率軸
subplot(2,1,2)
plot(f,S2/max(S2));
xlabel('f(單位：赫茲)'),ylabel('y(單位：伏)'),title('回波信號幅度譜');
grid on
axis tight