目录
1.问题说明
2.问题分析
3.代码分析
1.问题说明
这个是国赛的真题,我们这个里面只是浅谈,就是对于这个里面运用的过程仿真的思路进行说明,这个探测的波束问题实际上也是一个简单的过程仿真问题,也是需要去进行作图的,和之前的这个4个人相互追逐的问题属于一类问题--过程仿真问题;
实际赛题里面还考虑了重合率等等的相关信息,我们这个里面只是简单的说明一下仿真思想的体现,简单运用,编写过程的程序,所以不会考虑重合率等其他的问题;
2.问题分析
我们的这个游艇实际上就是来回的进行巡回的过程,我们需要计算经过几次可以全部检测完成这片的海域;
这个实际上就是一个数学问题,我们思路就是构建直线方程,联立方程求解点的坐标,根据点的坐标计算距离,知道这个水平行走距离大于L的时候,这个海域就会全部被检测完成;
这关键就是求解x1,y1和x2,y2点的坐标,这个时候就需要构建直线方程,就是简单的点斜式,利用已知的直线方程,这个里面的直线方程包括这个斜面的直线方程,还有就是x1y1这个点左边的直线方程,和x1y1右边的直线方程,左边的斜率使用的是kl表示的,kr表示的就是右边的斜率;
3.代码分析
clc;
clear;%海底的深度大小
L=10000;%%这个就是坡度的偏角
alfa=pi/100;%探测的张角大小
beta=pi/3;%%海水的深度
h=1000;%%海底平面的斜率
k=tan(alfa);%%左边的直线的斜率
kl=tan(pi/2-beta/2);%%右边的直线的斜率
kr=tan(beta/2-pi/2);%%左侧初始探测点的坐标
x2=0;
y2=0;%%初始点测量船的坐标
x1=(h+kl*x2-y2)/kl; %%这个就是直线方程联立求解的点的坐标
y1=h;xx1=x1; %%这四个式子就是把新的数据更新到我们的向量里面
xx2=x2; %%后面的作图的时候,需要用到这个向量
yy1=y1;
yy2=y2;for i=1:1000%%这个就是联立求解的点的坐标,上面已经出现过x2=(y1-kr*x1)/(k-kr);y2=k*x2;%%根据x2的坐标确定求解的x1=(h+kl*x2-y2)/kl;y1=h;%%不断地添加新的数据xx1=[xx1 x1];xx2=[xx2 x2];yy1=[yy1 y1];yy2=[yy2 y2];%%循环的终止的条件if x2>=Lbreak;end
end%%下面的这个就是绘制的斜面的图形,就只有一条线
plot([0 L],[0,L*tan(alfa)]);
hold on
plot([0 L],[h h],'k');%%下面是使用的循环语句绘制探测波束的图形
%%实际上就是一个迭代的过程,理解清楚就好了
for i=1:length(xx1)-1plot([xx1(i) xx2(i)],[yy1(i) yy2(i)],'r');plot([xx1(i) xx2(i+1)],[yy1(i) yy2(i+1)],'r');
end
相关的代码解释和注意说明已经作为注释写了进去,请注意查看理解;~~