如何查看窗函数的频谱

花如月

matlab产生的窗函数怎么查看频谱呢？都有哪些方法呢？最近看到一份资料上的方法如下:

1. n=51;
   
2. [h,w]=freqz(window,1);
   
3. subplot(1,2,1)
   
4. stem(window);
   
5. subplot(1,2,2);
   
6. plot(w/pi,20*log(abs(h)/abs(h(1))));

复制代码

使用的是freqz，觉得有些问题。请熟悉的朋友看看这样做是否理论上行得通呢？还有没有其他查看频谱特性的方法？先谢谢了

[ 本帖最后由 花如月 于 2008-1-8 21:39 编辑 ]

songzy41

原帖由 花如月 于 2008-1-8 21:38 发表
matlab产生的窗函数怎么查看频谱呢？都有哪些方法呢？最近看到一份资料上的方法如下:n=51;
[h,w]=freqz(window,1);
subplot(1,2,1)
stem(window);
subplot(1,2,2);
plot(w/pi,20*log(abs(h)/abs(h(1))));
使用的是freqz，觉得有些问题。请熟悉的朋友看看这样做是否理论上行得通呢？

上程序这样做求窗函数的频谱没有问题。要了解它为什么是正确的，首先要了解freqz函数。freqz主要用于求数字滤波器的频响曲线。freqz的调用一般是
[h,w]=freqz(b,a，N)
b是传递函数中分子的系数，a是传递函数中分母的系数，N是频谱单边时的频率个数（缺省为512），freqz说明给出计算的频率响应为：
                      jw                   -jw                     -jmw
          jw  B(e)    b(1) + b(2)e + .... + b(m+1)e
     H(e) = ---- = ------------------------------------
                   jw                      -jw                     -jnw
            A(e)        a(1) + a(2)e + .... + a(n+1)e
当a=1时，只留下分子的表达式，它相当于是对系数b作DFT或FFT。在楼主调用为freqz(window,1)，相当于FFT(b,1024)的傅里叶变换，即求了频谱。

zhangnan3509

谢谢宋老师的指导，是我的问题差点误导了如月啊:@L

witty01

问题理解的好深刻啊，又学到了，谢谢:lol

花如月

感谢宋老师的解答，收获颇多

hyl2323

matlab有个设计窗函数的GUI工具－－－wintool，我用过，好用的。