《交通管理问题》PPT幻灯片

1、1,数学建模实 验,王汝军 河西学院数学与统计学院,2,实验十 交通管理问题,王汝军 河西学院数学与统计学院,3,实验目的,1了解微分方程的一些基本概念。 2初步掌握微分方程模型建立、求解的基本方法和步骤。 3学习掌握用MATLAB软件中相关命令求解常微分方程的解析解。,4,实验内容,在城市道路的十字路口，都会设置红绿交通灯。为了让那些正行驶在交叉路口或离交叉路口太近而又无法停下的车辆通过路口，红绿灯转换中间还要亮起一段时间的黄灯。对于一名驶近交叉路口的驾驶员来说，万万不可处于这样进退两难的境地：要安全停车但又离路口太近；要想在红灯亮之前通过路口又觉得距离太远。那么，黄灯应亮多长时间才最为合理

2、呢？ 已知城市道路法定速度为v0 ，交叉路口的宽度为 I，典型的车身长度统一定为 L，一般情况下驾驶员的反应时间为T ，地面的磨擦系数为。(假设 I9 ， L4.5 ， 0.2， T1 s),5,实验准备,微分方程是研究函数变化过程中规律的有力工具，在科技、工程、经济管理、人口、交通、生态、环境等各个领域有着广泛的应用。如在研究牛顿力学、热量在介质中的传播、抛体运动、化学中液体浓度变化、人口增长预测、种群变化、交通流量控制等等过程中，作为研究对象的函数，常常要和函数自身的导数一起，用一个符合其内在规律的方程，即微分方程来加以描述。,6,1微分方程的基本概念 未知的函数以及它的某些阶的导数连同自

3、变量都由一已知方程联系在一起的方程称为微分方程。如果未知函数是一元函数，称为常微分方程。如果未知函数是多个变量的函数，称为偏微分方程。联系一些未知函数的多个微分方程称为微分方程组。微分方程中出现的未知函数的导数的最高阶数称为微分方程的阶。若方程中未知函数及其各阶导数都是一次的，称为线性常微分方程，一般表示为,7,(1),若（1）式中系数,（i=1,2,n）均与t 无关，称之为常系数,（或定常、自治、时不变）的。,8,建立微分方程模型要根据研究的问题作具体的分析。一般有以下三种方法： 根据规律建模：在数学、力学、物理、化学等学科中已有许多经过实践检验的规律和定律，如牛顿运动定律、基尔霍夫电流及电

4、压定律、物质的放射性规律、曲线的切线的性质等，这些都涉及某些函数的变化率。我们可以根据相应的规律，列出常微分方程。,9,微元法建模：利用微积分的分析法建立常微分方程模型，实际上是寻求一些微元之间的关系式，在建立这些关系式时也要用到已知的规律或定理。与第一种方法不同之处在于这里不是直接对未知函数及其导数应用规律和定理来求关系式，而是对某些微元来应用规律。,10,模拟近似法建模：在社会科学、生物学、医学、经济学等学科的实践中，常常要用模拟近似法来建立微分方程模型。这是因为，上述学科中的一些现象的规律性我们还不是很清楚，即使有所了解也并不全面，因此，要用数学模型进行研究只能在不同的假设下去模拟实际的

5、现象。如此模拟近似所建立的微分方程从数学上求解或分析解的性质，再去同实际情况作对比，观察这个模型能否模拟、近似某些实际的现象。,11,建立微分方程模型只是解决问题的第一步，通常需要求出方程的解来说明实际现象，并加以检验。,12,2微分方程通解的求解方法 （1）初等积分法 有些微分方程可直接通过积分来进行求解。例如，一阶常系数线性常微分方程 可化为,13,两边通过积分可得到通解y(t)为 其中 为任意的常数。有些常微分方程可用一些技巧（如分离变量法、积分因子法、常数变易法、降阶法等）化为可积分的方程而求得解析解。,14,（2）常系数线性微分方程求解 线性常微分方程的解满足叠加性原理，从而它的求解

6、可归结为求一个特解和相应齐次微分方程的解。一阶变系数线性常微分方程总可用这一思路来求得通解。高阶线性常系数微分方程可用特征根法求得相应齐次微分方程的基本解，再用常数变易法求特解。,15,例如，求 的通解。 解：特征方程为 在MATLAB命令框中输入命令 x=roots(1 0.2 3.92)%roots命令用来求多项式的根 求解得到一对共轭复根 x = -0.1000 + 1.9774i -0.1000 - 1.9774i,16,从而该微分方程的通解x(t)为 其中 A、B 为任意的常数。,17,一阶常微分方程组与高阶常微分方程可以互化，已给一个 n阶方程 (2) 设 (2)可化为一阶方程组,

