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

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） 设在任意时刻缉私船的坐标*，y，走私船到达Qat，c点，直线PQ与缉私船航线相切，切线PQ与y轴正方向夹角为（5） 设常数d=20，缉私船b=1.5d=30海里/小时，初始距离c=2d=40海里符号定义1、任意时刻缉私船的坐标*，y2、走私船到达Qat，c点3、切

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

7、关系可以得出，时间t时走私船的位置：2缉私船的速度，两船的初始位置c=2d。缉私船在*，y方向是速度分别为 得到微分方程：初始条件为：，3把与代入1中分别得出走私船的坐标模型求解 当 时当时模型三模型建立：根据我们的观测，方案一的走私船不能够到达p点又因为辑私船在*，y方向的速度分别为，既需求sink和cosk(1) 圆弧半径，缉私船的速度，两船初始距离(2)假设走私船和辑私船在t时在在坐标轴上的点为zm，n q*，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)-*和0.8*a+a*sin(p)-y，勾股定理得第三边为4编程时记*1=*, *(2)=y, *=(*(1),*(2),a=20方案一如下：缉私船在方向的速度分别为。在方案二如下：缉私船在方向的速度分别为。点的坐标为即为从两初值表中可看出，走私船按照方案一逃跑在点处就被抓到了，而按照方案二逃跑则是在点处被抓。所以按照方案一不可能到达点，按方案二可以到达点模型四走私船速度与时间的关系海里/小时缉私船以速度匀速追击 海里/小时模型建立建立判断模型：此模型用于确定走私船在何时改变方向利用数学软件matlab求得

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

10、s=0:0.05:2.8;>> *0=0,0;>> options=odeset('reltol',1e-6,'abstol',1e-9);>> t,*=ode45(js,ts,*0,options);>> d=20;b=1.5*d;c=2*d;>> plot(t,*(:,1),'-',t,*(:,2),'-'),grid>> ts=0:0.05:2.8;>> *0=0,0;>> options=odese

11、t('reltol',1e-6,'abstol',1e-9);>> t,*=ode45(js,ts,*0,options);>> d=20;b=1.5*d;c=2*d;>> z*=-5*cos(t/5)+20*t+5;>> zy=c*ones(length(t),1);plot(t,*(:,1),'-',t,*(:,2),'-'),gridpause,plot(z*,zy,*(:,1),*(:,2),grid 0  

12、0; 0      0   40.0000     0    0.0500    0.7500    0.0299   41.2990    2.3998    0.1000    

13、;1.5000    0.1194   42.5981    4.7980    0.1500    2.2500    0.2678   43.8971    7.1934    0.2000    3.0000

14、60;   0.4745   45.1962    9.5844    0.2500    3.7500    0.7385   46.4952   11.9698    0.3000    4.5000   

15、 1.0591   47.7942   14.3482    0.3500    5.2500    1.4351   49.0933   16.7185    0.4000    6.0000    1.8655 &#

16、160; 50.3923   19.0796    0.4500    6.7500    2.3488   51.6913   21.4303    0.5000    7.5000    2.8839   52.9904

17、60;  23.7699    0.5500    8.2500    3.4693   54.2894   26.0974    0.6000    9.0000    4.1034   55.5885   28.41

18、20    0.6500    9.7500    4.7847   56.8875   30.7132    0.7000   10.5000    5.5115   58.1865   33.0005    

19、;0.7500   11.2500    6.2821   59.4856   35.2734    0.8000   12.0000    7.0948   60.7846   37.5316    0.8500   12.75

20、00    7.9476   62.0836   39.7749    0.9000   13.5000    8.8388   63.3827   42.0033    0.9500   14.2500    9.76

21、64   64.6817   44.2167    1.0000   15.0000   10.7285   65.9808   46.4155    1.0500   15.7500   11.7230   67.2798  &

22、#160;48.5997    1.1000   16.5000   12.7479   68.5788   50.7698    1.1500   17.2500   13.8011   69.8779   52.9264    1.20

23、00   18.0000   14.8805   71.1769   55.0699    1.2500   18.7500   15.9839   72.4760   57.2021    1.3000   19.5000  &

24、#160;17.1090   73.7750   59.3212    1.3500   20.2500   18.2534   75.0740   61.4307    1.4000   21.0000   19.4149   76.3731

25、60;  63.5309    1.4500   21.7500   20.5906   77.6721   65.6232    1.5000   22.5000   21.7778   78.9711   77.7090   

26、 1.5500   23.2500   22.9733   80.2702   79.7901    1.6000   24.0000   24.0000   81.5692   81.5692    1.6500   24.7500

27、60;  24.7500   82.8683   82.8682    1.7000   25.5000   25.5000   84.1673   84.1673    1.7500   26.2500   26.2500   

28、85.4663   85.4663    1.8000   27.0000   27.0000   86.7654   86.7654    1.8500   27.7500   27.7500   88.0644   88.0644 &#

29、160;  1.9000   28.5000   28.5000   89.3634   89.3634    1.9500   29.2500   29.2500   90.6625   90.6624    2.0000   

30、30.0000   30.0001   91.9615   91.9615第二题function d*=js(t,*)a=30;b=1.6*a;c=a;d=5/36*pi;s=sqrt(a*t*cos(d)-*(1)2+(a*t*sin(d)+c-*(2)2);d*=b*(a*t*cos(d)-*(1)/s;b*(a*t*sin(d)+c-*(2)/s;>> ts=0:0.05:2;>> *0=0,0;>> options=odeset('reltol'

31、,1e-6,'abstol',1e-9);>> t,*=ode45(js,ts,*0,options);>> a=30;>> b=1.6*a;>> c=a;>> d=5/36*pi;>> z*=a*t*cos(d);>> zy=(c+a*t*sin(d).*ones(length(t),1);>> t,z*,*(:,1),zy,*(:,2)tzt0003000.051.35950.055130.63392.39920.12.71890.222831.26794.79310.154.0

32、7840.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.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.82982