Matlab 低通Butterworth滤波仿真m文件

Demyar

Matlab 低通Butterworth滤波仿真m文件

//////////////////////////////////////////////////////////////////////

clear;

P=bodeoptions;
7 C8 b* t" B6 P4 c$ h5 RP.Grid='on';
6 y7 a& T$ M( tP.XLim={[0 800]};
( E" d" D: j, v3 O# |: _+ ]P.XLimMode={'manual'};
P.FreqUnits='Hz';
P.MagUnits='abs';
P.FreqScale = 'linear';
n=2;f0=100;fs=6400;
[b,a]=butter(n,f0/(fs/2),'low')
. n2 i1 z  l. |/ @H=tf(b,a,1/fs);
( P0 y% p/ v+ z# ?& Y. nbode(H,P);grid on;     %绘制滤波器伯德图
! a: [9 W; L+ y4 H6 s" `" ~Ts=0.1;          % 分析的时间
- `  a5 F$ Y% j6 m5 B4 O' gt=0:1/fs:Ts;w=100*pi*2;     %w=2*pi*f0;  
u=800+100*sin(3*w*t)+500*[(t>=0.05)-(t>=0.05+1/fs)];      %input signal；直流800V+100Hz纹波+冲击
$ ~" Q+ J5 N9 l7 _0 `: vfor i=1:1:2
; q' _( X' @" j& k    y(i)=0;
end
# k9 m4 ^' H9 N8 f8 ]for i=3:1:Ts*fs+1
    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);
end
%figure;plot(t,u,t,y);grid on;   % 绘制滤波效果，100Hz衰减95%
figure;
, w" x4 _4 ^! c8 |plot(t,u);
grid on;   % 绘制滤波效果，100Hz衰减95%
2 @# z2 o. ^9 b! n( {1 V0 sfigure;
plot(t,y);
' j2 l& z% ?3 u& T4 qgrid on;   % 绘制滤波效果，100Hz衰减95%
6 `6 W( H$ j6 J/ x7 y  c- v/////////////////////////////////////////////////////////////////////////////

9 S0 C. B- A8 Y& ]& f' f- \& e

fanichicl

这资源都有啊 谢谢楼主了