怎么在python中使用optimize？相信很多沒有經(jīng)驗(yàn)的人對(duì)此束手無策，為此本文總結(jié)了問題出現(xiàn)的原因和解決方法，通過這篇文章希望你能解決這個(gè)問題。

靖遠(yuǎn)ssl適用于網(wǎng)站、小程序/APP、API接口等需要進(jìn)行數(shù)據(jù)傳輸應(yīng)用場(chǎng)景，ssl證書未來市場(chǎng)廣闊！成為創(chuàng)新互聯(lián)公司的ssl證書銷售渠道，可以享受市場(chǎng)價(jià)格4-6折優(yōu)惠！如果有意向歡迎電話聯(lián)系或者加微信：18982081108（備注：SSL證書合作）期待與您的合作！

Python的優(yōu)點(diǎn)有哪些

1、簡(jiǎn)單易用，與C/C++、Java、C# 等傳統(tǒng)語言相比，Python對(duì)代碼格式的要求沒有那么嚴(yán)格；2、Python屬于開源的，所有人都可以看到源代碼，并且可以被移植在許多平臺(tái)上使用；3、Python面向?qū)ο?，能夠支持面向過程編程,也支持面向?qū)ο缶幊蹋?、Python是一種解釋性語言，Python寫的程序不需要編譯成二進(jìn)制代碼，可以直接從源代碼運(yùn)行程序；5、Python功能強(qiáng)大，擁有的模塊眾多，基本能夠?qū)崿F(xiàn)所有的常見功能。

非線性方程組求解

SciPy中對(duì)非線性方程組求解是fslove()函數(shù)，它的調(diào)用形式一般為fslove(fun, x0)，fun是計(jì)算非線性方程組的誤差函數(shù)，它需要一個(gè)參數(shù)x，fun依靠x來計(jì)算線性方程組的每個(gè)方程的值（或者叫誤差），x0是x的一個(gè)初始值。

"""
計(jì)算非線性方程組：
  5x1+3 = 0
  4x0^2-2sin(x1x2)=0
  x1x2-1.5=0
"""
## 誤差函數(shù)
def fun(x):
  x0,x1,x2 = x.tolist()
  return[5*x1+3,4x0^2-2sin(x1x2),x1x2-1.5]

result = optimize.fsolve(fun,[1,1,1])
## result
[-0.70622057  -0.6  -2.5]

在計(jì)算非線性方程中的解時(shí)，比如像坐標(biāo)上升算法，其中需要用到未知數(shù)的導(dǎo)數(shù)，同樣，scipy的fslove()也提供了fprime參數(shù)傳遞未知數(shù)的雅各比矩陣從而加速計(jì)算，傳遞的雅各比矩陣每一行時(shí)某一方程對(duì)各個(gè)未知數(shù)的導(dǎo)數(shù)。對(duì)于上面的例子，我們可以寫下如下的雅各比矩陣傳入。

