CN104808659B - 船舶大圆航法航行航迹偏差的测定方法 - Google Patents
船舶大圆航法航行航迹偏差的测定方法 Download PDFInfo
- Publication number
- CN104808659B CN104808659B CN201510089844.2A CN201510089844A CN104808659B CN 104808659 B CN104808659 B CN 104808659B CN 201510089844 A CN201510089844 A CN 201510089844A CN 104808659 B CN104808659 B CN 104808659B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- point
- calculating
- starting point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000003556 assay Methods 0.000 title abstract 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title abstract 3
- 238000000034 method Methods 0.000 claims abstract description 70
- 238000012937 correction Methods 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 19
- 230000000694 effects Effects 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种船舶大圆航法航行航迹偏差的测定方法,该方法分为五个步骤:计算期望航线两点间最短距离;计算起点到当前位置间的航线最短距离;寻找期望航线上辅助参考点;计算辅助修正偏差角;计算航迹偏差。本发明采用在计划航线上取参考点的方法求取航迹偏差,有利于船舶快速修正航向,真正缩短航行时间;应用本航迹偏差的测定方法,可以更实时精确地计算出大圆航法航行时的航迹偏差,有利于提高船舶航迹控制精度,并以更经济更可靠的航行路线指导完成整个航段航行,达到节能减排的目的。
Description
技术领域
本发明属于船舶航迹与航向控制技术领域,具体涉及一种船舶大圆航法航行航迹偏差的测定方法。
背景技术
随着GPS定位系统的逐步商业化,航迹控制也逐步在自动舵中得到实现。如果船舶在航行过程中可以始终趋近期望航线航行,那么既可以减小航迹偏差提高航行效率,又可以节约燃料低碳环保。
求取大圆航法航迹偏差的传统方法一般采用分段方法,即将大圆航线等分为若干段后,将每段近似看做等航向线,再利用几何关系求解航迹偏差;该方法计算比较粗糙,计算的精度对分段个数的依赖性较强,分段数越多精确度越高;但分段数的选取并没有统一的标准,一般依靠经验选取,在航行距离很远时存在较大误差。如从美国东北部的波士顿到荷兰西部的鹿特丹,大圆航线距离约为1767海里,将大圆航线等分为若干段,取每段中间点切线的垂线距大圆航线500±20米的点作为船舶当前位置,计算航迹偏差。若将大圆航线等分4段,每段传统算法的航迹偏差约为14284米;若将大圆航线等分20段,每段传统算法的航迹偏差约为1040米;若将大圆航线等分50段,每段传统算法的航迹偏差约为605米。具体见附表1。
为了提高航迹控制精度、简化操作过程,迫切需要研究一种精度高、不依赖分段数的航迹偏差测定方法,便于船舶在大圆航法航行时的制导,减小航迹偏差,缩短航行距离,相应减少航行时间,提高航迹自动舵的性能。本航迹偏差测定方法已应用于某型船舶自动操舵仪,并取得了良好的航迹控制效果。
发明内容
本发明的目的在于提供一种船舶大圆航法航行航迹偏差的测定方法,旨在改善船舶的航迹控制效果,达到节约能源低碳环保的目的,使海洋运输更加安全可靠。
本发明是这样实现的,一种船舶大圆航法航行航迹偏差的测定方法,首先提出两个方法:方法1,已知两点经纬度,计算其大圆航法两点间的最短距离以及两点计划航向;方法2,已知起点经纬度、起点计划航向和大圆航法航行两点间最短距离,计算终点经纬度以及终点计划航向;其次结合上述两个方法计算大圆航法航迹偏差,具体包括以下几个步骤:
步骤一、计算期望航线两点间最短距离;
步骤二、计算起点到当前位置间的航线最短距离;
步骤三、寻找期望航线上辅助参考点;
步骤四、计算辅助修正偏差角;
步骤五、计算航迹偏差,并回到步骤二。
方法1,已知起点和终点的经纬度,计算大圆航法两点间的最短距离以及起点和终点的计划航向,具体方法为:
已知起点A和终点B的经纬度坐标分别为和计算起点A和终点B的归化纬度UA,UB,计算经度差λAB,计算大圆航法两点间的最短距离SAB以及起点A和终点B的计划航向ΨA,ΨB,某点的计划航向即该点所在弧切线方向与正北向的夹角;
1)、计算起点A和终点B的归化纬度:
其中和分别为起点A和终点B的经纬度坐标,f为地球扁率;
2)、计算经度差λAB:
λAB=ΔλAB+(1-M)fsinΨAB{θAB+MsinθAB[cos2θm+McosθAB(-1+2cos22θm)]};
其中经度差初值ΔλAB=λB-λA,θm为赤道到A,B所在弧线连线中点间的角距,ΨAB为起点A与终点B所在弧的延长线与赤道交点的航向角,θAB为起点A和终点B间的角距,M为引入的参数;
3)、计算大圆航法起点A和终点B间的最短距离SAB以及起点A和终点B的计划航向ΨA,ΨB:
SAB=bN(θAB-ΔθAB);
其中b为地球短半轴长,θAB为起点A和终点B间的角距,ΔθAB为起点A和终点B之间的角距偏差。
方法2,已知起点经纬度、起点计划航向和大圆航法航行两点间最短距离,计算终点经纬度以及终点计划航向,具体方法为:
已知起点A经纬度坐标起点A计划航向ΨA和大圆航法航行两点间最短距离SAB;计算起点A的归化纬度UA以及赤道到A点的角距θA,计算起点A和终点B间角距θAB,计算终点B的经度值λB,纬度值以及B点计划航向ΨB;
1)、计算起点A归化纬度UA以及赤道到A点的角距θA:
其中为起始点A的经纬度坐标,f为地球扁率,ΨA为A点计划航向;
2)、计算得到起点A和终点B间角距θAB为:
其中S为从A点出发的大圆航最短距离,ΔθAB为起点A和终点B之间的角距偏差;
3)、计算终点B经度值λB,纬度值以及B点计划航向ΨB:
λB=λA+ΔλAB;
步骤一,所述的计算期望航线两点间最短距离包括:
已知起点A和终点B经纬度坐标分别为计算起点A和终点B的计划航向ΨA和ΨB;起点A和终点B之间的最短距离SAB。
步骤二,所述的计算起点到当前位置间的航线最短距离包括:
已知船舶起点A和当前位置C的经纬度坐标和计算起点A到船舶当前位置C的最短距离SAC。
步骤三,所述的寻找期望航线上辅助参考点:
在计划航线上取一参考点使其与起点A的最短距离等于船舶当前位置到起点的距离;过D点做AB的切线,其计划航向为ΨD。
步骤四,所述的计算辅助修正偏差角的方法中直线CD与正北方向夹角为ΨC;则计算得到辅助修正偏差角β方法为:
β=ΨD-ΨC
步骤五,所述的计算航迹偏差中,计算CD长度为SCD,最终计算得到航迹偏差η:
η=SCD*sinβ
效果汇总
本发明的优点在于:
(1)、已知起点、终点经纬度值便可求出起点、终点计划航向,以及大圆航线两点间的最短长度;
(2)、已知起点经纬度、起点计划航向和大圆航线最短距离;便可求出终点经纬度值,以及终点计划航向;
(3)、采用在计划航线上取参考点的方法求取航迹偏差,有利于船舶快速修正航向,真正缩短航行时间;
(4)、应用在大圆航法航行时求航迹偏差的测定方法,可以更实时精确地计算出航迹偏差,有利于航迹的控制。
附图说明
图1为基于大圆航法航迹偏差测定方法流程图;
图2为椭球面上的参量表示图;
图3为航迹偏差解析计算图;
图4为航迹偏差解析计算局部放大图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1至图4所示,本发明是这样实现的,一种船舶大圆航法航行航迹偏差的测定方法,首先介绍两个方法:方法1,已知两点经纬度,计算其大圆航法两点间的最短距离以及两点计划航向;方法2,已知起点经纬度、起点计划航向和大圆航法航行两点间最短距离,计算终点经纬度以及终点计划航向;其次结合上述两个方法计算大圆航法航迹偏差,具体包括以下几个步骤:
S101:计算期望航线两点间最短距离;
S102:计算起点到当前位置间的航线最短距离;
S103:寻找期望航线上辅助参考点;
S104:计算辅助修正偏差角;
S105:计算航迹偏差,并回到步骤S102。
方法1,结合图2,已知起点A和终点B的经纬度和计算大圆航法起点A和终点B间的最短距离SAB以及起点A和终点B的计划航向ΨA,ΨB:
1)、计算起点A和终点B的归化纬度UA,UB:
式中和分别为起点A和终点B的经纬度坐标,f为扁率。
2)、计算经度差λAB:
a)、计算起点A和终点B间角距θAB:
其中,
F为符号值,可取±1。
b)、计算赤道与A,B连线中点线间的角距θm:
其中,
式中,ΨAB为起点与终点所在弧延长线与赤道交点的航向角。
c)、通过迭代法,计算经度差λAB:
λnew=ΔλAB+(1-M)fsinΨAB{θAB+MsinθAB[cos2θm+McosθAB(-1+2cos22θm)]} (6)
式中,起点A和终点B的经度差初值为ΔλAB=λB-λA,λAB为经度差,初始令λAB=ΔλAB
引入参数M:
则迭代误差errλ=λAB-λnew,当λAB的迭代误差errλ小于设定的误差界限(本发明中误差界限设定为1×10-12)时,计算结束,此时λnew即为所求的经度差λAB;否则当λ≥0时,令F=1;当λ<0时,令F=-1;将F对应值带入公式(3)中,重新迭代计算。
3)、计算大圆航法两点间的最短距离SAB以及起点A和终点B的计划航向ΨA和ΨB(某点计划航向即该点所在弧的切线与正北方向夹角):
a)、计算起点A和终点B间的最短距离SAB:
SAB=bN(θAB-ΔθAB) (8)
其中,ΔθAB为起点A和终点B间的角距偏差,具体计算方法如下:
分别引入参数N,P:
式中,a为地球长半轴长度,b为地球短半轴长度。
b)、计算起点A计划航向ΨA:
c)、计算终点B计划航向ΨB:
方法2,结合图2,已知起点A经纬度坐标起点A计划航向ΨA和大圆航法航行两点间最短距离SAB,计算终点B经纬度坐标以及终点计划航向ΨB:
1)、计算起点的归化纬度UA以及赤道到A点的角距θA:
式中,ΨA为起点A的计划航向。
2)、计算起点A和终点B间角距θAB:
分以下几个步骤求θAB:
a)、计算赤道与起点A和终点B连线中点线间的角距θm:
式中,θA为赤道到A点的角距。
令起点A和终点B间角距θAB的初值为SAB为从起点A出发到终点B的大圆航最短距离,b为地球短半轴长度。
为求θAB的初值,引入参数N:
式中,
式中,ΨAB为起点与终点所在弧延长线与赤道交点的航向角,ΨA为起点计划航向,UA为起点归化纬度。
b)、通过迭代法,计算起点A和终点B间角距θAB:
迭代误差为errθ=θnew-θAB,当θAB的迭代误差errθ小于设定的误差界限(本发明中误差界限设定为1×10-12)时,计算结束,此时θnew即为所求的起点A和终点B间角距θAB;否则令θAB=θnew带入式(16)继续计算。
其中,起点A和终点B间的角距偏差ΔθAB的计算方法为:
为求ΔθAB引入参数P:
式中,
式中,ΨAB为起点与终点所在弧延长线与赤道交点的航向角,ΨA为起点计划航向,UA为起点归化纬度。
3)、计算终点B的经度值λB,纬度值以及B点计划航向ΨB:
a)、计算终点B纬度值
式中,f为扁率,ΨAB为起点与终点所在弧延长线与赤道交点的航向角,ΨA为起点计划航向,UA为起点归化纬度,θAB为起点A和终点B间的角距。
b)、计算终点B的经度值λB,计划航向ΨB:
其中,起点A和终点B间的经度差初值ΔλAB的计算方法为:
λAB=λAB-(1-M)fsinΨAB{θAB+MsinθAB[cos2θm+McosθAB(-1+2cos22θm)]} (27)
为求ΔλAB引入参数M:
为求ΔλAB引入经度差λAB的计算方法为:
应用上述两个方法,按一下五个步骤实现对大圆航法航迹偏差的计算:
步骤一,计算期望航线两点间最短距离。
结合方法1,计算起点和终点的计划航向ΨA,ΨB,起点A和终点B之间的最短距离SAB:
SAB=bN(θAB-ΔθAB) (31)
步骤二,计算起点到当前位置间的最短距离。
结合方法1,计算船舶起点到当前位置间的最短距离SAC:
SAC=bN(θAC-ΔθAC) (32)
步骤三,寻找期望航线上辅助参考点:
结合图3,在计划航线上取一参考点使其与起点A的最短距离等于船舶当前位置到起点的距离;过D点做AB的切线,其计划航向为ΨD。
步骤四,计算辅助修正偏差角β:
1)、计算船舶当前位置的计划航向ΨC:
结合方法1,计算实际位置的计划航向ΨC为:
2)、计算参考点D的经纬度坐标和计划航向ΨD:
结合方法2计算出参考点D的经纬度以及D点计划航向ΨD;
其中,D点的经度为λD=λA+ΔλAD,D点的纬度以及D点计划航向ΨD的计算方法分别为:
其中ΔλAD为起点A与参考点D间的经度差初值。
3)、计算直线CD与正北方向夹角即当前位置计划航向ΨC:
由于本发明的算法可以实时计算航迹偏差,使船舶航迹不断趋于期望航线,所以CD是逐渐变短,C点在CD上的切线逐渐逼近直线CD,本发明中将C点在CD上的切线看做是直线CD。故结合方法1,可以计算出直线CD与正北方向夹角即当前位置计划航向ΨC:
4)、计算得到辅助修正偏差角β:
结合图3和图4,可以计算出辅助修正偏差角β:
β=ΨD-ΨC (37)
步骤五,计算航迹偏差η:
结合图3和图4,过C点做D点在AD上切线的垂线,记为η;η即为所要求的航迹偏差(由于本发明的算法可以实时计算航迹偏差,使船舶航迹不断趋于期望航线,所以CD逐渐变短,在计算航迹偏差时,本发明将CD看做是直线,便于计算)。
η=SCD*sinβ (38)
其中,结合方法1,CD的长度SCD的计算方法为:
SCD=bN(θCD-ΔθCD) (39)
其中θCD为实际位置C和参考点D间的角距,ΔθCD为实际位置C和参考点D间的角距偏差。
参考国际标准62065中Table G.4-Scenario 4中第7和第8航路点间的大圆航线,分别用本测定方法和传统算法进行航迹偏差的计算,计算结果如附表1所示。
附表1本测定方法和传统算法计算航迹偏差结果汇总表
本发明中用到的部分符号、物理意义及其单位汇总如附表2所示:
附表2变量汇总表
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性的劳动即可做出的各种修改或变形仍在本发明的保护范围之内。
Claims (1)
1.一种船舶大圆航法航行航迹偏差的测定方法,其特征在于,该船舶大圆航法航行航迹偏差的测定方法,结合大圆航法航线上已知两点的经纬度,计算两点间的最短距离和两点各自计划航向,以及已知起点经纬度、起点计划航向和大圆航法航行两点间最短距离,计算终点经纬度和终点计划航向,计算大圆航法航迹偏差,具体包括以下几个步骤:
步骤一、计算期望航线两点间最短距离;
步骤二、计算起点到当前位置间的航线最短距离;
步骤三、寻找期望航线上辅助参考点;
步骤四、计算辅助修正偏差角;
步骤五、计算航迹偏差,并回到步骤二;
其中,
1)已知起点和终点的经纬度,计算大圆航法两点间的最短距离以及起点和终点的计划航向,具体方法为:
已知起点A和终点B的经纬度坐标分别为和计算起点A和终点B的归化纬度UA,UB计算经度差λAB,计算大圆航法起点A和终点B间的最短距离SAB以及起点A和终点B的计划航向ψA,ψB,某点的计划航向即该点所在弧切线方向与正北向的夹角;
2)计算起点A和终点B的归化纬度UA和UB:
其中和分别为起点A和终点B的经纬度坐标,f为地球扁率;
3)计算经度差λAB:
λAB=ΔλAB-(1-M)fsinψAB{θAB+MsinθAB[cos2θm+McosθAB(-1+2cos22θm)]};
其中经度差初值ΔλAB=λB-λA,θm为赤道到A,B所在弧线连线中点间的角距,ψAB为起点A与终点B所在弧的延长线与赤道交点的航向角,θAB为起点A和终点B间的角距,M为引入的参数;
4)计算大圆航法起点A和终点B间的最短距离SAB以及起点A和终点B的计划航向ψA,ψB:
SAB=bN(θAB-ΔθAB);
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>N</mi>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mn>16384</mn>
</mfrac>
<mo>{</mo>
<mn>4096</mn>
<mo>+</mo>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>&lsqb;</mo>
<mo>-</mo>
<mn>768</mn>
<mo>+</mo>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>320</mn>
<mo>-</mo>
<mn>175</mn>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mn>1024</mn>
</mfrac>
<mo>{</mo>
<mn>256</mn>
<mo>+</mo>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>&lsqb;</mo>
<mo>-</mo>
<mn>128</mn>
<mo>+</mo>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>74</mn>
<mo>-</mo>
<mn>47</mn>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:
<mrow>
<msub>
<mi>&psi;</mi>
<mi>A</mi>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>cosU</mi>
<mi>B</mi>
</msub>
<msub>
<mi>sin&lambda;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>cosU</mi>
<mi>A</mi>
</msub>
<msub>
<mi>sinU</mi>
<mi>B</mi>
</msub>
<mo>-</mo>
<msub>
<mi>sinU</mi>
<mi>A</mi>
</msub>
<msub>
<mi>cosU</mi>
<mi>B</mi>
</msub>
<msub>
<mi>cos&lambda;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>cosU</mi>
<mi>A</mi>
</msub>
<msub>
<mi>sin&lambda;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<msub>
<mi>sinU</mi>
<mi>A</mi>
</msub>
<msub>
<mi>cosU</mi>
<mi>B</mi>
</msub>
<mo>+</mo>
<msub>
<mi>cosU</mi>
<mi>A</mi>
</msub>
<msub>
<mi>sinU</mi>
<mi>B</mi>
</msub>
<msub>
<mi>cos&lambda;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中a为地球长半轴长,b为地球短半轴长,θAB为起点A和终点B间的角距,ΔθAB为起点A和终点B之间的角距偏差;
5)已知起点经纬度、起点计划航向和大圆航法航行两点间最短距离,计算终点经纬度以及终点计划航向,具体方法为:
已知起点A经纬度坐标起点A计划航向ψA和大圆航法起点A和终点B间的最短距离SAB;计算起点A的归化纬度UA以及赤道到A点的角距θA,计算起点A和终点B间角距θAB,计算终点B的经度值λB,纬度值以及B点计划航向ψB;
6)计算起点A归化纬度UA以及赤道到A点的角距θA:
<mrow>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mfrac>
<mrow>
<msub>
<mi>tanU</mi>
<mi>A</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>cos&psi;</mi>
<mi>A</mi>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
其中为起始点A的经纬度坐标,f为地球扁率,ψA为A点计划航向;
7)计算起点A和终点B间角距θAB:
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>S</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
<mrow>
<mi>b</mi>
<mi>N</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>&Delta;&theta;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
其中SAB为大圆航法起点A和终点B间的最短距离,ΔθAB为起点A和终点B之间的角距偏差;b为地球短半轴长;N为待定参数;
8)计算终点B经度值λB,纬度值以及B点计划航向ψB:
λS=λA+ΔλAB;
<mrow>
<msub>
<mi>&psi;</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mfrac>
<mrow>
<msub>
<mi>sin&psi;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>-</mo>
<msub>
<mi>sinU</mi>
<mi>A</mi>
</msub>
<msub>
<mi>sin&theta;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>cosU</mi>
<mi>A</mi>
</msub>
<msub>
<mi>cos&theta;</mi>
<mrow>
<mi>A</mi>
<mi>B</mi>
</mrow>
</msub>
<msub>
<mi>cos&psi;</mi>
<mi>A</mi>
</msub>
</mrow>
</mfrac>
</mrow>
9)计算起点到当前位置间的航迹最短距离:
已知船舶起点A和当前位置C的经纬度坐标和计算起点A到船舶当前位置C的最短距离SAC;
10)找期望航线上辅助参考点:
在计划航线上取一参考点使其与起点A的最短距离等于船舶当前位置到起点的距离;过D点做AB的切线,其计划航向为ψD;
计算辅助修正偏差角的方法中,直线CD与正北方向夹角为ψC;则计算得到辅助修正偏差角β:
β=ψD-ψC
计算航迹偏差中,计算CD长度为SCD,最终计算得到航迹偏差η:
η=SCD*sinβ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510089844.2A CN104808659B (zh) | 2015-02-27 | 2015-02-27 | 船舶大圆航法航行航迹偏差的测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510089844.2A CN104808659B (zh) | 2015-02-27 | 2015-02-27 | 船舶大圆航法航行航迹偏差的测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104808659A CN104808659A (zh) | 2015-07-29 |
CN104808659B true CN104808659B (zh) | 2017-10-20 |
Family
ID=53693578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510089844.2A Expired - Fee Related CN104808659B (zh) | 2015-02-27 | 2015-02-27 | 船舶大圆航法航行航迹偏差的测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104808659B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106643732B (zh) * | 2016-12-29 | 2019-12-27 | 飞友科技有限公司 | 一种基于网络地图模拟描绘两地间航线的方法 |
CN106950955B (zh) * | 2017-03-17 | 2019-11-05 | 武汉理工大学 | 基于改进los引导算法的船舶航迹控制方法 |
CN106950959A (zh) * | 2017-04-01 | 2017-07-14 | 武汉理工大学 | 新型的船舶视距制导方法 |
CN106970628A (zh) * | 2017-05-19 | 2017-07-21 | 苏州寅初信息科技有限公司 | 一种智能无人船自动交易的控制方法及其无人船 |
CN107563541A (zh) * | 2017-07-28 | 2018-01-09 | 广东华风海洋信息系统服务有限公司 | 一种节省时间的船舶航线规划算法 |
CN107741232B (zh) * | 2017-10-20 | 2020-02-04 | 北京润科通用技术有限公司 | 一种测量船的航路规划方法及装置 |
CN108801262B (zh) * | 2018-05-31 | 2021-09-03 | 青岛中乌特种船舶研究设计院有限公司 | 一种船舶自动航行控制器航路规划与纠偏修正方法 |
CN111736582B (zh) * | 2019-03-19 | 2024-06-21 | 北京奇虎科技有限公司 | 路径处理方法、装置、电子设备及计算机可读存储介质 |
CN110163438B (zh) * | 2019-05-24 | 2023-07-28 | 上海海事大学 | 一种大圆航线距离效益的评估方法 |
CN110849376A (zh) * | 2019-12-04 | 2020-02-28 | 中国直升机设计研究所 | 基于公式的大圆航线二维地图显示方法 |
CN111782755B (zh) * | 2020-07-20 | 2021-05-25 | 中国人民解放军国防科技大学 | 基于虚拟网格字典的目标行进类意图识别方法和装置 |
CN112699531B (zh) * | 2020-12-08 | 2023-06-09 | 浩亚信息科技有限公司 | 一种通用航空低空飞行目视参考点建立方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090103720A1 (en) * | 2005-06-07 | 2009-04-23 | Manoj Karayil Thekkoott Narayanan | Method and system for secure and anti jamming wireless communication with high spectral efficiency |
US7860646B2 (en) * | 2007-04-16 | 2010-12-28 | The Boeing Company | Method and apparatus for routing ocean going vessels to avoid treacherous environments |
CN102306352A (zh) * | 2011-08-31 | 2012-01-04 | 上海海事大学 | 一种船舶调配方案自动化生成方法 |
CN103531043B (zh) * | 2013-09-27 | 2016-02-10 | 重庆大学 | 基于航线匹配的船舶点到点航行时间估算方法 |
-
2015
- 2015-02-27 CN CN201510089844.2A patent/CN104808659B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104808659A (zh) | 2015-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104808659B (zh) | 船舶大圆航法航行航迹偏差的测定方法 | |
CN104960520B (zh) | 基于Pure Pursuit算法的预瞄点确定方法 | |
CN105511464B (zh) | 一种船舶自主航行控制系统及控制方法 | |
CN103700286B (zh) | 一种舰载无人机自动着舰引导方法 | |
CN110609553B (zh) | 一种用于铺管船舶圆弧路径的los导引控制方法 | |
CN106950955B (zh) | 基于改进los引导算法的船舶航迹控制方法 | |
CN101872195B (zh) | 船舶海上航行航迹解析偏差生成方法 | |
JP4804032B2 (ja) | 船舶用自動航法援助システム | |
CN105843234A (zh) | 一种uuv对圆形障碍物几何绕行的二维航路规划方法 | |
CN103940420B (zh) | 等航程大圆航线设计方法 | |
CN107741232B (zh) | 一种测量船的航路规划方法及装置 | |
CN115562266A (zh) | 基于变参数视线法的无人船航迹控制方法及存储介质 | |
CN103376452A (zh) | 一种用单台声信标修正水下机器人位置误差的方法 | |
CN105892487A (zh) | 一种无人机8字形航迹控制方法 | |
CN108801262B (zh) | 一种船舶自动航行控制器航路规划与纠偏修正方法 | |
CN105549611A (zh) | 自主水下机器人的海洋环境自识别的航迹精确跟踪方法 | |
CN106950959A (zh) | 新型的船舶视距制导方法 | |
CN103557872A (zh) | 一种rnp中的综合系统误差实时计算方法 | |
CN111026116A (zh) | 船舶操纵控制方法、装置及电子设备 | |
CN117291020B (zh) | 一种欠驱动船舶的靠泊导引方法 | |
CN104794311A (zh) | 一种时间自由条件下的全局最优双脉冲轨道转移方法 | |
CN109631902B (zh) | 基于海洋测量的船上线在近点情况下最佳路径规划方法 | |
CN108489490B (zh) | 海测船上测线导航路径规划 | |
CN111158366B (zh) | 基于图搜索和几何曲线融合的路径规划方法 | |
US4768153A (en) | Automatic pilot system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171020 |
|
CF01 | Termination of patent right due to non-payment of annual fee |