11-5.m

matlab图像处理源码

%用ftrans2函数设计一个通带在0.1和0.6间（量化频率，这里1.0对应采样频率的一半或π弧度）的近似圆对称二维通带滤波器。
%ftrans2函数由一FIR滤波器创建二维滤波器
%首先用信号处理工具箱函数remez设计一个一维FIR通带滤波器
subplot(121)
colormap(jet(64))
b = remez(10,[0 0.05 0.15 0.55 0.65 1],[0 0 1 1 0 0]);
[H,w] = freqz(b,1,128,'whole');
plot(w/pi-1,fftshift(abs(H)));
%应用函数ftrans2和缺省 McClellan变换矩阵创建近似圆对称滤波器
subplot(122)
h = ftrans2(b);
freqz2(h)