| 查看: 641 | 回复: 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 ] |
回复此楼» 猜你喜欢
青基已中
已经有3人回复
-
散金币祈福
已经有186人回复
-
物理学I论文润色/翻译怎么收费?
已经有213人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有110人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有66人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请帮忙用matlab求矩阵特征值
已经有8人回复
- matlab 矩阵 求特征值，求标准正交向量，求对角矩阵
已经有3人回复
2楼2013-01-24 09:31:49
3楼2018-04-23 19:56:46
10