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

[3.编程] 请教程序错误

[复制链接]
发表于 2011-5-6 22:21:07 | 显示全部楼层 |阅读模式 来自 江苏南京
下面的代码希望在计算出J的某个特征向量q,然而命令似乎不能计算回代具体值指之后的J矩阵,如果加上人工计算后的J的表达式(代码中原本注释掉的一行),则计算正确。请问如何避免人工计算给J重新赋值?
restart;
with(linalg);

F[1] := A-(B+1)*X[1]+X[1]^2*X[2];
F[2] := B*X[1]-X[1]^2*X[2];

J := jacobian([F[1], F[2]], [X[1], X[2]]);

sol := solve({F[1] = 0, F[2] = 0, trace(J) = 0}, {B, X[1], X[2]});

assign(sol);
assume(A > 0);
# J := matrix(2, 2, [A^2, A^2, -1-A^2, -A^2])
ev := eigenvectors(J, 'radical');
q := ev[1][3][1];
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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

GMT+8, 2024-3-29 19:53 , Processed in 0.031017 second(s), 14 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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