11-8.m

这个确实有用 大家可以任意下载 完全免费的

%上面由单一的一维窗函数设计滤波器。
%这里示例来说明两个一维窗函数的情况。
[f1,f2] = freqspace(21,'meshgrid');
Hd = ones(21); 
r = sqrt(f1.^2 + f2.^2);
Hd((r<0.1)|(r>0.5)) = 0;
%Hd矩阵生成的方法和上面一样
%下面用hamming窗和blackman窗来设FIR滤波器
h = fwind1(Hd,hamming(21),blackman(21));
freqz2(h)
title('hamming窗和blackman窗');