EDA365电子论坛网

标题: Matlab 低通Butterworth滤波仿真m文件 [打印本页]

---

作者: Demyar    时间: 2018-12-29 07:00
标题: Matlab 低通Butterworth滤波仿真m文件
Matlab 低通Butterworth滤波仿真m文件
* X1 S% y0 q! _( G6 H
+ X: h2 T; l% j/ p6 f4 h" `

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

clear;

P=bodeoptions;
1 f- J+ m& X. e+ ?: K0 k" w# ?P.Grid='on';
P.XLim={[0 800]};
P.XLimMode={'manual'};
P.FreqUnits='Hz';
P.MagUnits='abs';
P.FreqScale = 'linear';
n=2;f0=100;fs=6400;
+ c9 q3 I; _0 w  c: d[b,a]=butter(n,f0/(fs/2),'low')
! H0 r5 Y5 k3 g9 m; N6 @$ b+ eH=tf(b,a,1/fs);
bode(H,P);grid on;     %绘制滤波器伯德图
# `6 ?6 y2 T3 XTs=0.1;          % 分析的时间
9 k: C2 t1 R$ `6 |6 R9 Ct=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纹波+冲击
for i=1:1:2
    y(i)=0;
$ S5 j6 M3 G: b" aend
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;
& F: n8 A# @7 h. iplot(t,u);
grid on;   % 绘制滤波效果，100Hz衰减95%
  K- v( ?" d% o0 A: ~0 ^, l2 ^figure;
! P: d$ i. Z; g, w1 C; y9 ^- y# N3 l& eplot(t,y);
grid on;   % 绘制滤波效果，100Hz衰减95%
/////////////////////////////////////////////////////////////////////////////

: t% N7 `& x3 h7 Y$ O/ ]' c

---

作者: fanichicl    时间: 2018-12-29 11:10
这资源都有啊 谢谢楼主了

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2