怎么用matlab编p-m滤波的程序

imread('eight.tif')

J = imnoise(I,'salt &pepper',0.02)

K = medfilt2(J)

imshow(J)title('岁基噪声干扰图像')

figure, imshow(K)title('medfilt2滤波图像')

X=Ja=2b=2

k=floor(a*b/2)+1

[M,N]=size(X)

uint8 Y=zeros(M,N)

funBox=zeros(a,b)

temp=zeros(a*b)

for i=1:M-a

for j=1:N-b

funBox=X(i:i+a,j:j+b)

temp=funBox(:)

tempSort=sort(temp)

Y(i,j)=tempSort(k)

end

end

figure, imshow(Y)title('自编程序滤波图像')

clear

close all

c=imread('123.png')%把彩色图片转化成灰度图片，256级

figure,imshow(c),title('原始图象')%显示原始图象

g=imnoise(c,'gaussian',0.1,0.002) %加入高斯噪声

figure,imshow(g),title('加入高斯噪声之后的图象') %显示加入高斯噪声之后的图象

%实验步骤二：用系统预定义滤波器进行均值滤波

n=1

A=fspecial('average',n)%生成系统预定义的3X3滤波器

Y=filter2(A,g)/255 %用生成的滤波器进行滤波,并归一化

figure,imshow(Y),title('系统函乎销谨数滤波图像')%显示滤波后的图象

a(1:n,1:n)=1 %a即n×n模板,元素全是1

p=size(g) %输入图像是p×q的,且p>n,q>n

x1=double(g)

x2=x1

%A(a:b,c:d)表示A矩阵的第a到b行,第c到d列的所有元素

for i=1:p(1)-n+1

for j=1:p(2)-n+1

c=x1(i:i+(n-1),j:j+(n-1)).*a %取出x1中从(i,j)开始的n行n列元素与模板相乘

s=sum(sum(c))%求c矩阵(即模板)中斗颂各元素之和

x2(i+(n-1)/2,j+(n-1)/2)=s/(n*n)%将模板各元素的均值赋给模板中心位置的元素

end

end

%未被赋值的元素取原值

d=uint8(x2)

%实验步骤三:用自己的编写的函数进行均值滤波

%调用自编函数进行均值滤波，n为模板大小

figure,imshow(d),title('自编程序滤波图像')%显示滤波后的图象

冲击响应不变法函数

[bz,az]=impinvar(b,a,Fs)

[bz,az]=impinvar(b,a)

例如：取采样频率f=1KHz,用双线性变换法设计五阶Butterworth低通数字滤波器，绘出模拟滤波器与数字滤波器的幅频与相频特性，MATLAB程序如下：

[z,p,k]=buttap(5) % 设计五阶Butterworth低通模拟滤波器原型

[zd,pd,kd]=bilinear(z,p,k,1000)%双线性变换得到低通数字滤波器

[b,a]=zp2tf(zd,pd,kd)%滤波器类型转换

w=128

freqs(b,a,w)

figure

freqz(b,a,w)

实例：1、设带通滤波镇拆器的滤波器中心频率为W0=2KHz,带宽为BW=100Hz, 取采样频率f=10kHZ,用脉冲相应不变法设计，设计五阶带通Butterworth数字滤波器,绘出数字滤波器的频谱特性

[z,p,k]=buttap(5)

[b,a]=zp2tf(z,p,k)

w=128

w0=2000

[bt,at]=lp2bp(b,a,w0,10000)

[bz,az]=impinvar(b,a,w)

freqz(bt,at,w)

2、直接设计五阶butterworth带通滤波器，绘出频谱图。（高端与低端截止频率分别为0.2和0.9）

figure

w=[0.2,0.9]

[b,a]=butter(5,w)

freqz(b,a)

3、设高通截止频率为w0=10000Hz, 取采样频率f=20000,用指正双线性变换法设计六阶高通Butterworth数字滤波器,绘出数字滤波器的频谱特性

[z,p,k]=besselap(6)

[b,a]=zp2tf(z,p,k)

w=128

w0=10000

[bt,at]=lp2hp(b,a,w0)

[bz,az]=bilinear(bt,at,20000)

freqz(bz,az,w)

4、取采样频率f=100Hz,用双线性变换法设计五阶Butterworth低通数字滤波器，绘出模拟滤波器与数字滤波器的幅频与相频特性

　御逗枣　[z,p,k]=buttap(5)

[zd,pd,kd]=bilinear(z,p,k,100)

[b,a]=zp2tf(zd,pd,kd)

w=128

freqs(b,a,w)

figure

freqz(b,a,w)

具体的得根据情况自己确定

---