|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
- syms m1 m2 c k1 k2 Y w den num real
- %这个real很重要的 如果在maple里面应该是 assume指令相当的
- % sym(m1,'real')
- % sym(m2,'real')
- % sym(c,'real')
- % sym(k1,'real')
- % sym(k2,'real')
- % sym(Y,'real')
- % sym(w,'real')
- % sym(den,'real')
- % sym(num,'real')
- % sym(,'real')
- % sym(,'real')
- M =[ m1 0; 0 m2]
- C =[ c -c ; -c c]
- K = [ k1,-k1;-k1,k1+k2]
- U0 = -w^2*M + i*w*C +K
- F=[ 0;k2*Y]
- U = inv(U0)*F
- U = inv(U0)*F
- U1=U(1,:)
- [num,den]=numden(det(U1))
- % collect(den,i)
- % U2=U(2,:)
- Imden=imag(-den)
- Imnum=imag(-num)
- imden = expand(Imden)
- imnum = expand(Imnum)
- realden = expand(real(-den))
复制代码
本题是 两自由度系统的一个例子的求解 |
|