From dfa85cca6464637aa6ad943da3f9ce8ccdfeca1f Mon Sep 17 00:00:00 2001
From: justanothercatgirl <sotov2070@gmail.com>
Date: Mon, 25 Nov 2024 19:19:25 +0300
Subject: cleanup

---
 101/plots.gp | 72 ++++++++++++++++++++++++++++++++++++++++++++++++++++++------
 1 file changed, 65 insertions(+), 7 deletions(-)

(limited to '101/plots.gp')

diff --git a/101/plots.gp b/101/plots.gp
index 702408a..ba43825 100644
--- a/101/plots.gp
+++ b/101/plots.gp
@@ -1,14 +1,72 @@
 set term pngcairo size 1000, 800
 set tmargin at screen 0.95
+set grid lw 2
 
 f1(x) = a1*x+b1
-fit f1(x) '.plot' using 1:2:3 yerr via a1, b1
+f2(x) = a2 + b2/(c2*x + d2)
+f3(x) = a3 * exp(b3*x - c3) + d3
 
-set output ''
-set label "" at graph 0.5, graph 1.025 center
-set xlabel ""
-set ylabel ""
+set label "График зависимости ускорений а от количества грузов N" at graph 0.5, graph 1.025 center
+set xlabel "N, ед."
+set ylabel "a, м/c^2"
+set xrange[0:11]
+set yrange[0:1.5]
+set output "ex1.png"
+fit f2(x) 'ex1_exp.plot' using 1:2:3 yerr via a2, b2, c2, d2
+set label "график a_{эксп} построен по мнк 'a + b/(cx + d)'" at 6, 0.18
+plot 	'ex1_calc.plot' using 1:2:3 with yerrorbars notitle lc 0 pt 1 lw 2,	\
+	'ex1_calc.plot' with lines lc rgb "red" title "a_{расч}(N)",		\
+	'ex1_exp.plot' using 1:2:3 with yerrorbars notitle lc 0 pt 1 lw 2,	\
+	f2(x) lc rgb "green" title "a_{эксп}(N)"
 
-plot 	'.plot' using 1:2:3 with yerrorbars notitle lc 0 pt 1 lw 2,	\
-	f1(x) title "" lc rgb "red",				\
+unset label
+set label "график зависимости a_{расч}, a_{эксп}, отклонения a от 1/m" at graph 0.5, graph 1.025 center
+set ytics nomirror
+set y2tics 
+unset xrange
+set yrange[0.05:0.18]
+set y2range[0:100]
+set y2label "delta a, %"
+set xlabel "1/M, кг^{-1}"
+set ylabel "a, м/c^2"
+fit f1(x) 'ex2_exp.plot' using 1:2 via a1, b1
+set output "ex2.png"
+plot 	'ex2_cal.plot'	axes x1y1 notitle lc 0 pt 1 lw 2,				\
+	'ex2_cal.plot'	with lines lc rgb "red" axes x1y1 title "a_{расч}, м/с^{2}", 	\
+	'ex2_exp.plot'	axes x1y1 notitle lc 0 pt 1 lw 2,				\
+	f1(x)		lc rgb "green" axes x1y1 title "a_{эксп}, м/с^{2}",		\
+	'ex2_dev.plot'	axes x1y2 notitle lc 0 pt 1 lw 2,				\
+	'ex2_dev.plot'	with lines lc rgb "blue" axes x1y2 title "delta a, %"
+
+unset label
+set label "график зависимости a_{расч}, a_{эксп}, отклонения a от delta m" at graph 0.5, graph 1.025 center
+unset xrange
+set yrange[0:2.8]
+set y2range[0:100]
+set y2label "delta a, %"
+set xlabel "delta m, кг"
+set ylabel "a, м/c^2"
+fit f1(x) 'ex3_exp.plot' using 1:2 via a1, b1
+fit f3(x) 'ex3_dev.plot' using 1:2 via a3, b3, c3, d3
+set output "ex3.png"
+plot 	'ex3_cal.plot'	axes x1y1 notitle lc 0 pt 1 lw 2,				\
+	'ex3_cal.plot'	with lines lc rgb "red" axes x1y1 title "a_{расч}, м/с^{2}", 	\
+	'ex3_exp.plot'	axes x1y1 notitle lc 0 pt 1 lw 2,				\
+	f1(x)		lc rgb "green" axes x1y1 title "a_{эксп}, м/с^{2}",		\
+	'ex3_dev.plot'	axes x1y2 notitle lc 0 pt 1 lw 2,				\
+	'ex3_dev.plot'	with lines lc rgb "blue" axes x1y2 title "delta a, %",		\
+	f3(x) axes x1y2 lt 2 lw 1 dashtype "--" notitle
+
+set output "ex3_mnk.png"
+unset label
+set label "график зависимости a_{расч}, a_{эксп}, отклонения a от delta m для МНК" at graph 0.5, graph 1.025 center
+set label sprintf("C = %.4f", a1) at 0.0255, 0.9
+set label sprintf("D = %.4f", b1) at 0.0255, 0.83
+unset xrange
+unset yrange
+unset y2tics
+unset y2label 
+fit f1(x) 'ex3_exp.plot' via a1, b1
+plot 	'ex3_exp.plot' notitle lc rgb "red" pt 1 lw 2,				\
+	f1(x) with lines lc 0 title "a(dm)"
 
-- 
cgit v1.2.3-70-g09d2