微分方程求解.doc

微分方程在科技、工程、经济管理以及生态、环境、人口、交通等各个领域中常用于建立数学模型，它是研究函数变化规律的有力工具.如在研究弹性物体的振动，电阻、电容、电感电路的瞬变，热量在介质中的传播，抛射体的轨迹，以及污染物浓度的变化，人口增长的预测，种群数量的演变，交通流量的控制等等过程中，作为研究对象的函数，要和函数的导数一起，用一个符合其内在规律的方程，即微分方程来描述. 建立微分方程只是解决问题的第一步，通常需要求出方程的解来说明实际现象，并加以检验.如果能得到解析解固然是便于分析和应用的，但是我们知道，只有线性常系数微分方程，并且自由项是某些特殊类型的函数时，才可以肯定得到这样的解，而绝大多数变系数方程、非线性方程都很难得到解析解.另外，有时即使能求出解析解，也会由于很难从解析解中计算函数的值而不实用.例如，容易求出初值问题 的解为 eed. 但是，对于给定的，要计算函数值还需要用数值积分的方法.于是对于用微分方程解决实际问题来说，数值解法就是一个十分重要的手段.为了将微分方程的数值解与解析解比较，本章首先简单介绍用MATLAB解微分方程的解析解（符号解）的方法，然后重点介绍用MATLAB程序计算微分方程的数值解法. 10.1 常微分方程（组）的MATLAB符号求解 在MATLAB系统中提供了常微分方程（组）符号解（Symbolic solution of ordinary di