example7_8.m

matlab图像处理源代码 大家多多分享

%首先，创建一个一维FIR低通滤波器，变换矩阵是可以根据需要定义的，这里定义了一个截
%断频
%率为0.6的十阶低通滤波器和一个截断频率为0.5的十四阶低通滤波器
b1=remez(10,[0 0.4 0.6 1],[1 1 0 0]);
[H1,w]=freqz(b1,1,64,'whole');
colormap(jet(64))
plot(w/pi-1,fftshift(abs(H1)))
b2=remez(14,[0 0.3 0.5 1],[1 1 0 0]);
[H2,w]=freqz(b2,1,64,'whole');
colormap(jet(64))
plot(w/pi-1,fftshift(abs(H2)))
%用于频率变换的两个一维滤波器的频率响应见图7-18和图7-19
%然后，用默认的McClellan变换生成期望的近似圆对称的滤波器
h1=ftrans2(b1);
figure,freqz2(h1,[32 32])
h2=ftrans2(b2);
figure,freqz2(h2,[32 32])