4 微分及方程建模.ppt

微分方程模 型 广州中医药大学数学模型课程 §3.1 微分方程的简单实例 §3.3 为什么要用三级火箭来发射人造卫星 §3.4 药物在体内的分布 §3.5 传染病模型 §3.6 糖尿病的诊断 §3.7 稳定性问题 §3.8 捕食系统的Volterra方程 §3.9 较一般的双种群生态系统 §3.10 分布参数法建模 例8 人口问题的偏微分方程模型 例9 交通流问题 作直线l1: x1+x2=K1及l2: x1+x2=K2， K1 K2，见图3-26。 dx1/dt0 dx2/dt0 图3-26 III II I k1 k2 dx1/dt0 dx2/dt0 dx1/dt0 dx2/dt0 有以下几个引理: 引理1 若初始点位于区域I中，则解 （x1(t)、x2(t)）从某一时刻起 必开此区域而进入区域II 引理2 若初始点（x1(0)、x2(0)）位于 区域II中，则（x1(t)，x2(t)）始 终位于II中，且： 引理3 若初始点位于区域III中,且对于 任意t ，（x1(t)，x2(t)）仍位于 III中，则当t→+∞时，（x1(t)， x2(t)）必以（K1,0）为极限点。 由引理1和引理2，初始点位于像限I和II的解必趋于平衡点（K1,0）。由引理3，初始点位于III且（x1(t)，x2(t)）始终位于III中的解最终必趋于平衡点（K1,0），而在某时刻进入区域II的解由引理最终也必趋于（K1,0）。易见只有上述三种可能，而在三种可能情况下（x1(t)，x2(t)）均以（K1,0）为极限，定理得证。 定理4的证明： 在研究实际课题时，数值解方法也许会用得更多。当解析解无法求得时，计算机作为强大的辅助工具发挥了它应起的作用。我校学生在研究1999年美国大学生数学建模竞赛题A（小行星撞击地球）时就遇到了一个棘手的问题：如何描述南极地区的生态系统，如何定量化地研究小行星撞击地球对南级生态环境的影响？在上网查阅了南极附近的海洋生态状况后，他们将南极附近的生物划分成三个部分：海藻、鳞虾和其他海洋生物。鳞虾吃海藻，其他海洋动物吃鳞虾，运用基本建模技巧建立了一个三房室系统模型。小行星的撞击会影响大气层的能见度，从而影响到海藻的生长（光合作用），进而影响到生物链中的其他生物。他们无法得到模型中的参数值（事实上，小行星撞击南极的事件并未发生过），就取了一系列不同的参数值，对不同参数值下模型的数值解进行了分析对比，研究了解对各参数变化的灵敏度，取得了十分有意义的结果并获得了当年国际竞赛的一等奖。 前面建立的模型都用了考察对象在系统中的均匀分布假设。这种方法建模被称为集中参数法。 考虑个体差异（或分布差异）的建模方法被称为分布参数法。分布参数法用于连续变量的问题时，得到的通常都是偏微分方程，无论建模还是求解都比较困难。仅举两个简单例子，来说明这种方法的应用。 即可令 t = 0时 G = G0, H = H0且 F1 ( G0,H0 ) = 0 F2 ( G0,H0 ) = 0 从而有 在测试过程中 G , H 均为变量，而我们关心的却只是它们的改变量，故令g = G – G0, h = H – H0 , 在（ 3.19 ）中将 展开，得到 其中 、 是g 和h 的高阶无穷小量。 很小时（即检测者至多为轻微病人时），为求解方 便，我们考察不包含它们的近似方程组 方程组（ 3.20 ）是一个非线性方程组，较难求解。当 、 首先，我们来确定右端各项的符号。从图 中可看出，当J(t)=0 时,若g ＞ 0 且 h = 0，则此人血糖浓度高于正常值，内分泌激素将促使组织吸收葡萄糖，并将其存储进肝 脏，此时有 ﹤ 0，从而应有： 0 其激素浓度将增加以抑制血糖浓度的增高，因而又有：: ＞ 0 反之，当J(t)=0而g=0且h＞0 时，此人激素浓度高于正常值，血糖浓度及激素浓度均将减少，从而必有 将方程组（ 3.20 ）改写成 其中 均为正常数。 （ 3.21 ）是关于 g、h的一阶常系数微分方程组，因激素浓度不易测得，对前式再次求导化为： 由于 故 或