一維搜索牛頓法_第1頁
一維搜索牛頓法_第2頁
一維搜索牛頓法_第3頁
一維搜索牛頓法_第4頁
一維搜索牛頓法_第5頁
已閱讀5頁，還剩3頁未讀， 繼續(xù)免費(fèi)閱讀

下載本文檔

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

文檔簡介

1、2013-2014 (1)專業(yè)課程實(shí)踐論文題目：一維搜索牛頓法一、算法理論牛頓法是一種函數(shù)逼近法，基本思想是：在極小點(diǎn)附近用 函數(shù)的二階泰勒多項(xiàng)式近似代替目標(biāo)函數(shù)，從而求得目標(biāo)函數(shù) 的極小點(diǎn)的近似值。對f (x)在xk點(diǎn)二階泰勒展開：f (x) = f(X ) + )(X f ) + 2 f (xk )(X - Xk )2 + 心-Xk )2) 略去高階項(xiàng)得f (x)牝 f (x)+*)(xx)+2Ek)(x一xk)2兩邊對 x 求導(dǎo)，f (x)機(jī) f (x ) + f (x )(x 一 x )令 f (x)=0 ，得到f x ) f (x ) kx=xk+1kf x )1f (x )k作為新

2、的迭代點(diǎn)，繼續(xù)迭代，直到達(dá)到精度，這樣就得到了函數(shù)f的一個駐點(diǎn)。以上過程即Newton法。二、算法框圖停，失敗X砰1七甘=妙TL) /S)三、算法程序#include#include using namespace std;double fun(double t)return (t*t*t-2*t+1);double dfun(double t)return (3*t*t-2);void NewtonIterative(double(*pfun)(double t),double (*pdfun)(double t) int maxflag=1000,k=1;double t0=1,err=0

3、.01,t;dot0=t;t=t0-pfun(t0)/pdfun(t0);k+;while(fabs(t-t0)err&k=maxflag)coutendl牛頓迭代失敗!迭代次數(shù)為k=kendl;elsecoutendl牛頓迭代成功!迭代次數(shù)為k=k-1endl;cout迭代結(jié)果：近似根為root=tendl;cout函數(shù)值近似為：f(root)=fun(t)endl;int main()NewtonIterative(fun,dfun);return 0;四、算法實(shí)現(xiàn)例1.試用Newton法求函數(shù) f (x) = x 4 - 4 x3 - 6 x 2 -16 x + 4的最優(yōu)解。Cx = 6,8 = 10-2 ) 0解：運(yùn)行程序 f(x) = 4x3 -12x2 -12x -16顯示結(jié)果運(yùn)行程序f(x) = 3 x 2 - 2 顯示結(jié)果。例3.試用 Newton法求函數(shù) f (x) = 4 x 3-12 x 2-12 x -16的最優(yōu)解。Cx = 6,

溫馨提示

  • 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)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論