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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
clear all;* i2 ]( B( V. p2 v: s9 q
0 r5 J: L6 S  h) K; j
close all;
- v9 z9 ^+ B  h; e1 f9 A4 ?
7 e" o0 k$ r5 P/ E# @$ E1 OP=input('inputimage data:');
/ }& O7 l2 j+ S& ?9 y8 z2 h
5 d  a2 U# k- ^d=size(P);
& n4 t" x$ ^2 f/ ~) g  u% X, P1 d- ~+ o/ e" }+ x( K' [9 x1 L+ g5 d
if(d(3)>1)
6 W- ]3 D1 b5 E* g) ?3 h6 l0 p  V  u/ K2 A5 O
    P=rgb2gray(P);$ }0 p5 a1 A! w4 H
6 o2 C0 t/ Q3 L/ ]0 N
end
0 D1 |$ h4 |' `0 ^; ^. ^3 q- _
, H9 d+ m  \3 l  |$ [" l5 u, xsubplot(131);imshow(P);
% B2 l( i% M! v1 K4 J- r/ u
" v) V8 ]' e9 l9 I& ztitle('原始图像');
- k, V6 t' m- r  h0 u5 F& i( B6 ]4 p; M  E% [; }
P1=imnoise(P,'salt& pepper',0.02);
5 w5 J! v% B3 d# D! G( j; J* s
7 r- l" {9 u5 _3 r6 Y/ ^subplot(132)
7 @0 M/ ^7 w7 _1 J* O8 V4 Z( Y2 s5 ]2 I- o/ ^% U( g8 X7 `
imshow (P1);* M! i; e$ C+ `- K) b
) `9 R2 J4 a( w' r0 {8 B
title('加入椒盐噪声后的图像');% I9 y. F8 k. B) s3 J$ ~

$ p5 h, _2 A* B8 `5 M. l: ff=double(P1);
/ L6 f# _) d* S, }/ c
, d0 [- e4 ~2 c# [0 U7 zg=fft2(f);%傅里叶变换" ^. I3 o3 z2 [, S1 R& g
5 D5 @9 X/ J1 g9 X% n' h; G
g=fftshift(g);%转换数据矩阵fftshift是针对频域的,将FFT的DC分量移到频谱中心, ~( c- S0 o0 Q, |
即对频域的图像,
, E0 H1 H! \& }" h' r5 a) o6 e/ U; c& Z( P/ _
[M,N]=size(g);9 Z2 L: `, _3 A3 O3 S7 L
; o% D, F0 T8 x8 {, M- S% I5 s1 l" u
D0=input('inputnonnegative dhreshold D0=');" s, j4 |3 z& k* C, A! U1 G5 @9 |1 ]+ N

+ B# |0 N, M; k+ ?n=input('input theorder of filtering n=');0 Y: V0 H! Y& e
9 M# m! z& ]# N' J
n1=fix(M/2);
4 O% H4 C9 J8 T: I( ~
0 t8 T5 \6 v, ?4 J* u2 z) Gn2= fix(N/2);( D% r% G. [" T+ R2 R; X6 Y0 c
% }& g; W9 i2 D+ D# P
for i=1:M* n' p. J/ B  f
+ {1 x! P5 `0 C
    for j=1:N  D3 u9 q6 {" C2 l- v' S/ j7 }

& A$ j2 Q# a$ @$ E' A        d=sqrt((i-n1)^2+(j-n2)^2);
& R& v% S: h/ @$ k. o, K. X2 `1 G* e3 f; j; _
h=1/(1+(d/D0)^(2*n));7 s; _- k5 T' s. C+ k; v) ]+ \

0 p9 f1 M5 L" v) ?        result(i,j)=h*g(i,j);- q' L* o' o! r- W

! \4 {2 I" ?; E) a0 K    end
/ D  ], x; X+ |5 L8 i# g& Q0 o4 p9 T8 b! _
end8 a- A* S) {% `) l' b0 j6 w: A
3 b) U" t& M! z% u
result=ifftshift(result);+ O$ }3 Z- l' H* `

+ o: u3 u7 h8 R8 l2 Z1 ^( C# HX2=ifft2(result);
: m: A! I% E$ w. }  k8 \' X
( ~  J0 b3 g5 K  v/ p$ ~% ~+ TX3=uint8(real(X2));
# y; O" {- }! A
& ~/ g6 u5 I2 A) ssubplot(133)
' y6 X. v; e8 m3 L
  i/ {" `/ r1 ^7 U% G3 A+ F4 c9 wimshow (X3);
7 c% ]( Q8 i8 n/ [1 h3 }. \1 k
4 _. L1 V- }$ J4 f1 B6 q& Otitle('Butterworth低通滤波所得图像');3 R9 \! _+ f4 Q8 v+ Q1 ~4 j

! D3 W6 `$ f# b* H3 P' H5 c  ?* l, `* R7 @
: D0 c/ |( S" H3 s1 S4 X

2 p' v1 B6 F* J& K# o0 k' r1 N$ @4 T- x8 R; v; ~

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

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

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

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

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