| 查看: 672 | 回复: 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的时候方程组的解。 |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有21人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
3楼2019-05-09 10:08:03