《matlab課程設(shè)計》word版-文庫吧
2025-01-01 07:28 本頁面
【正文】 %得公式k1 f2=h*feval(fun,x(k)+h/2,y(k,:)+f1/2,varargin{:})。 f2=f2(:)39。 %得公式k2 f3=h*feval(fun,x(k)+h/2,y(k,:)+f2/2,varargin{:})。 f3=f3(:)39。 %得公式k3 f4=h*feval(fun,x(k)+h,y(k,:)+f3,varargin{:})。 f4=f4(:)39。 %得公式k4 y(k+1,:)=y(k,:)+(f1+2*(f2+f3)+f4)/6。 %得y(n+1)end、實例求解源程序：%運行四階RK法clear, clc %清除內(nèi)存中的變量x0=0。xt=2。Num=100。h=(xtx0)/(Num1)。x=x0+[0:Num]*h。a=1。yt=1exp(a*x)。 %真值解fun=inline(39。y+139。,39。x39。,39。y39。)。 %用inline構(gòu)造函數(shù)f(x,y)y0=0。 %設(shè)定函數(shù)初值PointNum=5。 %設(shè)定取點數(shù)[x1,y1]=ode23(fun,[0,2],0)。[xr,yr]=MyRunge_Kutta(fun,x0,xt,y0,PointNum)。MyRunge_Kutta_x=xr39。MyRunge_Kutta_y=yr39。plot(x,yt,39。k39。,x1,y1,39。b39。,xr,yr,39。r39。)legend(39。真值39。,39。ode2339。,39。RungKutta法解39。,2)hold onplot(x1,y1,39。bo39。,xr,yr,39。r*39。)程序運行結(jié)果：MyRunge_Kutta_x = 0 MyRunge_Kutta_y = 0 二、變成解決以下科學(xué)計算問題：（
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1