《matlab課程設(shè)計》word版-文庫吧
2025-05-18 14:22 本頁面
【正文】 x0+[0:(PointNum1)]39。*h。 %得 x向量值 輸入求解的自變量范圍 求出待求簡單微分方程的真值解 用 MATLAB 自帶函數(shù) ode23求解待求微分方程 結(jié)束 用自編函數(shù)四階龍格－庫塔（ RK）方法求解待求微分方程 開始 for k=1:(PointNum) %迭代計算 f1=h*feval(fun,x(k),y(k,:),varargin{:})。 f1=f1(:)39。 %得公式 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。,
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1