From 1b055c87170fad129990d04ecbc3cb08b4f7017e Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Anton=20Luka=20=C5=A0ijanec?= <anton@sijanec.eu>
Date: Fri, 31 Mar 2023 12:23:05 +0200
Subject: z: zotksd, fizvaje10, infsemin

---
 fiz/vaje/10/.gitignore   |  1 +
 fiz/vaje/10/dokument.lyx | 49 +++++++++++++++++++++++++++++++++---------------
 fiz/vaje/10/tab.py       | 22 +++++++++++++++++++---
 3 files changed, 54 insertions(+), 18 deletions(-)

(limited to 'fiz')

diff --git a/fiz/vaje/10/.gitignore b/fiz/vaje/10/.gitignore
index 8c9c1c2..da2a8f9 100644
--- a/fiz/vaje/10/.gitignore
+++ b/fiz/vaje/10/.gitignore
@@ -1,3 +1,4 @@
 *.csv
 p.tsv
 podatki.lyx
+graf.tex
diff --git a/fiz/vaje/10/dokument.lyx b/fiz/vaje/10/dokument.lyx
index 3209653..bf0df5c 100644
--- a/fiz/vaje/10/dokument.lyx
+++ b/fiz/vaje/10/dokument.lyx
@@ -124,10 +124,10 @@ Teoretični model
 \begin_inset Formula $s_{0}\left(\nu\right)$
 \end_inset
 
- je
+ harmoničnega nihala s tremi parametri je
 \begin_inset Formula 
 \[
-s_{0}=A\frac{\omega}{\sqrt{\left(\omega^{2}+\omega_{0}^{2}\right)^{2}+4\beta^{2}\omega^{2}}}\text{{,}}
+s_{0}=\frac{A}{\sqrt{\left(B^{2}-\omega^{2}\right)^{2}+C^{2}\omega^{2}}}\text{{,}}
 \]
 
 \end_inset
@@ -145,15 +145,6 @@ kjer je
 \end_inset
 
 .
- Resonančna frekvenca, razbrana iz grafa — vrh, je 
-\begin_inset Formula $6,5$
-\end_inset
-
- radianov na sekundo, kar je 
-\begin_inset Formula $6,5\pi2$
-\end_inset
-
- nihajev na sekundo.
 \end_layout
 
 \begin_layout Standard
@@ -975,7 +966,8 @@ status open
 \begin_inset Caption Standard
 
 \begin_layout Plain Layout
-Izmerjeni in izračunani podatki
+Izmerjeni in izračunani podatki.
+ Podatek o amplitudi je natančen na milimeter.
 \end_layout
 
 \end_inset
@@ -986,6 +978,10 @@ Izmerjeni in izračunani podatki
 \end_inset
 
 
+\begin_inset Note Note
+status open
+
+\begin_layout Plain Layout
 \begin_inset Float figure
 placement H
 wide false
@@ -1044,14 +1040,14 @@ addlegendentry{meritve}
 
 
 \backslash
-addplot[blue] (x, {0.490752637*x/sqrt((x^2+6.6^2)^2+4*(-0.000179148908)^2*x^2)});
+addplot[blue] (x, {0.47583543/sqrt(-6.87050654^2-x^2)^2+1.2365817^2*x^2});
 \end_layout
 
 \begin_layout Plain Layout
 
 
 \backslash
-addlegendentry{Teoretični model}
+addlegendentry{$A=0.47583543$, $B=-6.87050654$, $C=1.2365817$}
 \end_layout
 
 \begin_layout Plain Layout
@@ -1222,10 +1218,33 @@ name "fig:graf-1"
 \end_inset
 
 
+\end_layout
+
+\end_inset
+
+
 \end_layout
 
 \begin_layout Standard
-Meritve se ne ujemajo z napovedano vrednostjo.
+\begin_inset CommandInset include
+LatexCommand input
+filename "graf.tex"
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Standard
+Iz grafa za harmonično nihanje razberemo 6,819 kot 
+\begin_inset Formula $\omega_{0}$
+\end_inset
+
+, kar pomeni, da je resonančna frekvenca matematičnega nihala 
+\begin_inset Formula $\SI{1,085}{\per\second}$
+\end_inset
+
+.
 \end_layout
 
 \end_body
diff --git a/fiz/vaje/10/tab.py b/fiz/vaje/10/tab.py
index c5c24db..2b0bbae 100755
--- a/fiz/vaje/10/tab.py
+++ b/fiz/vaje/10/tab.py
@@ -1,7 +1,11 @@
 #!/usr/bin/python3
 import pandas
-import sys
+from sys import argv, stderr
 from math import pi
+from numpy import sqrt, arange, argmax
+from scipy.optimize import curve_fit
+from matplotlib import pyplot
+from tikzplotlib import get_tikz_code
 l = 23.5
 df = pandas.read_csv("podatki.tsv", sep="\t")
 df.index.name = "idx"
@@ -9,6 +13,18 @@ df["x"] = (df["y"]/13*l)
 df["omega"] = ((9.81/df["l"])**(1/2))
 df["f"] = (df["omega"]*0.5/pi)
 df = df.astype(float).round(6);
-data = df.values
+def objective(x, A, B, C):
+	return A/sqrt((B**2-x**2)**2+C**2*x**2)
+popt, _ = curve_fit(objective, df.values[:, 3], df.values[:, 2])
+print(popt, file=stderr)
+pyplot.scatter(df.values[:, 3], df.values[:, 2])
+iksi = arange(min(df.values[:, 3]), max(df.values[:, 3]), 0.01)
+ipsiloni = objective(iksi, popt[0], popt[1], popt[2])
+print(f"omega nič je {iksi[argmax(ipsiloni)]}, kar je {iksi[argmax(ipsiloni)]/2/pi} Hz", file=stderr)
+pyplot.plot(iksi, ipsiloni, "--", color="red")
+pyplot.xlabel('$s_0 \\left[\\si{\\meter}\\right]$')
+pyplot.ylabel('$\\omega=2\\pi\\nu$')
+pyplot.title("$s_0\\left(\\omega\\right)$")
 if __name__ == "__main__":
-	df.to_csv("/dev/stdout", sep=",");
+	# df.to_csv("/dev/stdout", sep=argv[1][0]);
+	print(get_tikz_code(axis_width="0.75\\textwidth"))
-- 
cgit v1.2.3