Matlab 低通Butterworth滤波仿真m文件

Demyar

Matlab 低通Butterworth滤波仿真m文件
$ C; g$ N) g) L, p$ w2 ]0 N
9 }, K8 |$ u1 w1 k' f7 ]0 r& M$ T
$ U. V# U3 E7 g0 l( n; U. ~

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

clear;

P=bodeoptions;
, c( [, U: I* R2 v+ ~: ?+ I* iP.Grid='on';
P.XLim={[0 800]};
P.XLimMode={'manual'};
P.FreqUnits='Hz';
7 ^2 x8 m8 i1 ]8 G, ~9 AP.MagUnits='abs';
P.FreqScale = 'linear';
5 u- o4 [; u1 J8 U7 hn=2;f0=100;fs=6400;
[b,a]=butter(n,f0/(fs/2),'low')
H=tf(b,a,1/fs);
# F$ f4 _3 R+ V* M0 v5 v0 z# hbode(H,P);grid on;     %绘制滤波器伯德图
Ts=0.1;          % 分析的时间
/ U0 D2 I; y3 V( t! @; t# b; Vt=0:1/fs:Ts;w=100*pi*2;     %w=2*pi*f0;  
( F; A; K! x& d5 t% ^u=800+100*sin(3*w*t)+500*[(t>=0.05)-(t>=0.05+1/fs)];      %input signal；直流800V+100Hz纹波+冲击
for i=1:1:2
6 B0 X0 ~: r* d3 L) L- o: D    y(i)=0;
end
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
' {4 @# R# G! c( M& y%figure;plot(t,u,t,y);grid on;   % 绘制滤波效果，100Hz衰减95%
6 x6 T5 S2 t, y$ o* w0 gfigure;
plot(t,u);
% C8 C: P& P( ~. h9 M1 fgrid on;   % 绘制滤波效果，100Hz衰减95%
) `) X8 p5 V2 C' G4 A% t$ ^figure;
0 K  g4 w. R1 ]) rplot(t,y);
grid on;   % 绘制滤波效果，100Hz衰减95%
' g0 a7 u+ {6 w6 I/////////////////////////////////////////////////////////////////////////////

' }( W4 Y& S# b  Q8 g

fanichicl

这资源都有啊 谢谢楼主了