SciPy is an open source library of algorithms and mathematical tools for the Python programming language.
1 ответ

Как сделать Curve_Fit в Python

Мне нужно изогнуть набор данных, используя y = x / (a ​​+ x), где a - это параметр, который я должен получить из этого упражнения. import numpy as np import matplotlib.pyplot as plt from scipy.optimize import curve_fit x = [1, 2, 7, 10, 20, 70, 200,…
16 апр '14 в 15:08
1 ответ

Почему 6GB CSV-файл невозможно прочитать целиком в память (64 ГБ) в Numpy

У меня есть CSV-файл, размер которого составляет 6,8 ГБ, и я не могу прочитать его в память в массив массив, хотя у меня есть 64 ГБ ОЗУ CSV-файл содержит 10 миллионов строк, каждая строка содержит 131 запись (смесь int и float) Я пытался прочитать е…
31 янв '18 в 14:05
1 ответ

Решить для нелинейного уравнения, но минимизировать разницу по сравнению с начальным входом

Я ищу, чтобы найти значения для х и у, учитывая приведенную ниже формулу (риск портфеля). Код даёт мне результат, но у проблемы есть несколько решений (то есть более чем одна комбинация x а также y может решить уравнение). Меня интересует только реш…
16 авг '18 в 11:00
1 ответ

Попытка согласовать распределение вероятностей Pearson3 с помощью scipy

Я пытаюсь понять, как приспособить функцию распределения вероятностей, такую ​​как тип Пирсона 3, к набору данных (в частности, среднегодовое количество осадков в области). Я прочитал несколько вопросов по этому поводу, но я все еще что-то упускаю, …
2 ответа

Избегайте оценки функции с одним и тем же вводом несколько раз

Я пытаюсь использовать scipy.optimise.fsolve для решения функции. Я заметил, что функция оценивается с одним и тем же значением несколько раз в начале и в конце шагов итерации. Например, когда оценивается следующий код: from scipy.optimize import fs…
24 май '17 в 16:05
4 ответа

Расчет Z-баллов по среднему и st dev

Я хотел бы спросить, имеет ли какой-либо популярный пакет, такой как: numpy, scipy и т. Д., Встроенную функцию для вычисления Z-Score, если я уже знаю критическое значение, mean и st dev. Я делаю это обычно так: def Zscore(xcritical, mean, stdev): r…
2 ответа

Ошибка атрибута из объекта Minimizer, возвращенная из функции scipy.optimize.minimize()

Используя функцию scipy.optimize.minimize(), я показал разные результаты, используя разные методы для одной и той же целевой функции. Чтобы оценить пригодность к приему, я использую, чтобы рассматривать уменьшенный ци-квадрат в качестве первого крит…
04 июл '13 в 08:41
1 ответ

Подпись метода для якобиана функции наименьших квадратов в scipy

Может ли кто-нибудь привести пример предоставления якобиана функции наименьших квадратов в scipy? Я не могу понять, какую сигнатуру метода они хотят - они говорят, что это должна быть функция, но очень трудно определить, какие входные параметры и в …
19 окт '10 в 04:48
0 ответов

Выбор типа матриц scipy.sparse

Предположим, мне нужно выполнить вычисления A.dot(B).dot(A.T) A.dot(x) где A а также B оба scipy.sparse матрицы и x это плотный np.array, Как мне выбрать тип разреженных матриц для A а также B? я знаю это csc_matrix лучше для колонн-операций и csr_m…
21 янв '18 в 05:25
1 ответ

Почему мой 2D-интерполант генерирует матрицу с измененными осями в SciPy?

Я решаю дифференциальное уравнение с векторными входами y' = f(t,y), y(t_0) = y_0 где у0 = у (х) используя явный метод Эйлера, который говорит, что y_(i+1) = y_i + h*f(t_i, y_i) где t - вектор времени, h - размер шага, а f - правая часть дифференциа…
27 ноя '13 в 17:13
1 ответ

Многовариантный Нормальный PDF в Scipy

Пытаясь оценить Сципи multivariate_normal.pdf функционировать, но продолжать получать ошибки. MWE: import numpy as np from scipy.stats import multivariate_normal as mvnorm x = np.random.rand(5) mvnorm.pdf(x) дает TypeError: pdf() takes at least 4 ar…
13 июн '16 в 19:49
1 ответ

Ошибка импорта на scipy.misc.imrotate

Я пытаюсь создать свою собственную функцию преобразования радона, для которой мне нужно повернуть простое изображение, которое я создал. Согласно документации, функция находится в scipy.misc. Тем не мение, from scipy.misc import imrotate дает мне "н…
30 ноя '13 в 23:46
1 ответ

Подходит функция erf к данным

Так что у меня петля гистерезиса. Я хочу использовать функцию erf, чтобы соответствовать ее моим данным. Часть моего цикла показана черным цветом на нижнем графике. Я пытаюсь использовать scipy.optimize.curve_fit а также scipy.special.erf функция дл…
01 сен '16 в 08:17
1 ответ

Извлечение элементов из нескольких массивов в Numpy/Scipy

У меня есть три пустых массива: A, B и C. Я должен извлечь '1', который является общим (по пикселям) для любых двух из указанных массивов, и установить все остальные элементы как 0. import numpy as np A = np.array([[1,0,0,0,0], [1,0,0,0,0], [1,0,0,0…
19 янв '15 в 11:32
2 ответа

Использование scipy.optimize.fmin_slsqp

Я пытаюсь использовать пакет scipy.optimize, чтобы найти максимум моей функции стоимости. В данном конкретном случае: у меня есть список цен, которые меняются в течение дня. Для упрощения предположим, что день имеет 8 часов, а цена за каждый час выг…
13 фев '14 в 15:40
0 ответов

Могу ли я указать максимальную длину шага в scipy.optimize.minimize()?

Я хочу минимизировать скалярную функцию (энергию) одномерного массива переменных (атомные координаты). Функция и градиент вычисляются внешней программой. Моя проблема в том, что функция scipy делает большие шаги в X, что приводит к сбою f(X). Есть л…
29 ноя '17 в 18:43
1 ответ

Использование интерфейса Scipy cython_blas от Cython, не работающего на векторах Mx1 1xN

Это должно иметь дело с подобной проблемой здесь: вызов BLAS / LAPACK напрямую с использованием интерфейса SciPy и Cython, но отличается, потому что я использую фактический код в примере SciPy здесь _test_dgemm: https://github.com/scipy/scipy/blob/m…
07 июл '17 в 22:43
2 ответа

Как установить Scipy в виртуальной среде Python3.5 | Windows 10?

У меня есть Анаконда (Python 3.6) в моей Windows 10. Это включает в себя Scipy. Я также использую виртуальный Python 3.5 env для поддержки TensorFlow. Теперь проблема в том, что я не могу импортировать Scipy, пока я в этой виртуальной среде. Я попро…
11 май '17 в 09:59
1 ответ

Python Frequency filtering с, казалось бы, неправильными частотами

Сценарий ниже фильтрует частоты, обрезая все частоты больше 6. Однако вместо использования, казалось бы, правильной функции rfftfreq, fftfreq используется. Насколько я понимаю rfftfreq следует использовать вместе с rfft, Почему этот код работает, хо…
06 май '18 в 16:54
1 ответ

Как получить накопленные значения из массивов с метками времени в четных интервалах [PYTHON]

У меня есть массив с временными метками в секундах Например, этот массив с именем a: a = np.array(10,95,99,100,250) Тогда у меня есть массив b что дает мне количество акций, которые приходят в соответствии с отметками времени в: b = np.array(1,2,3,4…
25 июл '14 в 17:09