正文內容

數值微分計算方法實驗-資料下載頁

2025-12-28 06:50本頁面

　　

【正文】 E(1)=0。E(2)=abs(D(2)D(1))。 %計算導數近似值間隔序列的前兩個值 R(1)=1。R(2)=2*E(2)/(abs(D(1))+abs(D(2)))。 %輸出滿足精度要求的數值導數近似值 for i=3:max %計算導數近似值序列 h=h/10。 D(i)=(feval(f,x+2*h)+8*feval(f,x+h)8*feval(f,xh)+feval(f,x2*h))/(12*h)。 E(i)=abs(D(i)D(i1))。 R(i)=2*E(i)/(abs(D(i))+abs(D(i1)))。 if ((E(i)E(i1))amp。(R(i)=toler)) break。 end end n=i。 df=D(n)。 %輸出滿足精度要求的數值導數近似值 err=E(n)。 relerr=R(n)。 end 理查森外推法計算導數近似值的子函數 function df=Richarson(f,x,toler,delta,max,h0) %輸入輸出參數說明： %df 輸出的近似導數值 %f、 x 相應函數及其點 delta 相對誤差容忍上限， toler 誤差容忍上限， max 最大計算次數， h0 初始步長 h=h0。 D(1,1)=(feval(f,x+h)feval(f,xh))/(2*h)。 err=1。《數值方法》實驗報告 11 relerr=1。 j=1。 while((errdelta)amp。(relerrtoler)amp。(j=max))%進行是否終止判斷 h=h/2。 D(j+1,1)=(feval(f,x+h)feval(f,xh))/(2*h)。 for k=1:j %理查森外推公式計算高階精度的導數計算公式 D(j+1,k+1)=D(j+1,k)+(D(j+1,k)D(j+1,k))/(4^(k)1)。 end err=abs(D(j+1,j+1)D(j,j))。 relerr=2*err/(abs(D(j+1,j+1))+abs(D(j,j)))。 j=j+1。 end n=j df=D(n,n)。%輸出得到滿足得到精度要求的近似導數值 err relerr End 1 function df=difnew(X,Y) N=length(X)1。 A=Y。 for M=1:N+1 %對 X 序列重新排列得到新的插值點序列 if M==1 t=X。 else if M==N+1 for i=2:N+1 t(i)=X(i1)。 end else t(1)=X(M)。 for i=2:M t(i)=X(i1)。 end for i=M+1:N+1 t(i)=X(i)。 end end end for i=1:N+1 %依據重新排列的插值點序列，計算差商表以得到牛頓插《數值方法》實驗報告 12 值多項式 for i=j+1:N+1 A(i)=(A(i)A(i1))/(t(i)t(ij+1))。 end p=1。 df(M)=A(2)。 for k=2:N %累加求和求解插值點導數值 df(M)=df(M)+p*(t(1)t(k))*A(k+1)。 end end end end

點擊復制文檔內容

公司管理相關推薦

數值微分計算方法實驗-資料下載頁

數值計算方法(第4章)-資料下載頁

數值計算方法ppt課件(2)-資料下載頁

數值計算方法(第2章)-資料下載頁

數值計算方法(第3章)(1)-資料下載頁

計算方法課件劉師少第一章數值計算方法-資料下載頁

數值計算方法chap1誤差-資料下載頁

計算方法數值分析ppt課件-資料下載頁

數值計算方法課程設計(c語言)-資料下載頁

數值計算方法課程設計(c語言)-資料下載頁

線性方程組ax=b的數值計算方法實驗-資料下載頁

數值計算方法第七章-資料下載頁

微積分的數值計算方法romberg算法-資料下載頁

數值與計算方法第1章緒論-資料下載頁

數值計算方法(宋岱才版)課后答案-資料下載頁

數值分析計算方法課程介紹ok-資料下載頁

數值微分計算方法實驗-全文預覽

數值微分計算方法實驗-預覽頁

數值微分計算方法實驗-免費閱讀

數值微分計算方法實驗(存儲版)

數值微分計算方法實驗-文庫吧在線文庫