From c8548265806e368433eee27cde551a6ee13fd204 Mon Sep 17 00:00:00 2001
From: justanothercatgirl <sotov@twistea.su>
Date: Sat, 26 Oct 2024 22:28:42 +0300
Subject: Finished 117

---
 117/prak117.py | 47 ++++++++++++++++++++++++++++++++---------------
 1 file changed, 32 insertions(+), 15 deletions(-)

(limited to '117/prak117.py')

diff --git a/117/prak117.py b/117/prak117.py
index 68b5d6f..52d60c7 100644
--- a/117/prak117.py
+++ b/117/prak117.py
@@ -32,17 +32,15 @@ def task1():
         a = a_i(L[0], d)
         T = t/n
         data.append([d, t, a, T, T*T*a, a*a, sigma(a_i_sq, [L[0], d], [L[1], 0.001])])
-    print("%-16s" * 7 % ("d(mm)", "t(c)", "a(mm)", "T(c)", "x=T^2 (c^2*mm)", "y=a^2 (mm^2)", "S_y (mm^2)"))
+    print("%-16s" * 7 % ("d(cm)", "t(c)", "a(cm)", "T(c)", "x=T^2 (m*c^2)", "y=a^2 (m^2)", "S_y (cm^2) * 10^-3"))
     for i in data:
         l = i.copy()
-        l[I_D] *= 1000 # d
-        l[I_A] *= 1000 # a
-        l[I_X] *= 1000 # T^2 * a
-        l[I_Y] *= 1000000 # a^2
-        l[I_SY] *= 1000000 # S_{a^2}
+        l[I_D] *= 100 # d
+        l[I_A] *= 100 # a
+        l[I_SY] *= 1000 # S_{a^2}
         print("%-16.4f"  * 7 % tuple(l))
     with open("output_117_T_a.data", "w") as file:
-        for a, b in zip([i[I_T] for i in data], [i[I_A] for i in data]):
+        for a, b in zip([i[I_A] for i in data], [i[I_T] for i in data]):
             file.write(f"{a} {b}\n")
     X = np.array([i[I_X] for i in data])
     Y = np.array([i[I_Y] for i in data])
@@ -75,18 +73,37 @@ def g(l, T): return 4 * PI * PI * l / T / T
 
 def task2():
     n = 10
-    print('-' * 15, "УПРАЖНЕНИЕ 1", '-' * 15)
-    print("Так, ну поделить на 10 вручную легче, чем вводить это в комп а потом переписывать")
-    print("так что периоды посчитаете сами")
-    print("Введите t_01 и t_02 на одной строке через пробел")
-    T0 = avg(list(map(lambda x: float(x)/n, input().split())))
-    T_sigma = float(input("Введите погрешность одного измерения периода (sigma_суб.):\n"))
+    print('-' * 15, "УПРАЖНЕНИЕ 2", '-' * 15)
+    T_sigma = float(input("Введите погрешность одного измерения периода (sigma_суб.):\n")) / sqrt(6)
+    print("Так же, как в упражнении 1 вводите 'x, t_прямой, t_обратный' ")
+    print("чтобы посчитать среднее, пишите 'x, t1 t2 t3, t4 t5 t6' ")
+    print("в отличие от упражнения один, разделять запятыми обязательно")
+    xs = []; t_dir = []; t_rev = []; terr_dir = []; terr_rev = [];
+    while (inp := input()) != 'q':
+        inps = inp.split(',')
+        xs.append(float(inps[0]))
+        t_dir.append(avg(list(map(float, inps[1].split()))) / 10)
+        terr_dir.append(T_sigma / 10 / sqrt(len(inps[1].split())))
+        t_rev.append(avg(list(map(float, inps[2].split()))) / 10)
+        terr_rev.append(T_sigma / 10 / sqrt(len(inps[2].split())))
+    smallind = 0
+    smalldif = 0x7FFFFFFF
+    for i, (d, r) in enumerate(zip(t_dir, t_rev)):
+        if abs((r-d)) < smalldif:
+            smalldif = abs(d-r)
+            smallind = i
+    T0 = avg((t_dir[smallind], t_rev[smallind]))
     print("Введите L_пр и погрешность на одной строке (м)")
     l_pr = [*map(float, input().split())]
+    T_s = sqrt(terr_rev[smallind]**2 / 2 + terr_dir[smallind]**2 / 2)
     args = [l_pr[0], T0]
-    sigmas = [l_pr[1], T_sigma / sqrt(2) / sqrt(3)] # корни от 2 расчётов среднего
+    sigmas = [l_pr[1], T_s]
     g_ = (g(*args), sigma(g, args, sigmas))
-    print(f"T0 = {T0}±{sigmas[1]}\ng = {g_[0]}±{g_[1]}")
+    with open("output_117_T1.data", "w") as f1, open("output_117_T2.data", "w") as f2:
+        for x, t1, t2, e1, e2 in zip(xs, t_dir, t_rev, terr_dir, terr_rev):
+            f1.write(f"{x} {t1} {e1}\n")
+            f2.write(f"{x} {t2} {e2}\n")
+    print(f"T0[x={xs[smallind]}] = {T0}±{sigmas[1]}\ng = {g_[0]}±{g_[1]}")
 
 if __name__ == '__main__':
     task1()
-- 
cgit v1.2.3-70-g09d2