CN114166242B - 一种半球谐振陀螺检测信号不均匀性的校准方法及系统 - Google Patents
一种半球谐振陀螺检测信号不均匀性的校准方法及系统 Download PDFInfo
- Publication number
- CN114166242B CN114166242B CN202111280905.5A CN202111280905A CN114166242B CN 114166242 B CN114166242 B CN 114166242B CN 202111280905 A CN202111280905 A CN 202111280905A CN 114166242 B CN114166242 B CN 114166242B
- Authority
- CN
- China
- Prior art keywords
- compensation
- fitting
- hemispherical resonator
- parameters
- precession 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000001514 detection method Methods 0.000 title claims description 20
- 230000008859 change Effects 0.000 claims abstract description 14
- 238000013016 damping Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 10
- 238000009828 non-uniform distribution Methods 0.000 abstract description 6
- 230000006870 function Effects 0.000 description 18
- 238000004590 computer program Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000004075 alteration Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012067 mathematical method Methods 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 238000009987 spinning Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004070 electrodeposition Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process 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
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
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
本发明涉及一种半球谐振陀螺检测信号不均匀性的校准方法及系统,设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数;使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到振型的进动角θ顺时针和逆时针转动时,振型的进动角的角速度随振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到误差项Err关于振型的进动角θ的拟合函数,确定拟合函数中的拟合参数;基于各个补偿参数对两个方向的振动信号进行补偿校准,得到拟合参数随补偿系数变化的关系,根据拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对两个方向的振动信号进行补偿校准。使用最佳补偿系数,可以将标度因数随角度的不均匀分布抑制到最低。
Description
技术领域
本发明涉及半球谐振陀螺技术领域,尤其涉及一种半球谐振陀螺检测信号不均匀性的校准方法及系统。
背景技术
HRG(Hemispherical Resonator Gyro,半球谐振陀螺)是一种没有高速转子和活动支承的振动陀螺,具有精度高、质量小、体积小、启动时间短、高过载、高可靠性的特点。被誉为最具潜力的哥式振动陀螺。HRG的工作原理是基于半球壳谐振子绕中心轴旋转时产生的哥氏效应,而使其振型在环向相对壳体进动的物理机制。
半球谐振陀螺工作在二阶振动模态,常采用相隔45°分布的金属电极,通过静电力对HRG的谐振子进行振动驱动,并通过检测金属电极的静电电容获取谐振子的振幅、相位等振动信息。半球谐振陀螺的电极制作或者安装不均匀,探测电路的增益不均匀、串扰等都可能造成振动信号检测误差,进而造成角度解算误差,从而引入半球谐振陀螺的零偏或标度因数随角度分布不均匀。现有公开的文献或专利中对半球谐振陀螺电极增益误差的校准方法,对半球谐振陀螺的电极增益误差抑制效果有限。
发明内容
本发明针对现有技术中存在的技术问题,提供一种半球谐振陀螺检测信号不均匀性的校准方法及系统,使用最佳补偿系数,可以将标度因数随角度的不均匀分布抑制到最低,可以准确的校准半球谐振陀螺检测信号不均匀性
根据本发明的第一方面,提供了一种半球谐振陀螺检测信号不均匀性的校准方法,包括:步骤1,设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数;
步骤2,使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到所述振型的进动角θ顺时针和逆时针转动时,所述振型的进动角的角速度随所述振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到所述误差项Err关于所述振型的进动角θ的拟合函数,确定所述拟合函数中的拟合参数;
步骤3,基于各个所述补偿参数对所述两个方向的振动信号进行补偿校准,得到所述拟合参数随所述补偿系数变化的关系,根据所述拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对所述两个方向的振动信号进行补偿校准。
在上述技术方案的基础上,本发明还可以作出如下改进。
可选的,所述补偿参数包括补偿系数α和补偿系数β;
所述补偿系数α的对所述两个方向的振动信号x,y的补偿校准方法为:
x′=x;
y′=x-αy;
所述补偿系数β的对所述两个方向的振动信号x,y的补偿校准方法为:
x′=x;
y′=(1+β)y;
其中,x′,y′为通过所述补偿系数α和/或补偿系数β校准后的两通道检测信号。
可选的,所述步骤3中通过补偿系数α和补偿系数β中的一个或两个对所述两个方向的振动信号x,y进行一次或多次的补偿校准。
可选的,所述步骤2还包括:
通过幅度控制将半球谐振陀螺振幅控制为预设值,通过正交控制消除正交误差,通过角度控制实现改变或维持所述振型的进动角θ的角度。
可选的,所述两个方向的振动信号为:
x=a cosθcos(ωt)-qsinθsin(ωt);
y=asinθcos(ωt)+q cosθsin(ωt);
其中,ω为半球谐振陀螺二阶振动模态的圆频率,t为时间;确定a和q的值;
所述半球谐振陀螺的幅度控制动态方程为:
所述半球谐振陀螺的正交控制动态方程为:
所述半球谐振陀螺的角度控制动态方程为:
其中, 为阻尼系数;/>为阻尼不均匀系数;θτ为阻尼不均匀性的取向角度;Δω为半球谐振陀螺的频率裂解;θω为刚度不均匀的取向角;kWA为全角模式下半球谐振陀螺的标度因数;Ω为陀螺旋转的角速度;fas、fqc和fqs分别为控制力;
所述通过幅度控制将半球谐振陀螺振幅控制为预设值的方法为:通过对E进行PID控制获取控制力fas的控制信号,使E控制在预设的值;
所述通过正交控制消除正交误差的方法为:通过对Q进行PID控制获取控制力fqc的控制信号,使Q被控制到0;
所述通过角度控制实现改变或维持所述振型的进动角θ的角度的方法为:对力平衡模式,对振型的进动角θ进行PID控制,使其维持在固定位置,获得控制力fqs,通过测量控制力fqs获取输入角速度Ω信息;对全角模式,使控制力fqs为0,让陀螺的角度控制处于开环状态。
可选的,所述步骤2中所述振型的进动角的角速度随着所述振型的进动角θ改变时满足:
其中,Ω为陀螺旋转的角速度;k0为标度因数,k1,k2分别为标度因数的二次谐波余弦分量误差和正弦分量误差;B(θ)为半球谐振陀螺零偏。
可选的,所述步骤2中拟合得到所述误差项Err关于所述振型的进动角θ的拟合函数为:
Err=K0+K1 cos2θ+K2 sin2θ;
K0,K1,K2为所述拟合参数。
可选的,所述步骤3中当所述拟合参数K1,K2被抑制到0时,所对应的所述补偿参数为最佳的补偿参数。
根据本发明的第二方面,提供一种半球谐振陀螺检测信号不均匀性的校准系统,包括:补偿参数设置模块、拟合参数确定模块和补偿校准模块;
所述补偿参数设置模块,用于设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数;
所述拟合参数确定模块,用于使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到所述振型的进动角θ顺时针和逆时针转动时,所述振型的进动角的角速度随所述振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到所述误差项Err关于所述振型的进动角θ的拟合函数,确定所述拟合函数中的拟合参数;
所述补偿校准模块,用于基于各个所述补偿参数对所述两个方向的振动信号进行补偿校准,得到所述拟合参数随所述补偿系数变化的关系,根据所述拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对所述两个方向的振动信号进行补偿校准。
根据本发明的第三方面,提供了一种电子设备,包括存储器、处理器,所述处理器用于执行存储器中存储的计算机管理类程序时实现半球谐振陀螺检测信号不均匀性的校准方法的步骤。
根据本发明的第四方面,提供了一种计算机可读存储介质,其上存储有计算机管理类程序,所述计算机管理类程序被处理器执行时实现半球谐振陀螺检测信号不均匀性的校准方法的步骤。
本发明提供的一种半球谐振陀螺检测信号不均匀性的校准方法、系统、电子设备及存储介质,考虑存在着半球谐振陀螺的检测信号误差,造成半球谐振陀螺检测到的x,y信号是真实x,y信号互相混合而来,为了消除这种混合误差,引入补偿系数对x,y信号进行补偿,并计算各个补偿系数对x,y信号进行补偿后对应的误差项中关于振型的进动角的系数参数,确定系数参数为0时的补偿系数即为最佳补偿参数,可以将标度因数随角度的不均匀分布抑制到最低,可以准确的校准半球谐振陀螺检测信号不均匀性;补偿参数可以设置为不止一个的补偿系数,对一个信号进行正向和反向补偿,不同补偿系数进行累加或单独使用,更容易找到最佳补偿系数。
附图说明
图1为本发明提供的一种半球谐振陀螺检测信号不均匀性的校准方法的流程图;
图2为本发明一种半球谐振陀螺检测信号不均匀性的校准系统的结构图;
图3为本发明提供的一种可能的电子设备的硬件结构示意图;
图4为本发明提供的一种可能的计算机可读存储介质的硬件结构示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
图1为本发明提供的一种半球谐振陀螺检测信号不均匀性的校准方法的流程图,由图1可知,该校准方法包括:
步骤1,设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数;
在一种可能的实施例中,两个方向的振动信号为:
x=a cosθcos(ωt)-qsinθsin(ωt);
y=asinθcos(ωt)+q cosθsin(ωt);
其中,ω为半球谐振陀螺二阶振动模态的圆频率,t为时间;确定a和q的值。
补偿参数包括补偿系数α和补偿系数β。
补偿系数α的对两个方向的振动信号x,y的补偿校准方法为:
x′=x;
y′=x-αy;
补偿系数β的对两个方向的振动信号x,y的补偿校准方法为:
x′=x;
y′=(1+β)y;
其中,x′,y′为通过补偿系数α和/或补偿系数β校准后的两通道检测信号。
步骤2,使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到振型的进动角θ顺时针和逆时针转动时,振型的进动角的角速度随振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到误差项Err关于振型的进动角θ的拟合函数,确定拟合函数中的拟合参数。
Err=f1(θ)-f2(θ)。
在一种可能的实施例中,步骤2还包括:
通过幅度控制将半球谐振陀螺振幅控制为预设值,通过正交控制消除正交误差,通过角度控制实现改变或维持振型的进动角θ的角度。
具体的,可以通过半球谐振陀螺角度控制或者让半球谐振陀螺绕其敏感轴匀速转动,使振型的进动角θ匀速改变。
具体的,半球谐振陀螺的幅度控制动态方程为:
半球谐振陀螺的正交控制动态方程为:
半球谐振陀螺的角度控制动态方程为:
其中, 为阻尼系数;/>为阻尼不均匀系数;θτ为阻尼不均匀性的取向角度;Δω为半球谐振陀螺的频率裂解;θω为刚度不均匀的取向角;kWA为;Ω为陀螺旋转的角速度;fas、fqc和fqs分别为控制力。
通过幅度控制将半球谐振陀螺振幅控制为预设值的方法为:通过对E进行PID控制获取控制力fas的控制信号,使E控制在预设的值。
通过正交控制消除正交误差的方法为:通过对Q进行PID控制获取控制力fqc的控制信号,使Q被控制到0。
通过角度控制实现改变或维持振型的进动角θ的角度的方法为:对力平衡模式,对振型的进动角θ进行PID控制,使其维持在固定位置,获得控制力fqs,通过测量控制力fqs获取输入角速度Ω信息;对全角模式,使控制力fqs为0,让陀螺的角度控制处于开环状态。
在一种可能的实施例中,步骤2中振型的进动角的角速度随着振型的进动角θ改变时满足:
其中,Ω为陀螺旋转的角速度,其取决于角度控制的控制力大小,或者陀螺转动角速度大小;k0为标度因数,k1,k2分别为标度因数的二次谐波余弦分量误差和正弦分量误差;B(θ)为半球谐振陀螺零偏,其也随着θ的变化而变化,但不随输入角速度Ω变化。
在一种可能的实施例中,步骤2中拟合得到误差项Err关于振型的进动角θ的拟合函数为:
Err=K0+K1 cos2θ+K2 sin2θ;
K0,K1,K2为拟合参数。
具体实施中,可以通过合适的数学方法,例如快速傅里叶变换或者最小二乘法,得到拟合参数K0,K1,K2的测量拟合值。
步骤3,基于各个补偿参数对两个方向的振动信号进行补偿校准,得到拟合参数随补偿系数变化的关系,根据拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对两个方向的振动信号进行补偿校准。
在一种可能的实施例中,步骤3中通过补偿系数α和补偿系数β中的一个或两个对两个方向的振动信号x,y进行一次或多次的补偿校准。
两个补偿系数α、β分别表征两种检测增益误差源,步骤3中改变两个补偿系数α、β对半球谐振陀螺两个方向的振动信号进行校准时,补偿参数α,β的补偿顺序可以互换,也可以单独只对α,β之一进行补偿,也可以针对α,β或其中之一进行多次补偿。上述方法对α,β的值远小于1(特别是小于0.1)的情形特别有利。
步骤3中通过改变α和β的数值,测试测量拟合值K1,K2,可以得到K1,K2随α和β的变化关系。
在一种可能的实施例中,步骤3中当拟合参数K1,K2被抑制到0时,所对应的补偿参数为最佳的补偿参数。即所对应的补偿系数α和β即为最佳补偿参数,使用该最佳补偿参数,可以将标度因数随角度的不均匀分布抑制到最低。
提供一种半球谐振陀螺检测信号不均匀性的校准方法及系统,使用最佳补偿系数,可以将标度因数随角度的不均匀分布抑制到最低,可以准确的校准半球谐振陀螺检测信号不均匀性。
实施例1
本发明提供的实施例1为本发明提供的一种半球谐振陀螺检测信号不均匀性的校准方法的实施例,结合图1可知,该实施例中,半球谐振陀螺可用二维谐振子模型的数学物理方程描述,该方程可以等价为幅度控制动态方程、正交控制动态方程、角度控制动态方程、相位控制动态方程四个方程组,该校准方法的实施例包括:
半球谐振陀螺运动方程可以用如下二维谐振子模型表述:
其中x和y分别是半球谐振陀螺相差45°的两个方向的振动信号,而fx,fy为上述两个方向(或者其等效方向)上施加的控制力。k为全角模式标度因数,Ω为陀螺旋转的角速度,为阻尼系数,/>为阻尼不均匀系数,θτ为阻尼不均匀性的取向角度,ω为半球谐振陀螺二阶振动模态的圆频率,Δω为半球谐振陀螺的频率裂解,θω为刚度不均匀的取向角,θ为振型的进动角。x和y具有如下形式:
x=a cosθcos(ωt)-qsinθsin(ωt)
y=asinθcos(ωt)+q cosθsin(ωt)
在以x和y信号张成的空间中,上述方程是一个椭圆。a为椭圆的长轴,q为椭圆的短轴。
对x和y进行乘法解调,得到四个变量:
cx=LPF(x×cos(ωt+φ))
sx=LPF(x×sin(ωt+φ))
cy=LPF(y×cos(ωt+φ))
sy=LPF(y×sin(ωt+φ))
得到
构造五个中间变量:
振型的进动角θ可以用如下公式得到:
tan 2θ=S/R
由此可以得到幅度控制动态方程、正交控制动态方程、角度控制动态方程和相位控制动态方程分别为:
其中,和/>分别为E、Q、θ和φ对时间的导数,kWA为全角模式下半球谐振陀螺的标度因数。fac和fas为半球谐振陀螺振动模态椭圆长轴方向力的同相和正交分量,fqc和fqs为半球谐振陀螺振动模态椭圆短轴方向力的同相和正交分量;通过对E进行PID控制可以获取控制力fas的控制信号,使E控制在预设的值E0。通过对Q进行PID控制可以获取控制力fqc的控制信号,使Q被控制到0。对力平衡模式,对θ进行PID控制,使其维持在固定位置,可以获得控制力fqs,通过测量fqs获取输入角速度Ω信息;对全角模式,可以让fqs为0,让陀螺的角度控制处于开环状态,也可以认为施加fqs,控制进动角按照预设的方式改变。
在x和y电极上加的力为:
fx=fas cosθsin(ωt+φ)-fqcsinθcos(ωt+φ)-fqssinθsin(ωt+φ)
fy=fassinθsin(ωt+φ)+fqccosθcos(ωt+φ)+fqs cosθsin(ωt+φ)
依据上述模型进行控制时,存在着半球谐振陀螺的检测信号误差,例如半球谐振陀螺的电极位置由于电极制作工艺原因,正好相隔45°,由此会造成半球谐振陀螺检测到的x,y信号是真实x,y信号互相混合而来,为了消除这种混合误差,需要引入补偿系数α,其作用方法为:
x′=x
y′=x-αy
其中,x,y为经过α参数校准前的半球谐振陀螺相差45°的两个方向的振动信号,x′,y′为通过参数α校准后的半球谐振陀螺相差45°的两个方向的振动信号。
半球谐振陀螺电极间隙不均匀,或者半球谐振陀螺检测电路对x,y增益不一致,都会造成控制电路识别的x,y信号与真实的x,y信号不一致,为了消除这种增益误差,需要引入补偿系数β,其作用方法为:
x′=x
y′=(1+β)y
其中,x,y为经过β参数校准前的半球谐振陀螺相差45°的两个方向的振动信号,x′,y′为通过参数β校准后的半球谐振陀螺相差45°的两个方向的振动信号。
为了获取最佳α,β系数,本发明提供了如下方法:通过半球谐振陀螺角度控制,或者让半球谐振陀螺绕其敏感轴匀速转动,让半球谐振陀螺振型的进动角θ匀速改变。振型的进动角的角速度随着角度改变,其满足如下方程:
其中Ω为陀螺旋转的角速度,其取决于角度控制的控制力大小,或者陀螺转动角速度大小。k0为标度因数,k1,k2为标度因数的二次谐波余弦分量误差和正弦分量误差。B(θ)为半球谐振陀螺零偏,其也随着θ的变化而变化,但不随输入角速度Ω变化。当振型的进动角顺时针转动时,可以检测到随θ的变化关系为f1(θ),当振型的进动角逆时针转动时,可以检测到/>随θ的变化关系为f2(θ)。计算误差项:
Err=f1(θ)-f2(θ)
将测量得到的误差项可以写成如下形式
Err=K0+K1 cos2θ+K2 sin2θ
通过合适的数学方法,例如快速傅里叶变换或者最小二乘法,得到参数K0,K1,K2的测量拟合值。
改变α和β的数值,测试测量拟合值K1,K2,可以得到K1,K2随α和β的变化关系。当K1,K2被抑制到0时,所对应的α和β即为最佳补偿系数。使用最佳补偿系数,可以将标度因数随角度的不均匀分布抑制到最低。
当选取x信号对应的电极方向对应的振型的进动角规定θ=0°方向时,本方法将特别有利:当改变α时,K2将发生改变,而K1将保持不变,当K2被抑制为0时,所对应的α参数即为最佳值;当改变β时,K1将发生改变,而K2将保持不变,当K1被抑制为0时,所对应的β参数即为最佳值。
实施例2
本发明提供的实施例2为本发明提供的一种半球谐振陀螺检测信号不均匀性的校准系统的实施例,图2为本发明实施例提供的一种半球谐振陀螺检测信号不均匀性的校准系统结构图,结合图2可知,该实施例包括:补偿参数设置模块、拟合参数确定模块和补偿校准模块。
补偿参数设置模块,用于设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数。
拟合参数确定模块,用于使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到振型的进动角θ顺时针和逆时针转动时,振型的进动角的角速度随振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到误差项Err关于振型的进动角θ的拟合函数,确定拟合函数中的拟合参数。
补偿校准模块,用于基于各个补偿参数对两个方向的振动信号进行补偿校准,得到拟合参数随补偿系数变化的关系,根据拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对两个方向的振动信号进行补偿校准。
可以理解的是,本发明提供的一种半球谐振陀螺检测信号不均匀性的校准系统与前述各实施例提供的半球谐振陀螺检测信号不均匀性的校准方法相对应,半球谐振陀螺检测信号不均匀性的校准系统的相关技术特征可参考半球谐振陀螺检测信号不均匀性的校准方法的相关技术特征,在此不再赘述。
请参阅图3,图3为本发明实施例提供的电子设备的实施例示意图。如图3所示,本发明实施例提了一种电子设备,包括存储器1310、处理器1320及存储在存储器1310上并可在处理器1320上运行的计算机程序1311,处理器1320执行计算机程序1311时实现以下步骤:设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数;使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到振型的进动角θ顺时针和逆时针转动时,振型的进动角的角速度随振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到误差项Err关于振型的进动角θ的拟合函数,确定拟合函数中的拟合参数;基于各个补偿参数对两个方向的振动信号进行补偿校准,得到拟合参数随补偿系数变化的关系,根据拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对两个方向的振动信号进行补偿校准。
请参阅图4,图4为本发明提供的一种计算机可读存储介质的实施例示意图。如图4所示,本实施例提供了一种计算机可读存储介质1400,其上存储有计算机程序1411,该计算机程序1411被处理器执行时实现如下步骤:设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数;使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到振型的进动角θ顺时针和逆时针转动时,振型的进动角的角速度随振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到误差项Err关于振型的进动角θ的拟合函数,确定拟合函数中的拟合参数;基于各个补偿参数对两个方向的振动信号进行补偿校准,得到拟合参数随补偿系数变化的关系,根据拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对两个方向的振动信号进行补偿校准。
本发明实施例提供的一种半球谐振陀螺检测信号不均匀性的校准方法、系统及存储介质,考虑存在着半球谐振陀螺的检测信号误差,造成半球谐振陀螺检测到的x,y信号是真实x,y信号互相混合而来,为了消除这种混合误差,引入补偿系数对x,y信号进行补偿,并计算各个补偿系数对x,y信号进行补偿后对应的误差项中关于振型的进动角的系数参数,确定系数参数为0时的补偿系数即为最佳补偿参数,可以将标度因数随角度的不均匀分布抑制到最低,可以准确的校准半球谐振陀螺检测信号不均匀性;补偿参数可以设置为不止一个的补偿系数,对一个信号进行正向和反向补偿,不同补偿系数进行累加或单独使用,更容易找到最佳补偿系数。
需要说明的是,在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详细描述的部分,可以参见其它实施例的相关描述。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式计算机或者其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包括这些改动和变型在内。
Claims (6)
1.一种半球谐振陀螺检测信号不均匀性的校准方法,其特征在于,所述校准方法包括:
步骤1,设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数;
步骤2,使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到所述振型的进动角θ顺时针和逆时针转动时,所述振型的进动角的角速度随所述振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到所述误差项Err关于所述振型的进动角θ的拟合函数,确定所述拟合函数中的拟合参数;
步骤3,基于各个所述补偿参数对所述两个方向的振动信号进行补偿校准,得到所述拟合参数随补偿系数变化的关系,根据所述拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对所述两个方向的振动信号进行补偿校准;
所述补偿参数包括补偿系数α和补偿系数β;
所述补偿系数α的对所述两个方向的振动信号x,y的补偿校准方法为:
x′=x;
y′=x-αy;
所述补偿系数β的对所述两个方向的振动信号x,y的补偿校准方法为:
x′=x;
y′=(1+β)y;
其中,x′,y′为通过所述补偿系数α和/或补偿系数β校准后的两通道检测信号;
所述步骤2中所述振型的进动角的角速度随着所述振型的进动角θ改变时满足:
其中,Ω为陀螺旋转的角速度;k0为标度因数,k1,k2分别为标度因数的二次谐波余弦分量误差和正弦分量误差;B(θ)为半球谐振陀螺零偏;
所述步骤2中拟合得到所述误差项Err关于所述振型的进动角θ的拟合函数为:
Err=K0+K1 cos 2θ+K2 sin 2θ;
K0,K1,K2为所述拟合参数;
所述步骤3中当所述拟合参数K1,K2被抑制到0时,所对应的所述补偿参数为最佳的补偿参数。
2.根据权利要求1所述的校准方法,其特征在于,所述步骤3中通过补偿系数α和补偿系数β中的一个或两个对所述两个方向的振动信号x,y进行一次或多次的补偿校准。
3.根据权利要求1所述的校准方法,其特征在于,所述步骤2还包括:
通过幅度控制将半球谐振陀螺振幅控制为预设值,通过正交控制消除正交误差,通过角度控制实现改变或维持所述振型的进动角θ的角度。
4.根据权利要求1所述的校准方法,其特征在于,所述两个方向的振动信号为:
x=a cosθcos(ωt)-qsinθsin(ωt);
y=asinθcos(ωt)+q cosθsin(ωt);
其中,ω为半球谐振陀螺二阶振动模态的圆频率,t为时间;确定a和q的值;
所述半球谐振陀螺的幅度控制动态方程为:
所述半球谐振陀螺的正交控制动态方程为:
所述半球谐振陀螺的角度控制动态方程为:
其中, 为阻尼系数;/>为阻尼不均匀系数;θτ为阻尼不均匀性的取向角度;Δω为半球谐振陀螺的频率裂解;θω为刚度不均匀的取向角;kWA为全角模式下半球谐振陀螺的标度因数;Ω为陀螺旋转的角速度;fas、fqc和fqs分别为控制力;
所述通过幅度控制将半球谐振陀螺振幅控制为预设值的方法为:通过对E进行PID控制获取控制力fas的控制信号,使E控制在预设的值;
所述通过正交控制消除正交误差的方法为:通过对Q进行PID控制获取控制力fqc的控制信号,使Q被控制到0;
所述通过角度控制实现改变或维持所述振型的进动角θ的角度的方法为:对力平衡模式,对振型的进动角θ进行PID控制,使其维持在固定位置,获得控制力fqs,通过测量控制力fqs获取输入角速度Ω信息;对全角模式,使控制力fqs为0,让陀螺的角度控制处于开环状态。
5.一种半球谐振陀螺检测信号不均匀性的校准系统,其特征在于,包括:补偿参数设置模块、拟合参数确定模块和补偿校准模块;
所述补偿参数设置模块,用于设置半球谐振陀螺相差45°的两个方向的振动信号的补偿参数;
所述拟合参数确定模块,用于使半球谐振陀螺振型的进动角θ匀速改变,分别检测得到所述振型的进动角θ顺时针和逆时针转动时,所述振型的进动角的角速度随所述振型的进动角θ的变化关系f1(θ)和f2(θ),计算f1(θ)和f2(θ)的差为误差项Err,拟合得到所述误差项Err关于所述振型的进动角θ的拟合函数,确定所述拟合函数中的拟合参数;
所述补偿校准模块,用于基于各个所述补偿参数对所述两个方向的振动信号进行补偿校准,得到所述拟合参数随所述补偿系数变化的关系,根据所述拟合参数的值得到最佳的补偿参数,使用该最佳的补偿参数对所述两个方向的振动信号进行补偿校准;
所述拟合参数确定模块中所述振型的进动角的角速度随着所述振型的进动角θ改变时满足:
其中,Ω为陀螺旋转的角速度;k0为标度因数,k1,k2分别为标度因数的二次谐波余弦分量误差和正弦分量误差;B(θ)为半球谐振陀螺零偏;
所述拟合参数确定模块中拟合得到所述误差项Err关于所述振型的进动角θ的拟合函数为:
Err=K0+K1 cos 2θ+K2 sin 2θ;
K0,K1,K2为所述拟合参数;
所述补偿校准模块中当所述拟合参数K1,K2被抑制到0时,所对应的所述补偿参数为最佳的补偿参数。
6.一种计算机可读存储介质,其特征在于,其上存储有计算机管理类程序,所述计算机管理类程序被处理器执行时实现如权利要求1-4任一项所述的半球谐振陀螺检测信号不均匀性的校准方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111280905.5A CN114166242B (zh) | 2021-10-29 | 2021-10-29 | 一种半球谐振陀螺检测信号不均匀性的校准方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111280905.5A CN114166242B (zh) | 2021-10-29 | 2021-10-29 | 一种半球谐振陀螺检测信号不均匀性的校准方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114166242A CN114166242A (zh) | 2022-03-11 |
CN114166242B true CN114166242B (zh) | 2023-08-08 |
Family
ID=80477624
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111280905.5A Active CN114166242B (zh) | 2021-10-29 | 2021-10-29 | 一种半球谐振陀螺检测信号不均匀性的校准方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114166242B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114509057B (zh) * | 2022-03-14 | 2023-06-20 | 中国船舶重工集团公司第七0七研究所 | 一种谐振陀螺仪全角模式控制方法 |
CN115790665B (zh) * | 2023-01-05 | 2023-05-05 | 中国船舶集团有限公司第七〇七研究所 | 一种陀螺误差补偿方法、装置、电子设备及存储介质 |
CN116499497B (zh) * | 2023-06-28 | 2023-09-05 | 中国船舶集团有限公司第七〇七研究所 | 基于模态反转的半球谐振陀螺惯导系统误差自校准方法 |
CN116608890B (zh) * | 2023-07-21 | 2023-10-13 | 中国船舶集团有限公司第七〇七研究所 | 一种全角模式半球谐振陀螺的标度误差补偿方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6189382B1 (en) * | 1999-11-05 | 2001-02-20 | Litton Systems, Inc. | Vibratory sensor with self-calibration and low noise digital conversion |
CN108489512A (zh) * | 2018-02-28 | 2018-09-04 | 北京控制工程研究所 | 一种半球谐振陀螺标度因数的补偿标定方法及装置 |
CN109945894A (zh) * | 2019-04-08 | 2019-06-28 | 哈尔滨工业大学 | 半球谐振陀螺仪谐振子与外基座的安装误差补偿方法 |
CN111504295A (zh) * | 2020-04-09 | 2020-08-07 | 哈尔滨工程大学 | 一种克服速率积分振动陀螺低速自锁效应的方法 |
CN111878068A (zh) * | 2020-05-11 | 2020-11-03 | 中国科学院地质与地球物理研究所 | 一种高温固态谐振陀螺仪及其组成的钻井测量系统 |
CN112595303A (zh) * | 2020-12-10 | 2021-04-02 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种hrg全角模式下参数计算方法 |
CN112615251A (zh) * | 2020-12-15 | 2021-04-06 | 北京航天控制仪器研究所 | 一种用于原子陀螺仪的激光频率和功率的双重稳定方法及装置 |
CN113252019A (zh) * | 2021-05-13 | 2021-08-13 | 哈尔滨工业大学 | 一种前向放大系数不一致时的半球谐振陀螺振型角获取方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9574902B2 (en) * | 2011-10-14 | 2017-02-21 | Innalabs Limited | Calibration method for the scale factor of an axisymmetric vibratory gyroscope or gyrometer |
US9605964B2 (en) * | 2014-01-03 | 2017-03-28 | The Boeing Company | Gyro quadrature stabalization with demodulation phase error nulling |
-
2021
- 2021-10-29 CN CN202111280905.5A patent/CN114166242B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6189382B1 (en) * | 1999-11-05 | 2001-02-20 | Litton Systems, Inc. | Vibratory sensor with self-calibration and low noise digital conversion |
CN108489512A (zh) * | 2018-02-28 | 2018-09-04 | 北京控制工程研究所 | 一种半球谐振陀螺标度因数的补偿标定方法及装置 |
CN109945894A (zh) * | 2019-04-08 | 2019-06-28 | 哈尔滨工业大学 | 半球谐振陀螺仪谐振子与外基座的安装误差补偿方法 |
CN111504295A (zh) * | 2020-04-09 | 2020-08-07 | 哈尔滨工程大学 | 一种克服速率积分振动陀螺低速自锁效应的方法 |
CN111878068A (zh) * | 2020-05-11 | 2020-11-03 | 中国科学院地质与地球物理研究所 | 一种高温固态谐振陀螺仪及其组成的钻井测量系统 |
CN112595303A (zh) * | 2020-12-10 | 2021-04-02 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种hrg全角模式下参数计算方法 |
CN112615251A (zh) * | 2020-12-15 | 2021-04-06 | 北京航天控制仪器研究所 | 一种用于原子陀螺仪的激光频率和功率的双重稳定方法及装置 |
CN113252019A (zh) * | 2021-05-13 | 2021-08-13 | 哈尔滨工业大学 | 一种前向放大系数不一致时的半球谐振陀螺振型角获取方法 |
Non-Patent Citations (1)
Title |
---|
张熙等.载体运动条件下半球谐振陀螺的运动方程.惯性传感技术与应用研讨会文集.2020,460-465. * |
Also Published As
Publication number | Publication date |
---|---|
CN114166242A (zh) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114166242B (zh) | 一种半球谐振陀螺检测信号不均匀性的校准方法及系统 | |
JP5421693B2 (ja) | 自己較正型ジャイロスコープ・システム | |
RU2687101C2 (ru) | Способ калибровки вибрационного гироскопа | |
CN113899365B (zh) | 基于半球谐振陀螺检测通路非对称性参数辨识的驻波方位角测量方法 | |
CN101802552B (zh) | 确定轴向对称振动传感器转速的方法及实施其的惯性装置 | |
CN105783899B (zh) | 具有原位偏差自校准的高带宽科里奥利振动陀螺仪(cvg) | |
JP2007520716A (ja) | 振動構造ジャイロスコープにおけるバイアス誤差を低減する方法 | |
CN113551660B (zh) | 一种电极角度存在误差时的半球谐振陀螺振型角获取方法 | |
CN114858191B (zh) | 一种全角半球谐振陀螺检测电极误差标定方法 | |
WO2004046649A2 (en) | Method of calibrating bias drift with temperature for a vibrating structure gyroscope | |
CN109163734A (zh) | 一种基于双轴光纤旋转调制组合导航系统的自主标定方法 | |
CN114509057A (zh) | 一种谐振陀螺仪全角模式控制方法 | |
CN114370886B (zh) | 基于虚拟转动的全角模式振动陀螺测量误差自标定方法 | |
CN114858184A (zh) | 一种半球谐振子参数辨识方法 | |
US9927256B2 (en) | Calibration of vibrating gyroscope | |
Guo et al. | Adaptive compensation of damping asymmetry in whole-angle hemispherical resonator gyroscope | |
Wang et al. | Rate integrating hemispherical resonator gyroscope detection error analysis and compensation | |
RU2499225C1 (ru) | Способ определения курса путем поворота инерциального устройства | |
CN109143861B (zh) | 一种基于力矩器的变速倾侧动量轮系统主动振动抑制方法 | |
US9709399B2 (en) | Approach for control redistribution of coriolis vibratory gyroscope (CVG) for performance improvement | |
KR20180102487A (ko) | 자이로스코프 필드-내 예측 | |
US8997567B2 (en) | Method of angular measurement by means of a vibrating sensor to which modulated controls are applied | |
RU2526585C2 (ru) | Способ определения угла ориентации стоячей волны в твердотельном волновом гироскопе | |
Sun et al. | Investigation of cylindrical resonators’ damping asymmetry via analyzing q factor circumferential distribution | |
Song et al. | Parametric drive of a micro rate integrating gyroscope using discrete electrodes |
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 |