數(shù)值分析大作業(yè)-三、四、五、六、七(參考版)
2025-06-27 14:02本頁面
【正文】 程序執(zhí)行結果如下： adams4請輸入插值節(jié)點數(shù)n=3000。Adams預測校正公式的誤差39。)。plot(x,T,39。T=yf0。Adams預測校正公式法39。)。plot(x,y,39。)。title(39。k39。endsubplot(3,1,1)。 pre=pre1。 f2=f3。 %用局部截斷誤差進一步修正校正值 f=f1。 %用局部截斷誤差進一步修正預測值 mod1=y(j1)+h/24*(9*(l+2*cos(x(j)))+19*f5*f1+f2)。for j=k+1:n+1 pre1=y(j1)+h/24*(55*f59*f1+37*f29*f3)。 %fn3pre=0。 %fn1f2=y(k2)+2*cos(x(k2))。endf=y(k)+2*cos(x(k))。 k4=h*((y(k1)+k3)+2*cos(x(k1)+h))。 k2=h*((y(k1)+k1/2)+2*cos(x(k1)+h/2))。 %f0=cosx+sinx為精確值endy(1)=1。x=0:h:pi。)。程序執(zhí)行結果如下： classic_rk4請輸入插值節(jié)點數(shù)n=3000（2）用四階Adams預測校正算法求解，程序代碼如下：function adams4()n=input(39。經(jīng)典四階RK法的誤差39。)。 %計算經(jīng)典四階RK法的誤差plot(x,T,39。subplot(3,1,3)。經(jīng)典四階RK法39。)。plot(x,y,39。)。title(39。k39。 endsubplot(3,1,1)。 k4=(y(k)+h*k3)+2*cos(x(k)+h)。 k2=(y(k)+h*k1/2)+2*cos(x(k)+h/2)。for k=1:n f0(k+1)=cos(x(k))+sin(x(k))。k=2。 %f0=cosx+sinx為精確值h=pi/n。y(1)=1。請輸入插值節(jié)點數(shù)n=39。, R(2^(i1)))。endfprintf(39。while flag==1if abs(R(2^i)R(2^(i1)))/255eps flag=0。end%========計算滿足允許精度的Rn，并打印輸出========%i=1。end%========計算Cn========%function Cn=C(n)Cn=16/15*S(2*n)1/15*S(n)。endfor j=1:n Tn=Tn+h*(fx(x(j))+fx(x(j+1)))/2。x=zeros(1,n+1)。%========計算Tn========%function Tn=T(n)Tn=0。請輸入允許精度：eps=39。)。b=input(39。請輸入積分下限：a=39。end%Romberg求積法計算積分的通用程序function Romberg()clear。) %用藍色o線畫出x,b3圖像%試驗結果為：a1 = Columns 1 through 2 Columns 3 through 4 a2 = Columns 1 through 2 Columns 3 through 4 Columns 5 through 6 Columns 7 through 8 Columns 9 through 10 a3 = Columns 1 through 2 Columns 3 through 4 Columns 5 through 6 Columns 7 through 8 Columns 9 through 10 Columns 11 through 12 Columns 13 through 14 Columns 15 through 16 b1 = Columns 1 through 2 Columns 3 through 4 Columns 5 through 6 Columns 7 through 8 Columns 9 through 10 Columns 11 through 12 Columns 13 through 14 Columns 15 through 16 Columns 17 through 18 Columns 19 through 20 Columns 21 through 22 Columns 23 through 24 Columns 25 through 26 Columns 27 through 28 Columns 29 through 30 b2 = Columns 1 through 2 Columns 3 through 4 Columns 5 through 6 Columns 7 through 8 Columns 9 through 10 Columns 11 through 12 Columns 13 through 14 Columns 15 through 16 Columns 17 through 18 Columns 19 through 20 Columns 21 through 22 Columns 23 through 24 Columns 25 through 26
畢業(yè)設計相關推薦
鄂ICP備17016276號-1