找回密码
 注册
关于网站域名变更的通知
查看: 393|回复: 1
打印 上一主题 下一主题

Matlab 低通Butterworth滤波仿真m文件

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2018-12-29 07:00 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
Matlab 低通Butterworth滤波仿真m文件: l: a$ B1 p9 Q2 b+ z
; z% b" a! Q% X- w* o
+ n) L6 ^- n/ N6 D, C/ }6 E
8 i' c4 V+ m# d

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

clear;

P=bodeoptions;# @9 G" g' }8 x
P.Grid='on';
7 ^, C& Y0 y7 Z9 ^: S, ]P.XLim={[0 800]};% @) s: |8 `: o6 r1 E1 g
P.XLimMode={'manual'};: R, w& h4 |$ w/ x" p2 W, J8 D* ^
P.FreqUnits='Hz';7 U8 F* H& y. S
P.MagUnits='abs';
& p0 h! `3 ?$ a, H2 AP.FreqScale = 'linear';! p' f) J/ b6 M' o$ e% E; C- W
n=2;f0=100;fs=6400;5 d1 J" Z0 Z2 F* a
[b,a]=butter(n,f0/(fs/2),'low')1 V( N  n) Z6 L. M
H=tf(b,a,1/fs);
2 e: C. ?8 A# j2 C# e8 j! ^0 j0 rbode(H,P);grid on;     %绘制滤波器伯德图" v2 {$ S( w. v% h' X# \
Ts=0.1;          % 分析的时间
+ h. @& {. l; qt=0:1/fs:Ts;w=100*pi*2;     %w=2*pi*f0;  2 C5 e% J. g+ E, W6 h
u=800+100*sin(3*w*t)+500*[(t>=0.05)-(t>=0.05+1/fs)];      %input signal;直流800V+100Hz纹波+冲击% Q- x5 g, ~* x
for i=1:1:2
1 S- l7 P9 k3 |5 Z' A    y(i)=0;- O6 N# J) z' Z( t! s( H; N, u8 e+ T
end
4 t% M0 T" t. h# t% nfor i=3:1:Ts*fs+15 {1 _2 R% Y9 P) U: q
    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);
$ f- K7 e2 o6 S9 [$ f( _end
" }. N$ E0 }0 H%figure;plot(t,u,t,y);grid on;   % 绘制滤波效果,100Hz衰减95%
1 L3 Q$ J+ m& C: ofigure;
0 W/ p9 ^0 e- ?7 n  B9 [plot(t,u);
) D5 S2 k/ z! l& Ggrid on;   % 绘制滤波效果,100Hz衰减95%4 a) |# C  ~0 s: i" Z' @1 T
figure;
9 Y/ k8 u8 |# m. L( pplot(t,y);
& ~: ]6 s6 V& k/ G% j9 Jgrid on;   % 绘制滤波效果,100Hz衰减95%
) s5 N' l- a3 X2 ]- R6 W/////////////////////////////////////////////////////////////////////////////
% |: y5 k9 O* G+ N5 }  @
/ L  h- s0 A8 L1 n4 Z+ _0 Y" R/ a" _9 D6 m/ i

该用户从未签到

2#
发表于 2018-12-29 11:10 | 只看该作者
这资源都有啊 谢谢楼主了
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-7-11 05:56 , Processed in 0.109375 second(s), 23 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表