求助 一个简单的多元一次求和方程怎么编写LINGO

liwenhui

构造最小二乘解法；

1. SETS:
   
2.   VAR/ 1..6/: A, B, C;
   
3. ENDSETS
   
4. 
   
5. MIN = Z;
   
6. 
   
7. Z = @SUM( VAR( I): ( A( I)/ B( I) - C( I))^2);
   
8. 
   
9. @SUM( VAR( I): A( I)) = 57442053.49;
   
10. 
    
11. @SUM( VAR( I): B( I)) = 644527482.3;
    
12. 
    
13. DATA:
    
14. 
    
15. C = 0.05, 0.1, 0.2, 0.5, 0.8, 1;
    
16. 
    
17. ENDDATA
    
18. 
    
19. CALC:
    
20. 
    
21. @SET( 'GLOBAL', 1);
    
22. 
    
23. ENDCALC

复制代码

求解结果：

1. Global optimal solution found.
   
2.   Objective value:                             0.6000000E-09
   
3.   Objective bound:                             0.6000000E-09
   
4.   Infeasibilities:                             0.4348427E-07
   
5.   Extended solver steps:                               2
   
6.   Total solver iterations:                           634
   
7. 
   
8. 
   
9.                        Variable           Value        Reduced Cost
   
10.                               Z        0.000000            1.000000
    
11.                           A( 1)       0.3056445E+08        0.000000
    
12.                           A( 2)       0.9674469E-05        0.000000
    
13.                           A( 3)       0.5568388E-05        0.000000
    
14.                           A( 4)       0.2797364E-05        0.000000
    
15.                           A( 5)       0.2687761E+08        0.000000
    
16.                           A( 6)       0.3339973E-04        0.000000
    
17.                           B( 1)       0.6109298E+09        0.000000
    
18.                           B( 2)       0.9673430E-04        0.000000
    
19.                           B( 3)       0.2784192E-04        0.000000
    
20.                           B( 4)       0.5594727E-05        0.000000
    
21.                           B( 5)       0.3359770E+08        0.000000
    
22.                           B( 6)       0.3339973E-04        0.000000
    
23.                           C( 1)       0.5000000E-01        0.000000
    
24.                           C( 2)       0.1000000            0.000000
    
25.                           C( 3)       0.2000000            0.000000
    
26.                           C( 4)       0.5000000            0.000000
    
27.                           C( 5)       0.8000000            0.000000
    
28.                           C( 6)        1.000000            0.000000
    
29. 
    
30.                             Row    Slack or Surplus      Dual Price
    
31.                               1        0.000000           -1.000000
    
32.                               2      -0.1247501E-08        0.000000
    
33.                               3      -0.1219200E-08        0.000000
    
34.                               4      -0.4348427E-07        0.000000
    
35. 
    

复制代码