练习数值解

练习一：向前欧拉法解方程

y′(t)y(0)=−2y(t)+t=1

$$\begin{aligned} y'(t)&=-2y(t)+t\\ y(0)&=1 \end{aligned}$$
算法解释：对于一般的常微分方程
$$\begin{aligned} y'(t)&=f(t,y)\\ y(t_0)&=y_0 \end{aligned}$$
$\displaystyle \frac{y_{i+1}-y_i}{\Delta t}=f(t_i,y_i) \Rightarrow y_{i+1}=y_i+f(t_i,y_i)\Delta t$

clear;
%方程的系数
a=1.;
b=-2.;
c=1.;
%开始结束时间
tinit=0.;
tmax=5.;
%时间步数
maxt=3000;
dt=(tmax-tinit)/maxt;
%初始条件
y(1)=1.;
t(1)=tinit;
%时间循环（loop）
for j=1:maxt;
y(j+1)=y(j)+dt*(b*y(j)+c*t(j))/a;
t(j+1)=tinit+j*dt;
end;
%绘图
plot(t,y)
title('欧拉向前差分')
xlabel('t')
ylabel('y(t)')

练习二:中点法(精度比向前向后欧拉高一阶)

yi+1=yi−1+2f(ti,yi)Δt

$$y_{i+1}=y_{i-1}+2f(t_i,y_i)\Delta t$$
有点疑惑的是：初值需要给在$i=0,1$两处；可否用向后欧拉法$y_{i+1}=y_i+f(t_{i+1},y_{i+1})\Delta t$求出另外一个？比如$y_0=y_1-f(t_1,y_1)\Delta t$，但是这样的话就会有一个问题，向后欧拉法精度低一阶。