example7_7.m

This complete matlab for neural network

%首先，创建包含了期望的频率响应的矩阵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)