Matlab 低通Butterworth滤波仿真m文件

Demyar

Matlab 低通Butterworth滤波仿真m文件

* S$ Q$ B" A" z  F! w9 g# L
! C! b- @0 U3 f) O  L
* ]8 ~- u% ]2 @4 ^$ }* @+ f8 {

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

clear;

P=bodeoptions;
. u- A& e4 d6 eP.Grid='on';
P.XLim={[0 800]};
P.XLimMode={'manual'};
P.FreqUnits='Hz';
4 I8 ~, y8 ~% |) oP.MagUnits='abs';
( f) b8 A4 V: W, P2 N" CP.FreqScale = 'linear';
n=2;f0=100;fs=6400;
[b,a]=butter(n,f0/(fs/2),'low')
: u& X4 c) Z$ M2 b/ rH=tf(b,a,1/fs);
4 K. A/ O# h) H( R+ b+ V  e/ k. Y8 j# pbode(H,P);grid on;     %绘制滤波器伯德图
Ts=0.1;          % 分析的时间
t=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纹波+冲击
; C& |6 z0 v( [, o7 |7 Ifor i=1:1:2
8 {: T0 P8 E5 m" f7 O+ |    y(i)=0;
: O* m5 W4 O+ L' u+ Hend
for i=3:1:Ts*fs+1
# `* G" o; T4 f' O) a* j' I    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);
4 }- y* U; Q5 A# pend
4 _9 q) K4 }. o8 y%figure;plot(t,u,t,y);grid on;   % 绘制滤波效果，100Hz衰减95%
( i) Z& y" O0 L/ x5 d; Y- ifigure;
plot(t,u);
/ O* J% P$ H! p& s4 t! Jgrid on;   % 绘制滤波效果，100Hz衰减95%
  b5 V  @- @7 s# \& v* f& \figure;
plot(t,y);
grid on;   % 绘制滤波效果，100Hz衰减95%
/////////////////////////////////////////////////////////////////////////////

fanichicl

这资源都有啊 谢谢楼主了