Update interpolate script

ors-tx-interpolation/interpolate.py

@@ -37,37 +37,19 @@ def fit_xs_ys_spline(xs, ys, s):
   plt.plot(x, p(x))
   plt.show()
 
-b42_1w_data = [
-  (60,  3.2,  1.92),
-  (70,  13.7, 2.03),
-  (75,  18.9, 2.21),
-  (80,  24.8, 2.44),
-  (81,  26.3, 2.96),
-  (82,  27.6, 3.85),
-  (83,  28.3, 4.70),
-  (84,  29.3, 7.30),
-]
-b43_1w_data = [
-  (60, 3, 1.9),
-  (70, 13, 1.4),
-  (80, 23, 2.6),
-  (81, 24, 3.2),
-  (83, 26, 6.0),
-  (85, 27, 9.5),
-]
-b39_1w_data = [
-  (60, 1.3, 1.53),
-  (70, 11.3, 1.07),
-  (80, 20.4, 1.42),
-  (81, 21.2, 1.45),
-  (82, 22.1, 1.46),
-  (83, 23.0, 1.45),
-  (84, 24.0, 1.45),
-  (85, 24.9, 1.53),
-  (86, 25.9, 1.82),
-  (87, 26.7, 2.39),
-  (88, 27.5, 3.13),
-  (89, 28.4, 4.45),
+b28_1w_data = [
+  (60, 11.2, 1.70),
+  (70, 21.1, 1.72),
+  (71, 21.7, 1.73),
+  (72, 23.0, 1.72),
+  (73, 23.9, 1.74),
+  (74, 24.7, 1.82),
+  (75, 25.6, 1.92),
+  (76, 26.5, 2.21),
+  (77, 27.4, 2.88),
+  (78, 28.4, 3.94),
+  (79, 29.0, 5.28),
+  (80, 29.9, 7.31),
 ]
 b38_1w_data = [
   (60, 0, 1.81),
@@ -87,16 +69,73 @@ b38_1w_data = [
   (88, 27.0, 1.89),
   (89, 28.1, 3.00),
 ]
-b43_500mw_data = [
-  (60, 2.5, 1.77),
-  (70, 13.7, 1.82),
-  (73, 16.6, 1.97),
-  (74, 17.6, 2.07),
-  (75, 18.5, 2.51),
-  (76, 19.4, 3.77),
-  (77, 20.2, 6.17),
-  (78, 20.9, 8.81),
-  (79, 21.4, 10.6),
+b39_1w_data = [
+  (60, 1.3, 1.53),
+  (70, 11.3, 1.07),
+  (80, 20.4, 1.42),
+  (81, 21.2, 1.45),
+  (82, 22.1, 1.46),
+  (83, 23.0, 1.45),
+  (84, 24.0, 1.45),
+  (85, 24.9, 1.53),
+  (86, 25.9, 1.82),
+  (87, 26.7, 2.39),
+  (88, 27.5, 3.13),
+  (89, 28.4, 4.45),
+]
+b40_1w_data = [
+  (60, 1.5, 2.72),
+  (70, 10.7, 2.46),
+  (80, 20.2, 2.80),
+  (81, 20.8, 2.74),
+  (82, 22.0, 2.91),
+  (83, 22.5, 2.99),
+  (84, 23.6, 3.21),
+  (85, 24.4, 3.55),
+  (86, 25.3, 4.55),
+  (87, 26.0, 5.29),
+  (88, 27.1, 6.56),
+]
+b42_1w_data = [
+  (60,  3.2,  1.92),
+  (70,  13.7, 2.03),
+  (75,  18.9, 2.21),
+  (80,  24.8, 2.44),
+  (81,  26.3, 2.96),
+  (82,  27.6, 3.85),
+  (83,  28.3, 4.70),
+  (84,  29.3, 7.30),
+]
+b43_1w_data = [
+  (60, 3, 1.9),
+  (70, 13, 1.4),
+  (80, 23, 2.6),
+  (81, 24, 3.2),
+  (83, 26, 6.0),
+  (85, 27, 9.5),
+]
+n77_1w_data = [
+  (60, 5.0, 1.75),
+  (70, 14.9, 1.64),
+  (75, 19.8, 1.79),
+  (80, 24.8, 1.90),
+  (81, 25.7, 2.04),
+  (82, 26.5, 2.42),
+  (83, 27.3, 3.14),
+  (84, 28.2, 4.39),
+  (85, 29.0, 5.91),
+  (86, 29.8, 8.01),
+]
+n79_1w_data = [
+  (60, 6.1, 2.24),
+  (70, 16.2, 2.08),
+  (75, 21.2, 2.09),
+  (80, 26.2, 2.22),
+  (81, 27.0, 2.44),
+  (82, 27.8, 2.87),
+  (83, 28.6, 3.62),
+  (84, 29.4, 4.90),
+  (85, 30.0, 6.40),
 ]
 b38_500mw_data = [
   (60, 6.3, 1.29),
@@ -111,10 +150,25 @@ b38_500mw_data = [
   (80, 25.7, 7.87),
   (81, 26.2, 9.58),
 ]
+b43_500mw_data = [
+  (60, 2.5, 1.77),
+  (70, 13.7, 1.82),
+  (73, 16.6, 1.97),
+  (74, 17.6, 2.07),
+  (75, 18.5, 2.51),
+  (76, 19.4, 3.77),
+  (77, 20.2, 6.17),
+  (78, 20.9, 8.81),
+  (79, 21.4, 10.6),
+]
 
-fit_data(b42_1w_data, 2, "B42 2x1W")
-fit_data(b43_1w_data, 2, "B43 2x1W")
-fit_data(b39_1w_data, 2, "B39 2x1W")
-fit_data(b38_1w_data, 2, "B38 2x1W")
-fit_data(b43_500mw_data, 2, "B43 2x0.5W")
+fit_data(b28_1w_data,    2, "B28 2x1W")
+fit_data(b38_1w_data,    2, "B38 2x1W")
+fit_data(b39_1w_data,    2, "B39 2x1W")
+fit_data(b40_1w_data,    2, "B40 2x1W")
+fit_data(b42_1w_data,    2, "B42 2x1W")
+fit_data(b43_1w_data,    2, "B43 2x1W")
+fit_data(n77_1w_data,    2, "N77 2x1W")
+fit_data(n79_1w_data,    2, "N79 2x1W")
 fit_data(b38_500mw_data, 2, "B38 2x0.5W")
+fit_data(b43_500mw_data, 2, "B43 2x0.5W")