| 查看: 675 | 回复: 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的时候方程组的解。 |
回复此楼» 猜你喜欢
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有6人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 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