CN116430713A - 一种用于提高全角半球谐振陀螺的控制回路带宽的方法 - Google Patents

一种用于提高全角半球谐振陀螺的控制回路带宽的方法 Download PDF

Info

Publication number
CN116430713A
CN116430713A CN202310411145.XA CN202310411145A CN116430713A CN 116430713 A CN116430713 A CN 116430713A CN 202310411145 A CN202310411145 A CN 202310411145A CN 116430713 A CN116430713 A CN 116430713A
Authority
CN
China
Prior art keywords
full
control
hemispherical resonator
controller
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.)
Granted
Application number
CN202310411145.XA
Other languages
English (en)
Other versions
CN116430713B (zh
Inventor
郜中星
蒋丹
马燕
南方伯
崔浩宇
付雨鑫
于宗帅
赵聪宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao Harbin Engineering University Innovation Development Center
Original Assignee
Qingdao Harbin Engineering University Innovation Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qingdao Harbin Engineering University Innovation Development Center filed Critical Qingdao Harbin Engineering University Innovation Development Center
Priority to CN202310411145.XA priority Critical patent/CN116430713B/zh
Publication of CN116430713A publication Critical patent/CN116430713A/zh
Application granted granted Critical
Publication of CN116430713B publication Critical patent/CN116430713B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Gyroscopes (AREA)

Abstract

本发明公开了一种用于提高全角半球谐振陀螺的控制回路带宽的方法,包括以下步骤:构建输入‑输出控制系统模型,获取全角半球谐振陀螺控制系统被控对象;基于所述被控对象,通过双线性变换法对控制器进行离散化处理;对离散化处理后的所述控制器的若干参数进行标定,引入速度反馈提高全角半球谐振陀螺闭环控制回路的带宽。本发明基于全角半球谐振陀螺转动引起扰动频率与转速存在特定关系,通过抑制转动引起谐波干扰提高全角半球谐振陀螺的控制回路带宽。

Description

一种用于提高全角半球谐振陀螺的控制回路带宽的方法
技术领域
本发明属于智能化仪器仪表技术领域,尤其涉及一种用于提高全角半球谐振陀螺的控制回路带宽的方法。
背景技术
半球谐振陀螺是一种哥式振动陀螺,半球谐振陀螺具有两种工作模式,分别为力平衡模式与全角模式。力平衡模式需要通过施加驱动力抵消哥式力对半球谐振子进动状态的影响,受限于力平衡回路的饱和作用,半球谐振陀螺的动态范围难以得到提高;而在全角模式下,谐振振型始终处于自由进动的状态,陀螺可以获得无限大的带宽,但控制回路的带宽会受限于外围电路以及控制策略。半球谐振陀螺控制策略关乎谐振陀螺仪能否保持在谐振频率处以稳定的幅度振荡工作,这是半球谐振陀螺仪正常工作的前提,其控制效果决定了控制回路带宽大小和陀螺仪精度的上限。全角半球谐振陀螺在经典比例积分(Proportional-Integral,PI)控制下可以实现静态和低速转动状态下的高精度测量,但陀螺仪高速转动时,PI控制器的跟踪效果不理想,难以继续使全角半球谐振陀螺保持高精度的工作。这种由陀螺仪转速引起的干扰信号呈现出正弦规律的变化,为了有效抑制正弦干扰项,在原有的PI控制基础上并联谐振控制组成比例积分谐振控制(Proportional-Integral-Resonant,PIR),对特定频率的谐波形成有效抑制同时不影响低频段的控制效果。
发明内容
为解决上述技术问题,本发明提出了一种用于提高全角半球谐振陀螺的控制回路带宽的方法,克服全角半球谐振陀螺高速转动时经典PI控制器跟踪效果不好问题,提高跟踪效果。
为实现上述目的,本发明提供了一种用于提高全角半球谐振陀螺的控制回路带宽的方法,包括以下步骤:
构建输入-输出控制系统模型,获取全角半球谐振陀螺控制系统被控对象;
基于所述被控对象,通过双线性变换法对控制器进行离散化处理;
对离散化处理后的所述控制器的若干参数进行标定,引入速度反馈提高全角半球谐振陀螺闭环控制回路的带宽。
可选的,构建输入-输出控制系统模型具体包括:
控制对象半球谐振子可建模为一阶惯性环节,整体增益有包含D/A转换及高压放大增益的驱动增益、包含滤波放大及A/D转换增益的检测增益以及椭圆参数解算增益。
可选的,获取全角半球谐振陀螺控制系统被控对象的方法包括:
Figure BDA0004183218030000021
K=KdriveKHRGKsenseKpara
其中,G1(s)为全角半球谐振陀螺控制系统被控对象,K为陀螺控制系统总体增益,τ为时间常数,s为s域下的变量,Kdrive为驱动环节增益,KHRG为半球谐振子幅值增益,Ksense为检测环节增益,Kpara为椭圆参数解算环节增益。
可选的,所述控制器为经典PI控制器并联谐振控制组成PIR控制器。
可选的,基于所述控制器获取传递函数的方法包括:
Figure BDA0004183218030000022
其中,GPIR(s)为传递函数,kp、ki、kr分别为比例、积分、谐振环节增益,ω0为谐振控制器的谐振频率,ωc为谐振控制器截止频率。
可选的,通过双线性变换法对控制器进行离散化处理的方法包括:
Figure BDA0004183218030000031
其中,GPIR(z)为z域下的传递函数,GPIR(s)为s域下的传递函数,z为z域下的变量,s为s域下的变量,T为采样时间。
可选的,所述全角半球谐振陀螺控制系统包括:幅值控制、正交控制和频率控制。
可选的,对离散化处理后的所述控制器的若干参数进行标定,引入速率反馈提高全角半球谐振陀螺闭环控制回路的带宽的方法包括:
所述控制器的系数随着外部输入转速变化而变化,谐振控制中心频率始终跟踪外部转速的四倍频,有效抑制由转速引起的扰动,提高控制回路带宽。
本发明技术效果:本发明公开了一种用于提高全角半球谐振陀螺的控制回路带宽的方法,揭示了全角半球谐振陀螺转动引起扰动频率与转速存在特定关系,提出了通过抑制转动引起谐波干扰来提高控制回路带宽的方法。
附图说明
构成本申请的一部分的附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为本发明实施例全角半球谐振陀螺控制回路原理框图;
图2为本发明实施例转化为数学模型后的输入-输出控制系统结构框图;
图3为本发明实施例PI控制与PIR控制的伯德图;
图4为本发明实施例PI控制与PIR控制下的控制回路伯德图;
图5为本发明实施例速率反馈下的PIR控制系统结构框图;
图6为本发明整体实现的流程图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
如图1-6所示,本实施例中提供一种用于提高全角半球谐振陀螺的控制回路带宽的方法,包括:
建立输入-输出控制系统模型,确定陀螺转动频率与谐波干扰频率的数量关系,包括:全角半球谐振陀螺控制系统包括幅值控制、正交控制和频率控制。本发明主要研究幅值控制和正交控制,以幅值控制系统为例,如图1所示。图中:以预设幅值为输入信号R(s),以解算得到的幅值信号为输出信号C(s)。
将图1中各个模块进行建模后,其控制系统结构图如图2所示。得到被控对象的传递函数:
Figure BDA0004183218030000041
其中K=KdriveKHRGKsenseKpara
由于半球谐振子具有很高的Q值,其衰减时间常数τ达到百秒级别,属于大惯性环节。经典PI控制器可以实现低速下的有效控制,但半球谐振子高速转动时,PI控制跟踪速度慢,无法实现精准控制,进而导致系统不稳定。陀螺转速引起干扰项会呈现出正弦变化,其变化周期和转速存在固定关系。
陀螺控制系统中的扰动误差主要来源于半球谐振子极板间距不等、电路元件误差等带来的检测回路增益不匹配。因此,将x检测通道与y检测通道的增益系数kx和ky分别加入到解算回路中,图1半球谐振子解算过程中x,y方向的位移表达式变为:
Figure BDA0004183218030000051
以主波波腹a为例分析检测通道增益对椭圆参数解算的影响,利用图1中解调得到的cx、cy、sx和sy四个参数分别可以得到:
Figure BDA0004183218030000052
解算到的主波波腹变ac为:
Figure BDA0004183218030000053
由式(3)可以得出:当kx=ky时,ac=a,检测通道增益对信号解算没有影响;当kx≠ky时,若kx>ky,ac会随着角度的增加呈现正弦变化,而且ac变化频率为角度θ变化频率的四倍。
并引入速率反馈,实现变速情况下精准控制,包括:PIR控制是在传统的PI控制基础上并联谐振(R)控制,实现对正弦信号的无差控制,其传递函数为:
Figure BDA0004183218030000054
式中,kp、ki、kr分别为比例、积分、谐振环节增益,ω0为谐振控制器的谐振频率,ωc为谐振控制器截止频率。其中,截止频率ωc决定了谐振频率的范围,ωc越大,谐振区域的频率范围越大,一般取值为5-15rad/s,设计时选用5rad/s;谐振频率ω0决定谐振点的频率,与待消除的谐波频率有关,由式(3)得到的结论,谐波信号与陀螺角速率有关,最后谐振频率取角速率的四倍,并且通过速率反馈实时调整。
图3为PI控制和PIR控制的伯德图,从图中可知,加入谐振控制后,将谐振频率附近的增益提高,而不改变控制器其它频段的幅频特性。
将图2中的控制器替换为PIR控制器,则R(s)输入时的控制回路闭环传递函数ΦR(s)可以表示为:
Figure BDA0004183218030000061
如图4所示,根据式(16)绘制出PIR控制下的全角半球谐振陀螺控制回路伯德图,并与PI控制下的伯德图进行对比。由图4可知,选择合适的PIR参数,可以将特定频率点的幅频特性提高至-3dB以上。
变速情况下,引入速率反馈,其控制系统结构如图5所示,N(s)为陀螺转动引入的干扰项,
Figure BDA0004183218030000062
为陀螺解算出的角速率。PIR控制器的比例控制环节和积分控制环节参数保持不变,谐振控制的截止频率ωc保持为5rad/s,谐振频率ω0为角速率/>
Figure BDA0004183218030000063
的四倍。谐振环节增益kr不是定值,角速率越大,需要的增益越大,考虑到两者不存在线性关系,可以采用分段选择的方式确定kr,先在固定转速下调试不同转速范围内的kr最优值(每0.1Hz确定一个kr),变速情况下通过判断角速率所处范围,设定与其对应的kr,实现实时调整。
引入速率反馈后,实时改变谐振控制的谐振点位置并调整控制参数,提高谐振点的幅值裕度,进而提高系统控制回路带宽。
一种PIR控制器数字离散化的方法,包括:在实际的陀螺仪控制系统中处理的信号为数字信号,需要将PIR控制器s域的传递函数转换到z域,实现离散化。采用双线性变换法实现PIR控制器的离散化,离散公式为:
Figure BDA0004183218030000064
式中,T为采样时间,与数字控制系统时钟频率有关。
将式(4)通分整理,得:
Figure BDA0004183218030000071
将式(7)带入式(6)中,为了表达方便,记m3=kp,m2=2kpωc+ki+2krωc
Figure BDA0004183218030000072
n3=1,n2=2ωc,/>
Figure BDA0004183218030000073
得:
Figure BDA0004183218030000074
其中,分子、分母各系数为:
Figure BDA0004183218030000075
将离散域下的传函转化为差分方程。
先将式(8)化为:
(a3+a2z-1+a1z-2+a0z-3)Y(z)=(b3+b2z-1+b1z-2+b0z-3)X(z) (10)
将式(21)化为差分方程:
a3y(k)=b3x(k)+b2x(k-1)+b1x(k-2)+b0x(k-3)-a2y(k-1)-a1y(k-2)-a0y(k-3) (11)
之后将差分方程用C语言实现,完成PIR控制器在FPGA软核的实现。
如图6所示,基于上述PIR控制器设计方法以及离散化方法,引入速率反馈,实现控制器系数自适应控制,实现变速条件下控制系统的快速准确跟踪,提高控制系统带宽,具体包括:
步骤1:根据半球谐振陀螺控制系统模型确定PIR控制器各参数的大致范围。
步骤2:采用双线性变换法实现PIR控制器的离散化,并用C语言在FPGA软核实现。
步骤3:将半球谐振陀螺放置于隔振台上,以全角模式运行半球谐振陀螺。将kr置零,调节kp、ki至陀螺仪静态控制效果最佳,记录kp、ki
步骤4:将半球谐振陀螺固联在单轴转台上,并以全角模式运行半球谐振陀螺,转台依次设定固定转速为36°/s(0.1Hz)、72°/s(0.2Hz)、108°/s(0.3Hz)……720°/s(2Hz)。将kp、ki设置为步骤3中的记录值,ωc设置为5rad/s,ω0设置为转台转速的四倍,调整kr至陀螺仪控制效果最佳,依次记录kr
步骤5:在步骤4的基础上,引入速率反馈,将ω0、kr设置为变量,ω0等于陀螺仪角速度测量值的四倍,采用分段选择的方式调整kr,其它参数保持步骤3、4中的记录值。
以上,仅为本申请较佳的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应该以权利要求的保护范围为准。

Claims (8)

1.一种用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,包括以下步骤:
构建输入-输出控制系统模型,获取全角半球谐振陀螺控制系统被控对象;
基于所述被控对象,通过双线性变换法对控制器进行离散化处理;
对离散化处理后的所述控制器的若干参数进行标定,引入速度反馈提高全角半球谐振陀螺闭环控制回路的带宽。
2.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,构建输入-输出控制系统模型具体包括:
控制对象半球谐振子可建模为一阶惯性环节,整体增益有包含D/A转换及高压放大增益的驱动增益、包含滤波放大及A/D转换增益的检测增益以及椭圆参数解算增益。
3.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,获取全角半球谐振陀螺控制系统被控对象的方法包括:
Figure FDA0004183218010000011
K=KdriveKHRGKsenseKpara
其中,G1(s)为全角半球谐振陀螺控制系统被控对象,K为陀螺控制系统总体增益,τ为时间常数,s为s域下的变量,Kdrive为驱动环节增益,KHRG为半球谐振子幅值增益,Ksense为检测环节增益,Kpara为椭圆参数解算环节增益。
4.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,所述控制器为经典PI控制器并联谐振控制组成PIR控制器。
5.如权利要求4所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,基于所述控制器获取传递函数的方法包括:
Figure FDA0004183218010000021
其中,GPIR(s)为传递函数,kp、ki、kr分别为比例、积分、谐振环节增益,ω0为谐振控制器的谐振频率,ωc为谐振控制器截止频率。
6.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,通过双线性变换法对控制器进行离散化处理的方法包括:
Figure FDA0004183218010000022
其中,GPIR(z)为z域下的传递函数,GPIR(s)为s域下的传递函数,z为z域下的变量,s为s域下的变量,T为采样时间。
7.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,所述全角半球谐振陀螺控制系统包括:幅值控制、正交控制和频率控制。
8.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,对离散化处理后的所述控制器的若干参数进行标定,引入速率反馈提高全角半球谐振陀螺闭环控制回路的带宽的方法包括:
所述控制器的系数随着外部输入转速变化而变化,谐振控制中心频率始终跟踪外部转速的四倍频,有效抑制由转速引起的扰动,提高控制回路带宽。
CN202310411145.XA 2023-04-18 2023-04-18 一种用于提高全角半球谐振陀螺的控制回路带宽的方法 Active CN116430713B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310411145.XA CN116430713B (zh) 2023-04-18 2023-04-18 一种用于提高全角半球谐振陀螺的控制回路带宽的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310411145.XA CN116430713B (zh) 2023-04-18 2023-04-18 一种用于提高全角半球谐振陀螺的控制回路带宽的方法

Publications (2)

Publication Number Publication Date
CN116430713A true CN116430713A (zh) 2023-07-14
CN116430713B CN116430713B (zh) 2024-01-05

Family

ID=87085141

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310411145.XA Active CN116430713B (zh) 2023-04-18 2023-04-18 一种用于提高全角半球谐振陀螺的控制回路带宽的方法

Country Status (1)

Country Link
CN (1) CN116430713B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105203132A (zh) * 2015-10-21 2015-12-30 中国矿业大学(北京) 一种谐振式振动陀螺的输出频率检测方法
CN114955011A (zh) * 2022-04-12 2022-08-30 北京航空航天大学 一种dgvscmg飞轮模式下框架系统定角控制方法
CN115235444A (zh) * 2022-07-19 2022-10-25 青岛哈尔滨工程大学创新发展中心 一种用于测量全角半球谐振陀螺的控制回路带宽的方法
CN115574798A (zh) * 2022-09-30 2023-01-06 中国人民解放军火箭军工程大学 一种提高半球谐振陀螺综合性能的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105203132A (zh) * 2015-10-21 2015-12-30 中国矿业大学(北京) 一种谐振式振动陀螺的输出频率检测方法
CN114955011A (zh) * 2022-04-12 2022-08-30 北京航空航天大学 一种dgvscmg飞轮模式下框架系统定角控制方法
CN115235444A (zh) * 2022-07-19 2022-10-25 青岛哈尔滨工程大学创新发展中心 一种用于测量全角半球谐振陀螺的控制回路带宽的方法
CN115574798A (zh) * 2022-09-30 2023-01-06 中国人民解放军火箭军工程大学 一种提高半球谐振陀螺综合性能的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
XUKAI DING等: "Bias Modulation of Force-to-Rebalanced Micro Hemispherical Resonator Gyroscope Based on Mode-Rotation", 《IEEE》, vol. 22, no. 16, pages 15802 - 15816, XP011916874, DOI: 10.1109/JSEN.2022.3188287 *
南方伯等: "全角半球谐振陀螺控制回路的动态特性研究", 《仪器仪表学报》, vol. 43, no. 10, pages 42 - 49 *

Also Published As

Publication number Publication date
CN116430713B (zh) 2024-01-05

Similar Documents

Publication Publication Date Title
CN111578966B (zh) 一种基于lms算法的半球谐振子特征参数辨识方法
CN108459497A (zh) 一种基于adrc和nlpid的稳瞄伺服系统的控制方法
US9356607B2 (en) Controller unit and device for resetting an oscillator excited by a harmonic oscillation, and yaw rate sensor
CN108897218A (zh) 一种基于混合奇次重复控制器的磁悬浮转子奇次谐波振动抑制方法
CN115235444B (zh) 一种用于测量全角半球谐振陀螺的控制回路带宽的方法
CN106289208B (zh) 一种基于非线性自适应算法的磁轴承系统惯性轴辨识方法
CN114440933B (zh) 一种谐振陀螺仪旋转调制标度自校正系统
CN114440852A (zh) 一种半球谐振陀螺x/y通路驱动增益不对称性误差辨识方法
CN206378139U (zh) 一种正交偏置自消除的微机械陀螺仪
Deng et al. On vibration rejection of nonminimum-phase long-distance laser pointing system with compensatory disturbance observer
CN116430713B (zh) 一种用于提高全角半球谐振陀螺的控制回路带宽的方法
Raspopov et al. Tuning and calibration of a coriolis vibratory gyroscope with a metal resonator to operate in angular rate sensor mode
CN106842953A (zh) 一种无人直升机自适应低阶控制器
EP4339621A1 (en) Mems accelerometer systems
Li et al. Model analysis and resonance suppression of wide-bandwidth inertial reference system
CN114264318A (zh) 一种闭环光纤陀螺固有频率测试方法和装置
Zheng et al. Disturbance rejection in MEMS gyroscope: Problems and solutions
CN106094915A (zh) 能量注入式虚拟质量共振柱控制系统及控制方法
CN116300476B (zh) 基于转速环ladrc控制器的谐振抑制方法
CN117629247A (zh) 硅陀螺校正方法及装置
Tu et al. State estimation of a robot joint by a novel nonlinear tracking differentiator
CN118392216B (zh) 一种微机械陀螺标度因数温度非线性在线校准方法
CN118623857B (zh) 半球谐振陀螺双通道驱动增益不对称误差辨识方法、装置、设备及存储介质
CN116086449B (zh) 基于变结构控制的半球谐振陀螺稳定平台及其构建方法
Jia et al. In-run Self-calibration of Scale Factor Temperature Drifts for MEMS Gyroscope

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