机械社区

 找回密码
 注册会员

扫一扫,访问微社区

QQ登录

只需一步,快速开始

搜索
查看: 777|回复: 0

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

[复制链接]
发表于 2018-4-16 19:41:12 | 显示全部楼层 |阅读模式
% 无阻尼二自由系统固有频率的程序 20180401A
; R; w) U* n; t' O2 e# {m1=100;m2=100;k1=200;k2=400;% 5 w; ~* q* N2 i2 H) O
c1=0;c2=0;6 Y+ H- U/ y5 j. D% _( b: s0 x
x0=[1;0];xd0=[0;0];tf=50;dt=0.1;# F6 z, s; l" ^9 R
M=[m1,0;0,m2];k=[k1+k2,-k2;-k2,k2];+ S9 |6 H1 @, i* N+ V7 _
c=[c1+c2,-c2;-c2,c2];% R: ?3 j  I1 f& o
A=[zeros(2,2),eye(2);-M\k,-M\c];
- d+ ?( t* k2 _0 d7 q  iy0=[x0;xd0];4 ]" c6 S) K) Y+ w8 ^) F% P
for i=1:round(tf/dt)+1* f0 w. }! m5 E4 J; D2 a
    t(i)=dt*(i-1)
$ e/ c" N4 |# U( I    y(:,i)=expm(A*t(i))*y0;) T; G$ ~3 F" t5 ]( j8 Z
end" O5 A6 L3 W: E& K
subplot(2,1,1),plot(t,y(1,),grid  6 [) ^2 I7 C! L: s( o5 v7 O
subplot(2,1,2),plot(t,y(2,),grid% G1 r7 }5 E" w7 n$ b7 g+ ~. h
[p,lamda]=eig(A)
# ?9 y- ^5 u9 u. a3 D* P4 B& i5 c! m8 M& K( B. Y- x5 ]
. c" X7 @3 |/ _- R; J$ s

9 y. Z' H- ~) e) E0 o$ P; T7 Z( q" I
lamda =- b/ M# U5 i" z; G
  Columns 1 through 2
, f7 D8 s9 C- C, y/ M( j' q2 i# i% v  0.000000000000000 + 3.020447918044219i  0.000000000000000 + 0.000000000000000i
' ~1 A  b4 ~& p' _0 i- }* k  0.000000000000000 + 0.000000000000000i  0.000000000000000 - 3.020447918044219i! p8 z1 d0 e7 N
  0.000000000000000 + 0.000000000000000i  0.000000000000000 + 0.000000000000000i. o7 x% m8 S* d: l+ `6 z
  0.000000000000000 + 0.000000000000000i  0.000000000000000 + 0.000000000000000i
% k* `! \, e. b: b  Columns 3 through 4
" G0 H$ u1 s7 t( s6 i7 |3 j' a  0.000000000000000 + 0.000000000000000i  0.000000000000000 + 0.000000000000000i! f- a$ b( g* J' M; i; y  v
  0.000000000000000 + 0.000000000000000i  0.000000000000000 + 0.000000000000000i
+ ]( A) S6 n; i7 Z, r" H* f -0.000000000000000 + 0.936426384924271i  0.000000000000000 + 0.000000000000000i( E# D% \! Q! p' }* ~  c4 `# P
  0.000000000000000 + 0.000000000000000i -0.000000000000000 - 0.936426384924271i) \9 ~' b5 }+ ^1 {! a) m2 G" T

6 j: U; B8 c3 G& `0 X+ q  P8 J- W9 K8 I, [) [' X' B+ Z5 R

& T% d1 i* q, `2 ~/ I; E! l) `% y2 s

本帖子中包含更多资源

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

x
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2018-10-17 14:35 , Processed in 0.059884 second(s), 9 queries , Gzip On, MemCache On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

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