s105.m

简单的数字图像处理程序

f=[0 0.3 0.4 0.6 0.7 1];
m=[0 0 1 1 0 0];
%产生一个20阶的Parks-McClellan一维带通FIR滤波器
b=remez(20,f,m);
[h0,w]=freqz(b,1,512);
%freqz计算由矢量a(这里a=1)和b构成的数字滤波器：
figure(1),plot(f,m,'--',w/pi,abs(h0));
%理想（虚线）和实际（实线）的频响曲线
%频率变换法产生20阶的Parks-McClellan二维带通FIR滤波器
h=ftrans2(b);
figure(2),freqz2(h,[32 32]);    
%利用前面设计的滤波器进行图像分析
%读取原始图像
f=imread('3.jpg');
f=rgb2gray(f);
figure(3),imshow(f);           
fh=filter2(h,f);
figure(4),imshow(fh)