iirtest1.asv

黄金分割法

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,