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

《海上缉私模型论文数学建模MATLAB.doc》由会员分享，可在线阅读，更多相关《海上缉私模型论文数学建模MATLAB.doc（26页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、- .海上缉私问题建模 题目二组别：第五组组长：练佳翔组员：邵 力组员：庞雪梅海上缉私问题摘要针对海上缉私问题，要求出缉私船是否能追上走私船，或着是求缉私艇追上走私船的位置和时间，就需要知道走私船和缉私艇的位置坐标、大概的行驶路线、及二者的速度。对于走私船和缉私艇的位置坐标，可以由二者的行驶路线 、速度、行驶时间之间的关系得到。而走私船和缉私艇的位置坐标，可用三角函数、坐标关系、圆的位置关系求解。当缉私船追上走私船时，走私船和缉私艇的位置坐标一样，即二者的横坐标相等，纵坐标相等。在此期间，再加以MATLAB软件进展求解。关键字： 海上缉私 位置坐标 速度 MATLAB软件问题重述分别对以下情况

2、建立缉私船的位置和航线的数学模型，自己设定速度等参数，求数值解：(1) 走私船向正向非匀速直线行驶，其速度按正弦规律变化，如图1缉私船以速度匀速追击，为常数，两船初始距离图1(2) 两船速度大小都不变，走私船以速度沿着与正向成角的直线行驶，如图2缉私船的速度，两船初始距离取与，求数值解，并说明走私船按哪个角度逃跑较快？图2(3) 两船速度大小都不变，走私船以速度沿半径为的圆弧向点逃跑，现有两种方案，如图3问两种方案是否都能到达点？圆弧半径，缉私船的速度，两船初始距离方案1方案2图3(4)两船速度都大小不变，走私船以速度先向正向直线行驶，一段时间设尚未被缉私船追上后改变方向，沿着与正向成角的直线

3、行驶，如图4缉私船的速度，两船初始距离取，求数值解图4(4)(5) 开场两船速度大小都不变，走私船以速度向正向沿直线行驶，但当两船距离小于时，缉私船会发现被人追击，将沿正北方向以速度加速逃跑，如图5，缉私船的速度，两船初始距离，求数值解图5(6)实际在追击时，缉私船速度方向的改变并不连续，每隔时间变换一次角度，在两次变换之间，缉私船按直线运动假设两船速度大小都不变，走私船以速度向正向沿直线行驶，海里/小时，缉私船的速度海里/小时，两船初始距离海里，秒试画出缉私船的航线图，建立此时的追击模型，比拟与之前模型有何不同，并求数值解问题分析问题一：要确定缉私船追上走私船的位置及时间，就必须确定缉私船、

4、走私船的坐标。走私船的速度按正弦规律变化并向正东行驶，因此走私船的位移即向东行驶的距离可以用与坐标轴围成的面积表示。而缉私艇的坐标设为，再用缉私艇与走私船之间的位置关系，得出，的表达式。再根据两者坐标即可算出结果。问题二：走私船以速度沿着与正向成角的直线行驶。根据走私船的起始点与角的位置关系，及走私船的速度，可以算出走私船坐标。而缉私艇的坐标设为，再用缉私艇与走私船之间的位置关系，得出，的表达式。再根据两者坐标即可算出结果。问题三：由于走私船的行驶轨迹是圆弧，所以可以利用圆弧与所对圆心角的关系得出走私船的坐标的值。而缉私艇的坐标设为，再用缉私艇与走私船之间的位置关系，得出，的表达式。再根据两者

5、坐标即可算出结果。模型假设（1） 在整个追赶拦截的过程中，缉私船和走私船都不会有故障发生导致船不能正常行驶，在海上没有发生因为天气突然发生变化阻碍船前行的情况。（2） 走私船的行驶速度一直是呈正弦规律变化，缉私船一直是匀速前行。（3） 建立直角坐标系，在缉私船发现走私船时计时t=0，设此时走私船的位置在0，c，缉私船位置在0，0。（4） 设在任意时刻缉私船的坐标x，y，走私船到达Qat，c点，直线PQ与缉私船航线相切，切线PQ与y轴正方向夹角为（5） 设常数d=20，缉私船b=1.5d=30海里/小时，初始距离c=2d=40海里符号定义1、任意时刻缉私船的坐标x，y2、走私船到达Qat，c点3

