% 函数像 ode45 一样求解常微分方程。 % dydt = 一阶微分方程组% trange = 自变量的范围% yinit = 初始条件%tol =误差绝对公差,默认值= 1e-6; % 例子: % [t,y] = RK45(@(t,y)[y(2); (1-y(1)^2)*y(2)-y(1)],[0 20],[2 ; 0]); % 子图 (1,2,2); % plot(t,y,'-o') % title('van der Pol 方程 (\mu = 1) with RK45'); % xlabel('时间 t'); % ylabel('解决方案 y'); %图例('y_1','y_2'); % %[t2,y2] = ode45(@(t,y)[y(2);(1-y(1)^ 2)* y(2)-y(1)],[0 20],[2 ; 0]); % 子图 (1,2,1); % plot(t2,y2,
2021-10-30 23:11:08
2KB
matlab
1