用下列方法求f_第1頁
用下列方法求f_第2頁
用下列方法求f_第3頁
用下列方法求f_第4頁
用下列方法求f_第5頁
已閱讀5頁，還剩2頁未讀， 繼續(xù)免費閱讀

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

用下列方法求f(x)=x3-3x-1=0在x0=2附近的根，根的精確值x*=1.87938524......要求結果精確到四位有效數(shù)字。牛頓法+埃特金加速+斯蒂芬森弦截法+埃特金加速+斯蒂芬森二分法基本迭代法本次試驗我采用函數(shù)的方式來實現(xiàn)上式要求，函數(shù)程序如下：牛頓法function[x,k]=newton(f,df,x0,ep,N)%f為原函數(shù)%df為導函數(shù)ifnargin<5,N=500;endifnargin<4,ep=1e-3;endk=0;whilek<Nx=x0-feval(f,x0)/feval(df,x0);ifabs(x-x0)<epbreak;endx0=x;k=k+1;end弦切法function[x,k]=xianqie(f,x0,x1,ep,N)%f為原函數(shù)%x0,x1為兩個初值ifnargin<5,N=500;endifnargin<4,ep=1e-3;endk=0;whilek<Nx=x1-(x1-x0)*feval(f,x1)/(feval(f,x1)-feval(f,x0));ifabs(x-x1)<ep,break;endx0=x1;x1=x;k=k+1;end二分法function[x,k]=TwoDiv(f,a,b,ep)%f為原函數(shù)%a和b為x軸的區(qū)間x=(a+b)/2.0;k=0;whileabs(feval(f,x))>ep||(b-a>ep)iffeval(f,x)*feval(f,a)<0b=x;elsea=x;endx=(a+b)/2.0;k=k+1;end基本迭代法function[x,k]=diedai(phi,x0,ep,N)%phi為迭代式%ep為精度%N為迭代次數(shù)ifnargin<4,N=40;endifnargin<3,ep=1e-4;endk=0;whilek<Nx=feval(phi,x0);ifabs(x-x0)<epbreak;endx0=x;k=k+1;end計算結果分別為：、在命令行分別輸入以下語句：下面的語句表示二分法f=@(x)x^3-3*x-1;[x,k]=TwoDiv(f,1,2,1e-3)結果為分析結果：在命令行分別輸入以下語句：分別表示基本迭代法、迭代法的埃特金加速、迭代法的斯特芬森加速（a）、phi=@(x)(1+3*x)^(1/3);[x,k]=diedai(phi,2,1e-3,40)（b）、phi=@(x)x^3-3*x-1;>>[x,k]=aitken(phi,2,1e-3,20)（c）、p

溫馨提示

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

評論

0/150

提交評論