import math

V_1 = [235,300,325,350,375,400,450,500,250,200]
V_2 = [200,225,250,275,300,325,350,375,400,450]
V_3 = [300,350,400,450,500,550,600,650,700,525]

I_1 = [118.05, 149.66, 162.01, 175.5, 187.15, 199.4, 225.25, 250, 124.5,100.06]

I_2 = [67.55,75.81,84.45,92.55,100.89,109.46,118.01,126.06,134.75,151.32]
I_3 = [60.55,70.85,81.06,91.00,101.37,111.66,121.09,131.76,141.19,105.94]

V_max_1 = 500
I_max_1 = 250
R_max_1 = 2.0009

V_max_2 = 450
I_max_2 = 151.32
R_max_2 = 2.9806

V_max_3 = 700
I_max_3 = 141.19
R_max_3 = 4.9509

R_bridge_1 = 2.0326
R_bridge_2 = 3.0245
R_bridge_3 = 5.0645

R_V = 250
sig_V = 2.5 / 2
sig_I = 2.4
sig_l = 0.005
sig_d = 0.00001
s = 0.00105
d = 0.037

VI_sum_1 = 0
VI_sum_2 = 0
VI_sum_3 = 0

II_sum_1 = 0
II_sum_2 = 0
II_sum_3 = 0

RR_sum_1 = 0
RR_sum_2 = 0
RR_sum_3 = 0


for i in range (10):
  VI_sum_1 += V_1[i]*I_1[i]
  VI_sum_2 += V_2[i]*I_2[i]
  VI_sum_3 += V_3[i]*I_3[i]

  II_sum_1 += I_1[i]*I_1[i]
  II_sum_2 += I_2[i]*I_2[i]
  II_sum_3 += I_3[i]*I_3[i]

  RR_sum_1 += (V_1[i]*V_1[i]) / (I_1[i]*I_1[i])
  RR_sum_2 += (V_2[i]*V_2[i]) / (I_2[i]*I_2[i])
  RR_sum_3 += (V_3[i]*V_3[i]) / (I_3[i]*I_3[i])

VI_1 = VI_sum_1 / 10
VI_2 = VI_sum_2 / 10
VI_3 = VI_sum_3 / 10

II_1 = II_sum_1 / 10
II_2 = II_sum_2 / 10
II_3 = II_sum_3 / 10

RR_1 = RR_sum_1 / 10
RR_2 = RR_sum_2 / 10
RR_3 = RR_sum_3 / 10

R_1 = VI_1 / II_1
R_2 = VI_2 / II_2
R_3 = VI_3 / II_3

R_pract_1 = R_1 + (R_1*R_1 / R_V)
R_pract_2 = R_2 + (R_2*R_2 / R_V)
R_pract_3 = R_3 + (R_3*R_3 / R_V)

sig_rand_1 = (1/math.sqrt(10)) * math.sqrt(abs(RR_1 - (R_1*R_1)))
sig_rand_2 = (1/math.sqrt(10)) * math.sqrt(abs(RR_2 - (R_2*R_2)))
sig_rand_3 = (1/math.sqrt(10)) * math.sqrt(abs(RR_3 - (R_3*R_3)))

sig_sys_1 = R_1 * math.sqrt(((sig_V/V_max_1)*(sig_V/V_max_1))+((sig_I/I_max_1)*(sig_I/I_max_1)))
sig_sys_2 = R_2 * math.sqrt(((sig_V/V_max_2)*(sig_V/V_max_2))+((sig_I/I_max_2)*(sig_I/I_max_2)))
sig_sys_3 = R_3 * math.sqrt(((sig_V/V_max_3)*(sig_V/V_max_3))+((sig_I/I_max_3)*(sig_I/I_max_3)))

sig_R_1 = math.sqrt((sig_sys_1*sig_sys_1) + (sig_rand_1)*(sig_rand_1))
sig_R_2 = math.sqrt((sig_sys_2*sig_sys_2) + (sig_rand_2)*(sig_rand_2))
sig_R_3 = math.sqrt((sig_sys_3*sig_sys_3) + (sig_rand_3)*(sig_rand_3))

po_1 = (R_pract_1 * s) / 20
po_2 = (R_pract_2 * s) / 30
po_3 = (R_pract_3 * s) / 50

sig_po_1 = po_1*math.sqrt(((sig_R_1/R_max_1)*(sig_R_1/R_max_1))+2*((sig_d/d)*(sig_d/d))+((sig_l/20)*(sig_l/20)))

print("R по мосту P4833 (для l1)",R_bridge_1)
print("R по мосту P4833 (для l2)",R_bridge_2)
print("R по мосту P4833 (для l3)",R_bridge_3)

print("--------------------------------------")

print("R среднее (для l1)  = ",R_1)
print("R среднее (для l2)  = ",R_2)
print("R среднее (для l3)  = ",R_3)

print("--------------------------------------")

print("R проволоки (для l1)  = ",R_pract_1)
print("R проволоки (для l2)  = ",R_pract_2)
print("R проволоки (для l3)  = ",R_pract_3)

print("--------------------------------------")

print("Сигма случайная (для l1)",sig_rand_1)
print("Сигма случайная (для l2)",sig_rand_2)
print("Сигма случайная (для l3)",sig_rand_3)

print("--------------------------------------")

print("Сигма системная (для l1)",sig_sys_1)
print("Сигма системная (для l2)",sig_sys_1)
print("Сигма системная (для l3)",sig_sys_3)

print("--------------------------------------")

print("Сигма R (для l1)",sig_R_1)
print("Сигма R (для l2)",sig_R_2)
print("Сигма R (для l3)",sig_R_3)

print("--------------------------------------")

print("Ро (для l1)",po_1)
print("Ро (для l2)",po_2)
print("Ро (для l3)",po_3)

print("--------------------------------------")

print(sig_po_1)