精讲多练Matlab数学建模实例教程分析.ppt

综合实例 * * 综合实验案例 用 MATLAB 求解问题时，一般要经历建模和编程两个过程，只有在建模正确的前提下，方能得出正确的结果。 实例1：狼追击兔子的问题 狼追击兔子问题是欧洲文艺复兴时代的著名人物达.芬奇提出的一个数学问题。 当一个兔子正在它的洞穴南面60码处觅食时，一只恶狼出现在兔子正东的100码处。 当两只动物同时发现对方以后，兔子奔向自己的洞穴，狼以快于兔子一倍的速度紧追兔子不放。 狼是否会在兔子跑回洞穴之前追赶上兔子？ 问题分析 为了研究狼是否能够追上兔子，可以先考虑求出狼追兔子形成的追击曲线，然后根据曲线来确定狼是否能够追上兔子。 v1：兔子的速度（单位：码/秒） r：狼与兔子速度的倍数； v2：狼的速度（单位：码/秒），显然有v2=rv1 t：狼追击兔子的时刻（t=0时，表示狼开始追兔子的时刻） s2：在时刻t，狼跑过的路程（单位：码），s2=s2(t) s1：在时刻t，兔子跑过的路程（单位：码），s1=s1(t) Q(x1, y1) ：表示在时刻t时，兔子的坐标 P(x, y) ：表示在时刻t时，狼的坐标 变量说明 狼在追击过程中始终朝向兔子； 狼追击兔子的轨迹看作是一条光滑的曲线， 即将动点P(x,y)的轨迹看作一条曲线， 曲线方程表示为y=y(x)。 模型假设 以t＝0时，兔子的位置作为直角坐标原点， 兔子朝向狼的方向为x轴正向； 则显然有兔子位置的横坐标x1=0, y1=0。 对狼来说，当x＝100，y＝0， 即对于y=y(x)函数来说，有y|x=100=0 在t＝0刚开始追击时，狼的奔跑方向朝向兔子，此时即x轴负方向， 即有y’|x=100=0 （一）建模准备 1、追击方向的讨论 由于狼始终朝向兔子，则在狼在时刻t所在位置P(x,y)点在狼的轨迹处的切线方向在y轴上的截距为y1。 设切线上的动点坐标为（X，Y），则切线方程为 Y – y = y’(X-x) （1） 对于该方程是否有可利用的其它已知条件？ 它与y轴的交点即为(0,y1), 带入方程可得 y1=y-y’x （2） 且y1=v1t，即时刻t时兔子所跑过的路程。 （二）建立模型 在t时刻，兔子跑过的路程为 s1 = y1 = v1t （3） 由于狼的速度是兔子的r倍，则狼跑的路程为 s2 = rs1 = ry1 （4） 狼跑过的路程可以用对弧长的曲线积分知识得到 （5） 联立(2), (4), (5)得 （6） （6）两边对x求导得 （7） 微分方程（7）式的初始条件有： y|x=100=0, y’|x=100=0 2、狼与兔子速度关系的建模 要判定狼最终是否追上兔子， 可以通过（7）式 判定。 对（7）式， 当x＝0，如果计算求解得到y≥60 ，则视为没有追上； 当x＝0，如果计算求解得到y60 ，则视为兔子被追上； 3、是否追上的判断 由微分方程 得到其Matlab函数 function yy=odefunlt(x,y) %以狼在追击过程中的横坐标为自变量 yy(1,1)=y(2); yy(2,1)=sqrt(1+y(2).^2)./(2.*x); 模型求解 tspan=100:-0.1:0.1;%以狼的x坐标为自变量 y0=[0 0]; %下面只知道狼是否追上兔子，但是不易推得兔子刚刚到达窝边时，狼与兔之间的距离 [T,Y] = ode45(odefunlt,tspan,y0); n=size(Y,1); disp(狼的坐标(x=0.1)) disp(Y(n,1)) %通过追击曲线计算当狼的横坐标为0.1(即tspan=0.1)时，狼的纵坐标 主程序 运行结果： 狼的坐标(x=0.1) 62.1932 由运行结果可知，狼并没有追上兔子。 模型结果与分析 通过上面的结果已经知道狼并没有追上兔子。那么兔子跑回窝边时，狼与兔子之间的距离是多少？ 思考题 解决思路 可以对狼与兔子的追击过程通过计算机进行模拟，然后从模拟结果获取。 function sim_