def j(x):
  x0,x1,x2 = x.tolist()
  return[[0,5,0],[8*x0,-2*x2*cos(x1*x2],[0,x2,x1]]

result = optimize.fsolve(fun,[1,1,1],fprime=j)
#result
[-0.70622057  -0.6  -2.5]

scipy的內(nèi)部在實(shí)現(xiàn)fslove時(shí)應(yīng)該時(shí)應(yīng)該是利用了坐標(biāo)上升算法或者梯度相關(guān)優(yōu)化算法，但本人沒有考證，有興趣的可以看看源碼。

最小二乘擬合

關(guān)于最小二乘算法的理論這里并不想談，網(wǎng)上解釋的文章也挺多，在 optimize模塊中，可以使用leastsq()對(duì)數(shù)據(jù)進(jìn)行最小二乘擬合計(jì)算。 leastsq()的用法很簡(jiǎn)單，只需要將計(jì)箅誤差的函數(shù)和待確定參數(shù)的初始值傳遞給它即可。

x = np.array([8.19,2.72,6.39,8.71,4.7,2.66,3.78])
y = np.array([7.01,2.78,6.47,6.71,4.1,4.23,4.05])
def residual(p):
  k,b = p
  return y-(k*x+b)
r = optimize.leastsq(residual,[1,0])
k,b = r[0]
# print k
.613495349193
# print b
.79409254326

def func(x,p):
  """
    計(jì)算的正弦波 ：A*sin(2*pi*k*x+theta)
  """
  A,k,theta = p
  return A*sin(2*np.pi*k*x+theta)

def redis(p,y,x):
  return y-func(x,p)

x = np.linspace(0,2*np.pi,100)
A,k,theta = 10,0.34,np.pi/6
y0 = func(x,[A,k,theta])
# 加入噪聲
np.random.seed(0)
y1 = y0+2*np.random.randn(len(x))
p0 = [7,0.40,0]
# p0是A,k,theta的初始值，y1，x要擬合的數(shù)據(jù)
plsq = optimize.leastsq(redis, p0,args=(y1,x))
print [A,k,theta] #真是的參數(shù)值
print plsq[0]  #擬合后的參數(shù)值

對(duì)于像正弦波或者余弦波的曲線擬合，optimize提供curve_fit()函數(shù)，它的使用方式和leastq()稍有不同，它直接計(jì)算曲線的值，比如上面的擬合正弦波可以用cureve_fit()來寫。

def func2(x,p):
  """
    計(jì)算的正弦波 ：A*sin(2*pi*k*x+theta)
  """
  A,k,theta = p
  return A*sin(2*np.pi*k*x+theta)
ret,_=optimize.curve_fit(func2,x,y1,p0=p0)

該函數(shù)有一個(gè)缺點(diǎn)就是對(duì)于初始值敏感，如果初始頻率和真實(shí)頻率值差太多，會(huì)導(dǎo)致最后無法收斂到真是頻率。

局部最小值

optimize模塊還提供了常用的最小值算法如：Nelder-Mead、Powell、CG、BFGS、Newton-CG等，在這些最小值計(jì)算時(shí)，往往會(huì)傳入一階導(dǎo)數(shù)矩陣(雅各比矩陣)或者二階導(dǎo)數(shù)矩陣(黑塞矩陣)從而加速收斂，這些最優(yōu)化算法往往不能保證收斂到全局最小值，大部分會(huì)收斂到局部極小值。這些函數(shù)的調(diào)用方式為:

optimize.minimize(target_fun,init_val,method,jac,hess) 
target_fun：函數(shù)的表達(dá)式計(jì)算； 
init_val：初始值； 
method：最小化的算法； 
jac：雅各比矩陣 
hess：黑塞矩陣。

全局最小值算法

全局最小值使用optimize.basinhopping()來實(shí)現(xiàn)，這個(gè)函數(shù)首先要定義一個(gè)誤差計(jì)算方式，比如平方誤差函數(shù)，niter時(shí)迭代的次數(shù)，最后還需要一個(gè)局部極小值優(yōu)化方法，minimizer_kwargs傳入。比如上面的正弦函數(shù)擬合：

def func1(x,p):
  """
    計(jì)算的正弦波 ：A*sin(2*pi*k*x+theta)
  """
  A,k,theta = p
  return A*sin(2*np.pi*k*x+theta)
def func_error(p,y,x):
  return np.sum((y-func1(x,p)**2)
result = optimize.basinhopping(func_error,[1,1,1],niter=10,
              minimizer_kwargs={"method":"L-BFGS-B",
                        "args":(y1,x1)})
## [1,1,1]是傳入的初始值，args是需要擬合的數(shù)據(jù)

看完上述內(nèi)容，你們掌握怎么在python中使用optimize的方法了嗎？如果還想學(xué)到更多技能或想了解更多相關(guān)內(nèi)容，歡迎關(guān)注創(chuàng)新互聯(lián)成都網(wǎng)站設(shè)計(jì)公司行業(yè)資訊頻道，感謝各位的閱讀！

另外有需要云服務(wù)器可以了解下創(chuàng)新互聯(lián)scvps.cn，海內(nèi)外云服務(wù)器15元起步，三天無理由+7*72小時(shí)售后在線，公司持有idc許可證，提供“云服務(wù)器、裸金屬服務(wù)器、高防服務(wù)器、香港服務(wù)器、美國(guó)服務(wù)器、虛擬主機(jī)、免備案服務(wù)器”等云主機(jī)租用服務(wù)以及企業(yè)上云的綜合解決方案，具有“安全穩(wěn)定、簡(jiǎn)單易用、服務(wù)可用性高、性價(jià)比高”等特點(diǎn)與優(yōu)勢(shì)，專為企業(yè)上云打造定制，能夠滿足用戶豐富、多元化的應(yīng)用場(chǎng)景需求。

文章標(biāo)題：怎么在python中使用optimize-創(chuàng)新互聯(lián)
URL標(biāo)題：http://chinadenli.net/article38/dojgpp.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供網(wǎng)頁設(shè)計(jì)公司、網(wǎng)站建設(shè)、服務(wù)器托管、商城網(wǎng)站、網(wǎng)站導(dǎo)航、Google

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請(qǐng)盡快告知，我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng)，如需處理請(qǐng)聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時(shí)需注明來源： 創(chuàng)新互聯(lián)