Различия

Показаны различия между двумя версиями страницы.

--- calc:math [2021/03/18 14:36]
root
+++ calc:math [2023/03/14 15:39] (текущий)
root
@@ Строка 212: / Строка 212: @@
 <sxh python>
 from scipy.interpolate import RegularGridInterpolator
+import numpy as np
+import matplotlib.pyplot as plt
 def func(x, y):
@@ Строка 227: / Строка 229: @@
 #Нерегулярная сетка
-xs = (5 - (-5)) * np.random.random_sample(20) + (-5)
+xs = np.random.uniform(-5, 5, 50)
-ys = (2 - (-2)) * np.random.random_sample(20) + (-2)
+ys = np.random.uniform(-2, 2, 50)
 #Значение в точках
 values = func(xs, ys)
+# Линейная интерполяция на нерегулярной сетке
+baz = scipy.interpolate.LinearNDInterpolator(np.rollaxis(np.array([xs,ys]), 1), values)
+plt.imshow(baz((grid_x, grid_y)), extent=[-5,5,-2,2], origin='lower')
 # Значения на регулярной сетке
 data = scipy.interpolate.griddata(np.rollaxis(np.array([xs,ys]), 1), values, (grid_x, grid_y))
 # Интерполяция
 foo = RegularGridInterpolator((xi,yi),data)
-func(-0.5,1) # 1.0
+# точки
+plt.scatter(xs, ys)
+# исходное поле
+plt.imshow(func(grid_x, grid_y), extent=[-5,5,-2,2], origin='lower')
+# полученное поле на регулярной сетке
+plt.imshow(data, extent=[-5,5,-2,2], origin='lower')
-foo([-0.5,1])
+# результат интерполяции
+plt.imshow(foo((grid_x, grid_y)), extent=[-5,5,-2,2], origin='lower')
 </sxh>