时间戳

首页 > 解决方案 > 如何使用 scipy.optimize.minimize 修复推导条件的拟合

python - 如何使用 scipy.optimize.minimize 修复推导条件的拟合

问题描述

我试图在forcav=p(T,x)之后的推导始终为正的条件下x拟合函数。和的数据来自 excel 表格。是我想要得到的系数。cavconstant Tdp/dx (for constant T) > 0xTpz

我已经使用了这里的解决方案，以对派生 Python 的约束条件作为模板。这是我的代码，因为它现在提供错误消息：

import pandas as pd
import os
from scipy.optimize import minimize
import numpy as np

df = pd.read_excel(os.path.join(os.path.dirname(__file__), "./data.xlsx"))

T = np.array(df['T'], dtype=float)
x = np.array(df['x'], dtype=float)
p = np.array(df['p'], dtype=float)
p_s = 67

def cav(z,T,x):  #my function
    return x * p_s + x * (1 - x) * (z[0] + z[1] * T + z[2] * T ** 2 + z[3] * x + z[4] * x * T + z[5] * x * T ** 2) * p_s

def resid(p,T,x):
    return ((p-cav(T,x))**2).sum()

def constr(z):
    return np.gradient(cav(z,x,T))

con1 = {'type': 'ineq', 'fun': constr}

z0 = np.array([0,0,0,0,0,0], dtype=float)
res = minimize(resid,z0, args=(p,T,x), method='cobyla',options={'maxiter':50000}, constraints=con1)

和错误：

TypeError: resid() takes 3 positional arguments but 4 were given

我不明白我到底要为三个参数输入什么def。谢谢你的帮助！

标签： pythonpython-3.xscipyminimize

解决方案

resid该错误是因为除了初始猜测之外，您还传递了 3 个参数z0。

因此，必须改变的行是：

res = minimize(resid,z0, args=(T,x), method='cobyla',options={'maxiter':50000}, constraints=con1)

您的代码中的另一个问题是：

def resid(p,T,x):
    return ((p-cav(T,x))**2).sum()

你的方法cav需要三个参数，但你只传递了两个。所以这可能应该改为：

def resid(p,T,x):
    return ((p-cav(p,T,x))**2).sum()

推荐阅读