iirtest.m

黄金分割法

T=9.2e-6;
fs=73.33e6*2;
B=60.e6;
k=B/T;
N=round(T*fs);
n=-1*N/2:N/2;
f=exp(j*k*pi*(n/fs).^2);
subplot(3,1,1);
plot(real(f));
title('线性调频信号');
subplot(3,1,2);
plot(abs(fft(f)));
title('线性调频信号频谱')
f1=exp(-1*j*k*pi*((N/2-n)/fs).^2);
R=fft(f).*fft(f1);
r=ifft(R);
subplot(3,1,3);
plot(abs(r));
title('经匹配滤波后的信号');