kangtj 发表于 2008-4-22 20:12

求解非线性方程的问题

我用1stOpt编了个程序(源码如下),可是调试运行的时候没反应
不知道怎么回事。。有用1stOpt的高手帮我看看啊。。

Title "maxtrix-qiujie";
Parameters x(1:5)[-1,1];
//Variable ;
Constant F(1:10,1:5)=[

-123.401986666667,      15309.8526427374,      -1909789.96410041,      239549341.560198,      -30214602491.039,
-87.6375305,                7741.2805593795,      -689528.27440591,       61952373.811834,      -5616144306.013,
-59.06265104918,      3586.5719185062,      -223770.00498736,      14320890.485566,      -937932650.125,
-28.439255018957,      879.2160346562,                -29081.78717329,      1013450.252394,                -36742842.333,
2.815480861073,                87.4964063986,                713.65902354,                16697.447204,                196487.477,
33.027548857143,      1173.1652538497,      44394.49294689,                1769229.311412,                7343630.493 ,
64.045056222222,      4201.7966978657,      282266.01166872,      19393922.004712,      1360598257.263,
93.538933217391,      8819.304635927,                838503.24465969,      80419975.012864,      7783014155.159,
127.24891,                16253.0044120957,      2083630.02919474,      268092890.315688,      34616841265.119,
173.02841,               29985.3187910965,      5204399.56611482,      904686543.289774,      157501253705.876];

Constant Fr(1:5)=[
   1.05137753485,
   1957.02063533292 ,
   8214.58926666724 ,
   12077012.3703088,
   219286227.91329];
StartProgram;
Var   i,j,k: integer;
   sum,p,q,Z: double;
Begin
sum:=0 ;
i:= 1;
j:= 1;
k:= 1;

for k := 1 to 5 do
   q(k):=0;
   Z=0;
for i := 1 to 10 do

      for j:=1 to 5 do
         sum=F(i,j)*x(j)+sum;
      End;
      p(i)=exp(-sum);
      q(k)=p(i)*F(k,i)+q(k);
      Z=Z+exp(-sum);
End;
FunctionResult:= Fr(k)*Z-q(k);

End;
EndProgram;

dingd 发表于 2008-4-22 22:24

1:如果用的是网上流行的1.5PJ版,这种编程模式是没法运行的。
2:图片中的第一项公式中含有xk,这个下标k指的什么?好像不大对,是否有笔误?

kangtj 发表于 2008-4-23 00:20

公式里的fk就是F矩阵
Fr为Fi

我就是用网上的PJ版。不能调试运行的。
能帮我求解一下吗?
页: [1]
查看完整版本: 求解非线性方程的问题