| 查看: 665 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
1stopt求解变参数非线性方程组结果只有一组解 已有2人参与
|
||
|
用1stopt5.0运行的时候，非线性方程只有一组解，实际上应该有多组解 Constant A=2,gamma1=0.1,gamma2=0.2,k=500,w0=15,eps=0.1; LoopConstant omega=[10:1:20] ; Parameters x1,x2,x3,x4; Function -x2*omega+gamma1*x1+w0^2*x2/omega+3*k*x4*(x3^2+x4^2)/(4*omega^3)+gamma2*x3=A; x1*omega+x2*gamma1-w0^2*x1/omega-3*k*x3*(x3^2+x4^2)/(4*omega^3)+gamma2*x4=0; x1*gamma1-x4*omega+w0^2*x2/omega+(1+1/eps)*(3*k*x4*(x3^2+x4^2)/(4*omega^3)+gamma2*x3)=A; x2*gamma1+x3*omega-w0^2*x1/omega-(1+1/eps)*(3*k*x3*(x3^2+x4^2)/(4*omega^3)-gamma2*x4)=0; 结果显示： Constant A=2,gamma1=0.1,gamma2=0.2,k=500,w0=15,eps=0.1; LoopConstant omega=[10:1:20] ; Parameters x1,x2,x3,x4; Function -x2*omega+gamma1*x1+w0^2*x2/omega+3*k*x4*(x3^2+x4^2)/(4*omega^3)+gamma2*x3=A; x1*omega+x2*gamma1-w0^2*x1/omega-3*k*x3*(x3^2+x4^2)/(4*omega^3)+gamma2*x4=0; x1*gamma1-x4*omega+w0^2*x2/omega+(1+1/eps)*(3*k*x4*(x3^2+x4^2)/(4*omega^3)+gamma2*x3)=A; x2*gamma1+x3*omega-w0^2*x1/omega-(1+1/eps)*(3*k*x3*(x3^2+x4^2)/(4*omega^3)-gamma2*x4)=0; 这里我只放了omega=20的时候方程组的解。 |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有49人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
三无产品还有机会吗
已经有6人回复
3楼2019-05-09 10:08:03
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
版本问题吧，有11组解： 循环常数 omega 目标函数值 x1 x2 x3 x4 10 4.58635664500889E-26 0.00380294459721743 0.160294854647957 -0.0277043774026065 0.156228877092279 11 7.91721570819433E-27 0.00669308907034421 0.211859317209918 -0.0315743501340042 0.208493086819324 12 1.34137963790563E-24 0.0132338824083539 0.296344281124706 -0.0365113493760948 0.294971542948296 13 1.01936620072096E-23 0.0328237831973817 0.461292243502208 -0.040457396825942 0.468686530302756 14 9.15535540666567E-23 0.138909893810943 0.895184889404847 -0.00267827329179549 0.989325928727004 15 1.27077263082487E-20 0.746988968174331 1.13840265715959 0.822882483116559 2.42193671971171 16 1.52179093252339E-21 1.20677974111929 0.376059697385337 2.69161830257015 2.00602076273742 17 3.86982718524109E-24 0.042338001880337 -0.52754140284278 0.105065051557478 -0.52181906319924 18 1.78381813142587E-25 0.0196379178791081 -0.36165925952281 0.0595504858703478 -0.356709103497778 19 6.20181994642092E-26 0.011584651009512 -0.278277662864219 0.0405814648006044 -0.274614864703752 20 1.31099134650033E-27 0.00775740966522293 -0.227852327097322 0.0303034605957675 -0.225094128630944 |
2楼2019-04-16 10:45:05
sx644030092
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 117.7
- 散金: 40
- 红花: 1
- 帖子: 114
- 在线: 25.5小时
- 虫号: 4707607
- 注册: 2016-05-20
- 性别: MM
- 专业: 配位化学
4楼2019-05-23 10:22:05