Вопросы с меткой [scipy]

SciPy — библиотека для языка программирования Python с открытым исходным кодом, предназначенная для выполнения научных и инженерных расчётов.

Фильтрация
Сортировка
Метки
0 голосов
2 ответа
67 показов

Как ускорить вычисления вероятностей на python

Есть код, считает правильно, но долго - около 10 сек. Можно ли как-то ускорить вычисления? import scipy.stats as sps mo1, mo2, sig_h, sig_a, summ = 50, 50, 10, 10, 150 def find_prob(mo1, mo2, sig_h, ...
0 голосов
1 ответ
33 показа

Вычисление матричной экспоненты: ValueError: expected A to be a square matrix

Хочу посчитать: from jax.scipy.linalg import expm import jax.numpy as jnp from functools import lru_cache, reduce num_qubits=2 theta = jnp.asarray(np.pi*np.random.random((15,2,2,2,2,2,2,2,2))) def ...
0 голосов
1 ответ
74 показа

Построение огибающей по серединам отрезков гистограммы

Делаю курсовую, одним из заданий которой является построение огибающей для распределения. Как пример приведу биномиальное распределение, у меня есть следующий код: import matplotlib.pyplot as plt ...
0 голосов
0 ответов
38 показов

Разные результаты при использовании tuple и np.array в функции minimize библиотеки scipy.optimize

Выполняю минимизацию функции, и при использовании данных в формате tuple, и все корректно работает, но очень долго. Хотел оптимизировать код, перевел все на np.array, но теперь получаются другие ...
0 голосов
0 ответов
70 показов

Метод прогонки в scipy?

Коллеги, доброе время суток. Подскажите, пожалуйста, реализован ли метод прогонки для решения СЛАУ с 3-х диагональной матрицей в библиотеке scipy в "чистом" виде? Почитал про функцию ...
0 голосов
0 ответов
39 показов

ERROR: Dependency lookup for OpenBLAS with method 'pkgconfig' failed. Не удается найти cmake, даже если я его установил

Я пытался установить scipy из проекта в conda env, а затем в pip env, но scipy==1.7.0 выдает ошибку. Я использовал Python 3.8.18, а ему требуется Python >=3.7, <3.10. Я не понимаю, почему это ...
0 голосов
1 ответ
45 показов

Как правильно задать функцию для фитирования массива?

Имеется следующий код, читающий данные из текстового файла (файл можно скачать по ссылке): import matplotlib.pyplot as plt import numpy as np from scipy.optimize import curve_fit def ...
1 голос
1 ответ
1k показов

как поместить изображение на график?

скажите пожалуйсте как поместить картинку на график и задать ее положение, спасибо! Я использую pylab pylab.title(filename+' IP=539.1') pylab.plot(Omega,Sigmatot,color="black",linewidth=0.5) #Plots ...
0 голосов
1 ответ
127 показов

Ошибка KeyError: 4

Делаю курсовую работу по цифровой обработке сигнала, в открываемом файле содержатся точки сигнала, всё работало, в последнюю очередь нужно было расшифровать сообщение в сигнале. Написал вот такой вот ...
0 голосов
0 ответов
36 показов

Не могу подобрать функцию для метода curve_fit from scipy.optimize. График - какие-то периодические зубцы

import numpy as np import matplotlib.pyplot as plt from scipy.optimize import curve_fit from scipy.stats import sem x=[-10.0, -9.5, -9.0, -8.5, -8.0, -7.5, -7.0, -6.5, -6.0, -5.5, -5.0, -4.5, -4.0, -3....
0 голосов
1 ответ
64 показа

Как получить продолжение последовательности через обратное преобразования Фурье

У меня есть последовательность чисел N и я хочу получить её продолжение из частотного преобразования Фурье. from scipy.fft import fft, ifft import matplotlib.pyplot as plt N = np.arange(100) N_fft = ...
0 голосов
1 ответ
50 показов

Проблемы с t-test scipy

Есть такое задание. Нужно в python посчитать t-test по уже известным данным (оценка коэффициента и стандартное отклонение). Есть функция в библиотеке scipy, которая позволяет это сделать, но нужно ...
0 голосов
0 ответов
156 показов

В какой библиотеке можно постоить нелинейную регрессию?

У меня есть данные (df.head и .tail) датафрейм: date mentor_cnt mentee_cnt mentee_dynamic 0 2021-02-12 2 2 nan 1 2021-02-19 1 1 0,5 2 2021-03-06 1 1 1 3 2021-03-27 2 2 2 4 2021-04-10 4 4 2 5 ...
1 голос
1 ответ
121 показ

Есть ли библиотека для решения систем нелинейных уравнений для Java/Kotlin?

Нужно решить системы нелинейных уравнений на Java или Kotlin. И встал вопрос, а есть ли что-то похожее на fsolve из scipy.optimize для этих языков? Может кто знает и может подсказать?
0 голосов
2 ответа
568 показов

Подогнать функцию под кривую методом curve_fit модуля scipy

Всем привет, пытаюсь решить задачу, но не могу, то ли я делаю что то не так, то ли использую не те инструменты. Есть набор данных, далее представлены в коде. Есть формула "PM+(Q-P)(x)-(Q/M)*(x**2)...
0 голосов
1 ответ
66 показов

Для заданного y найти значения x

Помогите, пожалуйста, найти подход к решению следующей задачи: Пусть есть матрица X = (x1,…,xn), состоящая из временных рядов xi и вектор y. Суть задачи в том, чтобы научиться косвенно управлять ...
0 голосов
1 ответ
95 показов

scipy.stats.gamma.fit

В библиотеке scipy.stats согласно документации распределение gamma имеет только один параметр, но метод scipy.stats.gamma.fit(data) почему то возвращает кортеж из трех значений, а не одно! Почему так ...
0 голосов
1 ответ
73 показа

Эмпирическая функция распределения

В теме построение дискретной эмпирической функции распределения (cdf) описывается как построить график эмпирической функции распределения. А как вывести еë формулу? Чтобы потом по подобранной функции ...
0 голосов
1 ответ
621 показ

Интерполяция с помощью scipy

В целом задание звучит так: Існує набір даних speed, який є значенням швидкості, для транспортного засобу в певній точці спостереження. Очевидно, що дані має дискретний вид. Відомо, що спостереження ...
0 голосов
0 ответов
127 показов

sklearn бинарная классификация

Нужно предсказать, результат будет положительным или отрицательным. Есть обучающая выборка (.csv), тестовая - в них однозначные числа с дробью в 16 знаков, по типу : -5.072591732401996 -1....
0 голосов
2 ответа
217 показов

Восстановить матрицу

Есть матрица Y 20х20. Она бинарная. К ней применяют функцию "step" и получается другая бинарная матрица X 20х20. Задача: Мы знаем только матрицу Х. Мы знаем функцию, по которой получили ...
0 голосов
1 ответ
77 показов

Проблема со скоростью решения задачи SciPy.Minimaze (SLSQP)

Я пытаюсь написать решатель для производственной задачи, который должен распределять рабочее время по дням и рабочим местам, минимизируя переработку. В водных данных у меня 6000 переменных и 7000 ...
0 голосов
0 ответов
50 показов

Как корректно реализовать обработку нескольких одновременных условий при решении ОДУ scipy_ivp?

Решается задача описания работы двухстороннего пневматического цилиндра с нахождением положения поршней и величины давления в полостях. При описании работы пневмоцилиндра формируются следующие ...
2 голоса
1 ответ
597 показов

Поиск минимума функции с scipy.optimize.minimize

Прошу помочь с исправлением ошибки. Код оптимизации следующий: import numpy as np from scipy.optimize import Bounds, minimize m = 20 n = 10 np.random.seed(1) X = np.random.randn(m, n) y = np.random....
0 голосов
1 ответ
701 показ

Создать простой сплайн

Как сделать апроксимизацию сплайном в 1-2 строки кода и очень по простому? Допустим данные датафрейм, хотя можно и для списка - не так важно. function 1 240.9 2 241.18 3 241.166 4 241.2 5 ...
0 голосов
1 ответ
66 показов

Как можно кластеризовать пользователей по степени сходства их поведения?

Имеется pandas DataFrame: |'User1' |'User2' |'Correlation'| |--------|---------|-------------| |'Ivanov'|'Petrov' |0.9| |'Ivanov'|'Sidorov'|0.1| |'Petrov'|'Sidorov'|0.2| В итоге хотелось бы получить ...
1 голос
1 ответ
222 показа

Фон дендрограммы в Python заливается серым цветом

Мой код: credit_df = pd.read_csv('CC_GENERAL.csv') print('Dimension of array:{}'.format(credit_df.shape)) #Dimension of array:(8950, 18) from sklearn.decomposition import PCA pca = PCA(n_components=0....
3 голоса
2 ответа
1k показов

Как работает scipy.stats.binned_statistic_2d?

У меня есть список из примерно 1000 элементов. Каждый элемент содержит две координаты x и y, обе от 0 до 200. Мне необходимо разбить все точки на 25 секторов, каждый из которых имеет длину 40 на 40. ...
0 голосов
0 ответов
104 показа

Не получается минимизировать функцию потерь

Пытаюсь минимизировать функцию потерь радикальных оценок для оценки параметров нормального распределения Мой код выглядит следующим образом: import numpy as np from scipy import stats from scipy....
0 голосов
1 ответ
173 показа

не считается ковариация между двумя выборками из распределений из scipy.stats

Задача состоит в том чтобы посчитать ковариацию между двумя выборками из двух распределений (равномерного и распределения Максвелла). В задании предлагается использовать библиотеку numpy чтобы ...
0 голосов
1 ответ
67 показов

Ошибка при создании фрейма из матрицы svd

Мне нужно получить векторное представление слов с помощью tf-idf и понизить размерность векторов до 300, чтобы в дальнейшем я мог сложить этот вектор с предобученным word2vec. В качестве корпуса слов ...
1 голос
1 ответ
348 показов

Как методом интерполяции можно заполнить недостающие значения в колонке?

Как методом интерполяции можно заполнить недостающие значения в колонке "Production rate", какие есть варианты? df = pd.DataFrame({'Date': ['2010-01-10', '2010-01-09', '2010-01-08', ...
7 голосов
1 ответ
4k показов

Как ускорить цикл for для обхода массива?

Есть два массива. Двумя циклами for обхожу первый массив, на основании данных из которого, заполняю второй. При больших размерах массива программа выполняется очень долго. Подскажите методы ускорения ...
1 голос
1 ответ
746 показов

Как установить scipy с использованием pip?

Я уже пробовал методы обновления pip и wheel или же пробовать python3 -m pip install -U scipy, но все равно одна и та же ошибка и после чего pip пытается установить более старую версию scipy, в ...
0 голосов
2 ответа
257 показов

поиск минимума функции python [закрыт]

Прошу подсказать способ решения следующей задачи. Есть ли готовый метод в scipy Даны две функции g1, g2 g1 = c1*x1+c2*x2+c3*x3+c4*x4 g2 = m1*x1+m2*x2+m3*x3+m4*x4 необходимо найти такие целочисленные ...
0 голосов
2 ответа
496 показов

Как найти распределение вероятностей для непрерывной количественной переменной в python? [закрыт]

У меня есть список значений непрерывной количественной переменной(от 10 млрд до 80), всего 70 значений. Я хочу получить распределение вероятностей для каждого значения. P.S Я знаю про distplot у ...
4 голоса
1 ответ
115 показов

Как проинтегрировать уравнение из переменной?

Имеется код интегрирования некоторой функции f(x), но уравнение данной функции высчитывается программой заранее и записывается в переменную m_1. С не отформатированным уравнением (return (x)**2) ...
2 голоса
1 ответ
66 показов

Кластеризация scipy python

Есть вот такой 2D массив точек, делаю кластеризацию arr = [[-0.23849030299472848, 0.2863843214610768], [-0.23467813388970515, 0.2847352640189306], [-0.2317909064951051, 0.2841855891260814], [-0....
0 голосов
1 ответ
140 показов

Python scipy кластеризация

Столкнулся с такой проблемой, есть массив точек как на картинке, как я могу кластеризовать его и посчитать количество отдельных объектов в моём массиве? В данном случае хочу увидеть цифру 3 Раньше ...
0 голосов
1 ответ
542 показа

Python группировка списка

Как можно разделить такой массив на похожие элементы? list_coord = [[1, 1, 1], [2, 1, 1], [1, 1, 1], [1, 2, 2], [12, 11, 11], [11, 13, 13], [13, 12, 12], [7, 6, 6], [6, 6, ...
1 голос
1 ответ
2k показов

Построение аудио спектра в реальном времени

Всем привет. Мне нужен анализатор спектра в реальном времени как в звуковых плеерах. На языке python. Код для чтение wav файла и fft готов, но не понимаю как строить "столбики звука". ...
0 голосов
1 ответ
42 показа

Не работает программа с нейронной сетью

Сам код. import numpy as np from scipy.optimize import root def your_funcs(X): x, y = X # all RHS have to be 0 f = [x + y**2 - 4, np.exp(x) + x * y - 3] return f sol = ...
1 голос
1 ответ
285 показов

не получается корректно отобразить частотный спектр (Фурье)

Необходимо получить частоты присутствующих в сигнале гармоник. Изначально дан файл с самим сигналом и частота дискретизации 500 МГц. Сигнал выглядит следующим образом Для того, чтобы получить частоты ...
2 голоса
1 ответ
153 показа

Поиск циклов в данных

Пытаюсь найти циклы в данных представленных ниже. Как сделать это правильно и эффективно? import matplotlib.pyplot as plt from scipy.misc import electrocardiogram x = electrocardiogram()[2000:4000] ...
0 голосов
1 ответ
2k показов

Выбор метода интерполяции и экстраполяции

Есть задача выбора аккумуляторной батареи по требуемому времени работы и мощности. Первый шаг — рассчитать мощность по заданному времени исходя из таблицы разрядных характеристик. Для этого нужно ...
1 голос
1 ответ
2k показов

Аппроксимация функции с помощью функционала scipy.optimize.minimize

Пытаюсь аппроксимировать функцию с помощью полинома. Допустим реальные данные задаются следующей функцией: # Функция генерирует реальные значения def data_gen(x, alpha=1, beta=10, gama=5, per=1.75): ...
2 голоса
1 ответ
1k показов

Поиск коэффициентов уравнения

В уравнении вида y = alpha + k*x^n необходимо найти коэффициенты alpha, k, n. Показатель степени n от 0 до 1. При использовании функции curve_fit из scipy получаются совсем некорректные значения. ...
0 голосов
1 ответ
3k показов

Аппроксимация функции нескольких переменных

Для нахождения приближающей функции y = f(x) по методу наименьших квадратов и построения графиков я использовал следующие функции из пакета NumPy: t = np.polyfit(x, y, 3) f = np.poly1d(t) ax.plot(x, f(...
0 голосов
0 ответов
128 показов

Артефакты при изменении частоты дискретизации WAV файла

Реализую изменение частоты дискретизации WAV файла в кратное число раз на Python. Проблема появилась при реализации потоковой обработки файла. При обработке файла целиком такой проблемы не возникает, ...
1 голос
1 ответ
317 показов

SNR аудио файла

Мне нужно вычислить SNR аудио файла. Нашел информацию, что в есть готовая функция в пакете: scipy.stats.signaltonoise. Но не нашел примеров использования. Как эта функция собственно работает( ...