24小时热门版块排行榜

返回列表

hitcs

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 45.9
帖子: 47
在线: 18.9小时
虫号: 3591713
注册: 2014-12-12
性别: GG
专业: 金属材料的合金相、相变及

[求助] 求助，用1stopt 5.0解非线性方程组，两次得到的解不唯一？什么原因，大神帮忙已有1人参与

代码如下，9元非线性方程组，得到的结果不唯一：
Parameter x1=0.095843301[0,0.095843301],x2=0.034539474[0,0.034539474],x3=0.004934211[0,0.004934211],x4=0.044258373[0,0.044258373],x5=0.008672249[0,0.008672249],x6=0.046202153[0,0.046202153],x7=0.015550239[0,0.015550239],x8=0.5[0,0.5],x9=0.25[0,0.25];
Function
0=0.044258373*x1+0.044258373*((5.305*10^3/1.518)^0.5)*x1/x2^2*(1.518*x2^4*x8^2*x9/(2+2*((5.305*10^3/1.518)^0.5)*x1/x2^2)^4)^(1/6)-0.095843301*x4-2*0.095843301*(1.109*10^54*x4^4*x9^3)^(1/2);
0=0.044258373*x2+2*0.044258373*(1.518*x2^4*x8^2*x9/(2+2*((5.305*10^3/1.518)^0.5)*x1/x2^2)^4)^(1/6)-0.034539474*x4-2*0.034539474*(1.109*10^54*x4^4*x9^3)^(1/2);
0=0.044258373*x3+0.044258373*(2.012*x3)^(1/2)-0.004934211*x4-2*0.004934211*(1.109*10^54*x4^4*x9^3)^(1/2);
0=0.044258373*x5+2*0.044258373*(3.131*10^32*x5^4*x9^3)^(1/2)-0.008672249*x4-2*0.008672249*(1.109*10^54*x4^4*x9^3)^(1/2);
0=0.044258373*x6+0.044258373*(1.077*10^12*x6^2*x9)^(1/2)-0.046202153*x4-2*0.046202153*(1.109*10^54*x4^4*x9^3)^(1/2);
0=0.044258373*x7+0.044258373*(2.108*10^16*x7^2*x9^3)^(1/2)-0.015550239*x4-2*0.015550239*(1.109*10^54*x4^4*x9^3)^(1/2);
0=0.044258373*x8+0.5*0.044258373*((1.518*x2^4*x8^2*x9)^(1/6))*((2+2*(5.305*10^3/1.518)^0.5*x1/x2^2)^(1/3))-0.5*x4-2*0.5*(1.109*10^54*x4^4*x9^3)^(1/2);
0=0.044258373*x9+1.5*0.044258373*(1.109*10^54*x4^4*x9^3)^(1/2)+1.5*0.044258373*(3.131*10^32*x5^4*x9^3)^(1/2)+0.5*0.044258373*(1.077*10^12*x6^2*x9)^(1/2)+1.5*0.044258373*(2.108*10^16*x7^2*x9^3)^(1/2)+0.5*0.044258373*((5.305*10^3/1.518)^0.5)*x1/x2^2*(1.518*x2^4*x8^2*x9/(2+2*((5.305*10^3/1.518)^0.5)*x1/x2^2)^4)^(1/6)+0.5*0.044258373*(1.518*x2^4*x8^2*x9/(2+2*(5.305*10^3/1.518)^0.5*x1/x2^2)^4)^(1/6)-0.25*x4-2*0.25*(1.109*10^54*x4^4*x9^3)^(1/2);
0=x1+x2+x3+x4+x5+x6+x7+x8+x9+(1.109*10^54*x4^4*x9^3)^(1/2)+(3.131*10^32*x5^4*x9^3)^(1/2)+(1.077*10^12*x6^2*x9)^(1/2)+(2.108*10^16*x7^2*x9^3)^(1/2)+(5.305*10^3/1.518)^0.5*x1/x2^2*(1.518*x2^4*x8^2*x9/(2+2*(5.305*10^3/1.518)^0.5*x1/x2^2)^4)^(1/6)+(1.518*x2^4*x8^2*x9)^(1/6)*(2+2*(5.305*10^3/1.518)^0.5*x1/x2^2)^(1/3)+(1.518*x2^4*x8^2*x9/(2+2*(5.305*10^3/1.518)^0.5*x1/x2^2)^4)^(1/6)+2*(2.012*x3)^(1/2)-1;

回复此楼

» 猜你喜欢

大豆异黄酮分离已经有0人回复
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有79人回复
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名已经有1人回复
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生已经有5人回复
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究已经有12人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有1人回复
如何从铁电相到顺电相。已经有1人回复

1楼 2018-04-28 11:11:25

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

quanluhan

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3
帖子: 7
在线: 4.2小时
虫号: 3414828
注册: 2014-09-13
专业: 工程热物理相关交叉领域

楼主解决问题没？我也遇到了相同的问题，我也是解9元非线性方程组，感觉结果明显不对，但是代码也不知道哪里错了。。。请大神赐教啊

赞一下

回复此楼

2楼2018-06-08 14:59:25

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

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★
独孤神宇: 金币+1 2018-06-08 19:36:54

一楼的问题还是可以大概率得到唯一解的：

x1: 0.000415258421410393
x2: 0.0295645164484108
x3: 1.36712379593572E-5
x4: 2.29288172384266E-14
x5: 2.47606393479552E-9
x6: 1.35850617795293E-7
x7: 2.67957126216714E-9
x8: 0.427499066472031
x9: 0.121967698482867

赞一下(1人)

回复此楼

3楼2018-06-08 17:19:16

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