Tag: scipy

Умножение элементов в разреженном массиве с строками в матрице

Если у вас есть разреженная матрица X: >> X = csr_matrix([[0,2,0,2],[0,2,0,1]]) >> print type(X) >> print X.todense() <class 'scipy.sparse.csr.csr_matrix'> [[0 2 0 2] [0 2 0 1]] И матрица Y: >> print type(Y) >> print text_scores <class 'numpy.matrixlib.defmatrix.matrix'> [[8] [5]] … Как вы можете умножить каждый элемент X на строки Y. Например: [[0*8 2*8 0*8 […]

Проверьте верхнюю или нижнюю треугольную матрицу

Есть ли способ, используя numpy или scipy, проверить, является ли матрица нижней или верхней треугольной матрицей ?. Я знаю, как сделать функцию для проверки этого; но я бы хотел знать, имеют ли эти модули свои собственные функции. Я ищу в документации, но ничего не нашел.

Получение dtype массива результатов в numpy

Я хочу предварительно выделить память для вывода операции массива, и мне нужно знать, какой dtype это сделать. Ниже у меня есть функция, которая делает то, что я хочу, но это ужасно уродливо. import numpy as np def array_operation(arr1, arr2): out_shape = arr1.shape # Get the dtype of the output, these lines are the ones I […]

Глобальная минимальная кривая SciPy подходит

Я использую scipy.optimize.curve_fit , но я подозреваю, что он сближается с локальным минимумом, а не с глобальным минимумом. Я попытался использовать имитированный отжиг следующим образом: def fit(params): return np.sum((ydata – specf(xdata,*params))**2) p = scipy.optimize.anneal(fit,[1000,1E-10]) где specf – это кривая, которую я пытаюсь подобрать. Результаты в p , очевидно, хуже минимального значения, возвращаемого curve_fit даже когда […]

ImportError: невозможно импортировать имя linsolve

Я тестирую кусок кода Python, который содержит строку: from scipy import sparse, linsolve Когда я запускаю скрипт, я получаю сообщение об ошибке: from scipy import sparse, linsolve ImportError: cannot import name linsolve Быстрый поиск в Google показывает код для linsolve.py (размещен на Koders.com). Мой вопрос: Является ли linsolve частью scilab – или мне нужно скачать […]

Дамп разреженной матрицы в файл

У меня есть scipy.sparse.csr- матрица и вы захотите отправить ее в CSV-файл. Есть ли способ сохранить разреженность матрицы и записать ее в CSV?

Как решить рекуррентные отношения в python

Я пытаюсь написать код, чтобы дать числовые ответы на рекуррентное отношение. Сама связь проста и определяется следующим образом. Переменная x является целым числом p (i) = p (i + 2) / 2 + p (i-1) / 2, если i> 0 и i <x p (0) = p (2) / 2 p (i) = 1, если […]

Как построить интеграл сигнала с течением времени?

У меня есть зависящий от времени сигнал. Я хочу построить интеграцию с течением времени, когда ось x и значение интегрирования являются осью y. Есть ли способ Python сделать это? Чтобы быть более конкретным: У меня есть массив времени, time и сигнальный массив, signal . Они имеют одинаковую размерность. Мне нужно интегрировать signal течением time с […]

Как улучшить скорость odeint в Python?

Я использую Python и odeint из scipy пакета для решения большого количества (~ 10e6) связанных ODE. Систему уравнений можно сформулировать как сумму некоторых матричных умножений, и я использую numpy с поддержкой blas для этого. Моя проблема в том, что это занимает очень много времени. Когда я просматриваю код, я вижу, что большую часть времени он […]

(Python) Оценка значений доверительных интервалов регрессионного параметра с помощью scikits bootstrap

Я только что начал опробовать красивый загрузочный пакет, доступный через scikits: https://github.com/cgevans/scikits-bootstrap но я столкнулся с проблемой при попытке оценить доверительные интервалы для коэффициента корреляции от линейной регрессии. Полученные доверительные интервалы полностью выходят за пределы исходной статистики. Вот код: import numpy as np from scipy import stats import bootstrap as boot np.random.seed(0) x = np.arange(10) […]

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