用matlab设计低通滤波器来处理一张图片

Question

老师突然布置的一个作业，要求用matlab设计一个低通滤波器来处理随意的一张图片，要求有源程序，但matlab我不熟，所以求问有无人帮忙写下，或者说下哪本书的哪几页有此内容亦可。
好的会有加分。

Answer 1

%============IIR低通滤波器================
%============IIR低通滤波器================
Fp=1000;Fs=1200;Ft=1000;fp=Fp/Ft;fs=Fs/Ft;wp=2*tan(fp/2);ws=2*tan(fs/2);[n11,wn11]=buttord(wp,ws,1,100,'s'); %求低通滤波器的阶数和截止频率[b11,a11]=butter(n11,wn11,'s') %求S域的频率响应的参数[num11,den11]=bilinear(b11,a11,0.5);%双线性变换实现S域到Z域的变换[h,w]=freqz(num11,den11); %根据参数求出频率响应plot(w*2000/pi,abs(h));legend('用butter设计')%=================双线性变换法==============Fp=1000;Fs=1200;Ft=1000;fp=Fp/Ft;fs=Fs/Ft;wp=2*tan(fp/2);ws=2*tan(fs/2);[n11,wn11]=buttord(wp,ws,1,100,'s'); %求低通滤波器的阶数和截止频率[b11,a11]=butter(n11,wn11,'s') %求S域的频率响应的参数[num11,den11]=bilinear(b11,a11,0.5);%双线性变换实现S域到Z域的变换[h,w]=freqz(num11,den11); %根据参数求出频率响应plot(w*2000/pi,abs(h));legend('用butter设计')
%————————FIR低通————————%用窗函数设计低通滤波器的程序如下Fp=1000;Fs=1200;Ft=8000;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;rp=1;rs=50;p=1-10.^(-rp/20); %通带阻带波纹s=10.^(-rs/20);fpts=[wp ws];mag=[1 0];dev=[p s];[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev); %kaiserord 求阶数截止频率b21=fir1(n21,wn21,kaiser(n21+1,beta)); %由fir1设计滤波器[h,w]=freqz(b21,1); %得到频率响应plot(1250*w/pi,abs(h));title('FIR低通滤波器');

Answer 2

121212121