用czt函数进行频率细化的程序遇到问题，求达人

sanmuhaifeng

我最近想用CZT这个函数，来细化频率，结果出现问题，发现细化出来的幅值正确，可是频率不正确。
程序如下：

1. clc
   
2. clear
   
3. 
   
4. fs=1024;
   
5. N=1024;
   
6. n=0:N-1;
   
7. t=(n)/fs;
   
8. s=5*cos(2*pi*150*t)+10*cos(2*pi*200*t);
   
9. 
   
10. m=1024;   
    
11. f1=100;
    
12. f2=300;
    
13. w=exp(-j*2*pi*(f2-1)/(m*fs));  %%%表示为沿z平面的螺旋线点之间的比例，也表示为变换的步长
    
14. a=exp(j*2*pi*f1/fs);  %%%%表示为螺旋线开始点，即变换起点
    
15. 
    
16. y=fft(s,m);
    
17. fy=(0:length(s)-1)*fs/length(s);
    
18. 
    
19. z=czt(s,m,w,a);
    
20. fz=(0:length(s)-1)*(f2-f1)/length(s)+f1;
    
21. 
    
22. subplot(2,1,1)
    
23. plot(fy(1:512),abs(y(1:512))*2/m);
    
24. title('fft');
    
25. subplot(2,1,2)
    
26. plot(fz,abs(z)*2/m);
    
27. title('czt');
    

复制代码

结果如下：

sanmuhaifeng

问题解决了，我其中一句写错了。
w=exp(-j*2*pi*(f2-f1)/(m*fs));  %%%表示为沿z平面的螺旋线点之间的比例，也表示为变换的步长
应该是f2-f1， 不是f2-1。

tommytom

学习交流