| 查看: 1110 | 回复: 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; |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有123人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有2人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2018-06-08 17:19:16
2楼2018-06-08 14:59:25
状元糕
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 120.3
- 帖子: 61
- 在线: 41.4小时
- 虫号: 1675086
- 注册: 2012-03-08
- 性别: GG
- 专业: 金属材料表面科学与工程
4楼2018-07-04 22:54:51
