CN117516465B - 一种避免跳变的惯导***姿态角提取方法 - Google Patents
一种避免跳变的惯导***姿态角提取方法 Download PDFInfo
- Publication number
- CN117516465B CN117516465B CN202410010222.5A CN202410010222A CN117516465B CN 117516465 B CN117516465 B CN 117516465B CN 202410010222 A CN202410010222 A CN 202410010222A CN 117516465 B CN117516465 B CN 117516465B
- Authority
- CN
- China
- Prior art keywords
- angle
- steps
- following
- method comprises
- navigation system
- 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
- 238000000605 extraction Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 35
- 230000008569 process Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005295 random walk Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C1/00—Measuring angles
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
Abstract
本发明提供一种避免跳变的惯导***姿态角提取方法,根据惯导***解算姿态矩阵的元素得到俯仰角,判断俯仰角绝对值是否小于K°,如果是,更新航向角;如果否,不更新航向角,默认采用上周期提取的航向角;判断俯仰角是否大于零,如果是,通过横滚角与航向角之差提取横滚角;如果否,通过横滚角与航向角之和提取横滚角;最后对横滚角和航向角的取值范围进行规范化。本发明通过计算横滚角与航向角的和或差来提取横滚角,在全姿态范围内均不存在奇异性;当俯仰角处于±90°附近时,通过不更新航向角有效避免了姿态角的不确定跳变,本方法不改变欧拉角定义,不损失姿态角提取精度,有效避免了横滚角和航向角的无序跳变,满足用户和控制***的使用习惯。
Description
技术领域
本发明涉及测量测试技术领域,具体涉及一种避免跳变的惯导***姿态角提取方法。
背景技术
姿态角是载体坐标系与地理坐标系之间角位置关系最直观的描述方式,包含俯仰角、横滚角和航向角三个欧拉角。俯仰角表示载体纵轴与水平面之间的夹角,横滚角表示载体纵向对称面与纵轴所在铅垂平之间的夹角,航向角表示载体纵轴在水平面上的投影与当地子午线北向之间的夹角。欧拉角表述方式简单直观,容易理解,符合用户使用习惯,因此获得普遍应用并形成一系列经典的显示和控制方法。
从欧拉角定义来看,当俯仰角处于±90°时,飞行器的纵轴在水平面内的投影变成一个点,纵向铅垂面退化成纵向对称平面内的一条直线。所以从物理意义上看,尽管飞行器在空间中有确定的空间角位置, 但根本不存在有意义的航向角和横滚角, 这是姿态角定义上的奇异点, 反映了用姿态角描述刚体转动的局限性。
惯导***通常采用四元数或姿态矩阵微分方程进行姿态更新,导航算法本身没有奇异点,可以全姿态工作,姿态角提取仅仅是参数表达问题。在垂直条件下,航向角和横滚角都沿垂线方向定义,实际上相当于失去了一个自由度,航向角与横滚角的和或差可以确定,但无法提取唯一确定的航向角和横滚角信息。也就是说,虽然飞行器所处的空间角位置关系是唯一确定的,但可以有无数组欧拉角可以描述该位置,即欧拉角描述坐标转换关系存在多值性问题。
文献“惯性导航,科学出版社,2020,第三版,p260”给出了一种根据惯导***姿态矩阵提取姿态角的经典方法。方法根据姿态矩阵元素的反正弦和反正切函数求解俯仰角、横滚角和航向角,并通过真值表确定横滚角和航向角所在的象限。该方法能够在全范围内提取欧拉角信息,但当俯仰角处于±90°附近时,反正切函数内的数值出现分子分母均接近0的情况,横滚角和航向角的提取结果出现多值性,存在不确定跳变。横滚角和航向角的跳变本身是参数表达问题,不影响导航解算过程,但容易给用户或控制***造成不必要的困扰。
因此,需要一种可以解决全范围姿态角提取的角度跳变问题的方法。
发明内容
本发明是为了解决全范围姿态角提取的角度跳变问题,本发明提供一种避免跳变的惯导***姿态角提取方法,根据俯仰角计算结果判断提取流程,在俯仰角绝对值小于阈值时正常更新航向角,当俯仰角绝对值接近90度时默认采用上周期航向角,在已知航向角情况下,通过矩阵元素代数和计算横滚角与航向角的和或差来提取横滚角。该方法算法简洁,不存在原理性误差,可以有效解决垂直状态下航向角和横滚角的跳变问题。
本发明提供一种避免跳变的惯导***姿态角提取方法,包括以下步骤:
S1、根据惯导***解算姿态矩阵的元素得到俯仰角/>;
姿态矩阵为:
;
俯仰角为:/>;
其中,为姿态矩阵/>的第/>行第/>列元素,b为载体坐标系,n为导航坐标系;
S2、判断俯仰角绝对值是否小于K°,其中,K=90-k,k≤1,k为正数,如果是,更新航向角/>;如果否,不更新航向角/>,默认采用上周期提取的航向角/>;
S3、判断俯仰角是否大于零,如果是,进入步骤S4;如果否,进入步骤S5;
S4、通过横滚角与航向角之差提取横滚角/>,进入步骤S6;
S5、通过横滚角与航向角之和提取横滚角/>,进入步骤S6;
S6、对横滚角和航向角的取值范围进行规范化,得到规范化后横滚角和规范化后航向角/>;
规范化后航向角的值域为[0, />):/>;
规范化后横滚角的值域为[/>,/>):
;
一种避免跳变的惯导***姿态角提取方法完成。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,步骤S1中,的值域为[/>,/>)。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,步骤S2中,航向角为:
;
其中:atan2为标准库反正切函数,值域为[,/>)。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,步骤S2中,K为89.9992°。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,步骤S4中,横滚角为:
。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,横滚角与航向角之差为:
。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,对姿态矩阵的元素进行作差和求和处理得到:
。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,步骤S5中,横滚角为:
。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,横滚角与航向角之和为:
。
本发明所述的一种避免跳变的惯导***姿态角提取方法,作为优选方式,对姿态矩阵的元素进行作差和求和处理得到:
。
传统姿态提取方法中,由于载体的姿态角解算存在奇异点,惯导***在导航解算过程中,总是通过四元数或姿态矩阵方式进行姿态解算,然后再从姿态矩阵中提取姿态角供用户或控制***使用。
姿态矩阵的构成为:
;
传统方法直接根据姿态矩阵的元素计算俯仰角/>、横滚角/>和航向角/>:
; (1)
; (2)
; (3)
其中符号表示姿态矩阵/>的第/>行第/>列元素,atan2为标准库反正切函数,值域为[/>,/>)。
从数学公式上看,式(2)等价于:
;
区别仅在于atan2函数可以进一步确定值域的象限。可见,当俯仰角趋于±90°时,趋于0,在式(2)中横滚角/>计算式的分子/>趋于0,分母也趋于0。由于分子分母同时趋于0,分子分母的微小变化就会导致/>计算结果出现剧烈波动。同样,式(3)计算航向角/>也存在同样的情况。
从物理意义上看,在垂直条件下,航向角和横滚角都沿垂线方向定义,实际上相当于失去了一个自由度,航向角与横滚角的和或差可以确定,但无法提取唯一确定的航向角和横滚角信息。也就是说,虽然飞行器所处的空间角位置关系是唯一确定的,但可以有无数组欧拉角可以描述该位置,即欧拉角描述坐标转换关系存在多值性问题,因此正常提取过程会出现无序跳变。
本发明对姿态矩阵的元素进行加减并化简,得到:
; (4)
; (5)
当俯仰角为90°时,由式(4)可得:
; (6)
当俯仰角为-90°时,由式(5)可得:
; (7)
可见,当俯仰趋于±90°时,和/>仍然可以准确求解,这也与物理意义一致:在垂直条件下,航向角和横滚角都沿垂线方向定义,实际上相当于失去了一个自由度,航向角与横滚角的和或差可以确定,但无法提取唯一确定的航向角和横滚角信息。
根据式(4)可得:
;
根据式(5)可得:
;
因此,本发明使用式(1)计算俯仰角,不存在奇异性,但不再采用式(2)和(3)直接计算横滚角和航向角。
获得俯仰角后,当判断俯仰角不在±90°附近时,采用式(3)更新航向角没有奇异点,然后通过式(8)-(9)获得或/>后再求解横滚角/>。
当判断俯仰角在±90°附近时,式(3)计算航向角存在奇异,因此不再更新,相当于使用了前一周期解算的航向角,也就是在保证航向角不变的情况下计算横滚角,这样横滚角也不会发生跳变,从而实现了避免跳变的目的。
本发明具有以下优点:
本发明通过计算横滚角与航向角的和或差来提取横滚角,在全姿态范围内均不存在奇异性;当俯仰角处于±90°附近时,通过不更新航向角有效避免了姿态角的不确定跳变。该方法不改变欧拉角定义,不损失姿态角提取精度,且有效避免了横滚角和航向角的无序跳变,从而满足用户和控制***的使用习惯。
附图说明
图1为一种避免跳变的惯导***姿态角提取方法流程图;
图2a为现有技术通过导航解算姿态矩阵提取的俯仰角示意图;
图2b为现有技术通过导航解算姿态矩阵提取的横滚角示意图;
图2c为现有技术通过导航解算姿态矩阵提取的航向角示意图;
图3a为一种避免跳变的惯导***姿态角提取方法通过导航解算姿态矩阵提取的俯仰角示意图;
图3b为一种避免跳变的惯导***姿态角提取方法通过导航解算姿态矩阵提取的横滚角示意图;
图3c为一种避免跳变的惯导***姿态角提取方法通过导航解算姿态矩阵提取的航向角示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
实施例1
如图1所示,一种避免跳变的惯导***姿态角提取方法,包括以下步骤:
步骤1:参照图1,首先根据惯导***解算姿态矩阵的元素计算俯仰角/>。
姿态矩阵的构成为:
;
根据姿态矩阵的元素计算俯仰角/>:
;
其中,的值域为/>,符号/>表示姿态矩阵/>的第/>行第/>列元素。
步骤2:判断俯仰角绝对值<89.9992°时更新航向角/>;否则不更新,默认采用上周期提取的航向角。
根据姿态矩阵的元素更新航向角/>:
;
式中atan2为标准库反正切函数,值域为。
步骤3:判断俯仰角>0°时,通过矩阵元素代数和计算横滚角与航向角之差/>,并提取横滚角/>。
对姿态矩阵的元素进行处理得到:
;
横滚角与航向角之差为:
;
提取横滚角:
;
步骤4:判断俯仰角≤0°时,通过矩阵元素代数和计算横滚角与航向角之和/>,并提取横滚角/>。
对姿态矩阵的元素进行处理能够得到:
;
横滚角与航向角之和为:
;
提取横滚角:
:
步骤5:对横滚角和航向角的取值范围进行规范化。
计算规范化值域为的航向角/>:
;
计算规范化值域为的横滚角/>:
。
下面通过仿真试验对本发明的具体实施方式和实施效果做进一步说明。
仿真条件设置如下:设置载体初始俯仰角、横滚角和航向角为[89.999°; 5°;5°],以0.002°/s的俯仰角速率向上作抬头运动1s,最终俯仰角越过90°仍停留在89.999°。设置惯导***中陀螺零偏均方差为0.01°/h,随机游走噪声为0.001°/,导航解算采用陀螺角速率每10ms周期进行姿态矩阵更新。
图2a~2c为背景技术通过导航解算姿态矩阵提取的俯仰角、横滚角和航向角。在图2a中0.5s时俯仰角处于90°,受多值性和陀螺微小噪声影响,如图2b、2c所示,该时间段附近横滚角和航向角存在大幅度的不确定跳变。
图3a~3c为本发明技术通过导航解算姿态矩阵提取的俯仰角、横滚角和航向角。图3a中0.5s俯仰角处于90°时,图3b、3c中横滚角和航向角保持正常输出未发生跳变;在0.9s后,当图3a俯仰角越过90°且小于阈值89.9992°时,图3b、3c中横滚角和航向角出现正确的翻转。
实施例效果表明,本发明技术通过计算横滚角与航向角的和或差来提取横滚角,在全姿态范围内均不存在原理性误差;当俯仰角处于±90°附近时,通过不更新航向角有效避免了横滚角的不确定跳变。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种避免跳变的惯导***姿态角提取方法,其特征在于:包括以下步骤:
S1、根据惯导***解算姿态矩阵的元素得到俯仰角/>;
姿态矩阵为:
;
俯仰角为:/>;
其中,为姿态矩阵/>的第/>行第/>列元素,b为载体坐标系,n为导航坐标系;
S2、判断俯仰角绝对值是否小于K°,其中,K=90-k,k≤1,k为正数,如果是,更新航向角/>;如果否,不更新航向角/>,默认采用上周期提取的航向角/>;
S3、判断俯仰角是否大于零,如果是,进入步骤S4;如果否,进入步骤S5;
S4、通过横滚角与航向角之差提取横滚角/>,进入步骤S6;
S5、通过横滚角与航向角之和提取横滚角/>,进入步骤S6
S6、对横滚角和航向角的取值范围进行规范化,得到规范化后横滚角和规范化后航向角/>;
规范化后航向角的值域为/>:
;
规范化后横滚角的值域为/>:
;
一种避免跳变的惯导***姿态角提取方法完成。
2.根据权利要求1所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:步骤S1中,的值域为/>。
3.根据权利要求1所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:步骤S2中,航向角为:
;
其中:atan2为标准库反正切函数,值域为。
4.根据权利要求1所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:步骤S2中,K为89.9992°。
5.根据权利要求1所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:步骤S4中,横滚角为:
。
6.根据权利要求5所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:横滚角与航向角之差为:
。
7.根据权利要求6所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:对姿态矩阵的元素进行作差和求和处理得到:
。
8.根据权利要求1所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:步骤S5中,横滚角为:
。
9.根据权利要求8所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:横滚角与航向角之和为:
。
10.根据权利要求9所述的一种避免跳变的惯导***姿态角提取方法,其特征在于:对姿态矩阵的元素进行作差和求和处理得到:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410010222.5A CN117516465B (zh) | 2024-01-04 | 2024-01-04 | 一种避免跳变的惯导***姿态角提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410010222.5A CN117516465B (zh) | 2024-01-04 | 2024-01-04 | 一种避免跳变的惯导***姿态角提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117516465A CN117516465A (zh) | 2024-02-06 |
CN117516465B true CN117516465B (zh) | 2024-03-19 |
Family
ID=89749778
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410010222.5A Active CN117516465B (zh) | 2024-01-04 | 2024-01-04 | 一种避免跳变的惯导***姿态角提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117516465B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1851406A (zh) * | 2006-05-26 | 2006-10-25 | 南京航空航天大学 | 基于捷联惯性导航***的姿态估计和融合的方法 |
CN106342284B (zh) * | 2008-08-18 | 2011-11-23 | 西北工业大学 | 一种飞行载体姿态确定方法 |
CN104634348A (zh) * | 2015-03-12 | 2015-05-20 | 北京华力创通科技股份有限公司 | 组合导航中的姿态角计算方法 |
CN108534744A (zh) * | 2018-01-30 | 2018-09-14 | 歌尔科技有限公司 | 一种姿态角获取方法、装置和手柄 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5962771B2 (ja) * | 2012-12-12 | 2016-08-03 | 日産自動車株式会社 | 移動物***置姿勢角推定装置及び移動物***置姿勢角推定方法 |
CN111380501B (zh) * | 2020-03-25 | 2021-03-26 | 华中科技大学 | 基于双折射晶体视场效应的姿态角实时测量方法及装置 |
-
2024
- 2024-01-04 CN CN202410010222.5A patent/CN117516465B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1851406A (zh) * | 2006-05-26 | 2006-10-25 | 南京航空航天大学 | 基于捷联惯性导航***的姿态估计和融合的方法 |
CN106342284B (zh) * | 2008-08-18 | 2011-11-23 | 西北工业大学 | 一种飞行载体姿态确定方法 |
CN104634348A (zh) * | 2015-03-12 | 2015-05-20 | 北京华力创通科技股份有限公司 | 组合导航中的姿态角计算方法 |
CN108534744A (zh) * | 2018-01-30 | 2018-09-14 | 歌尔科技有限公司 | 一种姿态角获取方法、装置和手柄 |
Non-Patent Citations (2)
Title |
---|
基于空间机动平台的姿态算法研究;王延春;冯永新;潘成胜;;沈阳理工大学学报;20081215;第27卷(第06期);第30-34页 * |
陀螺测斜仪小角度井斜角测量的姿态提取方法;龙达峰;刘俊;张晓明;;传感技术学报;20130615;第26卷(第06期);第883-886页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117516465A (zh) | 2024-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104655131B (zh) | 基于istssrckf的惯性导航初始对准方法 | |
JP4615287B2 (ja) | 方位姿勢検出装置 | |
CN103900574B (zh) | 一种基于迭代容积卡尔曼滤波姿态估计方法 | |
CN109931955B (zh) | 基于状态相关李群滤波的捷联惯性导航***初始对准方法 | |
CN107063254B (zh) | 一种陀螺地磁组合的姿态解算方法 | |
CN105910606A (zh) | 一种基于角速度差值的方向修正方法 | |
CN108132053B (zh) | 一种行人轨迹构建方法、***及惯性测量装置 | |
CN110361683A (zh) | 基于双目标函数粒子群优化的磁力计校正方法 | |
CN108225370A (zh) | 一种运动姿态传感器的数据融合与解算方法 | |
WO2022160391A1 (zh) | 磁力计信息辅助的mems陀螺仪标定方法及标定*** | |
CN109631939B (zh) | 一种基于磁强计和加速度计的快速对准方法 | |
CN113587925B (zh) | 一种惯性导航***及其全姿态导航解算方法与装置 | |
CN108645404A (zh) | 一种小型多旋翼无人机姿态解算方法 | |
CN108592943B (zh) | 一种基于opreq方法的惯性系粗对准计算方法 | |
CN109470266A (zh) | 一种处理乘性噪声的星敏感器陀螺组合定姿方法 | |
CN110595503A (zh) | 基于李群最优估计的sins捷联惯性导航***晃动基座自对准方法 | |
CN113108790B (zh) | 机器人imu角度测量方法、装置、计算机设备及存储介质 | |
CN105004351A (zh) | 基于自适应upf的sins大方位失准角初始对准方法 | |
Chiang et al. | Constrained filtering method for attitude determination using GPS and gyro | |
US20170074689A1 (en) | Sensor Fusion Method for Determining Orientation of an Object | |
CN110737194B (zh) | 一种多mems传感器组合测姿方法 | |
CN117516465B (zh) | 一种避免跳变的惯导***姿态角提取方法 | |
CN113587926A (zh) | 一种航天器空间自主交会对接相对导航方法 | |
CN108801250A (zh) | 基于水下机器人的实时姿态获取方法及装置 | |
WO2016150312A1 (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 |