找回密码
 注册
关于网站域名变更的通知
查看: 482|回复: 1
打印 上一主题 下一主题

基于matlab的Butterworth低通滤波程序分享

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-7-8 15:39 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
clear all;, }. q) q) j5 n# F: q& f
7 \) @; M0 d7 m2 G& w1 b
close all;
% B9 l! |/ }* E$ v6 H/ e/ t$ y# p7 H" v  k
P=input('inputimage data:');
6 h6 A; r3 X' f- {+ P5 ]/ J" \1 k- W5 I- ]( K
d=size(P);; h9 }( O: |& t( C& j( J7 Q6 K6 C' K

' X& w  K$ |/ s5 M6 H/ qif(d(3)>1)" B% O' V* i% y% ]
: ^( K( j) r$ i2 O# u
    P=rgb2gray(P);
2 K  `3 l3 x9 w  t* F" }- E& U' a
/ O8 U. x# [% ]: N* @9 Cend8 [4 s5 U8 O: z* \9 N" u

# i# G% T3 G  Z( y5 J: A; v/ }subplot(131);imshow(P);
, o" S' u0 S' }- g( C" D* ?% I' k& {: r, P+ v* n3 ]
title('原始图像');
* _3 A1 |9 G' {' M) Z, z. K: w, G9 U% @: Q- f. e
P1=imnoise(P,'salt& pepper',0.02);; V8 k# }& P1 p) k) |7 w
9 V& H7 s7 u4 P" Y( i
subplot(132)( E! V$ [- x  A7 O9 B

; v) s& [% `8 g1 U( n, timshow (P1);" u; R$ j+ w8 r" }7 [( i
3 O5 h' g3 z7 G* ~/ o3 U
title('加入椒盐噪声后的图像');
( D) Z9 ^8 ?& u, S3 Q& v1 \3 j; x* a
0 N3 o$ @4 Y0 v  m" qf=double(P1);
7 C% F1 Z( m2 w1 p9 b& @; j2 O1 x" H! g4 Q# C) r6 W3 n$ T
g=fft2(f);%傅里叶变换3 f& ~0 p4 L9 n  L, n. ^
2 ^8 a, {2 Q4 s0 T) O8 l9 R0 R
g=fftshift(g);%转换数据矩阵fftshift是针对频域的,将FFT的DC分量移到频谱中心
) ^5 P6 a* O) S  j! X& H即对频域的图像,
0 A8 W& U3 B, n: W% n# r" h( v' W3 R
[M,N]=size(g);: _. E) R4 m+ V8 G8 R! K

3 f; v5 Q3 Z( u: O1 U7 Z; V% wD0=input('inputnonnegative dhreshold D0=');: m4 _- c7 z7 F0 m+ ]9 b5 R

5 L3 L$ `8 x$ H7 n0 j( u! w# cn=input('input theorder of filtering n=');8 M! j$ z: h% `+ B9 z; K, c) ?! E
/ {' L  Y5 g6 X$ u$ d5 g
n1=fix(M/2);
) r6 s8 [: s# A: C1 P3 M" _) y4 Q
n2= fix(N/2);
, C8 k+ m/ J9 `+ I* s0 t, Y* s2 B- J0 H* v0 a$ @( a) P1 G) H4 l
for i=1:M6 a/ Y; z* n* O. s
8 t: ?# c9 w4 G* A
    for j=1:N
. ^) V$ c) M6 Q! P( T
# h& V. j- |1 w" A$ a: ~9 @) D        d=sqrt((i-n1)^2+(j-n2)^2);( z4 K2 L+ h. B- V% w) Z! q
2 o* b; A# m. o# O
h=1/(1+(d/D0)^(2*n));
3 ], w+ C6 v6 l! h6 k' N" h& T1 G/ x- w) o# _4 ^# g
        result(i,j)=h*g(i,j);
, z3 ^0 {/ Y6 I% Z9 D& b+ N$ _6 r. W# v+ \: ~' M7 A% F
    end
& E, a" ]+ W1 a7 Z! Y5 i$ @; Q
' ]- L2 A7 Q; \# Y. g, Y4 ]end
: T6 X2 X% F& J6 x; q1 Q! H
  K9 p& F8 o* C  _2 xresult=ifftshift(result);7 C% Y$ C9 ]4 v9 y" }* p

' k# G/ Q5 G# @" M1 P+ YX2=ifft2(result);. k& a; Z$ N" G+ {

$ J, v" A2 E! T$ n" [X3=uint8(real(X2));& G% y3 l8 \3 m% C. q- P4 g

) k; T8 Z8 ^3 gsubplot(133)+ z  O! O% v( W7 q

. ~* |- p7 u9 X3 S; m, D. Gimshow (X3);
/ q6 m- e# k0 i8 }7 A9 X5 g, O: P4 S# m3 W; G, t
title('Butterworth低通滤波所得图像');
; v$ q9 D! `3 _- L+ d  w) M' w) U
7 X( X9 f$ d" G4 M/ y
9 }/ K1 @2 v( `, b+ {/ Q
, i/ z# K+ T! n1 V

: R5 P  }$ N/ Z8 i4 U

该用户从未签到

2#
发表于 2020-7-8 17:40 | 只看该作者
基于matlab的Butterworth低通滤波程序
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-24 17:01 , Processed in 0.171875 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表