具體內容如下
fromsympyimport*importnumpyasnpdefbacktracking_line_search(f,df,x,x_k,p_k,alpha0):rho=0.5c=10**-4alpha=alpha0replacements1=zip(x,x_k)replacements2=zip(x,x_k+alpha*p_k)f_k=f.subs(replacements1)df_p=np.dot([df_.subs(replacements1)fordf_indf],p_k)whilef.subs(replacements2)>f_k+c*alpha*df_p:alpha=rho*alphareplacements2=zip(x,x_k+alpha*p_k)returnalphadefstepest_line_search(f,x,x0,alpha0):df=[diff(f,x_)forx_inx]x_k=x0alpha=alpha0replacements=zip(x,x_k)len_df=sqrt(np.sum([df_.subs(replacements)**2fordf_indf]))whilelen_df>1e-6:p_k=-1*np.array([df_.subs(replacements)fordf_indf])alpha=backtracking_line_search(f,df,x,x_k,p_k,alpha)x_k=x_k+alpha*p_kreplacements=zip(x,x_k)len_df=np.sum([df_.subs(replacements)**2fordf_indf])returnx_kif__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)**2ans=stepest_line_search(f,x,np.array([1.2,1]),1)print"theminimalvalueinpoint:",ans分析:
這個采用的是backtracking line search來尋找alpha。
本文由 貴州做網站公司 整理發布，部分圖文來源于互聯網，如有侵權，請聯系我們刪除，謝謝！
c語言中正確的字符常量是用一對單引號將一個字符括起表示合法的字符常量。例如‘a’。數值包括整型、浮點型。整型可用十進制，八進制，十六進制。八進制前面要加0，后面...
2022年天津專場考試原定于3月19日舉行，受疫情影響確定延期，但目前延期后的考試時間推遲。 符合報名條件的考生，須在規定時間登錄招考資訊網(www.zha...
：喜歡聽，樂意看。指很受歡迎?！巴卣官Y料”喜聞樂見：[ xǐ wén lè jiàn ]詳細解釋1. 【解釋】：喜歡聽，樂意看。指很受歡迎。2. 【示例】：這是...
說了這么多次的房地產稅終于要來了，據最新消息稱，部分地區將開展房地產稅改革試點，那么具體會有哪些相關的措施推出呢?房地產稅一旦實施將會帶來哪些影響?下面我們來了解一下。房地產稅是政府對房地產征收的一種財產稅。它通常向業主或租戶等用戶征收。負責征收房地產稅的政府機構將對房地產價值進行估價，并以房地產價值的一定比例作為應繳稅款。為積極穩妥推進房地產稅立法改革，引導合理住房消費和節約集約利用土地資源，促...
在金融危機期間，失業率普遍會急劇.上升，因此保住工作則是應對金融危機最好的辦法。有了工作，才會有穩定的薪酬收入，因此可以支付日常生活中最重要的開銷。在金融危機時期，房地產市場、股票市場和期貨市場都處于歷史最低點，在這一階段擁有現金資產就是最具有保障的。以下是如何應對金融危機的相關內容。伴隨著如今鮑威爾已經承認美聯儲正在討論如何縮減購債規模，期權和債市交易市場已經發出了市場預警信號。自金融危機爆發以...
提起余額寶，相信對于很多人而言都不會感到陌生，全國至少有幾千萬人使用余額寶，余額寶的規模將近1.5萬億，直追國內四大銀行的存款規模。余額寶始終不是銀行，那么存進余額寶里面的錢會不會被盜走?下面來了解下。余額寶是貨幣基金，在實際生活中，貨幣基金的安全性還是非常高的，基本上不會出現本金虧損的情況。貨幣基金專門投向風險小的貨幣市場工具，區別于其他類型的開放式基金，具有高安全性、高流動性、穩定收益性，具有...
免責聲明：本站內容部分文字和圖片來源于網絡，版權屬于原作者所有，若有冒犯，請聯系我們進行刪除，謝謝。聯系郵箱：1051388971#qq.com （把#換成@）