iirtest1.m

黄金分割法

T=9.2e-6;
fs=73.33e6*2;
B=60e6;
k=B/T;
N=round(T*fs);
N=round(N/2)*2;
n=-1*N/2:N/2;
f=exp(j*k*pi*(n/fs).^2);
f=[zeros(1,N/2),f,zeros(1,N/2)];
subplot(3,1,1);
plot(n/fs*(1e+6),real(f(N/2:3*N/2)));
title('线性调频信号');xlabel('t(us)');ylabel('x');
subplot(3,1,2);
plot((1:2*N+1)*fs/(2*N+1)/(1.0e6),abs(fft(f)));
title('线性调频信号');xlabel('f(MHz)');ylabel('abs(F)');
f1=exp(-1*j*k*pi*((-1*n)/fs).^2);
f1=[zeros(1,N),f1];
R=fft(f).*fft(f1);
r=ifft(R);
subplot(3,1,3);
plot(log(abs(r(N/2-200:N/2+200))));
title('经匹配滤波后的信号');