skymyths 发表于 2007-9-19 19:36

请问如何对9×9的符号矩阵求特征值和特征向量

需要对一个9×9的符号矩阵求特征值和特征向量
符号定义及矩阵表示如下

syms k1 k2 k3 zta beta apha u v w H Y1 Y2 Y3 Y4 beta1 beta2 beta3 beta4
A=[0,k1,k2,k3,0,0,0,0,0;
-zta*u+k1*beta,k1*(1-apha)*u+zta,u*k2-v*apha*k1,u*k3-w*apha*k1,apha*k1,k1*beta1,k1*beta2,k1*beta3,k1*beta4;
-zta*v+k2*beta,v*k1-u*apha*k2,k2*(1-apha)*v+zta,v*k3-w*apha*k2,apha*k2,k2*beta1,k2*beta2,k2*beta3,k2*beta4;
-zta*w+k3*beta,w*k1-u*apha*k3,w*k2-v*apha*k3,k3*(1-apha)*w+zta,apha*k3,k3*beta1,k3*beta2,k3*beta3,k3*beta4;
(beta+H)*zta,-zta*u*apha+H*k1,-zta*v*apha+H*k2,-zta*w*apha+H*k3,(apha+1)*zta,zta*beta1,zta*beta2,zta*beta3,zta*beta4;
-Y1*zta,Y1*k1,Y1*k2,Y1*k3,0,zta,0,0,0;
-Y2*zta,Y2*k1,Y2*k2,Y2*k3,0,0,zta,0,0;
-Y3*zta,Y3*k1,Y3*k2,Y3*k3,0,0,0,zta,0;
-Y4*zta,Y4*k1,Y4*k2,Y4*k3,0,0,0,0,zta]

e=eig(A)

运行之后就没反应了,不知是什么地方出了问题,是否矩阵太大了?请高手指教,不胜感激!

[ 本帖最后由 eight 于 2007-9-19 23:04 编辑 ]

花如月 发表于 2007-9-19 20:56

唉,参数这么多,怕是matlab算不出来吧

咕噜噜 发表于 2007-9-20 13:27

如果是数值还可以,符号矩阵恐怕困难
页: [1]
查看完整版本: 请问如何对9×9的符号矩阵求特征值和特征向量