CN112987067B - 组合导航的洋流速度估计方法、装置和一种处理器 - Google Patents
组合导航的洋流速度估计方法、装置和一种处理器 Download PDFInfo
- Publication number
- CN112987067B CN112987067B CN202110204106.3A CN202110204106A CN112987067B CN 112987067 B CN112987067 B CN 112987067B CN 202110204106 A CN202110204106 A CN 202110204106A CN 112987067 B CN112987067 B CN 112987067B
- Authority
- CN
- China
- Prior art keywords
- ocean current
- auv
- coordinate system
- optimization function
- time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 claims abstract description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 19
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 230000009182 swimming Effects 0.000 claims 1
- 230000002349 favourable effect Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 28
- 230000008569 process Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/48—Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system
- G01S19/49—Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system whereby the further system is an inertial position system, e.g. loosely-coupled
-
- 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
- G01C21/165—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 combined with non-inertial navigation instruments
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
Abstract
本发明涉及一种组合导航的洋流速度估计方法、装置和一种处理器,所述一种组合导航的洋流速度估计方法,包括:建立最优化函数;所述最优化函数的参数包括:洋流速度;求解所述最优化函数,以估计出洋流速度。本发明所述的一种组合导航的洋流速度估计方法对于AUV工作在直线巡航状态下,也可完成对洋流速度的估计,且实时性高,本发明所述的洋流速度估计方法有利于提高导航精度,对于AUV在大深度海洋环境下的精准导航具有一定的实用价值。
Description
技术领域
本发明涉及惯性导航技术领域,具体涉及一种组合导航的洋流速度估计方法、装置和一种处理器。
背景技术
目前AUV水下导航多采用SINS/DVL组合导航形式,该模式下DVL有两种工作模式:探底模式(Bottom Tracking)和探水模式(Water Profiling),但是一般DVL的探底深度在百米量级,目前的SINS/DVL导航***中DVL的输出速度为相对速度,未能补偿洋流速度,造成大深度海洋环境中AUV导航精度非常低,有效估计出洋流并补偿,是提高AUV深海导航精度的关键。目前对于在深海环境下的估计洋流速度的方法主要有两种:第一种是将气象预报数据中的海流速度融合到SINS/DVL组合导航***中,第二种是采用高精度惯导,将洋流速度建模在SINS/DVL组合导航***中。
第一种估计洋流速度的方法一般需要从第三方机构购买,价格比较贵,一般在几万到几十万不等,并且数据存在较大延迟,一般延时在数小时到数天不等,这对于AUV的实时导航是不能接受的,数据精度差,一般分辨率在km量级,造成AUV的导航精度比较差;第二种估计洋流速度的方法理论上是可行的,但是经过实际水下导航验证的几乎没有,并且该方法需要考虑导航***的可观性,一般在转弯状态或者加减速状态下才能有效估计出洋流速度,但是AUV一般工作在直线巡航状态下,所以不具备普适性。现亟需提出一种可以在AUV直线巡航模式下估计出海流速度的方法,且要保证该方法具有较好的导航精度。
发明内容
有鉴于此,本发明的目的在于克服现有技术的不足,提供一种组合导航的洋流速度估计方法、装置和一种处理器。
为实现以上目的,本发明采用如下技术方案:一种组合导航的洋流速度估计方法,所述组合导航为SINS/DVL/GPS组合导航,包括:
建立最优化函数;所述最优化函数的参数包括:洋流速度;
求解所述最优化函数,以估计出洋流速度。
可选的,所述建立最优化函数,包括:
假设在i时刻,AUV在载体坐标系下的洋流速度为DVL探水模式下得到的DVL速度为/>DVL数据更新时间为/>SINS姿态角表示为(θp,θr,θy),(θp,θr,θy)分别表示俯仰角,滚转角和偏航角;
DVL坐标系和AUV坐标系重合,不存在转换关系,从运动学方面考虑,AUV在载体坐标系下的绝对速度满足以下公式:
其中,ENU代表东北天坐标系,AUV在东北天坐标系下解算出的位置和洋流速度呈线性关系,
时间T内,将GPS数据中的经纬度转换到东北天坐标系下的位置,表示为:
其中,i为i时刻,i时刻GPS的位置增量计算公式为:
最优化函数为:
条件二:GPS信号的精度大于预设精度值,且GPS信号的动态偏差值小于预设偏差值。
进一步的,当AUV每次上浮到水面停留时间短,以致每次只能采集到一个GPS观测数据时,所述最优化函数为:
可选的,所述求解所述最优化函数,包括:采用高斯牛顿方法求解所述最优化函数。
本发明还提供了一种组合导航的洋流速度估计装置,包括:
建模模块,用于建立最优化函数;所述最优化函数的参数包括:洋流速度;
求解模块,用于求解所述最优化函数,以估计出洋流速度。
本发明还提供了一种SINS/DVL/GPS组合导航***,包括:
如前面所述的一种组合导航的洋流速度估计装置。
本发明还提供了一种处理器,用于执行前面任一项所述的一种组合导航的洋流速度估计方法。
本发明采用以上技术方案,所述一种组合导航的洋流速度估计方法,所述组合导航为SINS/DVL/GPS组合导航,包括:建立最优化函数;所述最优化函数的参数包括:洋流速度;求解所述最优化函数,以估计出洋流速度。本发明所述的一种组合导航的洋流速度估计方法对于AUV工作在直线巡航状态下,也可完成对洋流速度的估计,且实时性高,本发明所述的洋流速度估计方法有利于提高导航精度,对于AUV在大深度海洋环境下的精准导航具有一定的实用价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明一种组合导航的洋流速度估计方法提供的流程示意图;
图2是本发明一种组合导航的洋流速度估计装置提供的结构示意图。
图中:1、建模模块;2、求解模块。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明的技术方案进行详细的描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施方式,都属于本发明所保护的范围。
图1是本发明一种组合导航的洋流速度估计方法提供的流程示意图。
如图1所示,本发明所述的一种组合导航的洋流速度估计方法,所述组合导航为SINS/DVL/GPS组合导航,包括:
S11:建立最优化函数;所述最优化函数的参数包括:洋流速度;
本发明所述的估计方法建立的假设如下:
假设在T时间内,洋流速度满足在各个方向上的大小是大致不变的,这是本模型构建的最为重要的假设。因为水下导航时没有GPS观测数据,无法实时估计洋流,为了保证实际的海洋环境满足这一假设,T在取值时不要过大。
进一步的,所述建立最优化函数,包括:
假设在i时刻,AUV在载体坐标系下的洋流速度为DVL探水模式下得到的DVL速度为/>DVL数据更新时间为/>SINS姿态角表示为(θp,θr,θy),(θp,θr,θy)分别表示俯仰角,滚转角和偏航角;
DVL坐标系和AUV坐标系重合,不存在转换关系,从运动学方面考虑,AUV在载体坐标系下的绝对速度满足以下公式:
其中,ENU代表东北天坐标系,AUV在东北天坐标系下解算出的位置和洋流速度呈线性关系,
时间T内,将GPS数据中的经纬度转换到东北天坐标系下的位置,表示为:
其中,i为i时刻,i时刻GPS的位置增量计算公式为:
最优化函数为:
其中,α为预设阈值,表示航向的直线度,表示AUV方位角的最大值,/>表示AUV方位角的最小值;考虑到实际不同方向上洋流速度是不一样的,对不同方向上的洋流做离散化处理,保证水下导航轨迹具有良好的线性度。
条件二:GPS信号的精度大于预设精度值,且GPS信号的动态偏差值小于预设偏差值。
具体的,可以通过计算GPS轨迹动态精度、水平精度因子、垂直精度因子、综合精度因子等因素进行综合评估。
在实际执行时,进行洋流速度估计过程需要AUV在水面定向游动一段时间,一般是几分钟,假设GPS的频率为1HZ,2分钟即可获取到120个点的SINS/DVL预测数据和GPS观测数据,其中,预测数据为SINS/DVL航位推算的输出位置,观测数据为GPS输出的位置。在满足上面提到的两个条件的前提下,获取到的120个点的SINS/DVL预测数据和GPS观测数据完全够用。
考虑到以下特殊情况,即假设AUV到水面只停留很短的时间,极端情况下每次上浮只有一个GPS观测数据,AUV就立即下潜。
在这种情况下,每次的GPS观测数据只有一个,但是下水下航行的SINS/DVL预测数据是很多的(与水下航行时间有关),此时的所述最优化函数为:
S12:求解所述最优化函数,以估计出洋流速度。
进一步的,所述求解所述最优化函数,包括:采用高斯牛顿方法求解所述最优化函数。
采用本发明所述的洋流速度估计方法进行洋流速度估计,组合导航***再结合估计出的洋流速度进行导航定位,表1是试验结果,实际的水下导航试验证明,加入上述洋流速度估计方法相比没有进行洋流速度估计可以大幅提高导航精度。
表1
本发明所述的一种组合导航的洋流速度估计方法对于AUV工作在直线巡航状态下,也可完成对洋流速度的估计,且实时性高,本发明所述的洋流速度估计方法有利于提高导航精度,对于AUV在大深度海洋环境下的精准导航具有一定的实用价值。
图2是本发明一种组合导航的洋流速度估计装置提供的结构示意图。
如图2所示,本发明所述的一种组合导航的洋流速度估计装置,包括:
建模模块1,用于建立最优化函数;所述最优化函数的参数包括:洋流速度;
求解模块2,用于求解所述最优化函数,以估计出洋流速度。
本发明所述一种组合导航的洋流速度估计装置的工作原理与上文所述的一种组合导航的洋流速度估计方法的工作原理相同,在此不再赘述。
本发明还提供了一种SINS/DVL/GPS组合导航***,包括:
如图2所述的一种组合导航的洋流速度估计装置。
此外,本发明还提供了一种处理器,用于执行图1所述的一种组合导航的洋流速度估计方法。
可以理解的是,上述各实施例中相同或相似部分可以相互参考,在一些实施例中未详细说明的内容可以参见其他实施例中相同或相似的内容。
需要说明的是,在本发明的描述中,术语“第一”、“第二”等仅用于描述目的,而不能理解为指示或暗示相对重要性。此外,在本发明的描述中,除非另有说明,“多个”的含义是指至少两个。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的实施例所属技术领域的技术人员所理解。
应当理解,本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行***执行的软件或固件来实现。例如,如果用硬件来实现,和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。
此外,在本发明各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器,磁盘或光盘等。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (5)
1.一种组合导航的洋流速度估计方法,所述组合导航为SINS/DVL/GPS组合导航,其特征在于,包括:
建立最优化函数;所述最优化函数的参数包括:洋流速度;
求解所述最优化函数,以估计出洋流速度;
其中,所述建立最优化函数,包括:
假设在i时刻,AUV在载体坐标系下的洋流速度为DVL探水模式下得到的DVL速度为/>DVL数据更新时间为/>SINS姿态角表示为(θp,θr,θy),(θp,θr,θy)分别表示俯仰角,滚转角和偏航角;
DVL坐标系和AUV坐标系重合,不存在转换关系,从运动学方面考虑,AUV在载体坐标系下的绝对速度满足以下公式:
其中,ENU代表东北天坐标系,AUV在东北天坐标系下解算出的位置和洋流速度呈线性关系,
时间T内,将GPS数据中的经纬度转换到东北天坐标系下的位置,表示为:
其中,i为i时刻,i时刻GPS的位置增量计算公式为:
最优化函数为:
条件二:GPS信号的精度大于预设精度值,且GPS信号的动态偏差值小于预设偏差值;
其中,当AUV每次上浮到水面停留时间短,以致每次只能采集到一个GPS观测数据时,所述最优化函数为:
2.根据权利要求1所述的洋流速度估计方法,其特征在于,所述求解所述最优化函数,包括:采用高斯牛顿方法求解所述最优化函数。
3.一种组合导航的洋流速度估计装置,其特征在于,执行如权利要求1所述的方法,所述装置包括:
建模模块,用于建立最优化函数;所述最优化函数的参数包括:洋流速度;
求解模块,用于求解所述最优化函数,以估计出洋流速度;
其中,所述建立最优化函数,包括:
假设在i时刻,AUV在载体坐标系下的洋流速度为DVL探水模式下得到的DVL速度为/>DVL数据更新时间为/>SINS姿态角表示为(θp,θr,θy),(θp,θr,θy)分别表示俯仰角,滚转角和偏航角;
DVL坐标系和AUV坐标系重合,不存在转换关系,从运动学方面考虑,AUV在载体坐标系下的绝对速度满足以下公式:
其中,ENU代表东北天坐标系,AUV在东北天坐标系下解算出的位置和洋流速度呈线性关系,
时间T内,将GPS数据中的经纬度转换到东北天坐标系下的位置,表示为:
其中,i为i时刻,i时刻GPS的位置增量计算公式为:
最优化函数为:
4.一种SINS/DVL/GPS组合导航***,其特征在于,包括:
如权利要求3所述的一种组合导航的洋流速度估计装置。
5.一种处理器,其特征在于,用于执行权利要求1至2任一项所述的一种组合导航的洋流速度估计方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110204106.3A CN112987067B (zh) | 2021-02-24 | 2021-02-24 | 组合导航的洋流速度估计方法、装置和一种处理器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110204106.3A CN112987067B (zh) | 2021-02-24 | 2021-02-24 | 组合导航的洋流速度估计方法、装置和一种处理器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112987067A CN112987067A (zh) | 2021-06-18 |
CN112987067B true CN112987067B (zh) | 2023-06-27 |
Family
ID=76349921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110204106.3A Active CN112987067B (zh) | 2021-02-24 | 2021-02-24 | 组合导航的洋流速度估计方法、装置和一种处理器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112987067B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113848959B (zh) * | 2021-09-13 | 2024-02-02 | 西北工业大学 | 一种周期更新估计洋流速度的水下滑翔机航向修正方法 |
CN114323001A (zh) * | 2021-12-15 | 2022-04-12 | 中国飞行试验研究院 | 一种基于洋流自主修正的多普勒测速精度试飞方法 |
CN114236173B (zh) * | 2021-12-16 | 2024-04-02 | 中国人民解放军国防科技大学 | 基于sins与gnss组合及dvl的洋流速度估计方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103017755B (zh) * | 2012-11-26 | 2016-01-27 | 哈尔滨工程大学 | 一种水下导航姿态测量方法 |
CN109253726B (zh) * | 2018-06-22 | 2020-05-05 | 东南大学 | 一种水下滑翔器导航定位***及上浮精度校正方法 |
CN110567454B (zh) * | 2019-09-08 | 2023-03-24 | 东南大学 | 一种复杂环境下sins/dvl紧组合导航方法 |
CN111829512B (zh) * | 2020-06-08 | 2024-04-09 | 中国航天空气动力技术研究院 | 一种基于多传感器数据融合的auv导航定位方法及*** |
CN111829511A (zh) * | 2020-06-08 | 2020-10-27 | 中国航天空气动力技术研究院 | 一种基于m估计的auv组合导航方法及*** |
-
2021
- 2021-02-24 CN CN202110204106.3A patent/CN112987067B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112987067A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112987067B (zh) | 组合导航的洋流速度估计方法、装置和一种处理器 | |
Specht et al. | Application of an autonomous/unmanned survey vessel (ASV/USV) in bathymetric measurements | |
CN105823480B (zh) | 基于单信标的水下移动目标定位算法 | |
CN109324330A (zh) | 基于混合无导数扩展卡尔曼滤波的usbl/sins紧组合导航定位方法 | |
CN106679662B (zh) | 一种基于tma技术的水下机器人单信标组合导航方法 | |
Turnbull et al. | Operational iceberg drift forecasting in Northwest Greenland | |
CN102323586B (zh) | 一种基于海流剖面的uuv辅助导航方法 | |
CN102213594A (zh) | 一种无人潜航器海流观测数据融合方法 | |
CN102749622A (zh) | 基于多波束测深的声速剖面及海底地形的联合反演方法 | |
CN107990891A (zh) | 基于长基线和信标在线标定的水下机器人组合导航方法 | |
CN109579850B (zh) | 基于对水速度辅助惯导的深水智能导航方法 | |
CN108871379B (zh) | 一种dvl测速误差在线标定方法 | |
CN110608744A (zh) | 一种具有动态避障功能的水质采样无人艇路径规划方法 | |
RU2467914C1 (ru) | Способ контроля мореходности судна и устройство для его осуществления | |
CN113686299A (zh) | 一种海上动态目标定位与移速预测方法 | |
JP6567665B2 (ja) | 船舶の航路(パス)の全地点における各々のドリフト(浮動)ベクトルを推定する方法 | |
CN111081071A (zh) | 一种船舶路径数据分析***及方法 | |
CN115792985A (zh) | 一种车辆定位方法、装置、电子设备、存储介质及车辆 | |
Popielarczyk | RTK water level determination in precise inland bathymetric measurements | |
US20180031376A1 (en) | Systems and methods for computing and displaying navigation indicators based on historic and real-time vessel data | |
CN113501114B (zh) | 一种基于深潜器无动力下潜实时信息的深海海流计算方法 | |
CN110136487A (zh) | 地图导航防触礁安全助航技术应用模型 | |
Matte et al. | A robust estimation method for correcting dynamic draft error in PPK GPS elevation using ADCP tilt data | |
CN117113798A (zh) | 漂浮物的漂移速度预测方法、装置、设备及存储介质 | |
CN111080788B (zh) | 一种海底地形绘制方法及装置 |
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 |