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

Matlab 低通Butterworth滤波仿真m文件

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
Matlab 低通Butterworth滤波仿真m文件' s$ G/ C; S( Q1 t* ^( i
# F/ a- g+ \6 p9 F& [
  w; d1 l4 _2 @% c

- q4 K3 x, B* l1 q1 @$ f; h. \' g( f

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

clear;

P=bodeoptions;$ Y* h, K* s1 J
P.Grid='on';8 {/ V& ~3 y1 o! k- v( D1 a
P.XLim={[0 800]};
6 ~0 w5 m' S$ j$ i9 C- R( |" GP.XLimMode={'manual'};; \6 ?0 D2 a: D- y
P.FreqUnits='Hz';
9 _* z$ Z8 m- MP.MagUnits='abs';
; V0 `$ O  A4 I, X9 tP.FreqScale = 'linear';, ~. U* r7 J& e4 ^) r
n=2;f0=100;fs=6400;6 e+ u5 E2 v% K' w7 J7 W4 C* U7 b. r  L* t
[b,a]=butter(n,f0/(fs/2),'low')
. n& ]1 R) A6 nH=tf(b,a,1/fs);
) _; F3 C3 |0 k% X: qbode(H,P);grid on;     %绘制滤波器伯德图
) i( f+ Q8 j% q! X: TTs=0.1;          % 分析的时间
8 t$ Z& s# w6 _/ }t=0:1/fs:Ts;w=100*pi*2;     %w=2*pi*f0;  
: n  U6 c  l) J5 a' d9 j- A8 `3 au=800+100*sin(3*w*t)+500*[(t>=0.05)-(t>=0.05+1/fs)];      %input signal;直流800V+100Hz纹波+冲击+ g. w! v6 j* |; l( K
for i=1:1:2# M: S& L; j& Z) r
    y(i)=0;
  ~0 D4 t" A; Oend
; S$ T6 S) ^! \1 R  W3 v" xfor i=3:1:Ts*fs+1
+ v& ?4 E* o% @, B! ~4 e* @' d    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);
- t% O6 r6 S9 N+ V" {* ]* Xend
! Q' D3 J/ d2 g! T1 x" v, o%figure;plot(t,u,t,y);grid on;   % 绘制滤波效果,100Hz衰减95%
. Y" }* f% L6 gfigure;
2 K. H3 E# q9 n4 R5 mplot(t,u); 5 U) N; N; `5 B) S* g
grid on;   % 绘制滤波效果,100Hz衰减95%) M% n/ l+ I6 s3 ^7 V. [1 h, G2 V: A
figure;
; `' A6 x8 c0 y: R  p! Yplot(t,y);1 J" m* S. m  \+ n9 ~+ u- Q
grid on;   % 绘制滤波效果,100Hz衰减95% 9 R! f) V) Z9 y- O" u! n
/////////////////////////////////////////////////////////////////////////////
5 C) j# T  ?* n: q0 Q; \0 [, \5 f5 F: A

- j+ J8 f9 ^- R

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-7-12 19:07 , Processed in 0.125000 second(s), 24 queries , Gzip On.

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

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

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