三阶、四阶龙格库塔函数matlab代码.doc
三阶龙格库塔法的计算公式为 三阶龙格库塔公式的Matlab程序代码function y DELGKT3_kutaf, h,a,b,y0,varvecat long;N b-a/h;y zerosN1,1;y1 y0;x ahb;var findsymf;for i2N1 K1 Funvalf,varvec,xi-1 yi-1; K2 Funvalf,varvec,xi-1h/2 yi-1K1*h/2; K3 Funvalf,varvec,xi-1h yi-1-h*K1K2*2*h; yi yi-1h*K14*K2K3/6;endat short;DELGKT3_kuta函数运行时需要调用下列函数 function fvFunvalf, varvec, varval var findsymf; if lengthvar4 if var1varvec1 fvsubsf,varvec1,varval1; else fvsubsf,varvec2,varval2; end else fvsubsf,varvec,varval; end三阶龙格库塔求解一阶常微分方程应用实例。用三阶龙格库塔法求下面常微分方程的数值解。 在编辑窗口输入下列程序段,然后执行该程序。syms x y;z2*x-3*y2;yyDELGKT3_kutaz,0.1,0,1,1,x y程序执行后得结果 四阶龙格库塔法的计算公式为 function y DELGKT4_lungkutaf, h,a,b,y0,varvecat long;N b-a/h;y zerosN1,1;y1 y0;x ahb;var findsymf;for i2N1 K1 Funvalf,varvec,xi-1 yi-1; K2 Funvalf,varvec,xi-1h/2 yi-1K1*h/2; K3 Funvalf,varvec,xi-1h/2 yi-1K2*h/2; K4 Funvalf,varvec,xi-1h yi-1h*K3; yi yi-1h*K12*K22*K3K4/6;endat short; 同理也要先调用Funval函数(见上)四阶龙格库塔求解一阶常微分方程应用实例。用四阶龙格库塔法求下面常微分方程的数值解。 在编辑窗口输入下列程序段,然后执行该程序。z1logx1yyDELGKT4_lungkutaz,0.1,0,1,1,x,y程序执行后得结果