让别人免费给你写程序基本上不可能,这个得花时间和精力。
clearclose
all
I
=
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('自编程序滤波图像')
%显示滤波后的图象
滤波器的理想频率响应函数为Hd(ejω),则其对应的单位脉冲响应为hd(n)=窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼hd(n)。由于hd(n)往往是无限长序列,且是非因果的迅搏裂,所以用窗函数。w(n)将hd(n)截断,并进行加权处理:
h(n)=hd(n)w(n)h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列,其频率响应函数H(ejω)为H(ejω)=用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N的取值。设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。
一般都选用Ⅰ型线性相位滤波器即滤波器阶数M为偶数,程序如下:
wp=ws=Ap=1As=100
dev=[Rp Rs]
[M,wc,beta,ftype]=kaiserord(f,a,dev)
M=mod(M,2)+M
plot(omega/pi,20*log10(abs(mag)))
运行程序可以得到滤波器的通阻带衰减,画出频率响应,若同阻带衰减不满足要求还可以使用滤波器的优化,一般使用的等波纹FIR进行优化。
扩展资料:
滤波器与机箱之间的一段连线银蠢会产生两种不良作用: 一个是机箱内部空间的电磁干扰会直接感应到这段线上,沿着电缆传出机箱,借助电缆辐射,使滤波器失效;另一个是外界干扰在被板上滤波器滤波之前,亩闭借助这段线产生辐射,或直接与线路板上的电路发生耦合,造成敏感度问题;
滤波阵列板、滤波连接器等面板滤波器一般都直接安装在屏蔽机箱的金属面板上。由于直接安装在金属面板上,滤波器的输入与输出之间完全隔离,接地良好,电缆上的干扰在机箱端口上被滤除,因此滤波效果相当理想。
参考资料来源:百度百科-滤波器
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)