信号基础知识----线性调频信号LFM
%关于线性调频信号（LFM）
%参考书目：声呐技术，第二章P33
SRE实战 互联网时代守护先锋，助力企业售后服务体系运筹帷幄！一键直达领取阿里云限量特价优惠。clc;close all;clear all;
%参数----------------------------------
f0=1000;
T=5;%时宽，时间总长
B=10;%带宽
fs=100;%采样频率
Ts=1/fs;%采样时间间隔
N=T/Ts;%采样点个数
k=B/T;%调频斜率
t=linspace(0,T,N);
y=exp(1i*(2*pi*f0*t+pi*k*t.^2)); %线性调频信号（LFM）的时间函数
%作图---------------------------------------------
subplot(211);
plot(t,real(y));title('LFM信号时域波形');xlabel('时间/s');ylabel('幅度');
Y=fftshift(fft(y));
f=linspace(0,fs,N);
subplot(212);
plot(f,abs(Y));title('LFM信号频谱');xlabel('频率/Hz');ylabel('幅度');
%fftshift的用处-------------------------------------
figure;
subplot(211);
plot(abs(fft(y)));
xlabel('\omega/\pi');ylabel('|e^j^\omega/|');title('fft变换');
f=linspace(0,fs,N);
subplot(212);
plot(f,abs(Y));title('LFM信号频谱');xlabel('频率/Hz');ylabel('幅度');title('利用了fftshift的命令');
%三维图像--------------------------------------------------------------------------------------------------------------
t=-T:0.05:T; %频移分辨力rou=0.88/T
v=-B:0.05:B; %时延分辨力rou=0.88/B
[t,v]=meshgrid(t,v); %401*201
r=pi*(k*t+v).*(T-abs(t));
X=sin(r)./r.*(T-abs(t)); %模糊函数，假设A=1
figure;
surf(t,v,abs(X))
关于fftshift还是好好研究一下：http://www.cnblogs.com/limanjihe/p/10014142.html
