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;

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有276人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复

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的回帖