數(shù)值分析第五版計算實習(xí)題(編輯修改稿)
2025-04-21 02:50 本頁面
【文章內(nèi)容簡介】 n(i+1)=yn。endXn=Yn。%以下根據(jù)經(jīng)典四階RK法公式求解for i=1:n xn=Xn(i)。 yn=Yn(i)。 k1=f(xn,yn)。 k2=f(xn+h/2,yn+h/2*k1)。 k3=f(xn+h/2,yn+h/2*k2)。 k4=f(xn+h,yn+h*k3)。 yn=yn+h/6*(k1+2*k2+2*k3+k4)。 Yn(i+1)=yn。enddisp(39。 改進(jìn)歐拉法 四階經(jīng)典RK法39。)。disp([Xn39。 Yn39。])結(jié)果如下： 改進(jìn)歐拉法 四階經(jīng)典RK法 1 1 1 （b）程序：clc。clear。a=0。b=1。%定義域H=[ ]。%步長y0=1/3。%初值f= @(x,y) 50*y+50*x^2+2*x。%微分函數(shù)xi=linspace(a,b,11)。Y=1/3*exp(50*xi)+xi.^2。%準(zhǔn)確解Ym=zeros(1,11)。for j=1:3 h=H(j)。 n=(ba)/h。 Xn=linspace(a,b,n+1)。%將定義域分為n等份 Yn=zeros(1,n)。%結(jié)果矩陣 Yn(1)=y0。%賦初值 for i=1:n xn=Xn(i)。 yn=Yn(i)。 k1=f(xn,yn)。 k2=f(xn+h/2,yn+h/2*k1)。 k3=f(xn+h/2,yn+h/2*k2)。 k4=f
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1