Tag: scipy

scipy.sparse dot чрезвычайно медленный в Python

Следующий код не будет завершен в моей системе: import numpy as np from scipy import sparse p = 100 n = 50 X = np.random.randn(p,n) L = sparse.eye(p,p, format='csc') XTdot(L).dot(X) Есть ли объяснение, почему это матричное умножение висит?

Чтение файла данных Matlab в Python, необходимость экспорта в CSV

Я прочитал файл Matlab, содержащий большое количество массивов в качестве набора данных, в Python, сохраняющий словарь Matlab под макетом имен переменных, используя команду: mat = loadmat('Sample Matlab Extract.mat') Есть ли способ, которым я могу использовать функцию Python для записи в функции csv, чтобы сохранить эту переменную словаря Matlab, которую я читал в Python, в виде […]

сохраняющая фигуру кусочно-кубическая интерполяция для 3D-кривой в питоне

У меня есть кривая в 3D-пространстве. Я хочу использовать кусочно-кубическую интерполяцию, сохраняющую фигуру, на ней аналогично pchip в matlab. Я исследовал функции, предоставляемые в scipy.interpolate, например interp2d, но функции работают для некоторых структур кривой, а не для точек данных, которые у меня есть. Любые идеи о том, как это сделать? Вот данные: x,y,z 0,0,0 0,0,98.43 […]

Как рассчитать вероятность скручивания кривой в scipy?

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

Ошибка при установке scipy на окнах

Я работаю над Windows-машиной и версией python версии 2.7. Я уже установил numpy, и теперь я пытаюсь установить scipy. Я попытался установить его, используя команду easy install, и загрузив zip-файл, а затем запустив файл setup.py. Когда я пытаюсь установить его с помощью easy_install c:\Python27\Scripts>easy_install scipy Searching for scipy Reading https://pypi.python.org/simple/scipy/ Best match: scipy 0.13.3 Downloading […]

Левые собственные векторы, не дающие правильной (марковской) стационарной вероятности в scipy

Учитывая следующую марковскую матрицу: import numpy, scipy.linalg A = numpy.array([[0.9, 0.1],[0.15, 0.85]]) Стационарная вероятность существует и равна [.6, .4] . Это легко проверить, если взять большую мощность матрицы: B = A.copy() for _ in xrange(10): B = numpy.dot(B,B) Здесь B[0] = [0.6, 0.4] . Все идет нормально. Согласно википедии : Стационарный вектор вероятности определяется как […]

Нужна помощь в решении нелинейного ODE второго порядка в python

Я не знаю, с чего начать с этой проблемы, поскольку у меня не было большого опыта в этом, но требуется решить эту часть проекта с помощью компьютера. У меня есть ODE 2-го порядка: m = 1220 k = 35600 g = 17.5 a = 450000 и b составляет от 1000 до 10000 с шагом 500. […]

Изменение размера 3D-изображения (и повторная выборка)

У меня есть трехмерное изображение мозга (давайте назовем его вспышкой), и на данный момент это 263 x 256 x 185. Я хочу изменить его размер другого изображения (назовите его whole_brain_bravo); 256 x 256 x 176 и (надеюсь) использовать интерполяцию lanczos для повторной выборки (Image.ANTIALIAS). Моя (неудачная) попытка: from scipy import ndimage as nd import nibabel […]

Как поднять массив numpy к власти? (что соответствует неоднократным матричным умножениям, а не элементарным)

Я хочу поднять 2-мерный array numpy, назовем его A , мощностью некоторого числа n , но до сих пор мне не удалось найти функцию или оператор для этого. Я знаю, что я мог бы применить его к типу matrix и использовать тот факт, что тогда (подобно тому, что было бы поведением в Matlab), A**n делает […]

Почему scipy.ndimage.io.imread возвращает PngImageFile, а не массив значений

У меня две разные машины с scipy 0.12 и PIL. На одной машине, когда я пытаюсь прочитать файл .png, он возвращает массив целых чисел с размером (wxhx 3): In[2]: from scipy.ndimage.io import imread In[3]: out = imread(png_file) In[4]: out.shape Out[4]: (750, 1000, 4) На другой машине, используя тот же файл изображения, возвращается объект PIL.PngImagePlugin.PngImageFile завернутый […]

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