24小时热门版块排行榜    

查看: 695  |  回复: 3

cooooldog

铁杆木虫 (著名写手)

ส็็็

[交流] 一个非线性方程组的求解 已有1人参与

这里有个非线性方程组, 具有如下特征:

1)雅克比矩阵的谱半径大于1,牛顿迭代?
2)有无穷多复数解
用mathematica的NSolve,不同版本给出不同结果;

已经发现了问题的原因:
1) 牛顿迭代收敛跟雅克比矩阵谱半径无关;
2) mathematica 8, 9 的NSolve有bug
谢谢
一个非线性方程组的求解
nonlinear system.png



[ Last edited by cooooldog on 2013-6-17 at 11:19 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

ส็็็็็็็็็็็็็็็็็็็็
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
方便他人,给出文本格式。
2楼2013-06-15 17:41:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

文本格式是这样的:

((300 (1920+8 x-21 y) (-80+y))/(7 (30+x)^3)+(2025 (208 x+5 (-3446+y)))/(52 (90+y)^2)+(300 (4500+80 x-21 z) (-100+z))/(7 (30+x)^3)-(1521 (15425-539 x+50 z))/(539 (39+z)^2)
-((300 (1920+8 x-21 y))/(7 (30+x)^2))-(2025 (-85+x) (208 x+5 (-3446+y)))/(52 (90+y)^3)+(2025 (49410+445 y-208 z) (-45+z))/(52 (90+y)^3)-(1521 (19680-539 y+90 z))/(539 (39+z)^2)
-((2025 (49410+445 y-208 z))/(52 (90+y)^2))-(300 (4500+80 x-21 z))/(7 (30+x)^2)-(1521 (-25+x) (539 x-25 (617+2 z)))/(539 (39+z)^3)-(1521 (-30+y) (539 y-30 (656+3 z)))/(539 (39+z)^3)
)
ส็็็็็็็็็็็็็็็็็็็็
3楼2014-08-04 17:57:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

这样更好一些:

{{(300 (1920+8 x-21 y) (-80+y))/(7 (30+x)^3)+(2025 (208 x+5 (-3446+y)))/(52 (90+y)^2)+(300 (4500+80 x-21 z) (-100+z))/(7 (30+x)^3)-(1521 (15425-539 x+50 z))/(539 (39+z)^2)},{-((300 (1920+8 x-21 y))/(7 (30+x)^2))-(2025 (-85+x) (208 x+5 (-3446+y)))/(52 (90+y)^3)+(2025 (49410+445 y-208 z) (-45+z))/(52 (90+y)^3)-(1521 (19680-539 y+90 z))/(539 (39+z)^2)},{-((2025 (49410+445 y-208 z))/(52 (90+y)^2))-(300 (4500+80 x-21 z))/(7 (30+x)^2)-(1521 (-25+x) (539 x-25 (617+2 z)))/(539 (39+z)^3)-(1521 (-30+y) (539 y-30 (656+3 z)))/(539 (39+z)^3)}}
ส็็็็็็็็็็็็็็็็็็็็
4楼2014-08-04 17:59:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cooooldog 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见