找回密码
 注册
Simdroid-非首页
查看: 106|回复: 2

复数矩阵求逆的问题

[复制链接]
发表于 2010-9-19 10:03:14 | 显示全部楼层 |阅读模式 来自 北京
a=1.65;
b=1.1;
L1=0.275;
m1=140;
m2=210;
m3=2150;
m4=3902;
m5=53;
c3=150;
c4=150;
c5=170;
k1=571420;
k2=571420;
k3=74400;
k4=121000;
k5=13000;
M=[m1,0,0,0,0;0,m2,0,0,0;0,0,m3,0,0;0,0,0,m4,0;0,0,0,0,m5]
C=[c3,0,-c3,a*c3,0;0,c4,-c4,-b*c4,0;-c3,-c4,c3+c4+c5,-a*c3-L1*c5+b*c4,-c5;...
        a*c3,-b*c4,-a*c3-L1*c5+b*c4,((a^2)*c3+(b^2)*c4+(L1^2)*c5),L1*c5;0,0,-c5,L1*c5,c5]
K=[k1+k3,0,-k3,a*k3,0;0,(k2+k4),-k4,-b*k4,0;-k3,-k4,(k3+k4+k5),(-a*k3+b*k4-L1*k5),-k5;...
        a*k3,-b*k4,(-a*k3+b*k4-L1*k5),((a^2)*k3+(L1^2)*k5+(b^2)*k4),L1*k5;0,0,-k5,L1*k5,k5]
Kt=[k1,0,0,0,0;0,k2,0,0,0]'
w=sym('w');
G=(-w^2)*M+w*i*C+K
G2=inv(G)
其中w为符号变量,为什么G矩阵求出来是正确的,但它的逆矩阵G2求出来数量级那么大,直接影响后面的积分工作?
发表于 2010-9-19 11:12:36 | 显示全部楼层 来自 河北廊坊
Simdroid开发平台
1# 水与火
对于用matlab做这种复杂的符号计算,我一般是不推荐的,建议如果可以的话,你将算式列出来,然后用数值计算去求解

评分

1

查看全部评分

回复 不支持

使用道具 举报

 楼主| 发表于 2010-9-19 18:31:08 | 显示全部楼层 来自 北京
我的目的是求附件中的H(w)矩阵。但一求G矩阵的逆就会产生很庞大的结果。还请高手帮帮忙。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-10-5 11:22 , Processed in 0.052671 second(s), 17 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表