花如月 发表于 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]
查看完整版本: 如何查看窗函数的频谱