matlab在数值分析中的应用 Runge-kutta幻灯片.pptVIP

求解： x0=[1.2; 0; 0; -1; tic, [t,y]=ode45(apolloeq,[0,20],x0); toc elapsed_time = 0.8310 length(t), plot(y(:,1),y(:,3)) ans = 689 得出的轨道不正确， 默认精度RelTol设置 得太大，从而导致的 误差传递，可减小该 值。 改变精度： options=odeset; options.RelTol=1e-6; tic, [t1,y1]=ode45(apolloeq,[0,20],x0,options); toc elapsed_time = 0.8110 length(t1), plot(y1(:,1),y1(:,3)), ans = 1873 min(diff(t1)) ans = 1.8927e-004 plot(t1(1:end-1),… diff(t1)) 例： x0=[1.2; 0; 0; -1; tic, [t1,y1]=rk_4(apolloeq,[0,20,0.01],x0); toc elapsed_time = 4.2570 plot(y1(:,1),y1(:,3)) % 绘制出轨迹曲线