Tag: scipy

Может ли какой-либо из стандартных пакетов Python реализовывать кубические и квартирные формулы закрытой формы?

Существуют точные (замкнутые) решения полиномиальных уравнений степени до и включительно 4. То есть существует квадратичная формула, кубическая формула и даже квартичная формула. Выполняют ли какие-либо из стандартных пакетов python (или даже некоторые нестандартные)? Я вижу, что Numpy предлагает возможность найти полиномиальные корни ( numpy.polynomial.polynomial.polyroots ), но эта процедура, по-видимому, находит корни численно с помощью метода […]

Найдите самый значительный угол скелета и сегментируйте скелет в этом углу

У меня есть изображения рудных швов, которые у меня сначала скелетонированы (средняя ось умножается на дистанционное преобразование), затем извлекаются углы (см. Зеленые точки). Это выглядит так: Проблема состоит в том, чтобы найти точку поворота, а затем отделить шов, разделив шов в точке поворота. Не все скелеты имеют точки поворота, некоторые из них довольно линейны, а […]

Поиск корней с матрицей SymPy

Редактирование: я нашел рабочее решение, но мне все равно понравилось бы больше объяснений тому, что здесь происходит: from scipy import optimize from sympy import lambdify, DeferredVector v = DeferredVector('v') f_expr = (v[0] ** 2 + v[1] ** 2) f = lambdify(v, f_expr, 'numpy') zero = optimize.root(f, x0=[0, 0], method='krylov') zero Оригинальный вопрос: Ниже мы имеем […]

Python: Как вычислить среднеквадратичную ошибку распределения?

Я поместил данные с GMM с данными, я хочу рассчитать среднеквадратичную ошибку модели, как я могу это сделать? Вот код для генерации данных import numpy as np import matplotlib.pyplot as plt from matplotlib.colors import LogNorm from sklearn import mixture import matplotlib as mpl from matplotlib.patches import Ellipse %matplotlib inline n_samples = 300 # generate random […]

Минимизация функции линейной комбинации данных со Scipy

Предположим, что у меня есть матрица X, где каждая строка представляет собой временной ряд. Например, X может быть матрицей размером 3 x 1000, что будет означать, что есть 3 временных ряда, каждый из которых состоит из 1000 временных точек. В дополнение к X, у меня есть один скаляр для каждого временного ряда в X. Я […]

индексы точки питона в KDTree

Учитывая список пунктов, как я могу получить их индексы в KDTree? from scipy import spatial import numpy as np #some data x, y = np.mgrid[0:3, 0:3] data = zip(x.ravel(), y.ravel()) points = [[0,1], [2,2]] #KDTree tree = spatial.cKDTree(data) # incices of points in tree should be [1,8] Я мог бы сделать что-то вроде: [tree.query_ball_point(i,r=0) for […]

ошибка переполнения

Я хочу заниматься музыкой, звуком и музыкой. До этого я создал образец сигнала с 10-секундной разверткой. У меня есть простой скрипт, который должен отображать некоторые сигналы. Первый сигнал – простой синус; второй – развертка; Оба с частотой чуть ниже частоты Найквиста, поэтому нет проблем. Код: #import import numpy as np import scipy.signal as sig import […]

AttributeError: объект 'Float' не имеет атрибута 'sqrt'

Я пытаюсь использовать функцию scipy.stats.linregress (). Я дважды проверил, что оба входа имеют одинаковый размер, и оба являются массивами. Однако я продолжаю получать следующую ошибку: Traceback (most recent call last): File "C:/Users/chianh/Desktop/Reliability Py Scripts/MedianRanks_Solver_Rev0.py", line 48, in <module> slope, intercept, r_value, p_value, std_err = stats.linregress(data, median_rank) File "C:\Users\chianh\AppData\Local\Programs\Python\Python35-32\lib\site-packages\scipy\stats\_stats_mstats_common.py", line 81, in linregress r_den = np.sqrt(ssxm […]

Решение целочисленного линейного программирования в Python SciPy

У меня есть Integer Linear Programming probem для решения. Итак, я решил o использовать python и SciPy spcecifically. У меня есть следующий короткий фрагмент. obj_function = running_instance_coeffs + new_instance_coeffs A = [[1]*len(running_instance_coeffs) + [1]*len(new_instance_coeffs), encoding_speed_running_instances + encoding_speed_new_instances] b = [MAX_INSTANCES, ((-1)*load/QoS)] bounds=[] for inst in running_instances: bounds.append([0,1]) for inst in instances: bounds.append([0,10]) res=linprog(obj_function, A_ub=A, b_ub=b, […]

Ошибка: исчерпана виртуальная память: не удается выделить память при установке SciPy

Я пытаюсь установить SciPy в моем Amazon AWS, который имеет 1 ГБ ОЗУ, и я закончил с ошибкой, говорящей, что «виртуальная память исчерпана: не может выделить память». В соответствии с этой ссылкой я сделал sudo /bin/dd if=/dev/zero of=/var/swap.1 bs=1M count=1024 sudo /sbin/mkswap /var/swap.1 sudo /sbin/swapon /var/swap.1 но это не сработало для меня. Кроме того, я […]

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