52 lines
1.4 KiB
Python
52 lines
1.4 KiB
Python
|
import sys
|
|||
|
|
import math
|
|||
|
|
|
|||
|
|
def avg(a):
|
|||
|
|
return sum(a) / len(a)
|
|||
|
|
|
|||
|
|
resistances = []
|
|||
|
|
|
|||
|
|
settings = {'r': 4, "R2": 0.0, "R3": 0.0}
|
|||
|
|
|
|||
|
|
for raw_i in open(sys.argv[1], "r"):
|
|||
|
|
i = raw_i.rstrip("\r\n")
|
|||
|
|
|
|||
|
|
if len(i) > 0:
|
|||
|
|
if i[0] == "#":
|
|||
|
|
cmd = i[1:].split()
|
|||
|
|
if cmd[0] in settings:
|
|||
|
|
settings[cmd[0]] = cmd[1]
|
|||
|
|
else:
|
|||
|
|
resistances.append([i.split()[0]] + [float(j.replace(',', '.')) for j in i.split()[1:]])
|
|||
|
|
|
|||
|
|
table_data = []
|
|||
|
|
|
|||
|
|
for i in resistances:
|
|||
|
|
Rxi = [round(r*float(settings["R2"])/float(settings["R3"]), int(settings['r'])) for r in i[1:]]
|
|||
|
|
R_avg = round(avg(Rxi), int(settings['r']))
|
|||
|
|
|
|||
|
|
table_data.append([i[0]] + Rxi + [R_avg])
|
|||
|
|
|
|||
|
|
from prettytable import PrettyTable
|
|||
|
|
|
|||
|
|
pt = PrettyTable()
|
|||
|
|
|
|||
|
|
pt.add_column("", list(range(1, len(table_data[0][1:-1]) + 1)) + ["❬Rx❭"])
|
|||
|
|
|
|||
|
|
for i in table_data:
|
|||
|
|
pt.add_column(i[0], i[1:])
|
|||
|
|
|
|||
|
|
print(pt)
|
|||
|
|
|
|||
|
|
#print(table_data)
|
|||
|
|
|
|||
|
|
print("Обчислюємо значення S для кожного стовпця:")
|
|||
|
|
for j in table_data:
|
|||
|
|
s = math.sqrt( sum([ (i - j[-1])**2 for i in j[1:-1] ]) / (len(j[1:-1])*(len(j[1:-1])-1)) )
|
|||
|
|
print(f"S ❬{j[0]}❭ = {round(s, int(settings['r']))}")
|
|||
|
|
|
|||
|
|
print("\nОбчислюємо значення σ для кожного стовпця:")
|
|||
|
|
for j in table_data:
|
|||
|
|
sigma = math.sqrt(avg([i**2 for i in j[1:-1]]) - avg(j[1:-1])**2)
|
|||
|
|
print(f"σ ❬{j[0]}❭ = {round(sigma, int(settings['r']))}")
|