花如月 发表于 2008-1-8 21:38

如何查看窗函数的频谱

matlab产生的窗函数怎么查看频谱呢?都有哪些方法呢?最近看到一份资料上的方法如下:n=51;
=freqz(window,1);
subplot(1,2,1)
stem(window);
subplot(1,2,2);
plot(w/pi,20*log(abs(h)/abs(h(1))));使用的是freqz,觉得有些问题。请熟悉的朋友看看这样做是否理论上行得通呢?还有没有其他查看频谱特性的方法?先谢谢了

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

songzy41 发表于 2008-1-9 08:39

本帖最后由 wdhd 于 2016-9-8 14:54 编辑

原帖由 花如月 于 2008-1-8 21:38 发表
matlab产生的窗函数怎么查看频谱呢?都有哪些方法呢?最近看到一份资料上的方法如下:n=51;
=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的调用一般是
=freqz(b,a,N)
b是传递函数中分子的系数,a是传递函数中分母的系数,N是频谱单边时的频率个数(缺省为512),freqz说明给出计算的频率响应为:
                      jw                   -jw                     -jmw
          jwB(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 发表于 2008-1-9 08:47

回复 2楼 的帖子

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

witty01 发表于 2008-1-9 14:17

回复 2楼 的帖子

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

花如月 发表于 2008-1-16 14:27

回复 2楼 的帖子

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

hyl2323 发表于 2008-1-18 19:31

matlab有个设计窗函数的GUI工具---wintool,我用过,好用的。
页: [1]
查看完整版本: 如何查看窗函数的频谱