| 查看: 539 | 回复: 2 | ||
[求助]
求含陀螺矩阵的特征值问题，求matlab方法
|
|
A为正定矩阵，B为反对称矩阵。求含陀螺力矩的特征值。下面图片为公式 求指导啊。 数据如下： E=205*10^9; m1=0.06;m2=0.015;m3=0.02;m4=0.07; l1=0.025;l2=0.03;l3=l1; d1=0.01;d2=0.008;d3=0.014; I1=pi*d1^4/64;I2=pi*d2^4/64;I3=pi*d3^4/64; k11=12*E*I1/l1^3;k12=6*E*I1/l1^2;k13=2*E*I1/l1; k21=12*E*I2/l2^3;k22=6*E*I2/l2^2;k23=2*E*I2/l2; k31=12*E*I3/l3^3;k32=6*E*I3/l3^2;k33=2*E*I3/l3; kb=1*10^0; J1=9*10^-6;J2=0.2*10^-6;J3=J2;J4=8*10^-6; Jp1=16*10^-6;Jp2=0.1*10^-6;Jp3=0.1*10^-6;Jp4=14*10^-6; K=[k11,k12,-k11,k12,0,0,0,0; k12,2*k13,-k12,k13,0,0,0,0; -k11,-k12,k11+k21+kb,-k12+k22,-k21,k22,0,0; k12,k13,-k12+k22,2*k13+2*k23,-k22,k23,0,0; 0,0,-k21,-k22,k21+k31+kb,-k22+k32,-k31,k32; 0,0,k22,k23,-k22+k32,2*k23+2*k33,-k32,k33; 0,0,0,0,-k31,-k32,k31,-k32; 0,0,0,0,k32,k33,-k32,2*k33]; M=[m1,0,0,0,0,0,0,0; 0,J1,0,0,0,0,0,0; 0,0,m2,0,0,0,0,0; 0,0,0,J2,0,0,0,0; 0,0,0,0,m3,0,0,0; 0,0,0,0,0,J3,0,0; 0,0,0,0,0,0,m4,0; 0,0,0,0,0,0,0,J4]; G=zeros(8);P=zeros(8); G(2,2)=-Jp1;G(4,4)=-Jp2;G(6,6)=-Jp3;G(8,8)=-Jp4; A=[M,P;P,K];B=[w*G,K;-K,P];  陀螺力矩.jpg [ Last edited by sugerbee on 2013-1-17 at 17:01 ] |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有138人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有25人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有29人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请帮忙用matlab求矩阵特征值
已经有8人回复
- matlab 矩阵 求特征值，求标准正交向量，求对角矩阵
已经有3人回复
2楼2013-01-24 09:31:49
3楼2018-04-23 19:56:46
