Я проверил доступный метод интерполяции в scipy, но не смог получить правильное решение для моего дела. предположим, что у меня есть 100 точек, координаты которых являются случайными, например, их положения x и y: x=np.random.rand(100)*100 y=np.random.rand(100)*100 z = f(x,y) #the point value calculated by certain function теперь я хочу получить точечное значение z новых равномерно отобранных […]
Я пытаюсь вычислить инверсию кумулятивного дважды нецентрального распределения f в python. Scipy имеет нецентральный f-дистрибутив (scipy.stats.ncf), но документация довольно бедна, и не кажется, что nc, параметр нецентральности, может быть списком (передать nc для числителя и знаменатель). Я знаю, что числитель и знаменатель – это нецентральные чи-квадратные распределения ( http://mathworld.wolfram.com/NoncentralF-Distribution.html ), но я не уверен, как […]
Я пытаюсь создать функцию calc(f,a,b) где x – это уравнение с переменной f, и я хочу поместить этот код в функцию. def calc(f, a, b): limits = [a, b] integral = odeint(lambda y, x : f, 0, limits) return integral[1] Эта функция получает интеграл, используя встроенную функцию odeint. Это то, что я пытаюсь сделать print […]
У Scipy есть превосходный пакет пространственного анализа, который включает K-мерное дерево. Я пытаюсь использовать функцию запроса, и она возвращает эту ошибку: ValueError: x должен состоять из векторов длины 6, но имеет форму (2,) Кто-нибудь знает, что эта ошибка имеет в виду? Из некоторых поисковых запросов Google я обнаружил, что он имеет этот общий формат: raise […]
Я использую griddata для повторной выборки 2-мерного массива на сетке. z.shape = (1000, 1000) x, y = np.arange(-5, 5, 0.01), np.arange(-5, 5, 0.01) newx, newy = np.arange(-2, 2, 0.1), np.arange(-2, 2, 0.1) griddata((x, y), z, (newx[None, :], newy[:, None])) Код должен: resample z (который представляет изображение) к новой более крупной или тонкой сетке новая сетка […]
Я использую Scipy SmoothBivariateSpline для создания кубического B-сплайна на двумерных данных. Теперь мне нужно написать кусочно-полиномиальное выражение для этой сплайновой кривой. Мой математический фон не очень силен, поэтому я не смог написать свой собственный алгоритм для преобразования из вывода t, c, k SmoothBivariateSpline в полиномиальное представление. Если это возможно, можете ли вы указать указатели на […]
У меня есть массив A экспортированный из excel, содержащий значения данных, как показано. 1-й столбец x и второй столбец y являются зависимыми переменными, а третий столбец z – независимыми переменными (выход). from xlrd import open_workbook Data = open_workbook("simple.xls") sheet = Data.sheet_by_name('Sheet1') A=[] # Read row by row for rownum in range(sheet.nrows): rowValues = sheet.row_values(rownum) A.append(rowValues) […]
Поскольку я понимаю LU-факторизацию, это означает, что матрицу A можно записать как A = LU для нижней треугольной матрицы L и верхнетреугольной матрицы U. Однако функции в scipy, относящиеся к LU-разложениям ( lu , lu_factor , lu_solve ), по-видимому, включают в себя третью матрицу P, такую, что A = PLU и P – матрица перестановок […]
Это очень простой вопрос, но он меня озадачивает. Я пытаюсь встроить некоторые scipy-процедуры в c-программу. Тем не менее, я не могу успешно завершить начальный шаг импорта любых scipy модулей. Я могу импортировать верхний уровень scipy, не получая нулевое возвращаемое значение, поэтому я уверен, что установка не проблема … PyObject *pckg_name, *pckg; pckg_name = PyString_FromString("scipy"); pckg […]
import numpy as np from scipy import signal data = np.array([[[13, 2, 1, np.nan, np.nan], [22, 1, 1, 4, 4], [4, 2, 3, 3, 4], [1, 1, 4, 1, 5], [2, 4, 5, 2, 1]], [[17, 7, 10, 6, np.nan], [np.nan, 7, 8, 6, 9], [6, 10, 9, 8, 10], [6, 8, 7, 10, 8], […]