为什么说chirp-z变换比fft变换更适合对窄带信号进行变换
本帖最后由 wdhd 于 2016-9-20 13:48 编辑程序
clear;clc
t=0:0.001:1;
h=3*sin(2*pi*120*t)+2*cos(2*pi*280*t);
Fs=1000;
f1=100;
f2=300;
m=length(h);
w=exp(-j*2*pi*(f2-f1)/(m*Fs));
a=exp(j*2*pi*f1/Fs);
y=fft(h,m);
z=czt(h,m,w,a);
fy=(0:m-1)*Fs/m;
fz=(0:m-1)*(f2-f1)/m+f1;
subplot(2,1,1)
plot(fy,abs(y))
title('fft')
subplot(2,1,2)
plot(fz,abs(z))
title('chirp z')
这个图能看出什么 请简单描述一下 答辩用
页:
[1]