Поскольку я ленив, я не хочу определять функцию для оптимизации таким образом: f = lambda x, a, b, …: a * x + c + … но я хочу сделать так: f = lambda x, p: p[0] * x + p[1] + … где p – массив начальных параметров. Проблема в том, что я не […]
Я ищу метод для решения системы линейных уравнений в Python. В частности, я ищу наименьший целочисленный вектор, который больше всех нулей и решает данное уравнение. Например, у меня есть следующее уравнение: и хотите решить , В этом случае наименьший целочисленный вектор, который решает это уравнение, равен , Однако как я могу определить это решение автоматически? […]
для проекта мне нужна эффективная функция в python, которая решает следующую задачу: Учитывая очень большой список X длинных разреженных векторов (=> большая разреженная матрица) и другая матрица Y, содержащая один вектор y, мне нужен список «расстояний», который у каждого элемента X. Таким образом, «расстояние "определяется следующим образом: Сравните каждый элемент двух векторов, всегда берете нижний […]
У меня есть текстовый корпус, содержащий более 1000 статей в отдельной строке. Я пытаюсь использовать иерархию иерархии с помощью Scipy в python для создания кластеров связанных статей. Это код, который я использовал для кластеризации # Agglomerative Clustering import matplotlib.pyplot as plt import scipy.cluster.hierarchy as hac tree = hac.linkage(X.toarray(), method="complete",metric="euclidean") plt.clf() hac.dendrogram(tree) plt.show() и я получил […]
В настоящее время я работаю над этим примером SciPy на Kernal Estimation. В частности, тот, который обозначен как «Одномерная оценка». В отличие от создания случайных данных, я использую возврат активов. Моя вторая оценка, хотя (и даже простая норма pdf, которую я создаю для сравнения) показывает плотность, которая достигает 20, что не имеет смысла … Мой […]
У меня есть следующий код: import scipy.optimize class demo(object): def get_square(self, var): return var ** 2 – 4 new = demo() scipy.optimize.fsolve(new.get_square(), 1) И я получил следующую ошибку: TypeError: get_square() missing 1 required positional argument: 'var' Но get_square () всегда должен иметь « self а « self не должен быть передан. В чем проблема?
Учитывая N наблюдаемых (x, y) точек данных с установленными стандартными отклонениями (sx, sy) , я сопоставляю линейную модель (y=a*x+b) с использованием пакета scipy.ODR . Обратите внимание, что значения (sx, sy) также учитываются при подгонке (через RealData ). Как только фитинг завершен, мне необходимо назначить правильные стандартные отклонения к установленным параметрам (a, b) . Пакет в […]
В настоящее время я пишу скрипт python, который должен сделать следующее: читайте аудиофайл, соответственно, wav-файл через scipy.io.wavfile.read (). вычислить спектрограмму данного wav-файла. записывайте данные из спектрограммы обратно в wav-файл. Когда я пытаюсь запустить скрипт, я получаю следующую ошибку: Traceback (most recent call last): File "demo2.py", line 61, in <module> F,T,S = scipy.signal.spectrogram(s_mono,rate,window,None,nOverlap,nFFT) AttributeError: 'module' object […]
Я пытаюсь использовать skimage.restoration.wiener , но у меня всегда получается изображение с кучей 1 (или -1), что я делаю неправильно? Оригинальное изображение происходит от Uni of Waterloo . import numpy as np from scipy.misc import imread from skimage import color, data, restoration from scipy.signal import convolve2d as conv2 def main(): image = imread("/Users/gsamaras/Downloads/boat.tif") psf = […]
Я использую интерполяцию 1d, используя scipy для временных рядов. Мои данные по оси x находятся в формате даты и времени, а ось y находится в режиме float, например: 3/15/2012 16:00:00 32.94 3/16/2012 16:00:00 32.95 3/19/2012 16:00:00 32.61 Теперь во время вычисления slope = (y_hi-y_lo) / (x_hi-x_lo) я получаю ошибку TypeError: unsupported operand type(s) for /: […]