黃金分割法、斐波那契法求極值_第1頁
黃金分割法、斐波那契法求極值_第2頁
黃金分割法、斐波那契法求極值_第3頁
全文預(yù)覽已結(jié)束

下載本文檔

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

文檔簡介

1、、用黃金分割法求極值:;%建立原方程函數(shù);function y=fx(x)if nargin=1 y=x+20/x;endend%a為區(qū)間下限,b為區(qū)間上限,e為精度;%fx(x)為原方程函數(shù);function xj,yj=huangjin(a,b,e)a=input(Please enter the value of a:);b=input(Please enter the value of b:);e=input(Please enter the value of e:);while b-ae x1=a+0.382*(b-a); x2=a+0.618*(b-a); if fx(x1)fx(

2、x2) b=x2; else a=x1; endend xj=(a+b)/2; yj=fx(xj);end此題中,a=0.2,b=1,e=0.01,程序運行結(jié)果為xj=0.9967,yj=21.0620。改變方程,將原函數(shù)方程改變。6、用斐波那契法求極值:;%建立原方程函數(shù);function y=fx(x)if nargin=1 y=3*x4-16*x3+30*x2-24*x+8;endend%建立斐波那契數(shù)列函數(shù);function fn=F(n)if nargin=1 y=0; y(1)=1; y(2)=1; for k=3:1:n y(k)=y(k-1)+y(k-2); end fn=y(n);endend%求解應(yīng)計算次數(shù)的函數(shù);%s為(b-a)/e的值,其中(a,b)為單峰區(qū)間,e為精度;function n=cishu(s)if nargin=1 n=1; while F(n)e if fx(x2)=fx(x1) b=x2; x2=x1; x1=a+b-x2; else a=x1; x1=x2; x2=a+b-x1; endendxj=(a+b)/2;yj=fx(xj);end此題中,a=-10,b=10,e=0.01,程序運行結(jié)果為:xj =

溫馨提示

  • 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)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論