复数矩阵求逆的问题
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=
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='
w=sym('w');
G=(-w^2)*M+w*i*C+K
G2=inv(G)
其中w为符号变量，为什么G矩阵求出来是正确的，但它的逆矩阵G2求出来数量级那么大，直接影响后面的积分工作？ 1# 水与火
对于用matlab做这种复杂的符号计算，我一般是不推荐的，建议如果可以的话，你将算式列出来，然后用数值计算去求解 我的目的是求附件中的H（w）矩阵。但一求G矩阵的逆就会产生很庞大的结果。还请高手帮帮忙。
页:
[1]