EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
Matlab 低通Butterworth滤波仿真m文件: l: a$ B1 p9 Q2 b+ z
; z% b" a! Q% X- w* o
+ n) L6 ^- n/ N6 D, C/ }6 E
8 i' c4 V+ m# d
////////////////////////////////////////////////////////////////////// clear; P=bodeoptions;# @9 G" g' }8 x
P.Grid='on';
7 ^, C& Y0 y7 Z9 ^: S, ]P.XLim={[0 800]};% @) s: |8 `: o6 r1 E1 g
P.XLimMode={'manual'};: R, w& h4 |$ w/ x" p2 W, J8 D* ^
P.FreqUnits='Hz';7 U8 F* H& y. S
P.MagUnits='abs';
& p0 h! `3 ?$ a, H2 AP.FreqScale = 'linear';! p' f) J/ b6 M' o$ e% E; C- W
n=2;f0=100;fs=6400;5 d1 J" Z0 Z2 F* a
[b,a]=butter(n,f0/(fs/2),'low')1 V( N n) Z6 L. M
H=tf(b,a,1/fs);
2 e: C. ?8 A# j2 C# e8 j! ^0 j0 rbode(H,P);grid on; %绘制滤波器伯德图" v2 {$ S( w. v% h' X# \
Ts=0.1; % 分析的时间
+ h. @& {. l; qt=0:1/fs:Ts;w=100*pi*2; %w=2*pi*f0; 2 C5 e% J. g+ E, W6 h
u=800+100*sin(3*w*t)+500*[(t>=0.05)-(t>=0.05+1/fs)]; %input signal;直流800V+100Hz纹波+冲击% Q- x5 g, ~* x
for i=1:1:2
1 S- l7 P9 k3 |5 Z' A y(i)=0;- O6 N# J) z' Z( t! s( H; N, u8 e+ T
end
4 t% M0 T" t. h# t% nfor i=3:1:Ts*fs+15 {1 _2 R% Y9 P) U: q
y(i)=-1*(a(2)*y(i-1)+a(3)*y(i-2))+b(1)*u(i)+b(2)*u(i-1)+b(3)*u(i-2);
$ f- K7 e2 o6 S9 [$ f( _end
" }. N$ E0 }0 H%figure;plot(t,u,t,y);grid on; % 绘制滤波效果,100Hz衰减95%
1 L3 Q$ J+ m& C: ofigure;
0 W/ p9 ^0 e- ?7 n B9 [plot(t,u);
) D5 S2 k/ z! l& Ggrid on; % 绘制滤波效果,100Hz衰减95%4 a) |# C ~0 s: i" Z' @1 T
figure;
9 Y/ k8 u8 |# m. L( pplot(t,y);
& ~: ]6 s6 V& k/ G% j9 Jgrid on; % 绘制滤波效果,100Hz衰减95%
) s5 N' l- a3 X2 ]- R6 W/////////////////////////////////////////////////////////////////////////////
% |: y5 k9 O* G+ N5 } @
/ L h- s0 A8 L1 n4 Z+ _0 Y" R/ a" _9 D6 m/ i
|