|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
需要对一个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 编辑 ] |
|