Вопросы по теме 'scipy-optimize'

Как получить оптимизированные параметры с помощью curve_fit для двух независимых переменных?
Я хочу подобрать кривую для двух независимых параметров x и y . И я хочу оптимизировать свои параметры a , b и c . Я попытался оптимизировать свои параметры, используя curve_fit в scipy. Но мои параметры не оптимизируются. Я использовал...
488 просмотров

Scipy optimize minimize всегда возвращает первоначальное предположение (SLSQP)
Как объясняется в названии, моя программа всегда возвращает первоначальное предположение. Для контекста программа пытается найти лучший способ разместить какой-либо продукт в нескольких магазинах. У каждого магазина есть прогноз ожидаемых продаж в...
1946 просмотров

Pandas векторизованный поиск корней
По сути, я пытаюсь использовать метод brentq row-wise в кадре данных pandas, чтобы получить корень функции, которая также принимает в качестве аргументов столбцы и константы. Что-то вроде следующего import pandas as pd import numpy as np...
154 просмотров

Вопрос по Scipy - Свернуть. Добавление дополнительных ограничений
Я пытаюсь использовать функцию scipy minimize для следующей оптимизации: V = np.matrix(pd.read_csv('V.csv'))` R = np.matrix(pd.read_csv('R.csv', index_col = 'Ticker'))` w0= list() for i in range(0, 84): w0.append(1/84) def...
136 просмотров
schedule 11.11.2022

Как работает алгоритм оптимизации SLSQP?
Я использую алгоритм SLSQP в openMDAO, но мне трудно понять, как он на самом деле работает. Я просто смотрю на обычный пример параболоида, который имеет 2 проектные переменные и направлен на минимизацию f без каких-либо ограничений. Распечатывая...
4492 просмотров

Подгонка функции f (x, y, z) квадратичным многочленом
Я пытаюсь подогнать функцию f (x, y, z) со следующим квадратичным многочленом: 3D-полином Некоторая искаженная сферическая поверхность в трех измерениях. Проблема связана с вычислением эффективных масс в физике твердого тела. Вот...
267 просмотров

Ошибка optim.fmin: IndexError: слишком много индексов для массива
Я пытаюсь оптимизировать функцию в python, используя optim.fmin из scipy. Функция должна оптимизировать вектор параметров, учитывая начальные условия и аргументы. Однако я продолжаю получать следующую ошибку, когда пытаюсь запустить оптимизацию,...
87 просмотров
schedule 01.11.2022

Проблема подбора кривой Python с пиковыми и плоскими (супер) гауссовыми сигналами
У меня есть стандартная гауссовская функция, которая выглядит так: def gauss_fnc(x, amp, cen, sigma): return amp * np.exp(-(x - cen) ** 2 / (2 * sigma ** 2)) И у меня есть функция fit_gaussian, которая использует scipy curve_fit, чтобы...
865 просмотров

разумная оптимизация с несколькими ограничениями
Как найти решение линейной регрессии с несколькими ограничениями на коэффициенты в Python?
57 просмотров

Есть ли функция numpy/scipy для расчета исходящего штрафа?
В моей задаче минимизации все ограниченные методы минимизации, такие как «L-BFGS-B», «TNC», не сходятся, а «Нелдер-Мид» сходятся очень хорошо. Поэтому я предпочитаю использовать «Nelder-Mead» с модифицированной функцией минимизации, например: def...
56 просмотров
schedule 27.09.2022

Как подогнать набор равномерно расположенных параллельных линий к сетке данных
Я хотел бы подогнать набор равномерно расположенных параллельных линий к сетке данных, используя метод минимизации наименьших квадратов. Я могу сопоставить независимые параллельные линии с моими данными. Однако, как и ожидалось, точки пересечения по...
146 просмотров

Программирование абсолютного отклонения как линейной программы
Я пытаюсь преобразовать сумму абсолютных отклонений в задачу линейного программирования, чтобы можно было использовать CPLEX (или другой решатель). Я застрял в том, как настраивать матрицы. Проблема в следующем: minimize abs(x1 - 5) + abs(x2 -...
207 просмотров

Scipy.optimize минимизировать несовместимые ограничения неравенства
Я пытаюсь минимизировать дисперсию портфеля из 100 ценных бумаг. def portvol(w, x): return np.dot(w.T, np.dot(x, w))*252 covmat = annreturn.cov() w0 = np.ones(len(covmat)) * (1 / len(covmat)) #equal weighting initially bounds = ((0,1),) *...
110 просмотров

Как узнать или найти вектор веса и удалить выбросы для всех точек данных, используя Constrained np.polyfit или curve_fit?
Я пытаюсь подогнать exponential decay к некоторым экспериментальным данным и использую polyfit в numpy. Кроме того, сама подгонка взвешена, то есть каждая точка данных имеет некоторый вес, связанный с ней, и polyfit function can find that optimal...
153 просмотров
schedule 13.04.2023

scipy.optimize.minimize не может найти решение методом наименьших квадратов
Мое намерение состояло в том, чтобы решить проблему подгонки ошибок l1 с помощью scipy.optimize, как это предлагается в Норма L1 вместо нормы L2 для функции стоимости в регрессионной модели , но я продолжаю получать неправильное решение, поэтому я...
88 просмотров

Как избежать ошибки значения в fsolve из scipy.optimize
Я пытаюсь решить длинный блок уравнений из реализации EES, используя scipy.optimze.fsolve. Но в этом блоке уравнений есть вызовы CoolProp, которые имеют диапазон проверки, и иногда он выдает ValueError. Я хочу знать, есть ли стратегия избежать...
60 просмотров
schedule 08.06.2022

Как вернуть веселье и jac в scipy.optimize.least_squares
Я использую scipy.optimize.least_squares для минимизации функции из 40 параметров, и для моей конкретной задачи я могу аналитически выразить якобиан моей функции. Однако функция, которую я минимизирую, и якобиан во многом используют одни и те же...
84 просмотров

Кривая оптимизации Scipy не соответствует реальным данным
Я пытаюсь подогнать убывающую экспоненциальную функцию к данным реального мира. У меня проблема с согласованием функции с фактическими данными. Вот мой код: def test_func(x, a, b, c): return a*np.exp(-b*x)*np.sin(c*x) my_time =...
68 просмотров

Есть ли способ инициализировать начальную точку scipy.optimize.linprog?
У меня есть последовательность линейных программ для решения. Каждый экземпляр отличается от предыдущего только тем, что A, границы и стоимость немного отличаются. Интуитивно, решения предыдущих проблем должны помочь. Как я могу реализовать это?...
104 просмотров
schedule 24.05.2024

Scipy curve_fit аварийно завершает работу при использовании Pyside2 QThread
Я использую pyside2 для создания приложения QT5, в котором я использую тысячи функций с помощью scipy.optimize.curve_fit. Чтобы избежать зависания графического интерфейса, я использую Qthread(), чтобы поместить эти тысячи вызовов curve_fit в...
71 просмотров
schedule 27.05.2023