110 lines
4.2 KiB
Python
110 lines
4.2 KiB
Python
|
#!/usr/bin/env python3
|
||
|
|
import numpy as np
|
||
|
|
import csv
|
||
|
|
import itertools
|
||
|
|
import math
|
||
|
|
|
||
|
|
column_label=0
|
||
|
|
column_capacity=3
|
||
|
|
column_resistance=4
|
||
|
|
|
||
|
|
|
||
|
|
config_parallel=2
|
||
|
|
config_cells=12
|
||
|
|
|
||
|
|
#config_parallel*config_cells*config_packs needed
|
||
|
|
batteriesneeded=config_parallel*config_cells
|
||
|
|
|
||
|
|
|
||
|
|
batteries=[]
|
||
|
|
batterylabels=[]
|
||
|
|
|
||
|
|
with open('12S_LiIon_Akkupack_selected.csv', 'r') as csvfile:
|
||
|
|
csvreader = csv.reader(csvfile, delimiter=';')
|
||
|
|
firstrow=True
|
||
|
|
for row in csvreader:
|
||
|
|
label=row[column_label]
|
||
|
|
capacity=row[column_capacity]
|
||
|
|
resistance=row[column_resistance]
|
||
|
|
if not firstrow:
|
||
|
|
capacity=float(capacity)
|
||
|
|
resistance=float(resistance)
|
||
|
|
|
||
|
|
batteries.append([label,capacity,resistance])
|
||
|
|
batterylabels.append(label)
|
||
|
|
firstrow=False
|
||
|
|
|
||
|
|
print(str(len(batteries))+" Batteries found")
|
||
|
|
if len(batteries)==batteriesneeded:
|
||
|
|
print("You have just enough batteries")
|
||
|
|
elif len(batteries)>batteriesneeded:
|
||
|
|
print("You have "+str(len(batteries)-batteriesneeded)+" batteries spare")
|
||
|
|
elif len(batteries)<batteriesneeded:
|
||
|
|
print("You need "+str(batteriesneeded-len(batteries))+" more batteries!")
|
||
|
|
exit()
|
||
|
|
|
||
|
|
#batteries=[["a",1],["b",2],["c",3],["d",4],["e",5],["f",6]]
|
||
|
|
|
||
|
|
bestcombination=[]
|
||
|
|
bestdist=100000000
|
||
|
|
|
||
|
|
combinations=[]
|
||
|
|
|
||
|
|
|
||
|
|
count=0
|
||
|
|
|
||
|
|
for comb_parallel1 in itertools.combinations(batteries,config_parallel):
|
||
|
|
rest1=[x for x in batteries if x not in comb_parallel1]
|
||
|
|
count+=1
|
||
|
|
print(count)
|
||
|
|
for comb_parallel2 in itertools.combinations(rest1,config_parallel):
|
||
|
|
rest2=[x for x in rest1 if x not in comb_parallel2]
|
||
|
|
for comb_parallel3 in itertools.combinations(rest2,config_parallel):
|
||
|
|
rest3=[x for x in rest2 if x not in comb_parallel3]
|
||
|
|
for comb_parallel4 in itertools.combinations(rest3,config_parallel):
|
||
|
|
rest4=[x for x in rest3 if x not in comb_parallel4]
|
||
|
|
for comb_parallel5 in itertools.combinations(rest4,config_parallel):
|
||
|
|
rest5=[x for x in rest4 if x not in comb_parallel5]
|
||
|
|
for comb_parallel6 in itertools.combinations(rest5,config_parallel):
|
||
|
|
rest6=[x for x in rest5 if x not in comb_parallel6]
|
||
|
|
for comb_parallel7 in itertools.combinations(rest6,config_parallel):
|
||
|
|
rest7=[x for x in rest6 if x not in comb_parallel7]
|
||
|
|
for comb_parallel8 in itertools.combinations(rest7,config_parallel):
|
||
|
|
rest8=[x for x in rest7 if x not in comb_parallel8]
|
||
|
|
for comb_parallel9 in itertools.combinations(rest8,config_parallel):
|
||
|
|
rest9=[x for x in rest8 if x not in comb_parallel9]
|
||
|
|
for comb_parallel10 in itertools.combinations(rest9,config_parallel):
|
||
|
|
rest10=[x for x in rest9 if x not in comb_parallel10]
|
||
|
|
for comb_parallel11 in itertools.combinations(rest10,config_parallel):
|
||
|
|
rest11=[x for x in rest10 if x not in comb_parallel11]
|
||
|
|
for comb_parallel12 in itertools.combinations(rest11,config_parallel):
|
||
|
|
combinations.append([comb_parallel1,comb_parallel2,comb_parallel3,comb_parallel4,comb_parallel5,comb_parallel6,comb_parallel7,comb_parallel8,comb_parallel9,comb_parallel10,comb_parallel11,comb_parallel12])
|
||
|
|
|
||
|
|
|
||
|
|
|
||
|
|
#evaluate current combination of comb_pack1 and comb_pack2
|
||
|
|
print(str(len(combinations))+" Combinations")
|
||
|
|
for i in combinations: #every combinations
|
||
|
|
|
||
|
|
seriescapacities=[]
|
||
|
|
for s in i: #all series
|
||
|
|
parallelcapacity=0
|
||
|
|
for p in s: #all parallels
|
||
|
|
parallelcapacity+=p[1]
|
||
|
|
seriescapacities.append(parallelcapacity)
|
||
|
|
|
||
|
|
meancapacity=np.mean(seriescapacities)
|
||
|
|
dist=0
|
||
|
|
for c in seriescapacities:
|
||
|
|
dist+=math.pow(abs(c-meancapacity),2)
|
||
|
|
|
||
|
|
if dist<bestdist:
|
||
|
|
print(str(dist)+" "+str(seriescapacities))
|
||
|
|
bestdist=dist
|
||
|
|
bestcombination=i
|
||
|
|
|
||
|
|
print("best combination with "+str(bestdist)+" distance:")
|
||
|
|
#print(bestcombination)
|
||
|
|
for i in bestcombination:
|
||
|
|
print(i)
|