Matlab 低通Butterworth滤波仿真m文件

Demyar

Matlab 低通Butterworth滤波仿真m文件

1 p1 S& B8 ]) ~' \* h
  w5 C, v0 G" C  J& [! {
* m: O  Z+ d# \: }6 G

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

clear;

P=bodeoptions;
7 u/ i& ~* h+ wP.Grid='on';
P.XLim={[0 800]};
P.XLimMode={'manual'};
P.FreqUnits='Hz';
P.MagUnits='abs';
P.FreqScale = 'linear';
n=2;f0=100;fs=6400;
" B% U" R3 n- Z" L5 U[b,a]=butter(n,f0/(fs/2),'low')
2 T3 C* j6 W+ G) QH=tf(b,a,1/fs);
bode(H,P);grid on;     %绘制滤波器伯德图
Ts=0.1;          % 分析的时间
5 L( L+ K5 @$ r/ d1 b% Zt=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纹波+冲击
* \; ]" X, s2 e( Ofor i=1:1:2
8 @: F" X) ^, U6 u; N    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);
9 b; [; V/ x" A8 Y; Mend
5 e- [2 K2 g% O) S( W  G2 w% O+ u%figure;plot(t,u,t,y);grid on;   % 绘制滤波效果，100Hz衰减95%
figure;
plot(t,u);
1 }" A, @5 i: f" h7 f" L+ Ngrid on;   % 绘制滤波效果，100Hz衰减95%
8 i6 p) S3 U4 afigure;
plot(t,y);
grid on;   % 绘制滤波效果，100Hz衰减95%
/////////////////////////////////////////////////////////////////////////////

( {# q# l# s4 T* c% u! l
/ F" P3 ]! ^# G' ]! u

fanichicl

这资源都有啊 谢谢楼主了