《渡轮航行路线的设计by王自伟_何庆明_张兴强》由会员分享，可在线阅读，更多相关《渡轮航行路线的设计by王自伟_何庆明_张兴强（17页珍藏版）》请在金锄头文库上搜索。
1、1 渡轮航行路线的设计摘 要我国是世界上的一个多河流国家，为了方便河流两岸人民的出行，渡轮成了交通的一个重要组成部分。本文以我国东部地区某市为例，采用数学分析方法建立了一系列数学模型来规划和设计渡轮的航行路线。对于问题一，考虑河水流速以及渡轮在静水中的速度在整个航行过程中保持不变，本文以北岸码头为原点、以南北岸码头连线为X轴、以河流的水流方向为Y轴建立恰当的直角坐标系，对渡轮航行的路线关于时间建立多元微分方程组，再运用三角函数的几何意义对方程组进行化简，得出了初步的微分方程组模型，并由此建立了模型一。使用 MATLAB 编程求解，在得不出解析解的情况下采用4 阶 Runge-Kutta 算法求
2、解微分方程组的数值解。 得出了渡轮的航行路线图 （下见图 4），以及在X轴方向、Y轴方向上的分位移关于时间的响应曲线（下见图5）。并分析数据后求得在该航行路线下所花费的时间666.7 ts，合 11分钟 6.7 秒。对于问题二，本文在问题一的基础上适当地改变坐标轴，以及各个码头的坐标；以北岸码头 2 为原点，以北岸码头 2 与南岸码头 2 的连线为X轴，以沿北岸码头 1 向北岸码头 2 方向的射线为Y轴，建立平面直角坐标系，实则是在模型一的基础上修改微分方程组的初值， 化简后得到新的微分方程组模型，并由此建立了模型二。 相同地，使用 MATLAB 编程，在得不出解析解的情况下采用4 阶 Run
3、ge-Kutta 算法求其的数值解。得出了渡轮的航行路线图（下见图6），以及在X轴方向、Y轴方向上的分位移关于时间的响应曲线（下见图7）。分析数据后可知在该路线下所花时间为578.6s，合 10 分钟 38.6 秒。对于问题三，本文在问题一和问题二的分析基础上，根据实际生活中河水流速不均匀分布的条件限制， 先使用 MATLAB 拟合出河水流速关于距北岸距离的曲线方程：2 10.0000030.0030210.733095vxxx。并由此修正模型一和模型二的微分方程组，得出模型三。同时修改微分方程组，并用MATLAB求数值解。根据求解结果，本文设计渡轮在该情形下的航行路线图（下见图10，图 12
4、）。并得出时间分别为735.5s和 608s。关键词：MATLAB 微分方程组Runge-Kutta 算法数据拟合小船过河2 一、问题重述某市有一条东西向的河流穿越市区，为方便两岸市民出行，设置了渡轮往返于南北岸间。如图 1 所示，南岸码头正对着北岸码头，河床宽度大约为1000d米。请你们通过建立数学模型，解决以下问题：1、假定河水流速1v ，渡轮在静水中的速度2v ，在整个航行过程中它们保持不变，试设计一条渡轮航行路线， 使得按照这条路线行驶从南岸码头出发正好能够到达北岸码头（或者从北岸码头出发正好能够到达南岸码头）。当11/vm s，22/vm s，时，画出航行路线图，并计算出按照这条路线
5、行驶所需的渡河时间。2、为缩短乘客航行时间，决定增设南岸码头和北岸码头各一个，如图2 所示。从北岸到南岸的乘客在北岸码头1登船，到南岸码头2 下船；而从南岸到北岸的乘客在南岸码头 1 登船，到北岸码头2下船。如果同一岸的两个码头间距设计为500 米，当11/vm s，22/vm s，时，画出航行路线图，并计算出按照这条路线行驶所需的渡河时间，与问题1 进行比较。3、下表 1 中是测出的水流速沿离北岸边距离的分布，在22/vm s的情况下，重新解决问题 1 和问题 2。图 1图 2表 1 水流速沿离北岸边距离的分布编号1 2 3 4 5 6 7 8 9 离北岸边距离（米）100 200 300
6、400 500 600 700 800 900 河水流速（米/秒）1.00 1.22 1.38 1.45 1.50 1.46 1.37 1.23 1.02 1000m水流方向北岸码头南岸码头北岸码头 2 南岸码头 2 1000m水流方向北岸码头1 南岸码头1 3 二、问题分析根据问题重述，可知这实质上是求曲线轨迹的微分方程问题。问题一对于问题一，为了方便我们更加清晰的剖析问题以及解决问题，我们对问题重述的图一做了适当旋转，以北岸码头为原点、以南北岸码头连线为X轴、以河流的水流方向为Y轴建立恰当的直角坐标系。由于船速与水速都是固定值，决定渡轮运动轨迹的只有船头指向。根据生活实际，我们做出合理假设
7、：渡轮在任意时刻的船头始终指向北岸码头。 由此得出渡轮运动时任一点的坐标，与其行驶方向同X轴的夹角的三角函数关系。然后通过对渡轮行驶途中的任意一时刻做运动分析，我们很容易就能得到一组任一点的分速度对时间的微分方程组。运用MATLAB 编程，便可得到渡轮运行的轨迹方程，并求出行驶时间。问题二对于问题二，我们在问题一的基础上做适当的修改，以北岸码头2 为原点，以北岸码头 2 与南岸码头 2 的连线为X轴，以沿北岸码头 1 向北岸码头 2 方向的射线为Y轴，重新建立平面直角坐标系。然后做和问题一相同的处理：即根据渡轮运动时任一点的坐标，与其行驶方向同X轴的夹角的三角函数关系， 以及任意一点分速度关于
8、时间的微分，得出新的微分方程组。相同地，使用MATLAB编程，便可得出渡轮的部分航行路线图，并得出航行时间。问题三对于问题三，考虑到实际生活中河水的流速并不是均匀分布的，我们可以先通过问题重述中给出的水流速度的测量数据做二次拟合，得到水速分布的曲线方程。再将得到的曲线方程看做一个整体分别替代问题一，问题二中的水速。然后用前面所述方法列出微分方程。并运用MATLAB 编程，并求解分析。便可得出渡轮行驶的路线以及所花费的时间。4 三、模型假设（1） 假设不计渡轮本身长度；（2） 假设渡轮速度恒为2v ；（3） 假设不计算渡轮启动终止时间；（4） 假设渡轮航行过程中畅通无阻；（5） 假设河水流速数据
9、符合某拟合函数；（6） 假设渡轮、河水流速不受风等天气因素影响；四、符号说明符号含义A南岸码头（ 1）B北岸码头（ 1）C南岸码头 2 D北岸码头 2 d南北岸的河床宽度S同一侧的两个码头的距离xX轴方向x t渡轮在x方向上的位移x t渡轮在x方向上的速度yY轴方向y t渡轮在y方向上的位移y t渡轮在y方向上的速度5 1v河水流速1vx河水流速关于x方向上位移函数2v渡轮在静水中的速度v合渡轮行驶时的合速度t渡轮行驶的某个时刻t渡轮在同一侧的两个码头间逆行所用时间t总渡轮航行所用的总时间渡轮船头朝向与 x轴夹角五、模型的建立和求解5.1 问题一模型的建立和求解5.1.1 模型一的建立以北岸码
10、头为原点，以沿北岸码头指向南岸码头的射线为x轴，沿北岸码头向东方向的射线为y轴，建立平面直角坐标系，如图3 所示：图 3设渡轮由南岸码头（ A 点）出发，驶向北岸码头（ B 点）。在 t 时刻，渡轮在x方6 向上的位移为 x t ，在y方向上的位移为 y t ；故渡轮在x方向上的速度为 x t ，在y方向上的速度为y t 。再设渡轮的船头朝向与水平x轴的夹角为，故将水速1v 与船速2v 分解后有：212cossinxtvytvv又因为渡轮运行时船头始终朝向北岸码头（B 点），所以有：2222cossinxxyyxy当0t时，渡轮还在南岸码头，故01000x，00y。综合上述条件，可得微分方程组
11、：2222 1220100000v xdxdtxyv ydyvdtxyxy将11 /vm s，22 /vm s带入微分方程求解即可。5.1.2 模型一的求解利用 MATLAB 编写程序，发现无法得出解析解， 故采取 4 阶 Runge-Kutta 算法求其数值解。在建立的boat1.m 文件，加入（ norm(x)1e-5）的限制条件，以保证渡轮在离北岸码头（ B 点）足够近的时候可以终止运算。求解后得出在水速11 /vm s、船速22 /vm s下，渡轮从南岸码头驶向北岸码头的运行路线（如图4），以及x，y关于时间 t 的响应曲线图（如图5）：7 图 4图 5由数值解法可知按照这条路线行驶所
12、需要的渡河时间666.7 ts，合 11 分钟 6.7秒。5.2 问题二模型的建立和求解5.2.1 模型二的建立与模型一不同，在模型二中我们以北岸码头2（D 点，下同）为原点，以沿北岸码头 2 指向南岸码头 2（C 点）的射线为x轴，沿北岸码头 1（B 点）向北岸码头 2 方向的射线为y轴，建立平面直角坐标系，如图6 所示：图 6 设渡轮由南岸码头1（A 点）出发，驶向北岸码头2（D 点）。在t时刻，渡轮在x方向上的位移为 x t ，在y方向上的位移为y t ；故渡轮在x方向上的速度为 x t ，8 在y方向上的速度为y t 。与模型一相似，渡轮的船头朝向与水平x轴的夹角为，渡轮运行时船头始终
13、朝向北岸码头2（D 点），故有与模型一相同的方程组：212cossinxtvytvv，2222cossinxxyyxy与模型一相比， A 点坐标的变化，使得初值条件发生改变：当0t时，渡轮还在南岸码头 1（A 点），故01000x，0500y。综合上述条件，可得微分方程组：2222 122010000500v xdxdtxyv ydyvdtxyxy将11 /vm s，22 /vm s带入微分方程求解即可。但是，值得考虑的是：当小船到达北岸码头2 后，需要逆水行驶到北岸码头1，才能实现二次利用，该过程中小船速度和水流速度在同一条水平线上，故它们的合速度 v合保持恒定，且21vvv合，时间Stv合
14、，即可求出。5.2.2 模型二的求解利用 MATLAB编程，与模型一相同，无法得出解析解，采取4 阶 Runge-Kutta算法求其数值解。在建立的boat2.m 文件中做与模型一一致的要求。求解后得出在水速11 /vm s、船速22 /vm s下，渡轮从南岸码头 1 出发，驶向北岸码头 2 的运行路线（如图 7），以及x，y关于时间的响应曲线图（如图8）：9 图 7 图 8 由数值解法可知按照这条路线行驶，乘客从南岸到北岸所需要的渡河时间578.7 ts。渡轮到达后， 从北岸码头 2 逆水行驶到北岸码头1，该过程中合速度 v合保持恒定，且212 1 /1 /vvvm sm s合，时间500=
15、500 1 /Smtsvm s合。综上可知，利用模型二方案，乘客从南岸码头1 到北岸码头2 需要花费时间578.6 ts，合 9 分钟 38.6秒。渡轮从南岸码头1 到北岸码头 1 总共需要花费的时间为= 1078.6 ttts总。5.2.3 模型一与模型二的比较模型一经过适当的修改， 例如 A点的坐标及微分方程的初值， 便可以得到模型二，两个模型的本质及算法是相同的。就两岸人民出行目的而言，模型二下渡轮的航行路线虽然长，但其所花费时间更少，更符合人们的实际生活需求；但渡轮在到达北岸码头 2 后，需逆水行驶到北岸码头1 才能执行下一趟任务，这样花费的时间和资源将会更多，而模型一就没有这个缺陷。
16、10 5.3 问题三模型的建立和求解5.3.1 河水流速的拟合在实际生活中，河水中水流的速度并不会均匀分布。根据表1 中测出的水流速度沿离北岸边距离的分布数据， 利用 MATLAB 拟合出水流速度1v 关于离北岸的距离x的关系曲线，如下图9：图 9 因为使用 MATLAB 拟合出的水速分布的曲线方程的二次项系数很小，故将数据表示改为长整型，并保留六位小数，使得二次项系数数据得以保留，曲线方程如下：2 10.0000030.0030210.733095vxxx5.3.2 模型三的建立针对问题一：在模型一的条件下，用1vx 替换1v ，其他条件保持不变，于是有：11 212cossinxtvytvxv,2222cossinxxyyxy当0t时，渡轮还在南岸码头，故01000x，00y。改变后的微分方程组：22222220.0000030.0030210.7330950100000v xdxdtxyv ydyxxdtxyxy将22 /vm s带入微分方程求解即可。针对问题二 : 在模型二的条件下，用1vx 替换1v ，其他条件保持不变，于是有：212cossinxtvytvxv,2222c
