b6.m

来自「MATLAB中的所有练习源程序,软件包形式,对学习MATLAB很有帮助」· M 代码 · 共 13 行

%首先，创建包含了期望的频率响应的矩阵Hd
[f1,f2] = freqspace(21,'meshgrid');
Hd = zeros(21); 
Hd(7:15,7:15) = 1;
Axis([-1 1 –1 1 0 1.2]),
colormap(jet(64))
mesh(f1,f2,Hd)

%期望的滤波器的频率响应见图7-16
%然后，设计通过此响应的滤波器
h = fsamp2(Hd);
figure,freqz2(h)