11-6.m

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

%用fsamp2函数设计一个近似对称的二维带通滤波器，通带范围[0.1 0.5]
%创建一个包含期望带通响应的矩阵Hd
%首先用函数freqspace创建频率范围向量f1和f2
subplot(121)
[f1,f2] = freqspace(21,'meshgrid');
Hd = ones(21); 
r = sqrt(f1.^2 + f2.^2);
Hd((r<0.1)|(r>0.5)) = 0;
colormap(jet(64))
mesh(f1,f2,Hd)
%设计一个通过这个响应的滤波器
subplot(122)
h = fsamp2(Hd);
freqz2(h)