CN116430713A - 一种用于提高全角半球谐振陀螺的控制回路带宽的方法 - Google Patents
一种用于提高全角半球谐振陀螺的控制回路带宽的方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000011426 transformation method Methods 0.000 claims abstract description 7
- 238000012546 transfer Methods 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 7
- 230000003321 amplification Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 238000011217 control strategy Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic 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转换增益的检测增益以及椭圆参数解算增益。
可选的,获取全角半球谐振陀螺控制系统被控对象的方法包括:
K=KdriveKHRGKsenseKpara
其中,G1(s)为全角半球谐振陀螺控制系统被控对象,K为陀螺控制系统总体增益,τ为时间常数,s为s域下的变量,Kdrive为驱动环节增益,KHRG为半球谐振子幅值增益,Ksense为检测环节增益,Kpara为椭圆参数解算环节增益。
可选的,所述控制器为经典PI控制器并联谐振控制组成PIR控制器。
可选的,基于所述控制器获取传递函数的方法包括:
其中,GPIR(s)为传递函数,kp、ki、kr分别为比例、积分、谐振环节增益,ω0为谐振控制器的谐振频率,ωc为谐振控制器截止频率。
可选的,通过双线性变换法对控制器进行离散化处理的方法包括:
其中,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)。
由于半球谐振子具有很高的Q值,其衰减时间常数τ达到百秒级别,属于大惯性环节。经典PI控制器可以实现低速下的有效控制,但半球谐振子高速转动时,PI控制跟踪速度慢,无法实现精准控制,进而导致系统不稳定。陀螺转速引起干扰项会呈现出正弦变化,其变化周期和转速存在固定关系。
陀螺控制系统中的扰动误差主要来源于半球谐振子极板间距不等、电路元件误差等带来的检测回路增益不匹配。因此,将x检测通道与y检测通道的增益系数kx和ky分别加入到解算回路中,图1半球谐振子解算过程中x,y方向的位移表达式变为:
以主波波腹a为例分析检测通道增益对椭圆参数解算的影响,利用图1中解调得到的cx、cy、sx和sy四个参数分别可以得到:
解算到的主波波腹变ac为:
由式(3)可以得出:当kx=ky时,ac=a,检测通道增益对信号解算没有影响;当kx≠ky时,若kx>ky,ac会随着角度的增加呈现正弦变化,而且ac变化频率为角度θ变化频率的四倍。
并引入速率反馈,实现变速情况下精准控制,包括:PIR控制是在传统的PI控制基础上并联谐振(R)控制,实现对正弦信号的无差控制,其传递函数为:
式中,kp、ki、kr分别为比例、积分、谐振环节增益,ω0为谐振控制器的谐振频率,ωc为谐振控制器截止频率。其中,截止频率ωc决定了谐振频率的范围,ωc越大,谐振区域的频率范围越大,一般取值为5-15rad/s,设计时选用5rad/s;谐振频率ω0决定谐振点的频率,与待消除的谐波频率有关,由式(3)得到的结论,谐波信号与陀螺角速率有关,最后谐振频率取角速率的四倍,并且通过速率反馈实时调整。
图3为PI控制和PIR控制的伯德图,从图中可知,加入谐振控制后,将谐振频率附近的增益提高,而不改变控制器其它频段的幅频特性。
将图2中的控制器替换为PIR控制器,则R(s)输入时的控制回路闭环传递函数ΦR(s)可以表示为:
如图4所示,根据式(16)绘制出PIR控制下的全角半球谐振陀螺控制回路伯德图,并与PI控制下的伯德图进行对比。由图4可知,选择合适的PIR参数,可以将特定频率点的幅频特性提高至-3dB以上。
变速情况下,引入速率反馈,其控制系统结构如图5所示,N(s)为陀螺转动引入的干扰项,为陀螺解算出的角速率。PIR控制器的比例控制环节和积分控制环节参数保持不变,谐振控制的截止频率ωc保持为5rad/s,谐振频率ω0为角速率/>的四倍。谐振环节增益kr不是定值,角速率越大,需要的增益越大,考虑到两者不存在线性关系,可以采用分段选择的方式确定kr,先在固定转速下调试不同转速范围内的kr最优值(每0.1Hz确定一个kr),变速情况下通过判断角速率所处范围,设定与其对应的kr,实现实时调整。
引入速率反馈后,实时改变谐振控制的谐振点位置并调整控制参数,提高谐振点的幅值裕度,进而提高系统控制回路带宽。
一种PIR控制器数字离散化的方法,包括:在实际的陀螺仪控制系统中处理的信号为数字信号,需要将PIR控制器s域的传递函数转换到z域,实现离散化。采用双线性变换法实现PIR控制器的离散化,离散公式为:
式中,T为采样时间,与数字控制系统时钟频率有关。
将式(4)通分整理,得:
其中,分子、分母各系数为:
将离散域下的传函转化为差分方程。
先将式(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转换增益的检测增益以及椭圆参数解算增益。
4.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,所述控制器为经典PI控制器并联谐振控制组成PIR控制器。
7.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,所述全角半球谐振陀螺控制系统包括:幅值控制、正交控制和频率控制。
8.如权利要求1所述的用于提高全角半球谐振陀螺的控制回路带宽的方法,其特征在于,对离散化处理后的所述控制器的若干参数进行标定,引入速率反馈提高全角半球谐振陀螺闭环控制回路的带宽的方法包括:
所述控制器的系数随着外部输入转速变化而变化,谐振控制中心频率始终跟踪外部转速的四倍频,有效抑制由转速引起的扰动,提高控制回路带宽。
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)
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 | 中国人民解放军火箭军工程大学 | 一种提高半球谐振陀螺综合性能的方法 |
-
2023
- 2023-04-18 CN CN202310411145.XA patent/CN116430713B/zh active Active
Patent Citations (4)
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)
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 |