1stopt求助

gouxu

用该软件求解一个34维的非线性方程组，版本是1.5专业版，遇到问题：遇到常数定义有误！
有关方程组的信息如下，望各位高手帮忙，为了这个方程组，我在这里卡壳几天了。。。

Title "node9";
//Parameters
x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34;
//Variable;
//Function
   x1-1=0;
   0.3479*x26*x4*cos(x17-x1)+0.3479*x26*x3*sin(x17-x1)-0.20304*x2+0.72654=0;
   0.1116*x5-0.268*x3+0.1564*x26*cos(x17-x1)=0;
   -5.1411*x4-1.9153*x26*sin(x17-x1)=0;
   -5*x5-5*x26+500=0;
   x7-1=0;
   0.6519*x27*x9*cos(x18-x6)+0.6519*x27*x8*sin(x18-x6)-0.19525*x7+0.3227=0;
   0.1667*x10-1.2465*x8+1.0798*x27*cos(x18-x6)=0;
   -13.4885*x9-11.6193*x27*sin(x18-x6)=0;
   -2.8571*x10-571.42*x27+628.562=0;
   x12-1=0;
   0.9162*x28*x14*cos(x19-x11)+0.9162*x28*x13*sin(x19-x11)+0.3468=0;
   0.1164*x15-0.8427*x13+0.7263*x28*cos(x19-x11)=0;
   -11.5631*x14-9.8964*x28*sin(x19-x11)=0;
   -2.8571*x15-285.71*x28+328.5665=0;
   -10000*x30-100*x16+10000=0;
   16.4474*x26*x3*sin(x1-x17)-16.4474*x26*x4*cos(x1-x17)+17.2414*x26*x29*sin(x17-x20)=0;
   16.4474*x26*x3*cos(x1-x17)+16.4474*x26*x4*sin(x1-x17)-16.4474*x26^2-17.2414*x26*x29*cos(x17-x20)=0;
   8.3472*x27*x8*sin(x6-x18)-8.3472*x27*x9*cos(x6-x18)+15.873*x27*x32*sin(x18-x23)=0;
   8.3472*x27*x8*cos(x6-x18)+8.3472*x27*x9*sin(x6-x18)-8.3472*x27^2-15.873*x27*x32*cos(x18-x23)=0;
   5.5157*x28*x13*sin(x11-x19)-5.5157*x28*x14*cos(x11-x19)+16.9492*x28*x34*sin(x19-x25)=0;
   5.5157*x28*x13*cos(x11-x19)-5.5157*x28*x14*sin(x11-x19)-5.5157*x28^2-16.9492*x28*x34*cos(x19-x25)=0;
   -17.2414*x26*x29*sin(x20-x17)+x29*x30*(1.3652*cos(x20-x21)-11.6041*sin(x20-x21))+x29*x31*(1.9422*cos(x20-x22)-10.5107*sin(x20-x22))=0;
   17.2414*x26*x29*cos(x20-x17)+x29*x30*(1.3652*sin(x20-x21)+11.6041*cos(x20-x21))+x29*x31*(1.9422*sin(x20-x22)+10.5107*cos(x20-x22))=0;
   x29*x30*(1.3652*cos(x21-x20)-11.6041*sin(x21-x20))+x30*x32*(1.1876*cos(x21-x23)-5.9751*sin(x21-x23))-1.25=0;
   x29*x30*(1.3652*sin(x21-x20)+11.6041*sin(x21-x20))+x30*x32*(1.1876*sin(x21-x23)+5.9751*cos(x21-x23))-0.5=0;
   x29*x31*(1.9422*cos(x22-x20)-10.5107*sin(x22-x20))+x31*x34*(1.282*cos(x22-x25)-5.5882*sin(x22-x25))-0.9=0;
   x29*x31*(1.9422*sin(x22-x20)+10.5107*cos(x22-x20))+x31*x34*(1.282*sin(x22-x25)+5.5882*sin(x22-x25))-0.3=0;
   -15.873*x27*x32*sin(x23-x18)+x30*x32*(1.1876*cos(x23-x21)-5.9751*sin(x23-x21))+x32*x33*(1.7094*cos(x23-x24)-13.6752*sin(x23-x24))=0;
   15.873*x27*x32*cos(x23-x18)+x30*x32*(1.1876*sin(x23-x21)+5.9751*cos(x23-x21))+x32*x33*(1.7094*sin(x23-x24)+13.6752*cos(x23-x24))=0;
   x32*x33*(1.7094*cos(x24-x23)-13.6752*sin(x24-x23))+x33*x34*(1.16*cos(x24-x25)-9.7632*sin(x24-x25))-1=0;
   x32*x33*(1.7094*sin(x24-x23)+13.6752*cos(x24-x23))+x33*x34*(1.16*sin(x24-x25)+9.7632*cos(x24-x25))-0.3=0;
   -16.9492*x28*x34*sin(x25-x19)+x31*x34*(1.282*cos(x25-x22)-5.5882*sin(x25-x22))+x33*x34*(1.16*cos(x25-x24)-9.7632*sin(x25-x24))=0;
   16.9492*x28*x34*cos(x25-x19)+x31*x34*(1.282*sin(x25-x22)+5.5882*cos(x25-x22))+x33*x34*(1.16*sin(x25-x24)+9.7632*cos(x25-x24))=0;

dingd

把其中的：
//Parameters
x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34;
//Variable;
//Function
改为：
Title "node9";
Parameters
x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34;
Function

另外建议找个使用手册看看，掌握基本的语法。还有将诸如“x1-1=0;”直接当常数就行了，减少不必要的计算量。

[ 本帖最后由 dingd 于 2009-2-25 09:32 编辑 ]

gouxu

谢谢版主！
官方的使用手册太简单了，有些地方写的有些敷衍了事。
这个方程组解了很多次，还是得不到理想的结果，方程（3） 0.1116*x5-0.268*x3+0.1564*x26*cos(x17-x1)=0;的值总是在11左右。该怎么处理呢？
另外，方程（1）x1-1=0应为x2-1=0。

gouxu

求得下面这组解
x=[1.00000000053084;1;41.6417935275771;1.67918155141273E-10;99.9999999408078;-1.15789143407611;
   1.00000032375681;1.08846539560608;0.00508615927325335;0.913894288129094;0.301315392179175;1.00000000042774;
   1.10898333813451;-0.010253319288071;4.11927280527527;0.414911604801637;-268.878371307429;-7.44648624437484;
   0.312118906165827;-14.4898148359848;3.55062671318044;553.242625189148;4.99311734397139;-1.71722030142607;
   -1.77497261528091;1.24135421545135E-11;1.09543052741996;1.10880726816146;-0.0533297747617175;
  0.995850883461644;1.11197385761069;-0.00368604660608643;-6.93352486500268;-2.66058550721011E-9;];
对应的函数值：
f=[0,0.5235,-0.0000,-0.0000,0.0000,0.0000,0.1269,-0.0216,0.0002,0.0000,0.0000,0.3486,0.3502,0.0000,0.0000,0.0000,-0.0000,0.0000,-0.0008,-0.0000,-0.0106,
0.0000,-0.0187,0.0016,-1.7683,-0.0000,-0.3414,-0.0004,-0.0919,0.2655,-0.8154,-0.0000,-0.0000,0.0000]
不知道这样一组解可不可行。因为x当中有些地方显然不合乎物理意义，但是求解了这么多次，也只有这一组比较好了。