数学建模社区-数学中国

标题: 用lingo解100*100的线性方程组怎么编程表述出方程组 [打印本页]

---

作者: 半度微笑    时间: 2015-7-20 12:36
标题: 用lingo解100*100的线性方程组怎么编程表述出方程组
大神们，100*100的线性方程组。我这里用9*9作为例子

sets:

r/1..9/:;

c/1..9/:;

link(r,c):a;

time/1..9/:t;

dose/1..9/:d;

endsets

data:

a=14 14 14 9 6 4 3 9 0

  14 14 14 14 9 6 4 1 2

  14 14 14 14 14 9 6 1 1

   9 14 14 14 14 14 9 1 1

   6 9 14 14 14 14 14 2 1

   4 6 9 14 14 14 14 3 2

   3 4 6 9 14 14 14 4 3

   0 0 1 1 1 2 3 14 14

   0 1 1 1 2 3 4 14 14；

d=500 500 500 500500 500 500 500 500;

enddata

min=@sum(link:a*t);

@for(dose(J):

@sum(r(I):t(J))<=dose(J));

@for(time(I):

);

End

---

作者: 森之张卫东    时间: 2015-7-20 20:31

1. sets:
   
2. 
   
3. column/1..9/:c,x;
   
4. 
   
5. row/1..9/:d;
   
6. 
   
7. link(row,column):a;
   
8. 
   
9. !time/1..9/:;
   
10. 
    
11. !dose/1..9/:;
    
12. 
    
13. endsets
    
14. 
    
15. data:
    
16. 
    
17. a=14 14 14 9 6 4 3 9 0 14 14 14 14 9 6 4 1 2 14 14 14 14 14 9 6 1 1
    
18. 
    
19.    9 14 14 14 14 14 9 1 1 6 9 14 14 14 14 14 2 1 4 6 9 14 14 14 14 3 2
    
20. 
    
21.    3 4 6 9 14 14 14 4 3 0 0 1 1 1 2 3 14 14 0 1 1 1 2 3 4 14 14;
    
22. 
    
23. d=500 500 500 500 500 500 500 500 500;
    
24. 
    
25. c=1 1 1 1 1 1 1 1 1 ;
    
26. 
    
27. enddata
    
28. 
    
29. min=@sum(link:a*x);
    
30. 
    
31. @for(row(i):@sum( column(j):a(i,j)*x(j) ) < d(i) );
    
32. 
    
33. 
    
34. End

复制代码

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5