机械社区

 找回密码
 注册会员

扫一扫,访问微社区

QQ登录

只需一步,快速开始

搜索
查看: 330|回复: 0

[matlab] 无阻尼二自由系统固有频率的程序

[复制链接]
发表于 2018-4-16 19:41:12 | 显示全部楼层 |阅读模式
% 无阻尼二自由系统固有频率的程序 20180401A
. \1 `/ {5 H' N: gm1=100;m2=100;k1=200;k2=400;% : Y/ P* r; Z# o* S  h& U5 Y
c1=0;c2=0;2 M) o  W4 v0 U$ t+ f4 S
x0=[1;0];xd0=[0;0];tf=50;dt=0.1;  J& U8 u- ~1 v$ m+ N) g
M=[m1,0;0,m2];k=[k1+k2,-k2;-k2,k2];
% I' j7 `0 D: q3 i' qc=[c1+c2,-c2;-c2,c2];
4 p: s" l9 [! p0 q/ OA=[zeros(2,2),eye(2);-M\k,-M\c];% U! y9 U$ {4 U: U  N. n
y0=[x0;xd0];! o' J- \' ]5 u& t/ N$ {9 G
for i=1:round(tf/dt)+1
1 O' C+ A. p5 C, D: T) G    t(i)=dt*(i-1)1 D. L1 n! _$ T# Z; h! j
    y(:,i)=expm(A*t(i))*y0;5 X4 h+ y( H' J3 P9 y- v
end
+ o$ J2 ^) z0 [4 ?: H2 `0 isubplot(2,1,1),plot(t,y(1,),grid  
3 s% F: W( b& i: Jsubplot(2,1,2),plot(t,y(2,),grid
' ]. P) [+ u; d[p,lamda]=eig(A)
+ }; q& C) _3 Y4 h% Q! P3 [2 @% @6 l

8 ~$ p1 }$ W5 e& Z+ b& y3 E- m1 R% u7 A1 L

" q7 }- `9 ?/ n# {lamda =& u, T# z/ ^4 p' C# h6 q
  Columns 1 through 2/ F4 @) L; ?4 F
  0.000000000000000 + 3.020447918044219i  0.000000000000000 + 0.000000000000000i
, W5 l5 W# Z8 l  0.000000000000000 + 0.000000000000000i  0.000000000000000 - 3.020447918044219i
$ M- }0 P, ]3 r! ?. X# v% w  0.000000000000000 + 0.000000000000000i  0.000000000000000 + 0.000000000000000i- q1 G5 `' B! D% \# o6 k- Z/ @% N
  0.000000000000000 + 0.000000000000000i  0.000000000000000 + 0.000000000000000i
. I; V# k1 \, M# C9 U  Columns 3 through 4
; Q& i: N% X8 B- m0 J  I* |8 a  0.000000000000000 + 0.000000000000000i  0.000000000000000 + 0.000000000000000i
" I) k/ S5 ^7 I& p  z% V. i' G, Y9 I  0.000000000000000 + 0.000000000000000i  0.000000000000000 + 0.000000000000000i
/ c' ]. k. n# T" v. A$ l% N+ ~ -0.000000000000000 + 0.936426384924271i  0.000000000000000 + 0.000000000000000i7 R6 @. u3 i& A3 J% h4 h, C
  0.000000000000000 + 0.000000000000000i -0.000000000000000 - 0.936426384924271i
* p! U$ k; ^! A
# K% S# X, V4 N% g1 g, r3 ^  P4 z
( k! [- j6 k, A( T: h% N
; }+ d7 T) r* b, T/ [" U5 o" Y% ^' h+ Q6 I

本帖子中包含更多资源

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

x
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册会员

本版积分规则

小黑屋|手机版|Archiver|中国机械社区 ( 京ICP备10217105号,京ICP证050210号,京公网安备11010802010176 )  

GMT+8, 2018-4-24 12:38 , Processed in 0.079904 second(s), 7 queries , Gzip On, MemCache On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

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