| 查看: 508 | 回复: 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 ] |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有180人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请帮忙用matlab求矩阵特征值
已经有8人回复
- matlab 矩阵 求特征值，求标准正交向量，求对角矩阵
已经有3人回复
2楼2013-01-24 09:31:49
3楼2018-04-23 19:56:46