| 查看: 556 | 回复: 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 ] |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有101人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请帮忙用matlab求矩阵特征值
已经有8人回复
- matlab 矩阵 求特征值，求标准正交向量，求对角矩阵
已经有3人回复
2楼2013-01-24 09:31:49
3楼2018-04-23 19:56:46
