求解非线性方程的问题
我用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; 1:如果用的是网上流行的1.5PJ版,这种编程模式是没法运行的。
2:图片中的第一项公式中含有xk,这个下标k指的什么?好像不大对,是否有笔误? 公式里的fk就是F矩阵
Fr为Fi
我就是用网上的PJ版。不能调试运行的。
能帮我求解一下吗?
页:
[1]