%% 定义计算的步长, 设置变量的初始值
clear;clc;close all
Delta=0.001; % 定义步长
t=0:Delta:20; % 定义自变量 t
n=length(t);
Y(:,1)=[2;0]; % 定义 x y 的初始值
%% 自定义欧拉法, 求解微分方程组
for k=1:n-1
Y(:,k+1)=Y(:,k)+Delta*f(t(k),Y(:,k));
end
x=Y(1,:);
y=Y(2,:);
%% 绘制 x y 的求解结果
figure
set(gcf,'units','normalized','position',[0.15 0.2 0.7 0.6]); % 设置 figure 窗口的位置和尺寸
subplot(1,2,1)
plot(t,x,'b')
xlabel('t')
ylabel('x')
subplot(1,2,2)
plot(t,y,'r')
xlabel('t')
ylabel('y')