Tag: scipy

Интеграция с использованием квадратурных

У меня возникло несколько вопросов, использующих квадратурную функцию в python 2.7 (часть модуля scipy.integrate). Уравнение, которое я пытаюсь интегрировать, просто: x/(d^2) – (x^2) Интеграция между пределами a и b. Однако мне нужно выполнить интеграцию при 40 различных значениях d и не уверен, как передать второй аргумент, чтобы зацикливать интеграцию по значениям d. Любая помощь будет […]

scipy curve_fit терпит неудачу при экспоненциальном подгонке

Когда я пытаюсь выполнить экспоненциальную подгонку с использованием curve_fit, scipy возвращает ошибку. Я делаю что-то неправильно? Удаление отрицательного знака из np.exp (-b * t) позволяет функции curve_fit работать, но возвращаемые значения не работают. #!/usr/bin/python import numpy as np import scipy as sp from scipy.optimize import curve_fit import scipy.optimize as opt import matplotlib.pyplot as plt x […]

Включает ли функция минимизации Scipy метод «COBYLA»?

Я использую алгоритм 'COBYLA' в функции optimize.minimize scipy .minimize (v.0.11 build для cygwin). Я заметил, что в этом случае bounds параметров не используются. Например, простой пример: from scipy.optimize import minimize def f(x): return -sum(x) minimize(f, x0=1, method='COBYLA', bounds=(-2,2)) возвращает: status: 2.0 nfev: 1000 maxcv: 0.0 success: False fun: -1000.0 x: array(1000.0) message: 'Maximum number of […]

индексирование гистограммы numpy

учитывая, что у меня есть трехмерная гистограмма или для простоты 3D-массив numpy формы (X, Y, Z) import numpy as np array = np.random.random((100,100,100)) Каков наилучший способ, используя numpy или scipy для получения значений показателей массива, которые удовлетворяют сферическим условиям? (index_x**2 + index_y**2 + index_z**2) <= radius**2 Obvioulsy, в более позднем состоянии, центр массива (0, 0, […]

Массив интеграции Python / Scipy

Я пытаюсь написать программу, которая выполняет следующие действия: Принимает значения V из массива Передает значения V в интеграл, который равен E Выходные интегральные результаты в массив I Участок I против V Уравнение выглядит довольно неприятно, но все является константой, отличной от V. Вот уравнение . Уравнение не очень важно. Как мне решить эту проблему? Моя […]

Как лучше всего масштабировать масштабную матрицу SciPy?

У меня есть разреженная матрица SciPy A , скажем, в формате CSR и вектор v соответствующей длины. Каков наилучший способ масштабирования строк A с помощью v , т. diag(v) * A Выполнения diag(v) * A ?

В чем разница между numpy array.tolist () и scipy.sparse tolist ()

import numpy as np from scipy.sparse import lil_matrix используя numpy, я получаю test_mat = (np.ones((4,6))) test_list = test_mat[0,:].tolist() дает test_list как список, который имеет 6 элементов. Однако я использую scipy.sparse test_mat = lil_matrix(np.ones((4,6))) test_list = test_mat[0,:].todense().tolist() дает test_list как список, который имеет один элемент, который, в свою очередь, имеет 6 элементов ( test_list[0] имеет 6 […]

Как получить сумму каждой строки и суммы каждого столбца в Scipy разреженных матрицах (csr_matrix и csc_matrix)?

У меня очень большая разреженная матрица Scipy (CSR_MATRIX). Я просто хочу знать, как я могу вычислить сумму значений для каждой строки, а также сумму значений для каждого столбца матрицы. У меня есть код, который выполняет ту же операцию, но использует CSC_MATRIX. Есть ли что-то другое между этими двумя для суммирования строк и столбцов? Я подумал, […]

SciPy шаг отклика, кажется, сломается для некоторых значений

Я использую SciPy вместо MATLAB в классе систем управления для построения шаговых ответов систем LTI. До сих пор он отлично работал, но у меня возникла проблема с очень конкретной системой. С помощью этого кода: from numpy import min from scipy import linspace from scipy.signal import lti, step from matplotlib import pyplot as p # Create […]

Оценка ошибки параметра с помощью optimize.fmin_l_bfgs_b

Есть ли способ оценить погрешность оцененных параметров оптимизации, рассчитанную с помощью scipy.optimize.fmin_l_bfgs_b (или любой из других ограниченных подпрограмм оптимизации в пакете scipy.optimize)?

Python - лучший язык программирования в мире.