擅长:python、mysql、java
<p>我认为处理边界的标准方法是使函数在参数超过边界时最小化(残差)。</p>
<pre><code>import scipy.optimize as optimize
def residuals(p,x,y):
if within_bounds(p):
return y - model(p,x)
else:
return 1e6
p,cov,infodict,mesg,ier = optimize.leastsq(
residuals,p_guess,args=(x,y),full_output=True,warning=True)
</code></pre>