python實(shí)現(xiàn)最速下降法
更新時(shí)間：2020年03月24日 14:22:49 作者：zcc_TPJH
這篇文章主要為大家詳細(xì)介紹了python實(shí)現(xiàn)最速下降法，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下
本文實(shí)例為大家分享了python實(shí)現(xiàn)最速下降法的具體代碼，供大家參考，具體內(nèi)容如下
代碼：
from sympy import *
import numpy as np
def backtracking_line_search(f,df,x,x_k,p_k,alpha0):
rho=0.5
c=10**-4
alpha=alpha0
replacements1=zip(x,x_k)
replacements2=zip(x,x_k+alpha*p_k)
f_k=f.subs(replacements1)
df_p=np.dot([df_.subs(replacements1) for df_ in df],p_k)
while f.subs(replacements2)>f_k+c*alpha*df_p:
alpha=rho*alpha
replacements2 = zip(x, x_k +alpha * p_k)
return alpha
def stepest_line_search(f,x,x0,alpha0):
df = [diff(f, x_) for x_ in x]
x_k=x0
alpha=alpha0
replacements=zip(x,x_k)
len_df = sqrt(np.sum([df_.subs(replacements) ** 2 for df_ in df]))
while len_df>1e-6:
p_k=-1*np.array([df_.subs(replacements) for df_ in df])
alpha = backtracking_line_search(f, df, x, x_k, p_k, alpha)
x_k=x_k+alpha*p_k
replacements = zip(x, x_k)
len_df=np.sum([df_.subs(replacements)**2 for df_ in df])
return x_k
if __name__=="__main__":
init_printing(use_unicode=True)
x1 = symbols("x1")
x2 = symbols("x2")
x = np.array([x1, x2])
f = 100 * (x2 - x1 ** 2)**2 + (1 - x1) ** 2
ans=stepest_line_search(f, x, np.array([1.2, 1]), 1)
print "the minimal value in point:",ans
分析：
這個(gè)采用的是backtracking line search來(lái)尋找alpha。
以上就是本文的全部?jī)?nèi)容，希望對(duì)大家的學(xué)習(xí)有所幫助，也希望大家多多支持腳本之家。
您可能感興趣的文章:
相關(guān)文章
python中sqllite插入numpy數(shù)組到數(shù)據(jù)庫(kù)的實(shí)現(xiàn)方法
本文給大家介紹python中sqllite插入numpy數(shù)組到數(shù)據(jù)庫(kù)的實(shí)現(xiàn)方法，在文章底部給大家提到了Python 操作sqlite數(shù)據(jù)庫(kù)及保存查詢numpy類型數(shù)據(jù)的實(shí)例代碼，需要的朋友參考下吧2021-06-06
這篇文章主要介紹了python實(shí)現(xiàn)ip查詢示例,需要的朋友可以參考下2014-03-03
Python的Scrapy爬蟲(chóng)框架簡(jiǎn)單學(xué)習(xí)筆記
這篇文章主要介紹了Python的Scrapy爬蟲(chóng)框架簡(jiǎn)單學(xué)習(xí)筆記,從基本的創(chuàng)建項(xiàng)目到CrawlSpider的使用等都有涉及,需要的朋友可以參考下2016-01-01
這篇文章主要介紹了python爬取12306驗(yàn)證碼，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2019-05-05
Python實(shí)現(xiàn)某論壇自動(dòng)簽到功能
這篇文章主要介紹了Python實(shí)現(xiàn)論壇自動(dòng)簽到功能，本文通過(guò)實(shí)例代碼給大家介紹的非常詳細(xì)，具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2019-08-08
這篇文章主要介紹了Python可視化神器pyecharts繪制餅狀圖，餅圖是用圓形及圓內(nèi)扇形的角度來(lái)表示數(shù)值大小的圖形，它主要用于表示一個(gè)樣本中各組成部分的數(shù)據(jù)占全部數(shù)據(jù)的比例2022-07-07
python利用thrift服務(wù)讀取hbase數(shù)據(jù)的方法
今天小編就為大家分享一篇python利用thrift服務(wù)讀取hbase數(shù)據(jù)的方法，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2018-12-12
昨晚我用python幫隔壁小姐姐P證件照然后發(fā)現(xiàn)
大家好,我是Lex 喜歡欺負(fù)超人那個(gè)Lex 建議大家收藏哦,以后幫小姐姐P自拍,證件照,調(diào)尺寸,背景,摳圖,直接10行代碼搞定,瞬間高大上2021-08-08
這篇文章主要為大家詳細(xì)介紹了python繪制多個(gè)曲線的折線圖，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2018-09-09
基于Python實(shí)現(xiàn)牛牛套圈小游戲的示例代碼
“幸運(yùn)牛牛套圈圈”套住歡樂(lè)，圈住幸福，等你來(lái)挑戰(zhàn)！這篇文章小編主要為大家介紹一款基于Python實(shí)現(xiàn)牛牛套圈小游戲，感興趣的小伙伴可以了解一下2023-02-02