MATLAB基于汉宁窗的FIR的低通滤波器的源代码及注释

很常见的设计题目

给你一个例程，只需要改一改参数就行了

clear all;

f=[0 019 02 03 031 059 06 08 081 1];

% 给定频率轴分点；

m=[0 0 1 1 0 0 1 1 0 0];

% 给定在这些频率分点上理想的幅频响应

N1=30;

N2=90;

% 取两种不同的滤波器长度；

b1=fir2(N1,f,m);

b2=fir2(N2,f,m);

% 得到两个滤波器；

subplot(311);

stem(b1,'');grid;

subplot(312);

stem(b2,'');grid;

M=128;

[h1,w]=freqz(b1,1,M,1);

[h2,w]=freqz(b2,1,M,1);

subplot(313);

plot(w,abs(h1),'b-',w,abs(h2),'g-');grid;

其中，f是归依化以后的频率 通过数字滤波器的采样频率算出来，根据通带和阻带算好f和m就行了

看一看help，这个函数应该有窗函数的选择 默认情况下是汉明窗

希望能够帮到你

假设是32点的高斯窗，用以下matlab语句就可以看到高斯窗时域和频域的情况：

N = 32;

wvtool(gausswin(N));

或者：

window=gausswin(n);

[h,w]=freqz(window,1);

subplot(1,2,1)

stem(window);

subplot(1,2,2);

plot(w/pi,20log(abs(h)/abs(h(1))));

哈明窗（Hamming）

哈明窗也是余弦窗的一种，又称改进的升余弦窗。哈明窗与汉宁窗都是余弦窗，只是加权系数不同。哈明窗加权的系数能使旁瓣达到更小。分析表明，哈明窗的第一旁瓣衰减为一42dB．哈明窗的频谱也是由3个矩形时窗的频谱合成，但其旁瓣衰减速度为20dB/（10oct），这比汉宁窗衰减速度慢。哈明窗与汉宁窗都是很有用的窗函数。