数学建模社区-数学中国

标题: 非线性方程组求解 [打印本页]
作者: shimian_2014    时间: 2014-9-2 18:38
标题: 非线性方程组求解
求助大神指点

非线性方程组求解问题.docx

27.64 KB, 下载次数: 7, 下载积分: 体力 -2 点

问题
作者: madio    时间: 2014-9-3 22:39
就是列方程求解吧，4个未知量，可以列5个方程，肯定可以解出来的
作者: 847717213    时间: 2014-9-3 23:53
还有Cz0 Cp0参数两个参数呢，你不会让求符号表达式吧！！！
（数模咨询http://shop108557885.taobao.com）
作者: 赵全    时间: 2014-9-4 21:25
不解啊，是否还有更好的方法
作者: mingtingqing    时间: 2014-10-6 16:00
将数值代进去，然后利用Lingo编程可以求解
作者: liwenhui    时间: 2017-6-5 17:59
LINGO代码如下：
  1. SETS:
  2. VAR/ 1..5/: T, F, CZ, CP;
  3. ENDSETS
  4. @FOR( VAR( I):
  5.   CZ( I) = -( KZ2 / KZ1) * T( I) - ( KZ2 / KZ1^2) + CZ0 * @EXP( KZ1 * T( I));
  6.   CP( I) = -( KP2 / KP1) * T( I) - ( KP2 / KP1^2) + CP0 * @EXP( KP1 * T( I));  
  7.   CZ( I) - CP( I) = F( I);  
  8. );
  9. DATA:
  10. F = 95.27, 62.85, 41.47, 27.36, 18.05;
  11. ENDDATA
  12. CALC:
  13. @SET( 'GLOBAL', 1);
  14. ENDCALC
复制代码
计算结果：
  1.   Feasible solution found.
  2.   Infeasibilities:                             0.3886242E-09
  3.   Extended solver steps:                               0
  4.   Total solver iterations:                          2827


  7.                                            Variable           Value
  8.                                                 KZ2        0.000000
  9.                                                 KZ1        99999.97
  10.                                                 CZ0        18.44234
  11.                                                 KP2        0.000000
  12.                                                 KP1        100000.1
  13.                                                 CP0       0.3923384
  14.                                               T( 1)       0.1663570E-04
  15.                                               T( 2)       0.1247606E-04
  16.                                               T( 3)       0.8318249E-05
  17.                                               T( 4)       0.4159366E-05
  18.                                               T( 5)        0.000000
  19.                                               F( 1)        95.27000
  20.                                               F( 2)        62.85000
  21.                                               F( 3)        41.47000
  22.                                               F( 4)        27.36000
  23.                                               F( 5)        18.05000
  24.                                              CZ( 1)        97.34081
  25.                                              CZ( 2)        64.21612
  26.                                              CZ( 3)        42.37140
  27.                                              CZ( 4)        27.95470
  28.                                              CZ( 5)        18.44234
  29.                                              CP( 1)        2.070812
  30.                                              CP( 2)        1.366122
  31.                                              CP( 3)       0.9014011
  32.                                              CP( 4)       0.5947027
  33.                                              CP( 5)       0.3923384

  35.                                                 Row    Slack or Surplus
  36.                                                   1        0.000000
  37.                                                   2        0.000000
  38.                                                   3        0.000000
  39.                                                   4        0.000000
  40.                                                   5        0.000000
  41.                                                   6        0.000000
  42.                                                   7        0.000000
  43.                                                   8        0.000000
  44.                                                   9        0.000000
  45.                                                  10        0.000000
  46.                                                  11        0.000000
  47.                                                  12        0.000000
  48.                                                  13        0.000000
  49.                                                  14        0.000000
  50.                                                  15        0.000000
复制代码





欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5