6、、切线PQ与y轴正方向夹角为4、初始位置c5、缉私船速度b6、走私船速度a7、走私船行驶t时间的路程k模型建立与求解模型一：模型建立：1根据题目所给走私船的速度的图，可以求解出速度的表达式：走私船的位移即向东行驶的距离可以用与坐标轴围成的面积表示如下：又因为走私船从处向正东行驶，因此走私船在时刻的位置在：2缉私船的速度，两船的初始位置c=2d。缉私船在x，y方向是速度分别为得到微分方程：初始条件为：，模型求解：用MATLAB求数值解，可得结果如下图。图1.1 xt，yt曲线图1.2 yx曲线模型二模型建立：( 1 )走私船以速度沿着与正向成角的直线行驶，时间t时走私船的行驶距离为根据角的三角函

7、数关系可以得出，时间t时走私船的位置：2缉私船的速度，两船的初始位置c=2d。缉私船在x，y方向是速度分别为 得到微分方程：初始条件为：，3把与代入1中分别得出走私船的坐标模型求解 当 时当时模型三模型建立：根据我们的观测，方案一的走私船不能够到达p点又因为辑私船在x，y方向的速度分别为，既需求sink和cosk(1) 圆弧半径，缉私船的速度，两船初始距离(2)假设走私船和辑私船在t时在在坐标轴上的点为zm，n qx，y，通过船经过的圆心角度p来表示点zP=a*t*360/(2*pi*r)r =a p= 180*t/pim=a-a*cos(p)n=0.8*a+a*sin(p)(3)zp连线与水

