微分方程详解.docVIP

第二章 微分方程 本章学习目的： 本章的主要目的在于：学习微分方程模型的建立、求解方法、分析结果及解决实际问题的全过程。 1．知道求解微分方程的解析法、数值解法以及图形表示解的方法； 2．理解求微分方程数值解的欧拉方法，了解龙格——库塔方法的思想； 3．熟练掌握使用MATLAB软件的函数求微分方程的解析解、数值解和图形解； 4．通过范例学习怎样建立微分方程模型和分析问题的思想。 §2.1 引例 在《大学物理》中，我们曾学习过简谐振动（如：弹簧振子、单摆），那就是一个典型的二阶常微分方程的模型。 这里我们讨论“倒葫芦形状容 器壁上的刻度问题”。 对于圆柱形状容器壁上的容积刻度，可以利用圆柱体体积公式：，其中容器的直径D为常数，体积V与相对于容器底部的任意高度H成正比，因此在容器壁上可以方便地标出容积刻度。 而对于几何形状不规则的容器，比如“倒葫芦形状”的容器壁上如何标出容积刻度呢？ 如图所示，建立坐标系，由微元法分析可知： ，其中x表示高度，直径是高度的函数，记为D（x） 如果该方程中的函数D(x)无解析表达式，只给出D(x)的部分测试数据，如何求解此微分方程呢？ h=0.2; d=[0.04,0.11,0.26,0.56,1.04,1.17]; x(1)=0;v(1)=0; for k=1:5 x(k+1)=x(k)+h; v(k+1)=v(k)+(h/2)*(pi/4)*(d