应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 求助，用1stopt 5.0解非线性方程组，两次得到的解不唯一？什么原因，大神帮忙
41/1返回列表
查看: 1252  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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;
回复此楼

» 猜你喜欢
1楼 2018-04-28 11:11:25
已阅   回复此楼   关注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
赞一下(1人) 回复此楼
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 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天基金会出结果吗？20260819 +9 kkkl_v 2026-08-19 9/450 2026-08-19 19:11 by 菜鸟遥感
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +53 医学老男孩 2026-08-13 118/5900 2026-08-19 18:14 by 医学老男孩
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 3/150 2026-08-19 18:04 by oCMhZAsgNFoO
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 PnssoZglMj3v 2026-08-19 4/200 2026-08-19 17:54 by oCMhZAsgNFoO
[基金申请] filecode，4个jtjc了 +10 ziyangfang 2026-08-19 11/550 2026-08-19 17:19 by wdc28
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +4 瞬息宇宙 2026-08-15 4/200 2026-08-19 16:23 by allmydream
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 时间戳变了，能看出什么问题？ +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
信息提示
关闭
请填处理意见
关闭