7、18,(3),19,反过来，在许多情况下，一阶微分方程组也可以化为高阶方程。所以一阶常微分方程组与高阶常微分方程的理论与方法在很多方面是相通的。一阶常系数线性微分方程组也可用特征根法进行求解。,20,3求微分方程（组）通解的MATLAB命令 求解微分方程（组）的解析解用函数dsolve。,21,实验方法与步骤,1dsolve命令的基本用法 下面以例题来予以说明： 例1求高阶方程 的通解. 输入命令： r=dsolve(D2y=cos(2*x)-y,y(0)=1,Dy(0)=0,x) 可得：,22,r = (1/2*sin(x)+1/6*sin(3*x)*sin(x)+(1/6*cos(3*x)

8、-1/2*cos(x)*cos(x)+4/3*cos(x) r=simple(r)%对r进行合并、分解化简 r = -1/3*cos(2*x)+4/3*cos(x),23,例2求下面微分方程组的通解,24,输入命令 x,y,z=dsolve(Dx=2*x-3*y+3*z,Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z); x=simple(x) x = -(-C1-C2*exp(-3*t)+C2-C3+C3*exp(-3*t)*exp(2*t),25, y=simple(y) y = -(C1*exp(-4*t)-C1-C2*exp(-4*t)-C2*exp(-3*t)+C2-C3

9、+C3*exp(-3*t)*exp(2*t) z=simple(z) z = (-C1+exp(4*t)*C1-C2*exp(4*t)+C2+exp(4*t)*C3)*exp(-2*t),26,2引例问题的分析与求解,首先，我们用模拟近似法对引例问题进行分析建模。 对于驶近交叉路口的驾驶员，在他看到黄色信号后要做出决定：是停车还是通过路口。如果他以法定速度（或低于法定速度）行驶，当决定停车时，他必须有足够的停车距离。当驾驶员决定通过路口时，必须有足够的时间让他能完全通过路口。这包括做出停车决定的反应时间以及通过停车所需的最短距离的驾驶时间，能够很快看到黄灯的驾驶员可以利用刹车距离将车停下来。,

10、27,于是，黄灯状态所应持续的时间包括驾驶员的反应时间，他通过交叉路口的时间以及通过刹车距离所需要的时间。 由题设可知城市道路法定速度为v0 ，交叉路口的宽度为I ，典型的车身长度统一定为L 。考虑到车通过路口实际上指的是车的尾部必须通过路口，因此，通过路口的时间为 .,28,现在我们来计算刹车距离：设w为汽车的重量，为磨擦系数，由牛顿力学知，地面对汽车的磨擦力为w ，其方向与汽车运动的方向相反。汽车在停车过程中，由牛顿第一动力定理有 fma 其中m为汽车质量（即 w/g，g 为重力加速度），a 为汽车的加速度， f是汽车所受的摩擦力。这里加速度a 是停车距离 关于时间t的二阶导数，所以行驶距

11、离 x与时间t 的关系可由下面的微分方程确定：,29,约去w，化简（4）式得,同时，我们知道，当 t0时，距离 x0，初速度是距离x在0时刻的一阶导数，于是可以给出方程（5）的初始条件,30,在MALAB命令框中输入命令 x=dsolve(D2x=-ug,x(0)=0,Dx(0)=v0,t) x = -1/2*ug*t2+v0*t 即得到停车距离x关于时间t 的解析式。,31,32,设黄灯闪烁的时间为A ，则A 的表达式为,33,结果分析,由假设知， I9 ， L4.5 ， T1s ，磨擦系数选取有代表性的 0.2，我们考虑当法定速度 v040、60、80km/h时，黄灯时间如表1所示，表1也

12、给出了与经验法黄灯时间的对比。,34,表1黄灯预测时间与经验法时间的对比,35,我们注意到，经验法的结果一律比我们预测的黄灯状态时间要短些，这使得我们联想起，许多城市交叉路口红、黄、绿灯的设计可能使得司机驾驶着的汽车在绿灯转变为红灯的时刻正处于交叉路口的位置。,36,练习与思考,1设一容积为 （单位：m3 ）的大湖受到某种化学废料的污染，污染物均匀地分布在湖中。若某时刻起污染源被切断，设湖水更新的速率是 （单位是： m3 天）。试建立求污染物的浓度下降至原来的5%所需时间的数学模型。美国密西根湖的容积为4871 103 （ m3 ），湖水的流量为3.663959132 103 （ m3 ），求污染中止后，污染物浓度下降到原来湖水污染浓度的3%所需要的时间。,37,练