常微分方程解法全解析：从理论基础到实践应用
在自然科学、工程技术及社会科学的诸多领域，常微分方程是描述动态过程演化规律的核心数学工具。无论是物理中的运动规律、化学中的反应动力学，还是工程中的控制系统、生态中的种群演化，都可以通过建立常微分方程模型来刻画。然而，建立微分方程仅仅是解决问题的第一步，要深入理解实际现象的演化趋势、验证模型的合理性，就必须求解方程并对解进行分析。
从解的形式来看，常微分方程的解可分为解析解和数值解。解析解是指能够用初等函数（如多项式、指数函数、三角函数等）及其有限次组合表示的解，这类解具有形式简洁、便于理论分析和应用的优点。但遗憾的是，解析解的存在具有极强的局限性——只有线性常系数微分方程，且自由项为指数函数、多项式、三角函数等特殊类型时，才能确保得到解析解。对于绝大多数变系数方程、非线性方程，即使形式看似简单（如 \( y' = y^2 + x^2 \)），也无法用初等函数表示其解，即所谓的“解不出来”。因此，对于大量实际问题中的常微分方程，数值解法成为不可或缺的核心手段。数值解法通过离散化处理，求解方程在若干离散点上的近似值，能够有效突破解析解法的限制，为实际问题提供可行的解决方案。本章将系统梳理常微分方程的数值解法体系，从离散化原理出发，依次介绍欧拉方法、改进欧拉方法、龙格-库塔方法、线性多步法等经典数值方法，同时涵盖一阶微分方程组与高阶微分方程的数值求解转化思路，并结合Matlab工具实现解法的工程化应用。
§ 1 常微分方程的离散化：数值解法的核心基础
数值解法的核心思想是“离散化”——将连续的微分方程转化为离散的代数方程（差分方程），通过求解差分方程得到原微分方程解在离散点上的近似值。本节以一阶常微分方程的初值问题为研究对象，奠定数值解法的理论基础。
