s107.m

简单的数字图像处理程序

f=[0 0.1 0.2 0.3 0.4 0.6 0.7 0.8 0.9 1];
m=[1 1 1 1 0 0 1 1 1 1];
%产生一个20阶的具有最小二乘线性相位的一维FIR带阻滤波器。
b=firls(20,f,m,'Hilbert');
[h0,w]=freqz(b,1,512,2);
plot(f,m,'--',w,abs(h0));%如图示 理想(虚线)和实际(实线)的频率响应曲线
%  在一维滤波器函数取理想带阻滤波器时，采用频率抽样法产生一个20阶的具有最小二乘线性相位的FIR带阻滤波器
h1=fsamp2(f,f,m,[5 5]);
figure(2),freqz2(h1,[32 32]);     %如图示  理想带阻滤波器
%  在一维滤波器函数取实际带阻滤波器时，采用频率抽样法产生一个20阶的具有最小二乘线性相位的FIR带阻滤波器
ff=1/512:1/512:1;
h2=fsamp2(ff,ff,h0,[5 5]);
figure(3),freqz2(h2,[32 32]);  % 如图示  实际带阻滤波器
%利用前面设计的滤波器进行图像分析
I=imread('3.jpg');
I1=rgb2gray(I);
fh1=filter2(h1,I1);  %利用h1进行滤波
figure(4),imshow(fh1);     %如图(b)所示
fh2=filter2(h2,I1); %利用h2进行滤波
figure(5),imshow(fh2);  %如图(c)所示