8、平的夹角设为k，构成一个三角形，需求三边，根据两座标，可得，两直角边为a-a*cos(p)-x和0.8*a+a*sin(p)-y，勾股定理得第三边为4编程时记x1=x, x(2)=y, x=(x(1),x(2),a=20方案一如下：缉私船在方向的速度分别为。在方案二如下：缉私船在方向的速度分别为。点的坐标为即为从两初值表中可看出，走私船按照方案一逃跑在点处就被抓到了，而按照方案二逃跑那么是在点处被抓。所以按照方案一不可能到达点，按方案二可以到达点模型四走私船速度与时间的关系海里/小时缉私船以速度匀速追击 海里/小时模型建立建立判断模型：此模型用于确定走私船在何时改变方向利用数学软件matlab

9、求得数值解知道见附件四1，假设是走私船一直沿着正向逃跑，缉私船需要4.1小时才可追上走私船。因此我们在建立航线模型时假设走私船在小时时改变航向建立航线模型缉私船在方向的速度分别为从数值解中，可以得知在该种情况下，缉私船在发现走私船开场追击走私船需要3.9小时才可以追上走私船。参考文献1 笑缘、国勇.数学建模.：中国财政经济，2011.1附件第一题functiondx=js(t,x)d=20;b=1.5*d;c=2*d;k=-5*cos(t/5)+20*t+5;s=sqrt(k-x(1)2+(c-x(2)2);dx=b*(k-x(1)/s;b*(c-x(2)/s; ts=0:0.05:2.8;

10、x0=0,0; options=odeset(reltol,1e-6,abstol,1e-9); t,x=ode45(js,ts,x0,options); d=20;b=1.5*d;c=2*d; plot(t,x(:,1),-,t,x(:,2),-),gridts=0:0.05:2.8;x0=0,0;options=odeset(reltol,1e-6,abstol,1e-9);t,x=ode45(js,ts,x0,options);d=20;b=1.5*d;c=2*d;zx=-5*cos(t/5)+20*t+5;zy=c*ones(length(t),1);plot(t,x(:,1),-,t

11、,x(:,2),-),gridpause,plot(zx,zy,x(:,1),x(:,2),grid 00040.000000.05000.75000.029941.29902.39980.10001.50000.119442.59814.79800.15002.25000.267843.89717.19340.20003.00000.474545.19629.58440.25003.75000.738546.495211.96980.30004.50001.059147.794214.34820.35005.25001.435149.093316.71850.40006.00001.8655

12、50.392319.07960.45006.75002.348851.691321.43030.50007.50002.883952.990423.76990.55008.25003.469354.289426.09740.60009.00004.103455.588528.41200.65009.75004.784756.887530.71320.700010.50005.511558.186533.00050.750011.25006.282159.485635.27340.800012.00007.094860.784637.53160.850012.75007.947662.08363

13、9.77490.900013.50008.838863.382742.00330.950014.25009.766464.681744.21671.000015.000010.728565.980846.41551.050015.750011.723067.279848.59971.100016.500012.747968.578850.76981.150017.250013.801169.877952.92641.200018.000014.880571.176955.06991.250018.750015.983972.476057.20121.300019.500017.109073.7

14、75059.32121.350020.250018.253475.074061.43071.400021.000019.414976.373163.53091.450021.750020.590677.672165.62321.500022.500021.777878.971177.70901.550023.250022.973380.270279.79011.600024.000024.000081.569281.56921.650024.750024.750082.868382.86821.700025.500025.500084.167384.16731.750026.250026.25

15、0085.466385.46631.800027.000027.000086.765486.76541.850027.750027.750088.064488.06441.900028.500028.500089.363489.36341.950029.250029.250090.662590.66242.000030.000030.000191.961591.9615第二题function dx=js(t,x)a=30;b=1.6*a;c=a;d=5/36*pi;s=sqrt(a*t*cos(d)-x(1)2+(a*t*sin(d)+c-x(2)2);dx=b*(a*t*cos(d)-x(1

16、)/s;b*(a*t*sin(d)+c-x(2)/s; ts=0:0.05:2; x0=0,0; options=odeset(reltol,1e-6,abstol,1e-9); t,x=ode45(js,ts,x0,options); a=30; b=1.6*a; c=a; d=5/36*pi; zx=a*t*cos(d); zy=(c+a*t*sin(d).*ones(length(t),1); t,zx,x(:,1),zy,x(:,2)tzt0003000.051.35950.055130.63392.39920.12.71890.222831.26794.79310.154.07840

17、.506731.90187.1760.25.43780.909932.53579.54160.256.79731.434833.169611.88320.38.15682.082933.803614.19370.359.51622.85534.437516.46580.410.87573.750935.071418.6920.4512.23524.76935.705320.8650.513.59465.907136.339322.97770.5514.95417.161436.973225.02350.616.31358.527337.607126.99650.6517.6739.999138

18、.241128.89180.719.032511.570438.87530.70560.7520.391913.233939.508932.43520.821.751414.981940.142834.07930.8523.110816.806640.776835.6380.924.470318.699841.410737.11280.9525.829820.653542.044638.5064127.189222.659942.678539.82311.0528.548724.711543.312541.06831.129.908226.80143.946442.24881.1531.267

19、628.921644.580343.37261.232.627131.066545.214344.44921.2533.986533.229245.848245.48971.335.34635.34646.482146.48191.3536.705536.705447.11647.1161.438.064938.064947.7547.74991.4539.424439.424448.383948.38381.540.783940.783849.017849.01781.5542.143342.143349.651749.65161.643.502843.502750.285750.28561

20、.6544.862244.862250.919650.91951.746.221746.221651.553551.55351.7547.581247.581152.187552.18731.848.940648.940652.821452.82141.8550.300150.353.455353.45531.951.659551.659554.089254.08921.9553.01953.01954.723254.7231254.378554.378355.357155.3572plot(t,x(:,1),-,t,x(:,2),-),gridpause,plot(zx,zy,x(:,1),x(:,2),grid第三题functiondx=js(t,x)a=20;b=1.4*a;c=0.8*a;p=180*t/pi;s=sqrt(a-a*cos(p)-x(1)2+(0.8*a+a*sin(p)-x(2)2)dx=b*(0.8*a+a*sin(p)-x(2)/s;b*(a-a*cos(p)-x(1)/s;ts=0:0.02:3;x0=0,0;options=odeset(reltol,1e-6,abstol,1e-9);t,x=ode45(js,ts,x0,options);. word.zl