CN117109639B - 半球谐振陀螺的温漂误差检测方法及系统 - Google Patents
半球谐振陀螺的温漂误差检测方法及系统 Download PDFInfo
- Publication number
- CN117109639B CN117109639B CN202311380657.0A CN202311380657A CN117109639B CN 117109639 B CN117109639 B CN 117109639B CN 202311380657 A CN202311380657 A CN 202311380657A CN 117109639 B CN117109639 B CN 117109639B
- Authority
- CN
- China
- Prior art keywords
- damping
- resonance frequency
- frequency coefficient
- circumferential
- angle
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 17
- 238000013016 damping Methods 0.000 claims abstract description 220
- 238000004364 calculation method Methods 0.000 claims abstract description 105
- 108010076504 Protein Sorting Signals Proteins 0.000 claims abstract description 93
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000012360 testing method Methods 0.000 claims description 92
- 230000008569 process Effects 0.000 claims description 18
- 238000013101 initial test Methods 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
-
- 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)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Gyroscopes (AREA)
Abstract
本发明涉及惯性导航技术领域,公开了一种半球谐振陀螺的温漂误差检测方法及系统,用于提高对温漂误差检测及补偿的准确率。方法包括:通过谐振频率信号序列对半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,对谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;对周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,对阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;对第一谐振频率系数集合以及第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过目标温度数据对半球谐振陀螺进行温漂补偿。
Description
技术领域
本发明涉及惯性导航技术领域,尤其涉及一种半球谐振陀螺的温漂误差检测方法及系统。
背景技术
半球谐振陀螺是一种具有惯导级性能的高精度固体振动陀螺, 其基于哥氏振动原理,利用半球唇壳的径向振动驻波进动效应感测基座旋转,其不含机械转子,具备测量精度高,工作寿命长,体积小,抗冲击、过载、辐射能力强,稳定性和可靠性高等优点,在惯导领域受到越来越多的关注和应用,对国家的工业智能化和武器装备现代化发展具有重要意义。
但是,半球谐振陀螺经常在复杂的温度环境下工作,包括温度变化快速、幅度较大的情况。这使得准确地捕捉和模型化温度引起的漂移变化变得更加困难。半球谐振陀螺的温度响应通常是非线性的,随着温度的变化,其性能和特性可能会出现复杂的变化。这就需要在漂移补偿算法中考虑非线性因素,增加了补偿算法的复杂性。在实际应用中,半球谐振陀螺的各个部分可能受到不同的温度影响,导致整个装置的温度分布不均匀。这可能会导致在不同部分出现不同程度的温度漂移,增加了补偿的难度。
发明内容
有鉴于此,本发明实施例提供了一种半球谐振陀螺的温漂误差检测方法及系统,用于提高对温漂误差检测及补偿的准确率。
本发明提供了一种半球谐振陀螺的温漂误差检测方法,包括:获取预设的测试参数集合,其中,所述测试参数集合包括:测试周期以及初始测试温度;通过所述测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集所述半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;基于所述振型角控制信号序列,通过所述谐振频率信号序列对所述半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对所述谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;对所述周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;对所述第一谐振频率系数集合以及所述第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过所述目标温度数据对所述半球谐振陀螺进行温漂补偿。
在本发明中,所述通过所述测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集所述半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列步骤,包括:对所述测试周期进行周期分割,得到稳态周期以及力反馈控制周期;通过所述稳态周期以及所述力反馈控制周期对所述半球谐振陀螺进行温度测试;在所述力反馈控制周期时,对所述半球谐振陀螺进行频率信号采集,得到对应的谐振频率信号序列,同时,对所述半球谐振陀螺进行振型角控制信号采集,得到对应的振型角控制信号序列。
在本发明中,所述基于所述振型角控制信号序列,通过所述谐振频率信号序列对所述半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对所述谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合步骤,包括:
通过预置的第一阻尼角计算公式,对所述振型角控制信号序列进行初始阻尼角计算,得到初始阻尼角数据,其中,所述第一阻尼角计算公式如下所示:
其中,为初始阻尼角,/>为第一振型角控制信号,/>为第二振型角控制信号,/>为第三振型角控制信号;
通过预置的第一周向阻尼不均计算公式,对所述振型角控制信号序列进行初始周向阻尼不均计算,得到初始周向阻尼不均数据,其中,所述第一周向阻尼不均计算公式如下所示:
其中,为初始周向阻尼不均数据;
对所述初始阻尼角数据以及所述初始周向阻尼不均数据进行稳态谐振频率分析,得到稳态谐振频率;
基于所述稳态谐振频率,对所述谐振频率信号序列进行周向阻尼不均以及阻尼角计算,得到所述周向阻尼不均数据集合以及所述阻尼角数据集合。
在本发明中,所述基于所述稳态谐振频率,对所述谐振频率信号序列进行周向阻尼不均以及阻尼角计算,得到所述周向阻尼不均数据集合以及所述阻尼角数据集合步骤,包括:
基于所述稳态谐振频率,通过预置的第二阻尼角计算公式,对所述谐振频率信号序列进行阻尼角转换,得到对应的阻尼角数据集合,其中,所述第二阻尼角计算公式如下所示:
其中,为阻尼角数据,/>为初始阻尼角谐振频率系数,/>为第一阻尼角谐振频率系数,/>为第二阻尼角谐振频率系数,/>为第三阻尼角谐振频率系数;/>为稳态谐振频率,/>为谐振频率信号;
基于所述稳态谐振频率,通过预置的第二周向阻尼不均计算公式,对所述谐振频率信号序列进行周向阻尼不均计算,得到对应的周向阻尼不均数据集合,其中,所述第二周向阻尼不均计算公式如下所示:
其中,为周向阻尼不均数据,/>为初始周向阻尼不均谐振频率系数,/>为第一周向阻尼不均谐振频率系数,/>为第二周向阻尼不均谐振频率系数,/>为第三周向阻尼不均谐振频率系数。
在本发明中,所述对所述周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合步骤,包括:
对所述周向阻尼不均数据集合进行序列转换,得到对应的周向阻尼不均序列;
对所述阻尼角数据集合进行序列转换,得到对应的阻尼角序列;
通过预置的第一谐振频率系数计算公式对所述周向阻尼不均序列进行第一谐振频率系数计算,得到第一谐振频率系数集合,其中,所述第一谐振频率系数计算公式如下所示:
其中,为第一谐振频率系数集合,/>为谐振频率信号序列,/>为周向阻尼不均序列;
对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合。
在本发明中,所述对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合步骤,包括:
通过预置的第二谐振频率系数计算公式对所述阻尼角序列进行第二谐振频率系数计算,得到第二谐振频率系数集合,其中,所述第二谐振频率系数计算公式如下所示:
其中,为第二谐振频率系数集合,/>阻尼角序列。
本发明还提供了一种半球谐振陀螺的温漂误差检测系统,包括:
获取模块,用于获取预设的测试参数集合,其中,所述测试参数集合包括:测试周期以及初始测试温度;
测试模块,用于通过所述测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集所述半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;
第一计算模块,用于基于所述振型角控制信号序列,通过所述谐振频率信号序列对所述半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对所述谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;
第二计算模块,用于对所述周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;
补偿模块,用于对所述第一谐振频率系数集合以及所述第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过所述目标温度数据对所述半球谐振陀螺进行温漂补偿。
本发明提供的技术方案中,获取预设的测试参数集合,其中,测试参数集合包括:测试周期以及初始测试温度;通过测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;基于振型角控制信号序列,通过谐振频率信号序列对半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;对周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;对所述第一谐振频率系数集合以及所述第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过所述目标温度数据对所述半球谐振陀螺进行温漂补偿。在本申请方案中,采用力反馈模式,施加控制力驱动陀螺振型进动到空间关系固定的三个位置,无需如旋转机构等其他设备进行控制,提高了实验精度;力反馈控制下,改变振型角方位后预留了缓冲时间,并以振型角停留在固定位置后一段时间内控制信号平均值作为输出,减小了振型角控制回路和噪声带来的干扰,通过获取预设的测试参数集合,包括测试周期和初始测试温度,可以对半球谐振陀螺进行温度测试。在测试过程中,采集半球谐振陀螺的谐振频率信号序列以及振型角控制信号序列。通过阻尼角数据集合,可以计算第二谐振频率系数,并根据第一谐振频率系数集合和第二谐振频率系数集合,对半球谐振陀螺进行温漂补偿,能够减小温度变化对谐振陀螺性能的影响,提高其稳定性和准确性。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中半球谐振陀螺的温漂误差检测方法的流程图。
图2为本发明实施例中通过测试参数集合对预置的半球谐振陀螺进行温度测试的流程图。
图3为本发明实施例中半球谐振陀螺的温漂误差检测系统的示意图。
附图标记:
301、获取模块;302、测试模块;303、第一计算模块;304、第二计算模块;305、补偿模块。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
为便于理解,下面对本发明实施例的具体流程进行描述,请参阅图1,图1是本发明实施例的一种半球谐振陀螺的温漂误差检测方法的流程图,如图1所示,包括以下步骤:
S101、获取预设的测试参数集合,其中,测试参数集合包括:测试周期以及初始测试温度;
S102、通过测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;
S103、基于振型角控制信号序列,通过谐振频率信号序列对半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;
S104、对周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;
S105、对第一谐振频率系数集合以及第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过目标温度数据对半球谐振陀螺进行温漂补偿。
需要说明的是,按温度谱进行温度试验,陀螺采用力反馈的工作模式,一个测试周期可以设置为900s,其中,初始环境温度为20℃,每经过一个周期提高0.5℃,试验时间t=15小时,录数频率r=100Hz(每秒录取100个数),试验时陀螺敏感轴指向地速为零点,此时陀螺输出振型角控制信号仅含常值漂移,前600s用于等待半球谐振陀螺达到稳态,后300s中另分3个100s的小周期,每个小周期的前20s通过力反馈控制将振型角分别控制到在空间中相对于电极轴夹角0°,22.5°和45°的位置,后80s通过信号输出模块录取陀螺的振型角控制信号序列,谐振频率信号序列。
由于谐振陀螺的振型角和谐振频率之间存在一定的关系,通过对振型角控制信号序列进行分析,可以计算出半球谐振陀螺在不同情况下的周向阻尼不均数据集合。同时,结合转换后的频率序列,还可以计算出阻尼角数据集合,这些数据在后续的温漂补偿中会用到。周向阻尼不均数据集合可以用于计算第一谐振频率系数,这是通过分析不同温度下的谐振频率和阻尼数据的关系得出的。同时,阻尼角数据集合可以用于计算第二谐振频率系数,这个系数同样会随温度变化而变化。温漂是指随着温度变化而引起的仪器性能变化。通过之前计算得到的第一和第二谐振频率系数集合,可以建立一个温漂补偿模型。这个模型可以在实际使用中,根据当前的谐振频率系数和温度变化,对半球谐振陀螺的性能进行补偿,以减小温度变化对谐振陀螺性能的影响。
需要说明的是,确定需要在半球谐振陀螺中实现的目标温度数据。使用温度-系数关系模型将目标温度数据转化为相应的第一谐振频率系数和第二谐振频率系数。这涉及将目标温度数据代入建立的模型中,以获得目标频率系数的估算值。将从目标温度数据匹配得到的第一谐振频率系数和第二谐振频率系数应用于半球谐振陀螺的控制系统。这意味着在不同温度下,根据实际工作温度与目标温度之间的差异,即时调整半球谐振陀螺的操作参数,以抵消温度引起的频率漂移。建立一个闭环控制系统,该系统使用实时监测的温度数据并将其与目标温度数据匹配的频率系数进行比较。在每个时间步中,控制系统将根据这些数据实时调整半球谐振陀螺的操作,以确保其在不同温度条件下保持稳定的谐振频率。
例如,假设有一个半球谐振陀螺用于精密导航,操作温度范围在-20℃至+40℃。首先采集了在不同温度条件下的第一谐振频率系数和第二谐振频率系数数据,以及相应的温度数据。然后,建立了一个温度-系数关系模型,可能是一个线性回归模型。接下来,确定了在各个温度点下所需的目标温度数据,以满足导航性能的要求。对于特定温度点,可以使用模型将目标温度数据匹配到对应的第一谐振频率系数和第二谐振频率系数。如果实际温度与目标温度不匹配,控制系统将根据这些差异即时调整谐振陀螺的操作参数,以确保它在不同温度下保持所需的性能。
通过执行上述步骤,获取预设的测试参数集合,其中,测试参数集合包括:测试周期以及初始测试温度;通过测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;基于振型角控制信号序列,通过谐振频率信号序列对半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;对周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;对第一谐振频率系数集合以及第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过目标温度数据对半球谐振陀螺进行温漂补偿。在本申请方案中,采用力反馈模式,施加控制力驱动陀螺振型进动到空间关系固定的三个位置,无需如旋转机构等其他设备进行控制,提高了实验精度;力反馈控制下,改变振型角方位后预留了缓冲时间,并以振型角停留在固定位置后一段时间内控制信号平均值作为输出,减小了振型角控制回路和噪声带来的干扰,通过获取预设的测试参数集合,包括测试周期和初始测试温度,可以对半球谐振陀螺进行温度测试。在测试过程中,采集半球谐振陀螺的谐振频率信号序列以及振型角控制信号序列。通过阻尼角数据集合,可以计算第二谐振频率系数,并根据第一谐振频率系数集合和第二谐振频率系数集合,对半球谐振陀螺进行温漂补偿,能够减小温度变化对谐振陀螺性能的影响,提高其稳定性和准确性。
在一具体实施例中,如图2所示,执行步骤S102的过程可以具体包括如下步骤:
S201、对测试周期进行周期分割,得到稳态周期以及力反馈控制周期;
S202、通过稳态周期以及力反馈控制周期对半球谐振陀螺进行温度测试;
S203、在力反馈控制周期时,对半球谐振陀螺进行频率信号采集,得到对应的谐振频率信号序列,同时,对半球谐振陀螺进行振型角控制信号采集,得到对应的振型角控制信号序列。
具体的,每个周期内包含不同阶段,包括稳态周期和力反馈控制周期。初始环境温度为20℃,在前600秒内,半球谐振陀螺逐渐达到稳态。这个阶段用于让陀螺达到所需的工作状态。在接下来的300秒内,将被分为3个小周期,每个小周期为100秒。在每个小周期内,前20秒将通过力反馈控制,将半球谐振陀螺的振型角控制到特定位置(0°、22.5°和45°)。剩余的80秒用于采集振型角控制信号和谐振频率信号。在稳态周期内,半球谐振陀螺将在初始环境温度下达到稳定状态,以确保温度对其性能的影响已经稳定。
在力反馈控制周期内,将逐渐提高温度,每个小周期提高0.5℃,并在每个小周期的最后80秒内进行数据采集。在每个小周期的后80秒内,进行信号采集:对于频率信号采集:使用录数频率r=100 Hz,意味着每秒采集100个数据点,将这些数据点作为谐振频率信号序列,在每个小周期的前20秒通过力反馈控制将半球谐振陀螺的振型角控制到特定位置。因此,在这20秒内,记录陀螺的振型角控制信号序列。
在一具体实施例中,执行步骤S103的过程可以具体包括如下步骤:
(1)通过预置的第一阻尼角计算公式,对振型角控制信号序列进行初始阻尼角计算,得到初始阻尼角数据,其中,第一阻尼角计算公式如下所示:
其中,为初始阻尼角,/>为第一振型角控制信号,/>为第二振型角控制信号,/>为第三振型角控制信号;
(2)通过预置的第一周向阻尼不均计算公式,对振型角控制信号序列进行初始周向阻尼不均计算,得到初始周向阻尼不均数据,其中,第一周向阻尼不均计算公式如下所示:
其中,为初始周向阻尼不均数据;
(3)对初始阻尼角数据以及初始周向阻尼不均数据进行稳态谐振频率分析,得到稳态谐振频率;
(4)基于稳态谐振频率,对谐振频率信号序列进行周向阻尼不均以及阻尼角计算,得到周向阻尼不均数据集合以及阻尼角数据集合。
需要说明的是,在对初始阻尼角数据以及初始周向阻尼不均数据进行稳态谐振频率分析,得到稳态谐振频率时,为了进行频率分析,建立半球谐振陀螺的数学模型。使用欧拉方程或其他相关的动力学方程来描述陀螺的运动。这个模型应该包括角速度、角度、阻尼等参数。在模型中设置初始条件,包括初始阻尼角数据和初始周向阻尼不均数据。这些数据将影响谐振频率的计算。稳态谐振频率是在系统达到稳定状态时的振动频率。在进行频率分析时可以使用数值求解方法,如数值积分或数值模拟,来模拟半球谐振陀螺的运动,并记录时间序列数据。获取时间序列数据后,可以使用傅立叶变换来将数据从时域转换为频域,最终确定稳态谐振频率。
在一具体实施例中,执行基于稳态谐振频率,对谐振频率信号序列进行周向阻尼不均以及阻尼角计算,得到周向阻尼不均数据集合以及阻尼角数据集合步骤的过程可以具体包括如下步骤:
(1)基于稳态谐振频率,通过预置的第二阻尼角计算公式,对谐振频率信号序列进行阻尼角转换,得到对应的阻尼角数据集合,其中,第二阻尼角计算公式如下所示:
其中,为阻尼角数据,/>为初始阻尼角谐振频率系数,/>为第一阻尼角谐振频率系数,/>为第二阻尼角谐振频率系数,/>为第三阻尼角谐振频率系数;/>为稳态谐振频率,/>为谐振频率信号;
(2)基于稳态谐振频率,通过预置的第二周向阻尼不均计算公式,对谐振频率信号序列进行周向阻尼不均计算,得到对应的周向阻尼不均数据集合,其中,第二周向阻尼不均计算公式如下所示:
其中,为周向阻尼不均数据,/>为初始周向阻尼不均谐振频率系数,/>为第一周向阻尼不均谐振频率系数,/>为第二周向阻尼不均谐振频率系数,/>为第三周向阻尼不均谐振频率系数。
在一具体实施例中,执行步骤S104步骤的过程可以具体包括如下步骤:
(1)对周向阻尼不均数据集合进行序列转换,得到对应的周向阻尼不均序列;
(2)对阻尼角数据集合进行序列转换,得到对应的阻尼角序列;
(3)通过预置的第一谐振频率系数计算公式对周向阻尼不均序列进行第一谐振频率系数计算,得到第一谐振频率系数集合,其中,第一谐振频率系数计算公式如下所示:
其中,为第一谐振频率系数集合,/>为谐振频率信号序列,/>为周向阻尼不均序列;
(4)对阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合。
在一具体实施例中,执行对阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合步骤的过程可以具体包括如下步骤:
(1)通过预置的第二谐振频率系数计算公式对阻尼角序列进行第二谐振频率系数计算,得到第二谐振频率系数集合,其中,第二谐振频率系数计算公式如下所示:
其中,为第二谐振频率系数集合,/>阻尼角序列。
本发明实施例还提供了一种半球谐振陀螺的温漂误差检测系统,如图3所示,该一种半球谐振陀螺的温漂误差检测系统具体包括:
获取模块301,用于获取预设的测试参数集合,其中,所述测试参数集合包括:测试周期以及初始测试温度;
测试模块302,用于通过所述测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集所述半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;
第一计算模块303,用于基于所述振型角控制信号序列,通过所述谐振频率信号序列对所述半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对所述谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;
第二计算模块304,用于对所述周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;
补偿模块305,用于对所述第一谐振频率系数集合以及所述第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过所述目标温度数据对所述半球谐振陀螺进行温漂补偿。
通过上述各个模块的协同工作,获取预设的测试参数集合,其中,测试参数集合包括:测试周期以及初始测试温度;通过测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;基于振型角控制信号序列,通过谐振频率信号序列对半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;对周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;对所述第一谐振频率系数集合以及所述第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过所述目标温度数据对所述半球谐振陀螺进行温漂补偿。在本申请方案中,采用力反馈模式,施加控制力驱动陀螺振型进动到空间关系固定的三个位置,无需如旋转机构等其他设备进行控制,提高了实验精度;力反馈控制下,改变振型角方位后预留了缓冲时间,并以振型角停留在固定位置后一段时间内控制信号平均值作为输出,减小了振型角控制回路和噪声带来的干扰,通过获取预设的测试参数集合,包括测试周期和初始测试温度,可以对半球谐振陀螺进行温度测试。在测试过程中,采集半球谐振陀螺的谐振频率信号序列以及振型角控制信号序列。通过阻尼角数据集合,可以计算第二谐振频率系数,并根据第一谐振频率系数集合和第二谐振频率系数集合,对半球谐振陀螺进行温漂补偿,能够减小温度变化对谐振陀螺性能的影响,提高其稳定性和准确性。
以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种半球谐振陀螺的温漂误差检测方法,其特征在于,包括:
获取预设的测试参数集合,其中,所述测试参数集合包括:测试周期以及初始测试温度;
通过所述测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集所述半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;
基于所述振型角控制信号序列,通过所述谐振频率信号序列对所述半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对所述谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;
对所述周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;
对所述第一谐振频率系数集合以及所述第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过所述目标温度数据对所述半球谐振陀螺进行温漂补偿。
2.根据权利要求1所述的半球谐振陀螺的温漂误差检测方法,其特征在于,所述通过所述测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集所述半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列步骤,包括:
对所述测试周期进行周期分割,得到稳态周期以及力反馈控制周期;
通过所述稳态周期以及所述力反馈控制周期对所述半球谐振陀螺进行温度测试;
在所述力反馈控制周期时,对所述半球谐振陀螺进行频率信号采集,得到对应的谐振频率信号序列,同时,对所述半球谐振陀螺进行振型角控制信号采集,得到对应的振型角控制信号序列。
3.根据权利要求1所述的半球谐振陀螺的温漂误差检测方法,其特征在于,所述基于所述振型角控制信号序列,通过所述谐振频率信号序列对所述半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对所述谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合步骤,包括:
通过预置的第一阻尼角计算公式,对所述振型角控制信号序列进行初始阻尼角计算,得到初始阻尼角数据,其中,所述第一阻尼角计算公式如下所示:
其中,为初始阻尼角,/>为第一振型角控制信号,/>为第二振型角控制信号,为第三振型角控制信号;
通过预置的第一周向阻尼不均计算公式,对所述振型角控制信号序列进行初始周向阻尼不均计算,得到初始周向阻尼不均数据,其中,所述第一周向阻尼不均计算公式如下所示:
其中,为初始周向阻尼不均数据;
对所述初始阻尼角数据以及所述初始周向阻尼不均数据进行稳态谐振频率分析,得到稳态谐振频率;
基于所述稳态谐振频率,对所述谐振频率信号序列进行周向阻尼不均以及阻尼角计算,得到所述周向阻尼不均数据集合以及所述阻尼角数据集合。
4.根据权利要求3所述的半球谐振陀螺的温漂误差检测方法,其特征在于,所述基于所述稳态谐振频率,对所述谐振频率信号序列进行周向阻尼不均以及阻尼角计算,得到所述周向阻尼不均数据集合以及所述阻尼角数据集合步骤,包括:
基于所述稳态谐振频率,通过预置的第二阻尼角计算公式,对所述谐振频率信号序列进行阻尼角转换,得到对应的阻尼角数据集合,其中,所述第二阻尼角计算公式如下所示:
其中,为阻尼角数据,/>为初始阻尼角谐振频率系数,/>为第一阻尼角谐振频率系数,/>为第二阻尼角谐振频率系数,/>为第三阻尼角谐振频率系数;/>为稳态谐振频率,/>为谐振频率信号;
基于所述稳态谐振频率,通过预置的第二周向阻尼不均计算公式,对所述谐振频率信号序列进行周向阻尼不均计算,得到对应的周向阻尼不均数据集合,其中,所述第二周向阻尼不均计算公式如下所示:
其中,为周向阻尼不均数据,/>为初始周向阻尼不均谐振频率系数,/>为第一周向阻尼不均谐振频率系数,/>为第二周向阻尼不均谐振频率系数,/>为第三周向阻尼不均谐振频率系数。
5.根据权利要求1所述的半球谐振陀螺的温漂误差检测方法,其特征在于,所述对所述周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合步骤,包括:
对所述周向阻尼不均数据集合进行序列转换,得到对应的周向阻尼不均序列;
对所述阻尼角数据集合进行序列转换,得到对应的阻尼角序列;
通过预置的第一谐振频率系数计算公式对所述周向阻尼不均序列进行第一谐振频率系数计算,得到第一谐振频率系数集合,其中,所述第一谐振频率系数计算公式如下所示:
其中,为第一谐振频率系数集合,/>为谐振频率信号序列,/>为周向阻尼不均序列;
对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合。
6.根据权利要求5所述的半球谐振陀螺的温漂误差检测方法,其特征在于,所述对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合步骤,包括:
通过预置的第二谐振频率系数计算公式对所述阻尼角序列进行第二谐振频率系数计算,得到第二谐振频率系数集合,其中,所述第二谐振频率系数计算公式如下所示:
其中,为第二谐振频率系数集合,/>阻尼角序列。
7.一种半球谐振陀螺的温漂误差检测系统,用以执行如权利要求1至6任一项所述的半球谐振陀螺的温漂误差检测方法,其特征在于,包括:
获取模块,用于获取预设的测试参数集合,其中,所述测试参数集合包括:测试周期以及初始测试温度;
测试模块,用于通过所述测试参数集合对预置的半球谐振陀螺进行温度测试,同时,采集所述半球谐振陀螺在温度测试过程中的谐振频率信号序列以及振型角控制信号序列;
第一计算模块,用于基于所述振型角控制信号序列,通过所述谐振频率信号序列对所述半球谐振陀螺进行周向阻尼不均数据计算,得到周向阻尼不均数据集合,同时,对所述谐振频率信号序列进行阻尼角计算,得到阻尼角数据集合;
第二计算模块,用于对所述周向阻尼不均数据集合进行第一谐振频率系数计算,得到第一谐振频率系数集合,同时,对所述阻尼角数据集合进行第二谐振频率系数计算,得到第二谐振频率系数集合;
补偿模块,用于对所述第一谐振频率系数集合以及所述第二谐振频率系数集合进行温度数据匹配,得到目标温度数据,并通过所述目标温度数据对所述半球谐振陀螺进行温漂补偿。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311380657.0A CN117109639B (zh) | 2023-10-24 | 2023-10-24 | 半球谐振陀螺的温漂误差检测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311380657.0A CN117109639B (zh) | 2023-10-24 | 2023-10-24 | 半球谐振陀螺的温漂误差检测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117109639A CN117109639A (zh) | 2023-11-24 |
CN117109639B true CN117109639B (zh) | 2023-12-22 |
Family
ID=88811413
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311380657.0A Active CN117109639B (zh) | 2023-10-24 | 2023-10-24 | 半球谐振陀螺的温漂误差检测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117109639B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005103618A1 (ja) * | 2004-04-20 | 2005-11-03 | Murata Manufacturing Co., Ltd. | ジャイロ装置 |
CN115790665A (zh) * | 2023-01-05 | 2023-03-14 | 中国船舶集团有限公司第七〇七研究所 | 一种陀螺误差补偿方法、装置、电子设备及存储介质 |
CN115824263A (zh) * | 2023-02-13 | 2023-03-21 | 中国船舶集团有限公司第七〇七研究所 | 基于半球谐振陀螺的阻尼修调方法及系统 |
CN115950451A (zh) * | 2023-01-05 | 2023-04-11 | 中国船舶集团有限公司第七0七研究所 | 一种全角半球谐振陀螺阻尼不均漂移补偿方法 |
CN116499497A (zh) * | 2023-06-28 | 2023-07-28 | 中国船舶集团有限公司第七〇七研究所 | 基于模态反转的半球谐振陀螺惯导系统误差自校准方法 |
CN116608890A (zh) * | 2023-07-21 | 2023-08-18 | 中国船舶集团有限公司第七〇七研究所 | 一种全角模式半球谐振陀螺的标度误差补偿方法 |
CN116698084A (zh) * | 2023-07-13 | 2023-09-05 | 中国人民解放军国防科技大学 | 一种对称式陀螺零偏漂移实时抑制与零偏稳定性提升方法 |
CN116817964A (zh) * | 2023-05-22 | 2023-09-29 | 厦门汇鸿嘉创科技有限公司 | 一种半球谐振陀螺仪阻尼不均匀误差补偿方法 |
CN116907541A (zh) * | 2023-07-14 | 2023-10-20 | 长安大学 | 一种半球谐振陀螺的误差补偿方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6473713B1 (en) * | 1999-09-20 | 2002-10-29 | American Gnc Corporation | Processing method for motion measurement |
US7437253B2 (en) * | 2004-07-29 | 2008-10-14 | The Boeing Company | Parametrically disciplined operation of a vibratory gyroscope |
-
2023
- 2023-10-24 CN CN202311380657.0A patent/CN117109639B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005103618A1 (ja) * | 2004-04-20 | 2005-11-03 | Murata Manufacturing Co., Ltd. | ジャイロ装置 |
CN115790665A (zh) * | 2023-01-05 | 2023-03-14 | 中国船舶集团有限公司第七〇七研究所 | 一种陀螺误差补偿方法、装置、电子设备及存储介质 |
CN115950451A (zh) * | 2023-01-05 | 2023-04-11 | 中国船舶集团有限公司第七0七研究所 | 一种全角半球谐振陀螺阻尼不均漂移补偿方法 |
CN115824263A (zh) * | 2023-02-13 | 2023-03-21 | 中国船舶集团有限公司第七〇七研究所 | 基于半球谐振陀螺的阻尼修调方法及系统 |
CN116817964A (zh) * | 2023-05-22 | 2023-09-29 | 厦门汇鸿嘉创科技有限公司 | 一种半球谐振陀螺仪阻尼不均匀误差补偿方法 |
CN116499497A (zh) * | 2023-06-28 | 2023-07-28 | 中国船舶集团有限公司第七〇七研究所 | 基于模态反转的半球谐振陀螺惯导系统误差自校准方法 |
CN116698084A (zh) * | 2023-07-13 | 2023-09-05 | 中国人民解放军国防科技大学 | 一种对称式陀螺零偏漂移实时抑制与零偏稳定性提升方法 |
CN116907541A (zh) * | 2023-07-14 | 2023-10-20 | 长安大学 | 一种半球谐振陀螺的误差补偿方法及系统 |
CN116608890A (zh) * | 2023-07-21 | 2023-08-18 | 中国船舶集团有限公司第七〇七研究所 | 一种全角模式半球谐振陀螺的标度误差补偿方法 |
Non-Patent Citations (6)
Title |
---|
An accurate thermoelastic model and thermal output error analysis of a hemispherical resonator gyroscope under varying temperatures;Zeyuan Xu , Guoxing Yi , Weidong Zhu;《Mechanical Systems and Signal Processing》;1-18 * |
Researched on the bias stability of the HRG affected by the temperature and the standing wave azimuth;Wanliang Zhao, Hao Yang, Lijun Song;《Measurement and Control》;1561-2158 * |
全角模式半球谐振陀螺阻尼误差补偿方法研究;郭锞琛 等;《自动化仪表》;6-9+13 * |
基于多项式模型的半球谐振陀螺温度漂移建模;李广胜;蒋英杰;孙志强;谢红卫;;测试技术学报(第06期);496-500 * |
张勇猛 ; 郭锞琛 ; 孙江坤 ; 余升 ; 肖定邦 ; 吴学忠.全角模式半球谐振陀螺的阻尼误差修调与补偿技术研究.《机械工程学报》.2021,145-152. * |
赵小明 ; 陈刚 ; 郑志鹏 ; 唐明浩 ; 贾晨凯.半球谐振陀螺惯导系统关键技术及试验验证.《2023年全国智能导航学术论文集》.2023,115-123 . * |
Also Published As
Publication number | Publication date |
---|---|
CN117109639A (zh) | 2023-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110608862A (zh) | 一种通过倾角传感器测量桥梁动态均值挠度方法 | |
CN110500238B (zh) | 确定风力发电机组的塔顶的振动位移的方法和设备 | |
CN113155114A (zh) | Mems惯性测量单元陀螺零位的温度补偿方法及装置 | |
CN108107233B (zh) | 加速度计标度因数的连续温度校正方法及系统 | |
CN102607598B (zh) | 陀螺仪硬件生成方法与装置 | |
CN110285832B (zh) | 光纤陀螺惯测产品的标定参数长期稳定性缺陷检测方法 | |
CN115435768A (zh) | 一种基于实时滑窗的半球谐振陀螺温度建模补偿方法 | |
CN111829503A (zh) | 一种光纤陀螺阈值测试方法及装置 | |
CN111624671A (zh) | 旋转加速度计重力梯度仪重力梯度解调相位角确定方法及装置 | |
Edamana et al. | Estimation with threshold sensing for gyroscope calibration using a piezoelectric microstage | |
CN117109639B (zh) | 半球谐振陀螺的温漂误差检测方法及系统 | |
Yan et al. | Improving calibration accuracy of a vibration sensor through a closed loop measurement system | |
CN110987016A (zh) | 一种动力调谐陀螺挠性组件动态特性校准系统 | |
CN108151764B (zh) | 一种激光陀螺高低温锁区测试装置及测试方法 | |
CN116821713B (zh) | 基于多变量动态时间规整算法的隔震效率评价方法及系统 | |
CN116499696B (zh) | 一种提高风洞模型试验模型姿态角动态精度的方法 | |
Ji et al. | Adaptive Kalman filter enhanced with spectrum analysis for wide-bandwidth angular velocity estimation fusion | |
Ren et al. | A new method for calibrating nonlinear coefficients of PIGA on linear vibrator | |
CN209623725U (zh) | 一种光纤线圈动态仿真测试装置 | |
CN108225377B (zh) | 一种通信铁塔倾角测量系统的性能测试方法 | |
CN111707288B (zh) | 一种无水平台的微型惯性传感器零位测量方法 | |
RU2515424C1 (ru) | Установка для измерения собственной частоты колебаний роторов силовых гироскопов | |
Arpaia et al. | An ultrasonic heading goniometer intrinsically robust to magnetic interference | |
CN117128956B (zh) | 基于角速度转换的动态倾角获取方法及应用该方法的设备 | |
RU2783189C1 (ru) | Способ контроля физических параметров резонатора твердотельного волнового гироскопа |
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 |