-
Notifications
You must be signed in to change notification settings - Fork 0
/
tst5.py
59 lines (52 loc) · 4.01 KB
/
tst5.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
import numpy as np
import matplotlib.pyplot as plt
macB,mocB,madB,modB,mohB = 1.47/3.3,0.6385/3.3,2.39/5.5,1.0372/5.5,0.9075/2.4
print(macB,mocB,madB,modB,mohB)
exit(5)
ub = np.array([18787.59321455486,18641.62006757802,15892.479356817481,15082.336542737812,12093.118453854811,11692.063942099503,11692.063942099503,10803.460397877916,10803.460397877916,10803.460397877916,7994.483304406775,7994.483304406775,7715.080702706282,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003,7671.669252414003])
lb = np.array([785.212222222205,2529.448333333302,2529.448333333302,2529.448333333302,2603.4536050771057,3505.7778597225515,4489.954303629127,4793.319581406923,5048.183664048219,5189.044605664496,6028.141590443057,6366.307650275304,6366.307650275305,6856.002243992102,6882.681832410486,6923.080443522095,6964.824649372798,6969.188125943316,6973.679327886701,6974.491032229499,6984.157415086197,6986.432287685,6993.823797942897,7003.490180799599,7263.036107297588,7302.416833252868,7369.53985207691,7520.933422979754,7520.933422979754,7552.024019150045,7595.775531357431,7639.543830558536,7642.003246593753,7643.1656830065,7671.669252414003])
f,a = plt.subplots()
t = len(lb)
t = np.arange(t)
a.plot(t,lb)
a.plot(t,ub)
plt.show()
exit(5)
for i in range(100):
if i % 10 == 1:
print(i)
exit(5)
lb = np.array([441.448,1140.18,1417.41,3893.65,4945.83,5148.36,7780.15,8004.81,8243.81,8630.21,8732.46,8857.74,9106.49,9407.55,9550.38,9634.65,9771.09,9777.15,10058.8,10067.3,10069.9,10087.9, 10237,10256.5,10290.6,10388.9,10591.4,10776.2,10808.7,10944.6,11325.8,11355.7,11397.5,11424.4,11431.8,11436.4,11446.7,11455.7,11493.4,11494.6,11501.5,11506.5,11509.7,11512.9,11518.6,11518.9,11536.3,11536.8,11539.1,11539.4,11555.8,11572.7,11589.9,11590.6,11602.9,11604.4,11613.3,11613.6,11624.8,11627.2,11628.3,11628.4,11631.5,11634.8,11634.6,11636.3,11640.2,11643.1,11645.5,11646.3,11646.6,11652.7,11655.3,11657.9,11663.7,11668.6,11674.7,11678.1,11678.1,11678.1])
ub = np.array([28581.3, 18057,14126.7,19140.5,24568.9,12489.9,21370.6,15703.1,14669.9,15728.2,17888.1,13346.5,11733.5,21297.6,14483.3,15709.7,13817.1,12066.7,11867.6,11739.2,11750.5,11719.8,16534.1,18719.5,15510.9,17678.7,14004.8,12872.9,13869.5, 11689,16580.5,13008.9,11705.9,12759.2,12402.1,12167.4,12157.3,12476.7,12449.4,12185.8,21564.7, 15778,12168.7,15730.7,11693.4,11695.5,12266.5,13950.5,11693.8, 11716, 11694,11990.9,11787.5,12007.6,11739.4,12060.8,11726.8,11750.1,12209.9,11976.8,12116.2,11940.5,11710.9,11711.2,11701.5,11756.6,11959.2,11713.2,11959.4,11793.2,11827.6, 11837,11686.5,11853.3,11850.2, 11869,11669.7,11678.1,11678.1,11678.1])
f,a = plt.subplots()
t = len(lb)
t = np.arange(t)
a.plot(t,lb)
a.plot(t,ub)
plt.show()
# t = 0
# purchase P_h_Q = |0.152206| 0.24| 0|
# chp0: P_h_Q = | 0.4| 0| 0|
# chp0: P_h_Q = |0.342265| 0| 0.1|
# caes c/d/h = |0.494471| 0| 0|
# load: P_h_Q = | 0.4| 0.24| 0.1|
# t(0)| 95| 96|t = 1
# purchase P_h_Q = |0.123207|0.193812| 0|
# chp0: P_h_Q = |0.342265|0.046188| 0.1|
# chp0: P_h_Q = | 0.4| 0| 0|
# caes c/d/h = |0.465472| 0| 0|
# load: P_h_Q = | 0.4| 0.24| 0.1|
# t(1)| 124| 124|t = 2
# purchase P_h_Q = | 0|0.369858| 0|
# chp0: P_h_Q = | 0.25| 0.12| 0.25|
# chp0: P_h_Q = |0.262322|0.110142| 0|
# caes c/d/h = | 0|0.487678| 0|
# load: P_h_Q = | 1| 0.6| 0.25|
# t(2)| 95| 94|t = 3
# purchase P_h_Q = |0.535684|0.775248| 0|
# chp0: P_h_Q = | 0.31906|0.0647521|0.140192|
# chp0: P_h_Q = | 0.25| 0.12|0.259808|
# caes c/d/h = | 0|0.495256| 0|
# load: P_h_Q = | 1.6| 0.96| 0.4|
# penalty = | 0| 0|
# t(3)| 64| 64|