声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1975|回复: 2

[1stopt] 求解非线性方程的问题

[复制链接]
发表于 2008-4-22 20:12 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我用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;

待解方程

待解方程
回复
分享到:

使用道具 举报

发表于 2008-4-22 22:24 | 显示全部楼层
1:如果用的是网上流行的1.5PJ版,这种编程模式是没法运行的。
2:图片中的第一项公式中含有xk,这个下标k指的什么?好像不大对,是否有笔误?
 楼主| 发表于 2008-4-23 00:20 | 显示全部楼层
公式里的fk就是F矩阵
Fr为Fi

我就是用网上的PJ版。不能调试运行的。
能帮我求解一下吗?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-20 07:04 , Processed in 0.051038 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表