эффект Поккельса

.idea/.gitignore

@@ -0,0 +1,8 @@
+# Default ignored files
+/shelf/
+/workspace.xml
+# Datasource local storage ignored files
+/dataSources/
+/dataSources.local.xml
+# Editor-based HTTP Client requests
+/httpRequests/

.idea/inspectionProfiles/Project_Default.xml

@@ -0,0 +1,32 @@
+<component name="InspectionProjectProfileManager">
+  <profile version="1.0">
+    <option name="myName" value="Project Default" />
+    <inspection_tool class="DuplicatedCode" enabled="true" level="WEAK WARNING" enabled_by_default="true">
+      <Languages>
+        <language minSize="61" name="Python" />
+      </Languages>
+    </inspection_tool>
+    <inspection_tool class="Eslint" enabled="true" level="WARNING" enabled_by_default="true" />
+    <inspection_tool class="PyPep8Inspection" enabled="true" level="WEAK WARNING" enabled_by_default="true">
+      <option name="ignoredErrors">
+        <list>
+          <option value="E722" />
+        </list>
+      </option>
+    </inspection_tool>
+    <inspection_tool class="PyPep8NamingInspection" enabled="true" level="WEAK WARNING" enabled_by_default="true">
+      <option name="ignoredErrors">
+        <list>
+          <option value="N806" />
+        </list>
+      </option>
+    </inspection_tool>
+    <inspection_tool class="PyUnresolvedReferencesInspection" enabled="true" level="WARNING" enabled_by_default="true">
+      <option name="ignoredIdentifiers">
+        <list>
+          <option value="tuple.a" />
+        </list>
+      </option>
+    </inspection_tool>
+  </profile>
+</component>

.idea/inspectionProfiles/profiles_settings.xml

@@ -0,0 +1,6 @@
+<component name="InspectionProjectProfileManager">
+  <settings>
+    <option name="USE_PROJECT_PROFILE" value="false" />
+    <version value="1.0" />
+  </settings>
+</component>

.idea/misc.xml

@@ -0,0 +1,4 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<project version="4">
+  <component name="ProjectRootManager" version="2" project-jdk-name="Python 3.9" project-jdk-type="Python SDK" />
+</project>

.idea/modules.xml

@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<project version="4">
+  <component name="ProjectModuleManager">
+    <modules>
+      <module fileurl="file://$PROJECT_DIR$/.idea/Эффект Поккельса.iml" filepath="$PROJECT_DIR$/.idea/Эффект Поккельса.iml" />
+    </modules>
+  </component>
+</project>

.idea/Эффект Поккельса.iml

@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<module type="PYTHON_MODULE" version="4">
+  <component name="NewModuleRootManager">
+    <content url="file://$MODULE_DIR$" />
+    <orderEntry type="inheritedJdk" />
+    <orderEntry type="sourceFolder" forTests="false" />
+  </component>
+</module>

data.tex

@@ -0,0 +1,32 @@
+\chapter{Данные}
+\label{cha:data}
+
+\begin{table}[h!]\label{tab:res}	
+\centering	
+\begin{tabular}{|c|c|c|c|}
+\hline
+m & $r_m$, мм & $r_m^2$, мм$^2$ & $\Delta r_m^2$, мм$^2$\\
+\hline
+1 &9&81& 54\\
+\hline
+1 &9&81& 54\\
+\hline
+2&15&225&60\\
+\hline
+2&12&144&48\\
+\hline
+3&19&361&76\\
+\hline
+3&17&289&68\\
+\hline
+4&23&529&69\\
+\hline
+4&20.5&420&62\\
+\hline
+5&26.5&702&80\\
+\hline
+5&23.5&552&71\\
+\hline
+\end{tabular}
+\caption{Измерения тёмных колец интерференционной картины при поперечном выходном поляризаторе}
+\end{table}

main.py

@@ -0,0 +1,56 @@
+from matplotlib import pyplot as plt
+import numpy as np
+
+fig = plt.figure(figsize=(7, 4))
+
+plt.grid(linestyle='--')
+
+plt.xlabel('m', fontsize=15)
+plt.ylabel('r^2, mm^2', fontsize=15)
+
+r = [9, 15, 19, 23, 26.5]
+r2 = [9, 12, 17, 20.5, 23.5]
+r3 = []
+err = [3, 3, 2, 2, 2, 2, 1.5, 1.5, 1.5, 1.5]
+
+for i in range(5):
+    r3.append(r[i]**2)
+    r3.append(r2[i]**2)
+
+for i in range(10):
+    print(r3[i])
+
+err2 = []
+for i in range(10):
+    err2.append(2*np.sqrt(r3[i])*err[i])
+print(err2)
+
+m = [1, 1, 2, 2, 3, 3, 4, 4, 5, 5]
+
+plt.errorbar(m, r3, xerr=0, yerr=err2, fmt='o', linewidth=1)
+
+def MNK(m, r2):
+    #r2 = r**2
+    up = 0
+    down = 0
+    up2 = 0
+    down2 = 0
+    for i in range(10):
+        up+=m[i]*r2[i]
+        down+=m[i]**2
+        up2+=r2[i]**2
+        down2+=m[i]**2
+    k=up/down
+    print(k, 1/np.sqrt(10)*np.sqrt(up2/down2-k**2))
+
+
+MNK(m, r3)
+x = np.linspace(1, 5, 10)
+y = 117.45*x
+
+plt.plot(x, y, linewidth=2, label='k = 117$\pm 5$')
+
+
+plt.title('МНК приближение')
+plt.legend(loc='upper center', fontsize=15)
+plt.show()