diff options
| -rw-r--r-- | .gitignore | 1 | ||||
| -rw-r--r-- | 117/prak117.data | 18 | ||||
| -rw-r--r-- | 117/prak117.gp | 31 | ||||
| -rw-r--r-- | 117/prak117.py | 47 |
4 files changed, 81 insertions, 16 deletions
@@ -2,3 +2,4 @@ output* venv/ __pycache__/ *.png +*.log diff --git a/117/prak117.data b/117/prak117.data index c258fcd..ade21ae 100644 --- a/117/prak117.data +++ b/117/prak117.data @@ -22,6 +22,22 @@ 450 21.063 500 25.883 q -17.222 17.205 0.1 +5, 15.475, 16.429 +15, 15.667, 16.511 +25, 15.795, 16.547 +35, 15.956, 16.639 +45, 16.120, 16.727 +55, 16.372, 16.814 +65, 16.508, 16.895 +75, 16.695, 16.986 +85, 16.906, 17.080 +95, 17.131, 17.224 +97.5, 17.173, 17.194 +98, 17.222, 17.205 +98.5, 17.283, 17.207 +100, 17.262, 17.218 +105, 17.354, 17.243 +115, 17.607, 17.356 +q 0.735 0.001 diff --git a/117/prak117.gp b/117/prak117.gp new file mode 100644 index 0000000..7896f31 --- /dev/null +++ b/117/prak117.gp @@ -0,0 +1,31 @@ +set term png +set output 'output_117_T_a.png' +set xlabel "a(m)" +set ylabel "T(c)" +set grid +plot 'output_117_T_a.data' using 1:2 with lines title "T(a)", \ + 'output_117_T_a.data' linecolor rgb '#FF0000' pointtype 5 notitle + +set output 'output_117_2.png' +set xlabel "d(mm)" +set ylabel "T(c)" +f(x) = a*x+b +g(x) = c*x+d +fit f(x) 'output_117_T1.data' using 1:2:3 via a,b +fit g(x) 'output_117_T2.data' using 1:2:3 via c,d +plot f(x) title "T_{прям.}(d)",\ + 'output_117_T1.data' using 1:2:3 pointtype 6 notitle with yerrorbars , \ + g(x) title "Т_{обр.}(d)", \ + 'output_117_T2.data' using 1:2:3 pointtype 6 notitle with yerrorbars + +set output 'output_117_2_zoomed.png' +set xrange[90:110] +replot + +unset xrange +set xlabel "T^2⋅a, c^2⋅m" +set ylabel "a^2, m^2" +set output 'output_117_mnk.png' +h(x) = m*x+l +fit h(x) 'output_117_mnk.data' using 1:2:3 via m,l +plot h(x) title "y(x)", 'output_117_mnk.data' using 1:2:3 pointtype 6 notitle with yerrorbars 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() |