CN108507570B - 舰载分布式局部基准状态检测与姿态信息重构方法 - Google Patents
舰载分布式局部基准状态检测与姿态信息重构方法 Download PDFInfo
- Publication number
- CN108507570B CN108507570B CN201810182388.XA CN201810182388A CN108507570B CN 108507570 B CN108507570 B CN 108507570B CN 201810182388 A CN201810182388 A CN 201810182388A CN 108507570 B CN108507570 B CN 108507570B
- Authority
- CN
- China
- Prior art keywords
- attitude
- local reference
- state detection
- distributed
- ship
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
- G01C21/203—Specially adapted for sailing ships
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种舰载分布式局部基准状态检测与姿态信息重构方法,利用船舶二维水弹性理论中的模态分析法,建立姿态预测模型,利用分布式姿态基准的姿态测量信息完成模型参数估计和姿态预测,实现对待监测局部基准的状态检测,如果局部基准出现故障,给出姿态信息重构结果。本发明提出的姿态预测模型,既考虑了刚体运动又考虑了船体角变形的影响,相对传统基于主惯导基准信息直接引射的方法,具有更高的姿态重构精度。
Description
技术领域
本发明涉及惯性导航领域,特别是一种舰载分布式局部基准状态检测与姿态信息重构方法。
背景技术
大型舰船上雷达、火炮、导弹等武器装备在工作时需要高精度姿态基准信息的支持。惯性导航设备由于可以提供实时、连续、可靠的高精度姿态信息,成为现代舰载武器系统不可或缺的一部分。舰船一般装备高精度主惯导姿态基准,如激光陀螺单轴/双轴旋转惯导系统,各战位点的坐标系都是以主惯导的坐标系为参考。但是,舰船的结构是挠性的,在日照、载荷变化和海浪作用下将产生角变形,如日照引起的静态角变形一天之内最大可达度量级,海浪引起的动态角变形也最大也可达到十几角分。由于安装差异和船体角变形的影响,直接利用高精度主惯导提供的姿态基准信息不能满足广泛布置于全舰武器系统的使用需求,各主要作战单元不得不配置与其作战需求相适应的惯性导航系统,用以确定各作战单元的坐标系,称为局部坐标基准或局部基准。
中国专利201410076442.4公开的“基于加速度计输出增量的捷联惯性导航方法”就是一种给出舰船作战单元姿态和位置信息的局部基准。局部基准依靠加速度计和陀螺仪等惯性器件测量出载体的线运动和角运动,由惯性导航方程得到载体相对于当地地理坐标系的姿态信息。但是,目前的局部基准通常由不同厂家生产和配备。根据局部基准所用的惯性元件不同,分为挠性陀螺局部基准、光纤陀螺局部基准、激光陀螺局部基准等。这些局部基准由于体制、工作方式和接口协议不同,相互之间独立工作,无法形成信息互通。当某一局部基准的性能下降或出现故障时,无有效的检测手段,且受船体变形的影响,直接复制主惯导的姿态信息将产生较大的误差。因此,基于局部基准独立工作的应用模式,对武器系统来说是非常脆弱的。
中国专利201710242262.2公开的“基于分布式惯性网络的高精度传递对准方法”,介绍一种基于分布式惯性网络的局部基准故障检测与隔离方法。该发明将所有节点局部基准测量的姿态信息传递到某一节点,采用广义似然比法,检测惯导的状态,当惯导出现故障时,将主惯导姿态补偿固定安装角后传递给故障节点。这种方法可以实现紧急情况下的姿态传递,但是忽略了船体角变形的影响,无法满足高精度的应用需求。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种舰载分布式局部基准状态检测与姿态信息重构方法,考虑船体角变形的影响,满足高精度的应用需求。
为解决上述技术问题,本发明所采用的技术方案是:一种舰载分布式局部基准状态检测与姿态信息重构方法,包括以下步骤:
1)构建分布式姿态基准网络,启动姿态基准设备,各姿态基准设备通过网络通信系统向信息处理计算机发送姿态测量结果;
2)根据船舶力学中的二维水弹性理论建立全舰姿态预测模型,利用所述全舰姿态预测模型和所述姿态基准设备的姿态测量结果,实时估计全舰姿态预测模型参数;利用所述姿态基准设备的姿态测量结果,对待监测局部基准的姿态进行预测;
3)如果待监测局部基准系统无姿态信息输出,或者姿态测量值与第2) 步预测得到的姿态值之差超过设定的阈值,则判定为局部基准状态异常;
4)根据第3)步的局部基准状态检测结果,如果局部基准出现状态异常,将第2)步得到的姿态预测结果发送至出现故障的局部基准战位点,实现姿态信息重构。
步骤1)中,所述分布式姿态基准网络包括依次双向连接的姿态基准系统、网络交换机和信息处理计算机。
步骤1)中,姿态测量结果表示为:
Θj(x1,t),Θj(x2,t),…,Θj(xm,t);
其中,Θj(xi,t),i=1,2,…,m表示第i个姿态基准系统的姿态测量结果,下标j=x,y,z分别表示横摇角、纵摇角和航向角,m为姿态基准系统的数量。
步骤2)中,全舰姿态预测模型表达式为:
Θj(xk,t)=a1jΘj(x1,t)+a2jΘj(x2,t)+…+amjΘj(xm,t);其中, Aj=[a1j a2j …amj]T, Ej=[θ1j(xk) θ2j(xk) … θ lj(xk)]T,1≤k≤m,θrj(x)为角振动模态,r表示模态的阶数,r=1表示刚体的角振动模态,r=2,…,l表示挠曲变形的角振动模态。
Zj=[Θj(xm,t0) Θj(xm,t0+h) … Θj(xm,t1)]T
Xj=[a1j a2j … a(m-1)j]T
步骤2)中,t时刻待监测局部基准的姿态预测值表示为:
步骤3)中,设定的阈值为9角分。
与现有技术相比,本发明所具有的有益效果为:本发明提出了一种局部姿态基准状态检测方法,可利用主惯导基准和分布式局部基准的姿态输出,预测待监测部位局部基准的姿态,为局部基准的状态检测提供了一种有效技术途径;本发明提出的姿态预测模型,既考虑了刚体运动又考虑了船体角变形的影响,相对传统基于主惯导基准信息直接引射的方法,具有更高的姿态重构精度;本发明充分利用舰船现有的主惯导基准和局部基准设备以及通信网络,只需增加一个信息处理计算机,具有成本低,实施简便的特点。
附图说明
图1船体坐标系示意图;
图2局部基准状态检测与姿态重构流程图;
图3分布式姿态基准状态监测网络结构图;
图4姿态基准构成图;
图5姿态基准在舰船上的安装布置示意图。
具体实施方式
本发明包括以下步骤:
第一步:构建分布式姿态基准网络。利用舰船的主惯导基准和分布式局部基准以及舰船上的网络通信系统,构建全舰姿态预测网络。
第二步:启动姿态基准设备。完成主惯导基准和局部基准设备启动,各姿态基准设备通过网络通信系统向信息处理计算机发送姿态测量结果。
第三步:姿态预测。首先,根据船舶力学中的二维水弹性理论建立全舰姿态预测模型。其次,利用建立的姿态预测模型和第二步姿态基准设备的姿态测量值,实时估计模型参数。最后,利用第二步姿态基准设备的姿态测量结果,对待监测局部基准的姿态进行预测。
第四步:局部基准状态检测。如果待监测局部基准无姿态信息输出,或者姿态测量值与第三步预测得到的姿态值之差超过设定的阈值,则判定为局部基准状态异常。
第五步:姿态信息重构。根据第四步的状态检测结果,如果局部基准出现状态异常,将第三步得到的姿态预测结果发送至出现故障的局部基准战位点,实现姿态信息重构。
为了便于描述,本发明建立了船体坐标系,如图1所示。定义为:坐标原点 o位于船尾,x轴沿艏艉线方向,指向船艏,y轴指向左舷,z轴垂直甲板面向上,x-y-z符合右手定则。定义绕ox轴的角运动为横摇角,绕oy轴的角运动为纵摇角,绕oz轴的运动为航向角。
本发明是一种舰载分布式局部基准状态检测与姿态信息重构方法,利用船舶二维水弹性理论中的模态分析法,建立姿态预测模型,利用分布式姿态基准的姿态测量信息完成模型参数估计和姿态预测,实现对待监测局部基准的状态检测,如果局部基准出现故障,给出姿态信息重构结果,流程如图2所示,包括以下几个步骤:
步骤一:构建分布式姿态基准网络。图3所示为分布式姿态基准状态监测网络结构图,包括姿态基准(1)(Attitude Measurement System,AMS)、网络交换机 (2)和信息处理计算机(3)。姿态基准(1)包括主惯导基准(11)和分布式局部基准(12),如图4所示。为了便于描述,根据姿态基准在ox轴上的安装距离xi,i=1,2,…,m 的不同,分别称为AMS1、AMS2、…、AMSm。各姿态基准在舰船上的分布如图 5所示。信息处理计算机用于姿态信息处理、状态检测和姿态信息重构。网络交换机用于姿态基准与信息处理计算机之间的通信。
步骤二:启动姿态基准设备。姿态基准完成启动,并将姿态测量结果通过网络发送至信息处理计算机。在t时刻,姿态基准AMS1、AMS2、…、AMSm的姿态测量结果为:
Θj(x1,t),Θj(x2,t),…,Θj(xm,t) (1)
式中,Θj(xi,t),i=1,2,…,m表示AMSi姿态测量值,下标j=x,y,z分别表示横摇角、纵摇角和航向角。
步骤三:姿态预测。信息处理计算机接收到姿态基准设备发送的姿态信息,利用全舰姿态预测模型实现待监测局部基准的姿态预测。姿态预测步骤如下:
(3.1)建立姿态预测模型。应用船舶二维水弹性理论中的模态分析法(见董艳秋著的《船舶波浪外荷和水弹性》第81-88页)可得船体在波浪中的姿态运动表示为角振动模态与主坐标叠加的形式,即:
其中,θrj(x)为角振动模态,由船体的结构参数确定,为常数;prj(t)为主坐标,表示船体在外界驱动力作用下的运动响应,是随时间变化的参数;r表示模态的阶数,r=1表示刚体的角振动模态,r=2,…,l表示挠曲变形的角振动模态;l表示总的角振动模态数;下标j=x,y,z分别表示横摇、纵摇和航向方向。
根据(2)式,AMS1、AMS2、…、AMSm处的姿态角可表示为:
将(3)式写为矩阵形式:
将上式简记为:
Yj=BjPj (5)
式中,
Yj=[Θj(x1,t) Θj(x2,t) … Θj(xm,t)]T
Pj=[p1j(t) p2j(t) … p lj(t)]T
式中,符号“T”表示矩阵转置。
根据(5)式,当姿态基准设备的数量大于等于待估计角振动模态的阶数,即 m≥l时,方程有解。利用最小二乘法可求得姿态运动的主坐标参数:
式中,“-1”表示矩阵求逆,Yj为AMS1、AMS2、…、AMSm测量的姿态信息,Bj为系数矩阵,由角振动模态参数确定。
求得主坐标参数,可得舰船上任意位置xk处的姿态预测值为:
Θj(xk,t)=θ1j(xk)p1j(t)+θ2j(xk)p2j(t)+…+θ lj(xk)p lj(t) (8)
式中,1≤k≤m。将(8)式写为矩阵形式:
Θj(xk,t)=Ej TPj (9)
式中,
Ej=[θ1j(xk) θ2j(xk) … θ lj(xk)]T (10)
将(7)式代入(9)式,可得:
令
(11)式可以转换为:
由于船体的角振动模态由船体结构参数确定,所以Aj为常量。令 Aj=[a1j a2j …amj]T,aij,i=1,2,…,m为矩阵Aj的元素。(13)式可以展开为:
Θj(xk,t)=a1jΘj(x1,t)+a2jΘj(x2,t)+…+amjΘj(xm,t) (14)
公式(14)即为姿态预测模型,可见,任意位置的姿态值可以表示为其它位置姿态测量值的线性组合。只要确定模型参数aij,i=1,2,…,m,就可以根据AMS1、 AMS2、…、AMSm测量得到的姿态输出,预测全舰任意点的姿态。由于考虑了挠曲变形的影响,相对传统的基于主惯导基准直接传递的方法,具有更高的姿态预测精度。
(3.2)模型参数估计。利用分布式姿态基准的姿态输出,对其中一套姿态基准的状态进行监测,并假设待监测的局部基准为AMSm(可任意选择)。在t0~t1时段,假设所有的姿态基准均正常工作。根据(14)式,利用t0~t1时段 AMSi,i=1,2,…,m-1的姿态输出Θj(xi,t),i=1,2,…,m-1和AMSm的姿态输出Θj(xm,t)构建参数估计矩阵:
式中,h为采样时间间隔。上式可简记为:
Zj=ΓjXj (16)
其中,
Zj=[Θj(xm,t0) Θj(xm,t0+h) … Θj(xm,t1)]T
Xj=[a1j a2j … a(m-1)j]T
利用最小二乘法可求得预测模型参数为:
(3.3)预测姿态。在t1~t2时段,信息处理计算机利用 AMSi,i=1,2,…,m-1的姿态输出和(3.2)步估计得到的参数,根据(14)式预测AMSm安装位置的姿态:
步骤四:局部基准状态检测。在t时刻,如果AMSm的无信息输出或姿态测量值Θj(xm,t)与(3.3)预测得到的姿态值之差超过设定的阈值Tj (如局部基准的姿态精度为3角分(1σ),按照3σ准则,Tj可设定为9角分),则判定AMSm出现故障。状态检测公式如下:
Claims (6)
1.一种舰载分布式局部基准状态检测与姿态信息重构方法,其特征在于,包括以下步骤:
1)构建分布式姿态基准网络,启动姿态基准设备,各姿态基准设备通过网络通信系统向信息处理计算机发送姿态测量结果;
2)根据船舶力学中的二维水弹性理论建立全舰姿态预测模型,利用所述全舰姿态预测模型和所述姿态基准设备的姿态测量结果,实时估计全舰姿态预测模型参数;利用所述姿态基准设备的姿态测量结果,对待监测局部基准的姿态进行预测;全舰姿态预测模型表达式为:
Θj(xk,t)=a1jΘj(x1,t)+a2jΘj(x2,t)+…+amjΘj(xm,t);其中,
Ej=[θ1j(xk) θ2j(xk) … θ lj(xk)]T,1≤k≤m,
θrj(x)为角振动模态,r表示模态的阶数,r=1表示刚体的角振动模态,r=2,...,l表示挠曲变形的角振动模态;Θj(xi,t),i=1,2,…,m表示第i个姿态基准系统的姿态测量结果,下标j=x,y,z分别表示横摇角、纵摇角和航向角,m为姿态基准系统的数量;
3)如果待监测局部基准系统无姿态信息输出,或者姿态测量值与第2)步预测得到的姿态值之差超过设定的阈值,则判定为局部基准状态异常;
4)根据第3)步的局部基准状态检测结果,如果局部基准出现状态异常,将第2)步得到的姿态预测结果发送至出现故障的局部基准战位点,实现姿态信息重构。
2.根据权利要求1所述的舰载分布式局部基准状态检测与姿态信息重构方法,其特征在于,步骤1)中,所述分布式姿态基准网络包括依次双向连接的姿态基准系统、网络交换机和信息处理计算机。
3.根据权利要求2所述的舰载分布式局部基准状态检测与姿态信息重构方法,其特征在于,步骤1)中,姿态测量结果表示为:
Θj(x1,t),Θj(x2,t),...,Θj(xm,t);
其中,Θj(xi,t),i=1,2,…,m表示第i个姿态基准系统的姿态测量结果,下标j=x,y,z分别表示横摇角、纵摇角和航向角,m为姿态基准系统的数量。
6.根据权利要求1所述的舰载分布式局部基准状态检测与姿态信息重构方法,其特征在于,步骤3)中,设定的阈值为9角分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810182388.XA CN108507570B (zh) | 2018-03-06 | 2018-03-06 | 舰载分布式局部基准状态检测与姿态信息重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810182388.XA CN108507570B (zh) | 2018-03-06 | 2018-03-06 | 舰载分布式局部基准状态检测与姿态信息重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108507570A CN108507570A (zh) | 2018-09-07 |
CN108507570B true CN108507570B (zh) | 2020-10-27 |
Family
ID=63376124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810182388.XA Active CN108507570B (zh) | 2018-03-06 | 2018-03-06 | 舰载分布式局部基准状态检测与姿态信息重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108507570B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110895149B (zh) * | 2019-12-04 | 2021-07-27 | 中国人民解放军国防科技大学 | 局部基准传递对准精度内场测试系统及测试方法 |
CN112762935B (zh) * | 2020-12-23 | 2022-08-12 | 赫星(厦门)电子有限公司 | 一种基于船体姿态监测的定位方法及系统 |
CN114786128B (zh) * | 2022-06-20 | 2022-09-13 | 中国船舶重工集团公司第七0七研究所 | 基于数据质量控制的船舶惯导传递对准自适应切换方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5640325A (en) * | 1988-09-13 | 1997-06-17 | Litton Systems, Inc. | Sensor array dynamic position and orientation determination system |
US7321305B2 (en) * | 2005-07-05 | 2008-01-22 | Pinc Solutions | Systems and methods for determining a location of an object |
CN102322873B (zh) * | 2011-08-23 | 2014-02-26 | 北京航空航天大学 | 一种分布式pos地面演示验证系统 |
CN103256928B (zh) * | 2013-04-28 | 2015-05-20 | 南京航空航天大学 | 一种分布式惯性导航系统及其姿态传递对准方法 |
CN106979781B (zh) * | 2017-04-12 | 2020-11-06 | 南京航空航天大学 | 基于分布式惯性网络的高精度传递对准方法 |
-
2018
- 2018-03-06 CN CN201810182388.XA patent/CN108507570B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108507570A (zh) | 2018-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108507570B (zh) | 舰载分布式局部基准状态检测与姿态信息重构方法 | |
Stokey et al. | Development of the REMUS 600 autonomous underwater vehicle | |
US20090043436A1 (en) | Automatic Vessel Position Holding Control Method and Controller | |
CN106979781B (zh) | 基于分布式惯性网络的高精度传递对准方法 | |
CN110083057B (zh) | 基于水翼运动姿态的pid控制方法 | |
CN110132308A (zh) | 一种基于姿态确定的usbl安装误差角标定方法 | |
Herrero et al. | Improving parameter estimation efficiency of a non linear manoeuvring model of an underwater vehicle based on model basin data | |
CN110333369B (zh) | 基于水面gps校正的uuv的dvl测速系统及自适应去噪方法 | |
CN112925332B (zh) | 无人艇与水下无人潜器联合系统的协同交会对接控制方法 | |
CN112241589B (zh) | 一种基于系统可观测度的多水下无人艇协同运动编队设计方法 | |
CN113156825A (zh) | 一种自适应反步滑模控制的舰载光电跟踪系统设计方法 | |
CN110908404B (zh) | 一种基于数据驱动的auv智能观测运动方法 | |
CN109975760B (zh) | 基于双矢量重构的基阵安装误差角高精度标定方法和装置 | |
CN103335650A (zh) | 一种基于惯性测量单元的坐标失配测量方法 | |
Sasaki et al. | Feasibility study for telexistence on a ship-measurement of delay time of satellite communication | |
KR102077621B1 (ko) | 잠수함의 심도 산출 장치 및 그 방법 | |
Velasco et al. | Measurements of Hydrodynamic Parameters and control of an underwater torpedo-shaped vehicle | |
CN111649744A (zh) | 一种基于动力学模型辅助的组合导航定位方法 | |
Souza et al. | Development of an asv trimaran for testing and implementation of auv navigation algorithms | |
Le et al. | Estimation of Ship Hydrodynamic Coefficients in Harbor Maneuvers and Its Applications | |
Allotta et al. | An IMU and USBL-aided buoy for underwater localization | |
Söder et al. | Motion-based monitoring of racking stresses in ro-ro ships | |
Aaltonen et al. | Implementation of a low-cost autonomous underwater vehicle using open source ROS components with consumer class sonar technologies | |
CN109425328A (zh) | 一种测量垂直剖面温度、盐度、声速的小型无人船设计 | |
JP7420358B1 (ja) | 位置保持システム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |