CN111578966A - 一种基于lms算法的半球谐振子特征参数辨识方法 - Google Patents
一种基于lms算法的半球谐振子特征参数辨识方法 Download PDFInfo
- Publication number
- CN111578966A CN111578966A CN202010273032.4A CN202010273032A CN111578966A CN 111578966 A CN111578966 A CN 111578966A CN 202010273032 A CN202010273032 A CN 202010273032A CN 111578966 A CN111578966 A CN 111578966A
- Authority
- CN
- China
- Prior art keywords
- lms algorithm
- parameter identification
- harmonic oscillator
- signal
- following
- 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
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Geometry (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Manufacturing & Machinery (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Gyroscopes (AREA)
Abstract
本发明提供一种基于LMS算法的半球谐振子特征参数辨识方法,包括:采集振动陀螺的检测信号进行数据处理,包括将检测信号与参考信号相乘,解算得到椭圆坐标系下的参数方程。通过幅度控制、正交控制和频率控制,使参数逐渐收敛至稳定状态。根据解算得到的椭圆坐标系下的参数方程,构建LMS算法参数辨识模型。根据LMS算法参数辨识模型的输出,求解得到待求的特征参数,实现特征参数辨识功能。本发明具有简洁性,只需要将半球谐振陀螺接入驱动检测装置中,其余解算均由处理器完成。该方法具有适用性,驱动检测装置可以适用于绝大多数半球谐振子的测定。该方法具有实时性和准确性,可以在陀螺装配完成后进行测试,且测试结果准确,可以直接使用。
Description
技术领域
本发明涉及一种基于LMS算法的半球谐振子特征参数辨识方法,属于智能化仪器仪表领域。
背景技术
由于实际工艺水平的限制,制备得到的非理想半球谐振子往往会存在许多误差。因此,为了准确评估半球谐振子的性能指标,判断其是否达到合格标准,就需要对以下半球谐振子特征参数进行测试:刚度轴各向异性Δω;刚度轴失准角θω;阻尼轴各向异性Δ(1/τ);阻尼轴失准角θτ。
目前,半球谐振子特征参数的辨识,大多是通过精确测量半球谐振子来得到较为准确的特征参数值。这种方法操作和运算较为复杂,且需要十分精密的测量仪器。为此,本专利提出一种更为简洁且不需要精密测量仪器的辨识方法。
发明内容
本发明的目的是为了克服参数辨识困难的问题而提出的一种基于LMS算法的半球谐振子特征参数辨识方法,用以获取精确的特征参数:Δω、Δ(1/τ)、θω、θτ,为半球谐振子的性能评估提供一种便捷有效的方法。
本发明的目的是这样实现的:步骤如下:
步骤1:非理想半球谐振子在椭圆轨道坐标系下的运动微分方程为:
式中:a的数值反映了陀螺的振动幅度大小;q的数值反映了陀螺偏离理想简正模态的程度;θ为角度信号;是振动信号和参考信号之间出现微小相差;Δω是刚度轴各向异性;θω是刚度轴失准角;Δ(1/τ)是阻尼轴各向异性;θτ是阻尼轴失准角;
步骤2:全角模式时,对陀螺施加幅度控制、正交控制、频率控制和相位控制;
步骤3:构建LMS算法滤波器,
步骤5:求解刚度轴各向异性Δω;
步骤6:求解刚度轴失准角θω。
本发明还包括这样一些结构特征:
1.步骤2具体是:
由幅度控制,维持a=a_setting;
由正交控制,维持q的常数项为0但残余正弦变化的部分,且q<<a;因此,可以得到下列公式:
展开其中的第二、第三项公式得到:
2.步骤3具体包括:
(1)对于原始信号经过延迟处理后得到对应不同延迟的输入信号x(k);
(2)对权向量信号w(k)进行初始化设置;
(3)根据(1)中得到x(k)与其对应权向量信号w(k)的转置wT(k)相乘,得到该时刻的输出信号y(k);
(4)将参考信号d(k)与根据(3)中得到该时刻的输出信号y(k)作差,得到误差信号e(k);
(5)权向量信号更新方法为:
w(k+1)=w(k)+μe(k)x(k)
(6)整体算法运算过程为:
则有:
与现有技术相比,本发明的有益效果是:与现有技术相比,本发明的有益效果是:现有技术需要使用精密仪器对半球谐振子实物进行测量,通过一系列计算后得出半球谐振子的特征参数,操作复杂,且在谐振子装配后不可测量。而本发明则具有以下优点:首先,该方法具有简洁性,只需要将半球谐振陀螺接入驱动检测装置中,其余解算均由处理器完成。其次,该方法具有适用性,驱动检测装置可以适用于绝大多数半球谐振子的测定。最后,该方法具有实时性和准确性,可以在陀螺装配完成后进行测试,且测试结果准确,可以直接使用。
此外,由于本专利中所涉及的方法需要将半球谐振子的振型旋转起来,旋转越快则辨识速率越快。因此,我们可以将振动陀螺固联在转台上,借助转台的快速旋转来提高特征参数的辨识速率。
附图说明
图1是基于LMS算法的半球谐振子特征参数辨识方法流程图。
图2是LMS算法滤波器结构图。
图3是基于LMS算法的半球谐振子特征参数辨识方法具体实施的结构框图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
本发明的步骤为:
步骤1:非理想半球谐振子在椭圆轨道坐标系下按照以下微分方程进行运动:
步骤2:全角模式时,对陀螺施加幅度控制、正交控制、频率控制和相位控制。
由幅度控制,维持a=a_setting;
由正交控制,维持q的常数项为0但残余正弦变化的部分,且q<<a;因此,可以得到下列公式:
展开其中的第二、第三项公式得到:
步骤3:构建LMS算法滤波器,包括以下步骤:
(1)对于原始信号经过延迟处理后得到对应不同延迟的输入信号x(k);
(2)对权向量信号w(k)进行初始化设置;
(3)根据(1)中得到x(k)与其对应权向量信号w(k)的转置wT(k)相乘,得到该时刻的输出信号y(k);
(4)将参考信号d(k)与根据(3)中得到该时刻的输出信号y(k)作差,得到误差信号e(k);
(5)权向量信号更新方法如下:
w(k+1)=w(k)+μe(k)x(k)
(6)整体算法运算过程如下:
步骤5:求解刚度轴各向异性Δω。
由于LMS算法滤波器的特性可知:
步骤6:求解刚度轴失准角θω。
结合附图对本发明的进行详细的描述:
如图1所述,本发明提供了一种基于LMS算法的半球谐振子特征参数辨识方法,图1是基于LMS算法的半球谐振子特征参数辨识方法流程图,包括以下步骤:
步骤1:采集振动陀螺的检测信号进行数据处理,包括将检测信号与参考信号相乘,解算得到椭圆坐标系下的参数方程。通过幅度控制、正交控制和频率控制,使参数逐渐收敛至稳定状态。
步骤2:根据步骤1中解算得到的椭圆坐标系下的参数方程,构建LMS算法参数辨识模型。
步骤3:根据LMS算法参数辨识模型的输出,求解得到待求的特征参数,实现特征参数辨识功能。
步骤1包括:
步骤1-1:非理想半球谐振子在椭圆轨道坐标系下按照以下微分方程进行运动:
步骤1-2:全角模式时,对陀螺施加幅度控制、正交控制、频率控制和相位控制。
由幅度控制,维持a=a_setting;由正交控制,维持q的常数项为0但残余正弦变化的部分,且q<<a;可以得到下列公式:
展开其中的第二、第三项公式得到:
步骤3:将步骤2中LMS算法参数辨识模型的输出进行解算,求解得到待求的特征参数Δω和θω,实现特征参数辨识功能。
步骤3-1:求解刚度轴非对称度Δω。
由于LMS算法滤波器的特性可知:
步骤3-2:求解刚度轴失准角θω。
下面以对于半球谐振子特征参数Δ(1/τ)和θτ的求解,验证本发明所提出的一种基于LMS算法的半球谐振子特征参数辨识方法的正确性和有效性。
此后计算方法均与前文一致,运算结果为:
w5(k)的输出为:-γΩ;
图3是基于LMS算法的半球谐振子特征参数辨识方法具体实施的结构框图,如图3所示:经过驱动、检测、信号处理合成以及解算,最终可以得到非理想半球谐振子在椭圆轨道坐标系下的运动微分方程,并将其作为LMS算法参数辨识模型的输入信号。LMS算法参数辨识模型的输出经过解算得到了半球谐振子特征参数Δω、θω、Δ(1/τ)和θτ,证实了本方法的正确性和有效性。
本发明提供了一种基于LMS算法的半球谐振子特征参数辨识方法,具体实现该技术方案的方法和途径很多,以上所述即是本发明的优选实施方案,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (6)
1.一种基于LMS算法的半球谐振子特征参数辨识方法,其特征在于:步骤如下:
步骤1:非理想半球谐振子在椭圆轨道坐标系下的运动微分方程为:
式中:a的数值反映了陀螺的振动幅度大小;q的数值反映了陀螺偏离理想简正模态的程度;θ为角度信号;是振动信号和参考信号之间出现微小相差;Δω是刚度轴各向异性;θω是刚度轴失准角;Δ(1/τ)是阻尼轴各向异性;θτ是阻尼轴失准角;
步骤2:全角模式时,对陀螺施加幅度控制、正交控制、频率控制和相位控制;
步骤3:构建LMS算法滤波器,
步骤5:求解刚度轴各向异性Δω;
步骤6:求解刚度轴失准角θω。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010273032.4A CN111578966B (zh) | 2020-04-09 | 2020-04-09 | 一种基于lms算法的半球谐振子特征参数辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010273032.4A CN111578966B (zh) | 2020-04-09 | 2020-04-09 | 一种基于lms算法的半球谐振子特征参数辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111578966A true CN111578966A (zh) | 2020-08-25 |
CN111578966B CN111578966B (zh) | 2021-07-16 |
Family
ID=72089910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010273032.4A Active CN111578966B (zh) | 2020-04-09 | 2020-04-09 | 一种基于lms算法的半球谐振子特征参数辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111578966B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112506040A (zh) * | 2020-11-16 | 2021-03-16 | 哈尔滨工程大学 | 一种全角半球谐振陀螺单通道控制系统及方法 |
CN112595303A (zh) * | 2020-12-10 | 2021-04-02 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种hrg全角模式下参数计算方法 |
CN115031713A (zh) * | 2022-04-07 | 2022-09-09 | 哈尔滨工程大学 | 一种用于自校准半球谐振陀螺检测信号非线性的方法 |
CN115077561A (zh) * | 2022-06-15 | 2022-09-20 | 青岛哈尔滨工程大学创新发展中心 | 一种自适应补偿半球谐振子阻尼各向异性的方法及系统 |
CN117387589A (zh) * | 2023-12-11 | 2024-01-12 | 四川图林科技有限责任公司 | 一种半球谐振陀螺的谐振子检测电极输出信号滤波方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1164356A1 (en) * | 2000-06-13 | 2001-12-19 | Ford Global Technologies, Inc. | System and method for estimating sensor errors |
CN107631738A (zh) * | 2017-09-20 | 2018-01-26 | 哈尔滨工业大学 | 一种陀螺飞轮偶不平衡辨识方法 |
CN109084741A (zh) * | 2018-07-04 | 2018-12-25 | 中国人民解放军国防科技大学 | 一种圆柱壳体振动陀螺谐振子频率裂解方法及系统 |
CN109945894A (zh) * | 2019-04-08 | 2019-06-28 | 哈尔滨工业大学 | 半球谐振陀螺仪谐振子与外基座的安装误差补偿方法 |
CN110260850A (zh) * | 2019-06-14 | 2019-09-20 | 哈尔滨工程大学 | 一种半球谐振子的高压稳定偏置电路 |
CN110968961A (zh) * | 2019-12-17 | 2020-04-07 | 哈尔滨理工大学 | 一种连续回转电液伺服马达参数辨识方法 |
-
2020
- 2020-04-09 CN CN202010273032.4A patent/CN111578966B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1164356A1 (en) * | 2000-06-13 | 2001-12-19 | Ford Global Technologies, Inc. | System and method for estimating sensor errors |
CN107631738A (zh) * | 2017-09-20 | 2018-01-26 | 哈尔滨工业大学 | 一种陀螺飞轮偶不平衡辨识方法 |
CN109084741A (zh) * | 2018-07-04 | 2018-12-25 | 中国人民解放军国防科技大学 | 一种圆柱壳体振动陀螺谐振子频率裂解方法及系统 |
CN109945894A (zh) * | 2019-04-08 | 2019-06-28 | 哈尔滨工业大学 | 半球谐振陀螺仪谐振子与外基座的安装误差补偿方法 |
CN110260850A (zh) * | 2019-06-14 | 2019-09-20 | 哈尔滨工程大学 | 一种半球谐振子的高压稳定偏置电路 |
CN110968961A (zh) * | 2019-12-17 | 2020-04-07 | 哈尔滨理工大学 | 一种连续回转电液伺服马达参数辨识方法 |
Non-Patent Citations (5)
Title |
---|
BOUJO E.ETC: "Robust identification of harmonic oscillator parameters using the adjoint Fokker-Planck equation", 《PROCEEDINGS OF THE ROYAL SOCIETY A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES》 * |
JI YANG,ETC: "Scale factor error compensation of FOG based on fractional tap-length LMS method", 《PROCEEDINGS OF THE 2016 IEEE 11TH CONFERENCE ON INDUSTRIAL ELECTRONICS AND APPLICATIONS (ICIEA)》 * |
叶中付: "《统计信息处理》", 31 January 2009 * |
李巍等: "基于D-最优的半球谐振陀螺仪误差参数辨识研究", 《哈尔滨理工大学学报》 * |
王旭: "半球谐振陀螺误差建模补偿与力平衡控制方法研究", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112506040A (zh) * | 2020-11-16 | 2021-03-16 | 哈尔滨工程大学 | 一种全角半球谐振陀螺单通道控制系统及方法 |
CN112506040B (zh) * | 2020-11-16 | 2024-01-30 | 哈尔滨工程大学 | 一种全角半球谐振陀螺单通道控制系统及方法 |
CN112595303A (zh) * | 2020-12-10 | 2021-04-02 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种hrg全角模式下参数计算方法 |
CN112595303B (zh) * | 2020-12-10 | 2022-09-23 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种hrg全角模式下参数计算方法 |
CN115031713A (zh) * | 2022-04-07 | 2022-09-09 | 哈尔滨工程大学 | 一种用于自校准半球谐振陀螺检测信号非线性的方法 |
CN115031713B (zh) * | 2022-04-07 | 2023-02-10 | 哈尔滨工程大学 | 一种用于自校准半球谐振陀螺检测信号非线性的方法 |
CN115077561A (zh) * | 2022-06-15 | 2022-09-20 | 青岛哈尔滨工程大学创新发展中心 | 一种自适应补偿半球谐振子阻尼各向异性的方法及系统 |
CN115077561B (zh) * | 2022-06-15 | 2023-03-10 | 青岛哈尔滨工程大学创新发展中心 | 一种自适应补偿半球谐振子阻尼各向异性的方法及系统 |
CN117387589A (zh) * | 2023-12-11 | 2024-01-12 | 四川图林科技有限责任公司 | 一种半球谐振陀螺的谐振子检测电极输出信号滤波方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111578966B (zh) | 2021-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111578966B (zh) | 一种基于lms算法的半球谐振子特征参数辨识方法 | |
CN110346720B (zh) | 一种马达非线性参数的测试方法及装置 | |
CN113252019B (zh) | 一种前向放大系数不一致时的半球谐振陀螺振型角获取方法 | |
CN109084741B (zh) | 一种圆柱壳体振动陀螺谐振子频率裂解方法及系统 | |
CN104697553B (zh) | 光纤陀螺捷联惯导系统加速度计内杆臂标定方法 | |
CN108489512B (zh) | 一种半球谐振陀螺标度因数的补偿标定方法及装置 | |
CN114440852B (zh) | 一种半球谐振陀螺x/y通路驱动增益不对称性误差辨识方法 | |
CN113155114A (zh) | Mems惯性测量单元陀螺零位的温度补偿方法及装置 | |
CN106932125B (zh) | 一种硅谐振压力传感器的补偿方法 | |
CN115127533B (zh) | 一种基于谐振陀螺振型旋转调制下的误差拟合补偿方法 | |
CN103808349A (zh) | 矢量传感器的误差校正方法和装置 | |
CN111089576A (zh) | 光纤陀螺实际输出值确定方法、光纤陀螺阀值测试方法 | |
CN111272193A (zh) | 一种基于噪声功率谱估计的mems陀螺仪在线频差识别方法 | |
CN114858184A (zh) | 一种半球谐振子参数辨识方法 | |
JP2012508867A (ja) | コリオリジャイロの動作状態のシミュレーション方法 | |
CN104457792B (zh) | 一种在无机械转动条件下测量光纤陀螺标度因数的方法 | |
CN106199270A (zh) | 一种石英晶体谐振器等效电路参数的测量方法 | |
CN116124279A (zh) | 一种确定谐振频率的方法、装置、设备及存储介质 | |
CN114264318A (zh) | 一种闭环光纤陀螺固有频率测试方法和装置 | |
CN112097895B (zh) | 一种传感器频响的拓宽方法 | |
Tang et al. | An improved scale factor calibration model of MEMS gyroscopes | |
CN103196450B (zh) | 基于模拟电路形式的卡尔曼滤波方法及模拟电路 | |
CN112131768A (zh) | 一种基于模态和频率的谐振式加速度计优化方法 | |
CN114088077B (zh) | 一种改进的半球谐振陀螺信号去噪方法 | |
CN109579973A (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 |