用czt函数进行频率细化的程序遇到问题，求达人
我最近想用CZT这个函数，来细化频率，结果出现问题，发现细化出来的幅值正确，可是频率不正确。程序如下：
clc
clear
fs=1024;
N=1024;
n=0:N-1;
t=(n)/fs;
s=5*cos(2*pi*150*t)+10*cos(2*pi*200*t);
m=1024;
f1=100;
f2=300;
w=exp(-j*2*pi*(f2-1)/(m*fs));%%%表示为沿z平面的螺旋线点之间的比例，也表示为变换的步长
a=exp(j*2*pi*f1/fs);%%%%表示为螺旋线开始点，即变换起点
y=fft(s,m);
fy=(0:length(s)-1)*fs/length(s);
z=czt(s,m,w,a);
fz=(0:length(s)-1)*(f2-f1)/length(s)+f1;
subplot(2,1,1)
plot(fy(1:512),abs(y(1:512))*2/m);
title('fft');
subplot(2,1,2)
plot(fz,abs(z)*2/m);
title('czt');
结果如下：
问题解决了，我其中一句写错了。
w=exp(-j*2*pi*(f2-f1)/(m*fs));%%%表示为沿z平面的螺旋线点之间的比例，也表示为变换的步长
应该是f2-f1， 不是f2-1。 学习交流
页:
[1]