lfm_signal.m

开发系统的程序

% 演示线性调频信号及其频谱图
tao=1*10^(-4);    
f0=0;
B=1*10^6;
u=B/tao;
t=-tao/2:tao/1023:tao/2;
s=exp(j*2*pi*(f0*t+0.5*u*t.^2));         %信号表达式
subplot(1,2,1)
plot(t,s,'r'),axis([-tao/2 tao/2 -1 1]); 
grid on
S0=fft(s,1024);
S1=abs(S0);
S2=fftshift(S1);         %信号频谱
ll=length(S2);
f1=(1023/tao)*(0:511)/1024;
f=[-1023/(tao*2),-fliplr(f1),f1(2:length(f1))];      %频率轴
subplot(1,2,2)
plot(f,S2/max(S2));
grid on
axis([-B B 0 1]);