常微分方程描述自然现象。物体冷却过程符合牛顿冷却定律。一杯热水放在房间里。热水温度逐渐下降。房间温度保持稳定。温度变化率与温差成正比。温差大温度下降快。温差小温度下降慢。这是一个典型的一阶常微分方程。

建立数学模型。设热水温度为T。房间温度为T0。温度变化率是dT/dt。根据牛顿冷却定律得到方程dT/dt=-k(T-T0)。k是一个正常数。负号表示温度下降。这个方程包含未知函数T及其导数。这是一个常微分方程。

求解这个方程。分离变量法适合这个方程。将方程写成dT/(T-T0)=-kdt。两边同时积分。左边积分得到ln|T-T0|。右边积分得到-kt C。C是积分常数。于是有ln|T-T0|=-kt C。取指数运算得到T-T0=e^{-kt C}。令A=e^C。得到通解T=T0 Ae^{-kt}。

确定常数A。初始时刻t=0热水温度为T1。代入方程得到T1=T0 A。所以A=T1-T0。最终解为T=T0 (T1-T0)e^{-kt}。这个解描述温度变化过程。温度从T1开始。指数衰减到T0。k值决定衰减速度。k大衰减快。k小衰减慢。

这个模型可以验证。实际测量热水温度。记录不同时刻温度值。与理论计算比较。确定k的数值。实验结果与理论吻合。证明模型正确。这个模型有很多应用。法医推断死亡时间。尸体温度随时间下降。测量尸体温度。估计死亡时间。工业生产需要控制温度。这个模型提供理论依据。

种群增长也符合微分方程。兔子在草原上繁殖。没有天敌食物充足。兔子数量增长很快。数量变化率与现有数量成正比。设兔子数量为N。变化率是dN/dt。建立方程dN/dt=rN。r是增长率。这个方程很简单。解这个方程。分离变量得到dN/N=rdt。两边积分得到lnN=rt C。取指数得到N=e^{rt C}。令A=e^C。得到N=Ae^{rt}。初始数量为N0。代入得到A=N0。最终解为N=N0e^{rt}。这个解表示指数增长。实际增长不会无限。环境资源有限。种群竞争加剧。增长会变慢。

改进这个模型。考虑环境容纳量K。当N接近K时增长变慢。变化率与N和(K-N)成正比。建立方程dN/dt=rN(1-N/K)。这个方程称为逻辑斯蒂方程。这个方程非线性。求解比较复杂。分离变量得到dN/[N(1-N/K)]=rdt。左边可以分解为(1/N 1/(K-N))dN。两边积分得到ln|N|-ln|K-N|=rt C。整理得到ln|N/(K-N)|=rt C。取指数得到N/(K-N)=e^{rt C}。令A=e^C。得到N/(K-N)=Ae^{rt}。解出N得到N=K/(1 Be^{-rt})。B是常数。初始条件确定B。这个解描述S形曲线。开始增长慢。然后加快。最后变慢。接近环境容纳量。这个模型更符合实际。可以预测种群数量。管理野生动物。控制害虫数量。

弹簧振动也用到微分方程。弹簧挂着一个物体。拉一下物体振动。建立坐标系。平衡位置为原点。物体位置用x表示。根据牛顿第二定律。物体质量m乘加速度等于合力。弹簧力与位移成正比。方向相反。弹簧力是-kx。k是弹性系数。空气阻力与速度成正比。阻力是-cdx/dt。c是阻力系数。得到方程md²x/dt²=-kx-cdx/dt。整理得到md²x/dt² cdx/dt kx=0。这是一个二阶常微分方程。

求解这个方程。假设解为x=e^{rt}。代入方程得到特征方程mr² cr k=0。解这个二次方程。根的情况决定振动形式。根为复数时系统振动。根为实数时系统不振动。阻力很小的时候系统振动。阻力很大的时候系统慢慢回到平衡位置。

无阻力情况比较简单。c=0方程变为md²x/dt² kx=0。特征方程是mr² k=0。根为r=±i√(k/m)。令ω=√(k/m)。得到通解x=Acosωt Bsinωt。A和B是常数。初始位置和速度确定常数。这个解表示简谐振动。振幅不变。频率固定。实际总有阻力。振动逐渐停止。

有阻力情况更复杂。特征方程mr² cr k=0。判别式Δ=c²-4mk。Δ小于0系统振动。根为复数实部为负。振幅指数衰减。Δ等于0临界阻尼。最快回到平衡位置。Δ大于0过阻尼。慢慢回到平衡位置。这些情况都有应用。汽车减震器需要临界阻尼。尽快停止振动。钟摆需要小阻力。振动时间长。

电路系统也类似微分方程。RLC串联电路。电源电压U。电阻R电感L电容C。电流I。根据基尔霍夫电压定律。电阻电压IR。电感电压LdI/dt。电容电压Q/C。Q是电荷。电流是电荷变化率I=dQ/dt。方程是Ld²Q/dt² RdQ/dt Q/C=U。这个方程与弹簧振动方程形式相同。数学本质一样。解法相同。不同领域相同数学模型。

数值方法求解微分方程。很多方程没有解析解。计算机求数值解。欧拉方法最简单。考虑方程dy/dx=f(x,y)。初始条件y(x0)=y0。取步长h。计算下一个点y1=y0 hf(x0,y0)。x1=x0 h。继续这个过程。得到一系列点。近似解曲线。欧拉方法简单。精度不高。

改进欧拉方法。预估校正方法。先用欧拉方法预估。再用梯形公式校正。提高精度。龙格库塔方法更精确。常用四阶龙格库塔。计算四个斜率。加权平均。精度高。应用广泛。实际计算常用这些方法。

微分方程稳定性分析很重要。生态系统平衡状态。微小扰动后系统行为。平衡点满足方程右边为零。线性化方程。特征根实部符号决定稳定性。实部负稳定。实部正不稳定。稳定平衡能保持。不稳定平衡会破坏。这些分析很重要。预测系统长期行为。

传染病模型用微分方程。人群分三类。易感者S。感染者I。康复者R。接触传染。建立方程。dS/dt=-βSI。dI/dt=βSI-γI。dR/dt=γI。β是传染率。γ是恢复率。这个模型预测疫情发展。初始感染者少。指数增长。达到峰值。然后下降。最终停止。防控措施改变参数。减少接触降低β。医疗提高恢复率γ。模型指导疫情防控。

微分方程应用广泛物理学生物学经济学工程学。描述变化规律。通过建立方程。求解方程。分析解的性质。理解现象预测未来。学习微分方程很重要。掌握基本方法。解决实际问题。