로드 중...
검색 중...
일치하는 항목 없음
millco.py
1
17
18if__name__ =="__main__":
19
20 # [1]
21 from슬롯 사이트magic가져오기GamsInteractive
22슬롯 사이트 = GamsInteractive()
23
24 # [2]
25 가져오기numpyasnp
26m = 슬롯 사이트exchange_container
27
28 # 모델 세트 및 매개변수
29사이트 = m.addSet('사이트', 레코드 = ["1", "2", "3"])
30밀스 = m.addSet('밀스', 레코드 = ["밀 A","밀 B", "밀 C"])
31dist = m.addParameter('거리', [사이트, 공장], 레코드 = np.array([[ 8,15,50],
32 [10,17,20],
33 [30,26,15]]))
34공급 = m.addParameter('공급', [사이트], 레코드 = np.array([20,30,45]))
35수요 = m.addParameter('수요', [밀스], 레코드 = np.array([30,35,30]))
36cost_per_haul = m.addParameter('cost_per_haul', 레코드 = 4)
37
38 # 모델 변수
39m.addAlias(',사이트)
40m.addAlias('', 밀스)
41배송 = m.addVariable('배', '긍정적', [현장, 공장])
42obj = m.addVariable('obj', '무료')
43
44 # [3]
45슬롯 사이트슬롯 사이트('''
46# 슬롯 사이트 모델
47방정식 defcost; defcost.. cost_per_haul*sum((s,m), ship(s,m)*dist(s,m)) =e= obj;
48방정식 공급 부족; defsupply(들).. sum(m, ship(s,m)) =e= 공급(들);
49방정식 수요; defdemand(m).. sum(s, ship(s,m)) =e= 수요(m);
50
51모델 밀코 / 모두 /;''')
52
53 # [4]
54슬롯 사이트슬롯 사이트('lp를 사용하여 millco min obj를 해결하세요;')
55
56 # [5]
57인쇄(m["배"].pivot())
58인쇄(f'총 비용은 obj.records["level"][0]'이 될 것입니다)
59
60 # [6]
61슬롯 사이트create('젠모델')
62m = 슬롯 사이트exchange_container
63
64 # 일반 모델 세트 및 매개변수
65i = m.addSet('나', 레코드 = ['나'+str(i)foriin범위(6)], 설명 ='방정식 색인')
66j = m.addSet('j', 레코드 = ['j'+str(j)forj에범위(9)], 설명 ='변수 인덱스')
67A = m.addParameter('A', [i,j], 레코드 = np.array([[1,1,1,0,0,0,0,0,0],
68 [0,0,0,1,1,1,0,0,0],
69 [0,0,0,0,0,0,1,1,1],
70 [1,0,0,1,0,0,1,0,0],
71 [0,1,0,0,1,0,0,1,0],
72 [0,0,1,0,0,1,0,0,1]]))
73b = m.addParameter('b', [i], 레코드 = np.array([20,30,45,30,35,30]))
74c = m.addParameter('c', [j], 레코드 = np.array([8,15,50,10,17,20,30,26,15]))
75x = m.addVariable('x', '긍정적', [j])
76obj = m.addVariable('obj', '무료')
77
78 # [7]
79슬롯 사이트슬롯 사이트('''
80# 일반 슬롯 사이트 모델
81수식 defobj; defobj.. 4*sum(j, c(j)*x(j)) =e= obj;
82방정식 e(i); e(i).. sum(j, A(i,j)*x(j)) =e= b(i);
83
84모델 세대 / 모두 /;''')
85
86 # [8]
87슬롯 사이트슬롯 사이트('lp를 사용하여 gen min obj를 해결하세요;')
88
89 # [9]
90인쇄(x.records)
91인쇄(f'총 비용은 obj.records["level"][0]'이 될 것입니다)
92
93 # [10]
94슬롯 사이트gams_cleanup(closedown=사실입니다)