matlab最優(yōu)化牛頓法_第1頁
matlab最優(yōu)化牛頓法_第2頁
matlab最優(yōu)化牛頓法_第3頁
matlab最優(yōu)化牛頓法_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進行舉報或認(rèn)領(lǐng)

文檔簡介

1、0實驗報告日期：2oii年12月2日、實驗概述實驗名稱】：牛頓法實驗性質(zhì)】：驗證性實驗?zāi)康募耙蟆浚号浜险n堂教學(xué)，培養(yǎng)學(xué)生動手能力，根據(jù)牛頓法求極小值的思想設(shè)計程序?；驹怼浚号ｎD法的迭代公式：x(k+1)，x(k)-2f(xk)-1f(xk)其中2kf(Xk)是f(x)在Xk處的Hesse矩陣。實施環(huán)境】：MATLAB7.7.0二、實驗內(nèi)容【項目內(nèi)容及要求】用牛頓法求解以下問題min(x1-1)4+x22三、實驗過程【實驗操作步驟】functionx1k=newton(x1,j)%x1為初始點x1=8,8;j=1e-10;hs=inline(x-1)A4+yA2);ezcontour(hs

2、,-1010-1010);holdon;symsxyf=(x-1)A4+yA2;gradl=jacobian(f,x,y);%求梯度grad2=jacobian(gradl,x,y);%求Hesse矩陣k=0;whilelgradlz=subs(subs(gradl,x,xl(l),y,xl(2);%求梯度值grad2z=subs(subs(grad2,x,xl(l),y,xl(2);%求Hesse矩陣x2=xl-inv(grad2z)*(gradlz);%牛頓迭代公式ifnorm(xl-x2)jbreak;elseplot(xl(l),x2(l),xl(2),x2(2),-r*);k=k+l;xl=x2;endend【實驗結(jié)果】xl=8,8;j=le-l0;xlk=newton(xl,j)xl=l.0000指導(dǎo)

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論