非线性方程组求解
求助大神指点就是列方程求解吧,4个未知量,可以列5个方程,肯定可以解出来的 还有Cz0 Cp0参数两个参数呢,你不会让求符号表达式吧!!!
(数模咨询http://shop108557885.taobao.com) 不解啊,是否还有更好的方法 将数值代进去,然后利用Lingo编程可以求解 LINGO代码如下:SETS:
VAR/ 1..5/: T, F, CZ, CP;
ENDSETS
@FOR( VAR( I):
CZ( I) = -( KZ2 / KZ1) * T( I) - ( KZ2 / KZ1^2) + CZ0 * @EXP( KZ1 * T( I));
CP( I) = -( KP2 / KP1) * T( I) - ( KP2 / KP1^2) + CP0 * @EXP( KP1 * T( I));
CZ( I) - CP( I) = F( I);
);
DATA:
F = 95.27, 62.85, 41.47, 27.36, 18.05;
ENDDATA
CALC:
@SET( 'GLOBAL', 1);
ENDCALC计算结果: Feasible solution found.
Infeasibilities: 0.3886242E-09
Extended solver steps: 0
Total solver iterations: 2827
Variable Value
KZ2 0.000000
KZ1 99999.97
CZ0 18.44234
KP2 0.000000
KP1 100000.1
CP0 0.3923384
T( 1) 0.1663570E-04
T( 2) 0.1247606E-04
T( 3) 0.8318249E-05
T( 4) 0.4159366E-05
T( 5) 0.000000
F( 1) 95.27000
F( 2) 62.85000
F( 3) 41.47000
F( 4) 27.36000
F( 5) 18.05000
CZ( 1) 97.34081
CZ( 2) 64.21612
CZ( 3) 42.37140
CZ( 4) 27.95470
CZ( 5) 18.44234
CP( 1) 2.070812
CP( 2) 1.366122
CP( 3) 0.9014011
CP( 4) 0.5947027
CP( 5) 0.3923384
Row Slack or Surplus
1 0.000000
2 0.000000
3 0.000000
4 0.000000
5 0.000000
6 0.000000
7 0.000000
8 0.000000
9 0.000000
10 0.000000
11 0.000000
12 0.000000
13 0.000000
14 0.000000
15 0.000000
页:
[1]