ch10_2_1.m

来自「常用图像处理源码合集」· M 代码 · 共 8 行

%利用频率变换法生成一个2-D滤波器
b = firpm(10,[0 0.4 0.6 1],[1 1 0 0]);
h = ftrans2(b);
[H,w] = freqz(b,1,64,'whole');%得到b的1-D频率响应
colormap(jet(64))%图形颜色映射
plot(w/pi-1,fftshift(abs(H)))
figure, freqz2(h,[32 32])