24小时热门版块排行榜    

查看: 1098  |  回复: 3

hitcs

铜虫 (初入文坛)

[求助] 求助,用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;
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quanluhan

新虫 (初入文坛)

楼主 解决问题没?我也遇到了相同的问题,我也是解9元非线性方程组,感觉结果明显不对,但是代码也不知道哪里错了。。。请大神赐教啊
2楼2018-06-08 14:59:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


独孤神宇: 金币+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
3楼2018-06-08 17:19:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

状元糕

铁虫 (小有名气)

大神,能不能分享一份1stopt5.o给我,非常谢谢!2207929801@qq.com
潜力无限
4楼2018-07-04 22:54:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hitcs 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 能源材料化学课题组招收硕士研究生8-10名 +3 脱颖而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[考研] 281求调剂(0805) +3 烟汐忆海 2026-03-16 7/350 2026-03-17 20:16 by peike
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 312求调剂 +4 陌宸希 2026-03-16 5/250 2026-03-17 17:09 by ruiyingmiao
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 8/400 2026-03-17 16:57 by ruiyingmiao
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 9/450 2026-03-17 10:29 by xujiaoszu
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
信息提示
请填处理意见