chirp.m

chirp信号matlab仿真

%----------------------------------------
%    不加窗进行旁瓣抑制
%----------------------------------------
clear all;
close all;
fs=100e6;
T=5e-6;
B=10e6;
k=B/T;
t=[-T/2:1/fs:T/2-1/fs];
Ns=length(t);
f=[-fs/2:fs/Ns:fs/2-fs/Ns];
y=exp(j*pi*k*t.^2);
Y=fft(y);
H=conj(Y);
h=ifft(H);
Y1=fftshift(fft(y));
% delay=10;
Sig=Y.*H;
sig=ifftshift(ifft(Sig));
figure,subplot(2,1,1),plot(t,real(y)),
title('输入信号实部');
subplot(2,1,2),plot(t,imag(y)),
title('输入信号虚部');
figure,subplot(2,1,1),plot(f,abs(Y1)),xlabel('f(单位：Hz)'),
ylabel('y(单位：伏)'),title('Chirp信号频谱图)');axis([-5e7,5e7,0,80]);
subplot(2,1,2),plot(f,angle(Y1)),title('输入信号相位谱');
figure,plot(t,abs(sig)),xlabel('单位t(秒)'),ylabel('单位v（伏）'),
title('输出信号');axis([-2.5e-6,2.5e-6,0,800]);
figure,subplot(2,1,1),plot(t,real(sig)),title('输出信号实部'),
subplot(2,1,2),plot(t,imag(sig)),title('输出信号虚部');
figure,subplot(2,1,1),plot(f,abs(fftshift(Sig))),xlabel('单位f(Hz)'),ylabel('|H（f）|'),title('输出信号频谱');
subplot(2,1,2),plot(f,angle(fftshift(Sig))),xlabel('单位f(Hz)'),ylabel('Phase(f)'),title('输出信号相位谱');
figure,subplot(2,1,1),plot(real(h)),title('滤波器信号实部');
subplot(2,1,2),plot(imag(h)),title('滤波器信号虚部');
figure,subplot(2,1,1),plot(f,abs(fftshift(H))),title('滤波器频谱');
% subplot(2,1,2)