最小二乘拟合，令人困惑的 python scipy 赋值查询

Question

我是新手，已经看了几个小时，将其与其他最小二乘拟合示例进行比较，但示例代码似乎没有点击。

从 http://docs.scipy.org/doc/scipy/reference/tutorial/optimize.html 代码是

>>> from numpy import *
>>> x = arange(0,6e-2,6e-2/30)
>>> A,k,theta = 10, 1.0/3e-2, pi/6
>>> y_true = A*sin(2*pi*k*x+theta)
>>> y_meas = y_true + 2*random.randn(len(x))
>>> def residuals(p, y, x):
...     A,k,theta = p
...     err = y-A*sin(2*pi*k*x+theta)
...     return err
>>> def peval(x, p):
...     return p[0]*sin(2*pi*p[1]*x+p[2])
>>> p0 = [8, 1/2.3e-2, pi/3]
>>> print(array(p0))
[  8.      43.4783   1.0472]
>>> from scipy.optimize import leastsq
>>> plsq = leastsq(residuals, p0, args=(y_meas, x))
>>> print(plsq[0])
[ 10.9437  33.3605   0.5834]
>>> print(array([A, k, theta]))
[ 10.      33.3333   0.5236]
>>> import matplotlib.pyplot as plt
>>> plt.plot(x,peval(x,plsq[0]),x,y_meas,'o',x,y_true)
>>> plt.title('Least-squares fit to noisy data')
>>> plt.legend(['Fit', 'Noisy', 'True'])
>>> plt.show()

在残差函数中，为什么可以将 3 个对象 A，k、theta 分配给 p？任何帮助不胜感激

Answer 1

残差函数利用元组解包来检索序列p的元素。元组解包的一个更简单的例子可能是：

x = ['One', 'Two', 'Three']
a,b,c = x
print(a) # One

在残差函数的情况下，它使其更易于阅读（分别使用 A, k, theta 而不是p[0], p[1], p[2]）。