用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的时候方程组的解。 |