Python曲线拟合多变量

我试图用多个变量拟合一个函数，我的fit_函数返回两个值，我需要找到适合这两个值的最佳参数

下面是示例代码

import numpy as np
from scipy.optimize import curve_fit

# Fit function returns two values
def func(X, a, b, c):
    x,y = X
    val1 = np.log(a) + b*np.log(x) + c*np.log(y)
    val2 = np.log(a)-4*val1/3
    return (val1,val2)

# some artificially noisy data to fit
x = np.linspace(0.1,1.1,101)
y = np.linspace(1.,2., 101)
a, b, c = 10., 4., 6.
z ,v = func((x,y), a, b, c) * 1 + np.random.random(101) / 100

# initial guesses for a,b,c:
p0 = 8., 2., 7.
   
curve_fit(func, (x,y), (z,v), p0)

它适用于一个返回值的函数，但不适用于两个。它给出：N=3不得超过M=2误差

if n > m:
     raise TypeError('Improper input: N=%s must not exceed M=%s' % (n, m))    
Improper input: N=3 must not exceed M=2

我需要同时找到使

val1-z

和

val2-v

之间的残差最小化的参数

我错过了什么

这就是我的输入数据的样子

---

我需要同时适合z/x和v/x的参数。

scipy.optimize.curve_fit是通过比较func参数列表（a，b，c）的长度（3）和因变量（z，v）的长度（2）来检查是否有至少与拟合参数一样多的数据点。是的，z和v都有三个以上的数据点，但（z，v）的长度是两个。

正如其他人所观察到的，您的函数需要返回带有输入数据形状的内容，因此您需要更改错误函数的输出形状。由于

scipy

执行最小二乘函数，这是通过使函数返回

np.sqrt（val1**2+val2**2）

实现的

然而，对于这种类型的问题，我更喜欢在我编写的

scipy

周围使用包装器，以简化处理多个组件的过程，称为

在

symfit

中，此示例问题将按如下方式解决：

from symfit import parameters, variables, log, Fit, Model
import numpy as np
import matplotlib.pyplot as plt

x, y, z1, z2 = variables('x, y, z1, z2')
a, b, c = parameters('a, b, c')

z1_component = log(a) + b * log(x) + c * log(y)
model_dict = {
    z1: z1_component,
    z2: log(a) - 4 * z1_component/3
}
model = Model(model_dict)
print(model)

# Make example data
xdata = np.linspace(0.1, 1.1, 101)
ydata = np.linspace(1.0, 2.0, 101)
z1data, z2data = model(x=xdata, y=ydata, a=10., b=4., c=6.) + np.random.random(101)

# Define a Fit object for this model and data. Demand a > 0.
a.min = 0.0
fit = Fit(model_dict, x=xdata, y=ydata, z1=z1data, z2=z2data)
fit_result = fit.execute()
print(fit_result)

# Make a plot of the result
plt.scatter(xdata, z1data, s=1, color='blue')
plt.scatter(xdata, z2data, s=1, color='green')
plt.plot(xdata, model(x=xdata, y=ydata, **fit_result.params).z1, color='blue')
plt.plot(xdata, model(x=xdata, y=ydata, **fit_result.params).z2, color='green')

输出：

z1(x, y; a, b, c) = b*log(x) + c*log(y) + log(a)
z2(x, y; a, b, c) = -4*b*log(x)/3 - 4*c*log(y)/3 - log(a)/3

Parameter Value        Standard Deviation
a         2.859766e+01 1.274881e+00
b         4.322182e+00 2.252947e-02
c         5.008192e+00 5.497656e-02
Fitting status message: b'CONVERGENCE: REL_REDUCTION_OF_F_<=_FACTR*EPSMCH'
Number of iterations:   23
Regression Coefficient: 0.9961974241602712

z1（x，y；a，b，c）=b*log（x）+c*log（y）+log（a）
z2（x，y；a，b，c）=-4*b*log（x）/3-4*c*log（y）/3-log（a）/3
参数值标准差
a 2.859766e+01 1.274881e+00
b 4.322182e+00 2.252947e-02
c 5.008192e+00 5.497656e-02
拟合状态消息：b'收敛：根据定义，如果给定一个
x
值，函数应返回一个
y
值。您似乎返回了两个
y
值，我不确定您在解决什么问题。如果你用手把方程写进拟合函数，它会是什么样子？我加上图片我认为你不能得到相同的参数来拟合两条曲线，因为它们不是相同的曲线。左边看起来是单调递增的，而右边的趋势则是平缓的。据我所知，你一次只能装一条线。现在，如果你想要一个适合两个模型（log对正弦和线性），那么考虑SCIKIT学习。