海上缉私模型论文数学建模MATLAB..docVIP

PAGE PAGE 21 实用标准 精彩文档 海上缉私问题 建模 题目二 组别：第五组 组长：练佳翔 组员：邵 力 组员：庞雪梅 海上缉私问题 摘要 针对海上缉私问题，要求出缉私船是否能追上走私船，或着是求缉私艇追上走私船的位置和时间，就需要知道走私船和缉私艇的位置坐标、大概的行驶路线、及二者的速度。对于走私船和缉私艇的位置坐标，可以由二者的行驶路线 、速度、行驶时间之间的关系得到。而走私船和缉私艇的位置坐标，可用三角函数、坐标关系、圆的位置关系求解。当缉私船追上走私船时，走私船和缉私艇的位置坐标相同，即二者的横坐标相等，纵坐标相等。在此期间，再加以MATLAB软件进行求解。 关键字： 海上缉私 位置坐标 速度 MATLAB软件 问题重述 分别对以下情况建立缉私船的位置和航线的数学模型，自己设定速度等参数，求数值解： (1) 走私船向正东方向非匀速直线行驶，其速度按正弦规律变化，如图1．已知缉私船以速度匀速追击，（为常数），两船初始距离．  图1 (2) 两船速度大小都不变，走私船以速度沿着与正东方向成角的直线行驶，如图2．已知缉私船的速度，两船初始距离．取与，求数值解，并说明走私船按哪个角度逃跑较快？ 图2 (3) 两船速度大小都不变，走私船以速度沿半径为的圆弧向点逃跑，现有两种方案，如图3．问两种方案是否都能到达点？已知圆弧半径，缉私船的速度，两船初始距离． 方案1　　　　　　　　　　　　　　　　方案2 图3 (4)两船速度都大小不变，走私船以速度先向正东方向直线行驶，一段时间（设尚未被缉私船追上）后改变方向，沿着与正东方向成角的直线行驶，如图4．已知缉私船的速度，两船初始距离．取，求数值解． 图4 (4) (5) 开始两船速度大小都不变，走私船以速度向正东方向沿直线行驶，但当两船距离小于时，缉私船会发现被人追击，将沿正北方向以速度加速逃跑，如图5．已知，，缉私船的速度，两船初始距离，求数值解． 图5 (6) 实际在追击时，缉私船速度方向的改变并不连续，每隔时间变换一次角度，在两次变换之间，缉私船按直线运动．若两船速度大小都不变，走私船以速度向正东方向沿直线行驶，（海里/小时），缉私船的速度（海里/小时），两船初始距离（海里），（秒）．试画出缉私船的航线图，建立此时的追击模型，比较与之前模型有何不同，并求数值解． 问题分析 问题一：要确定缉私船追上走私船的位置及时间，就必须确定缉私船、走私船的坐标。走私船的速度按正弦规律变化并向正东行驶，因此走私船的位移即向东行驶的距离可以用与坐标轴围成的面积表示。而缉私艇的坐标设为，再用缉私艇与走私船之间的位置关系，得出，的表达式。再根据两者坐标即可算出结果。 问题二：走私船以速度沿着与正东方向成角的直线行驶。根据走私船的起始点与角的位置关系，及走私船的速度，可以算出走私船坐标。而缉私艇的坐标设为，再用缉私艇与走私船之间的位置关系，得出，的表达式。再根据两者坐标即可算出结果。 问题三：由于走私船的行驶轨迹是圆弧，所以可以利用圆弧与所对圆心角的关系得出走私船的坐标的值。而缉私艇的坐标设为，再用缉私艇与走私船之间的位置关系，得出，的表达式。再根据两者坐标即可算出结果。 模型假设 在整个追赶拦截的过程中，缉私船和走私船都不会有故障发生导致船不能正常行驶，在海上没有发生因为天气突然发生变化阻碍船前行的情况。 走私船的行驶速度一直是呈正弦规律变化，缉私船一直是匀速前行。 建立直角坐标系，在缉私船发现走私船时计时（t=0），设此时走私船的位置在（0，c），缉私船位置在（0，0）。 设在任意时刻缉私船的坐标（x，y），走私船到达Q（at，c）点，直线PQ与缉私船航线相切，切线PQ与y轴正方向夹角为 设常数d=20，缉私船b=1.5d=30海里/小时，初始距离c=2d=40海里 符号定义 1、任意时刻缉私船的坐标（x，y） 2、走私船到达Q（at，c）点 3、切线PQ与y轴正方向夹角为 4、初始位置c 5、缉私船速度b 6、走私船速度a 7、走私船行驶t时间的路程k 模型建立与求解 模型一： 模型建立： （1）根据题目所给走私船的速度的图，可以求解出速度的表达式： 走私船的位移即向东行驶的距离可以用与坐标轴围成的面积表示如下： 又因为走私船从处向正东行驶，因此走私船在时刻的位置在： （2）缉私船的速度，两船的初始位置c=2d。缉私船在x，y方向是速度分别为 得到微分方程： 初始条件为： ， 模型求解： 用MATLAB求数值解，可得结果如图所示。 图1.1 x（t），y（t）曲线 图1.2 y（x）曲线 模型二 模型建立： ( 1 )走私船以速度沿着与正东方向成角的直线行驶，时间t时走私船的行驶距离为 根据角的三角函数关系可以得出，时间t时走私船的位置