24小时热门版块排行榜    

查看: 640  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sugerbee

新虫 (初入文坛)

[求助] 求含陀螺矩阵的特征值问题,求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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

猜猜

新虫 (初入文坛)

3楼2018-04-23 19:56:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

sugerbee

新虫 (初入文坛)

这个已经做出来了。
2楼2013-01-24 09:31:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:38 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见