CN104729531B - 一种mems陀螺静电平衡电压值的确定方法 - Google Patents
一种mems陀螺静电平衡电压值的确定方法 Download PDFInfo
- Publication number
- CN104729531B CN104729531B CN201410833945.1A CN201410833945A CN104729531B CN 104729531 B CN104729531 B CN 104729531B CN 201410833945 A CN201410833945 A CN 201410833945A CN 104729531 B CN104729531 B CN 104729531B
- Authority
- CN
- China
- Prior art keywords
- mode
- input
- data
- value
- group
- 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
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
Abstract
一种MEMS陀螺静电平衡电压值的确定方法,改进了常规的穷举算法来确定调整电路的输入值,先以m步长搜索,初选出符合给定判断依据范围的值,缩小搜索范围,然后以一步长搜索,最终找到理想的调整电压值;其判断依据为两个模态频率相近、两模态间的耦合值小和零偏值小,通过简单的软件方式实现了MEMS陀螺静电平衡电压值的确定,搜索效率高,操作简单,最大程度上满足了环形MEMS陀螺电路非理想结构调整的需求。
Description
技术领域
本发明涉及一种静电平衡电压值的确定方法,特别是一种MEMS陀螺静电平衡电压值的确定方法,适用于环形MEMS陀螺电路非理想结构的调整。
背景技术
近年来,随着MEMS技术的发展,微机械陀螺仪以其体积小、功耗低、精度高等优点越来越受到人们的关注,在民用消费领域、航空航天以及现代国防领域具有广泛的应用前景。
然而,在微机械加工过程中,误差的存在不可避免。加工误差会引起结构的非理想性(主要包括质量和弹性的非理想性),造成两个模态的频率不相等,进而导致信噪比降低,并产生较大的零点漂移,影响精度无法满足系统测试要求。因此,研究静电平衡调整方法具有重要的实际意义。传统调整方法有激光平衡法,即通过对硅谐振环上不均匀处进行激光打平和粘附来修复这种非理想性,但该方法成本高,最小单位材料的消除和粘附都会引起谐振频率很大程度的改变,所以调整精度不是很高,并且容易造成不可恢复的损坏。
目前主要是增加调整电极,通过静电力的方法改变谐振环的有效弹性系数来弥补结构上的非理想性。虽然已有方法给出了用静电力调整陀螺仪的不等弹力,也给出了较详尽的分析调整过程,通过频率响应函数的奈奎斯特图来求得结构非理想参数,并给出了电压的表达式,然而该方法实验环节较薄弱,理论指导公式只适合调整电极特定的排列顺序,同时对测试设备要求较高,要用到网络动态信号分析仪来测试系统的奈奎斯特图,成本较高,并且该推导公式适用于表头,未必适用于将表头、ASIC以及单片机为一体的陀螺电路,使用范围较窄。
发明内容
本发明的技术解决问题是:克服现有方法的不足之处,提供了一种MEMS(MicroElectronic Mechanical System,微电子机械系统)陀螺静电平衡电压值的确定方法,首先采用m步长搜索平衡电压的潜在解,最后通过1步长搜索确定最终解,通过简单的软件方式实现了MEMS陀螺静电平衡电压值的确定,搜索效率高,操作简单,最大程度上满足了环形MEMS陀螺电路非理想结构调整的需求。
本发明的技术解决方案是:一种MEMS陀螺静电平衡电压值的确定方法,所述确定方法由位于MEMS陀螺振动环外侧的平衡调整电路实现,所述平衡调整电路包括第一平衡调整电路和第二平衡调整电路,第一平衡调整电路包括第一输入端、第一输出端和第二输出端,第二平衡调整电路包括第二输入端、第三输出端和第四输出端,步骤如下:
(1)根据第一输入端输入值b取值范围和第二输入端输入值b′取值范围确定平衡调整电路的潜在解集S(b,b′),一组b和b′的值为一个潜在解,所述b和b′均为8位二进制数;
(2)第一输入端输入值b和第二输入端输入值b′均以m个A/D最小分辨率,即m个LSB为步长,在S(b,b′)中进行取值,获得m步长筛选后的潜在解集为S1(b,b′),S1(b,b′)中第一输入端输入值b和第二输入端输入值b′的取值均为m的整数倍;
(3)将步骤(2)中确定的潜在解集S1(b,b′)中第一个值b1和b′1分别作用于平衡调整电路的第一输入端和第二输入端,首先使第一模态实现恒幅谐振,记录第一模态的固有谐振频率值L11、第一模态对第二模态的耦合量D11和第二模态的输出零偏值E11,然后使第二模态实现恒幅谐振,记录第二模态固有谐振频率值L21、第二模态对第一模态的耦合量D21和第一模态的输出零偏值E21;获得一组数据P1=(b1,b′1,L11,D11,E11,L21,D21,E21);
(4)按顺序取出潜在解集S1(b,b′)中的第一输入端输入值b和第二输入端输入值b′输入值,重复步骤(3),取得j组数据,其中第i组数据记作Pi=(bi,b′i,L1i,D1i,E1i,L2i,D2i,E2i);所述j为S1(b,b′)中潜在解的个数,i∈[1,j];;(5)计算步骤(4)中确定的j组数据中每组数据的频差L1i-L2i,利用每组数据的频差对该组数据进行筛选,获得频差筛选后的数据集合;具体为:若该组数据的频差处于预先设定的频差范围内,则保留该组数据,否则删除该组数据;
(6)在频差筛选后的数据集合中,按照预先设定的筛选条件选择两组数据,筛选获得的第一组数据中第一输入端的输入值为b11,第二输入端的输入值为b′11,第二组数据中第一输入端的输入值为b22,第二输入端的输入值为b′22,所述|b11-b22|<=16;|b′11-b′22|<=16;
(7)利用步骤(6)中得到的两组数据输入端的输入值计算1步长筛选第一输入端输入值取值区间[bb1,bb2]和第二输入端输入值取值区间[bb′1,bb′2];
(8)第一输入端输入值b和第二输入端输入值b′分别以1个A/D最小分辨率为步长在[bb1,bb2]和[[bb′1,bb′2]中取值,获得1步长筛选后的潜在解集为S2(b,b′);
(9)用步骤(8)获得的1步长筛选后的潜在解集S2(b,b′)更新步骤(2)中的S1(b,b′),重复步骤(3)~步骤(4),选择第一模态对第二模态的耦合量D1i、第二模态的输出零偏值E1i、第二模态对第一模态的耦合量D2i和第一模态的输出零偏值E2i绝对值的和最小的一组数据为最终数据,所述最终数据为:Pbest=(bbest,b′best,L1best,D1best,E1best,L2best,D2best,E2best);
(10)利用步骤(9)确定的bbest和b′best,计算出四组电极加的电压值大小,具体为:
若bbest的最高位为1,则V1为28V,V3可调,若bbest的最高位为0,则V3为28V,V1可调;所述调整电压的大小由bbest的低七位决定,转换公式为:电压值=低七位的十进制数/127×28;
若b′best的最高位为1,则V2为28V,V4可调,若b′best的最高位为0,V4为28V,V2可调,所述调整电压的大小由b′best的低七位决定,转换公式为:电压值=低七位的十进制数/127×28。
所述步骤(2)中m的取值范围为:8~16。
所述步骤(3)中L11,D11,E11,L21,D21,E21均通过n次试验求平均获得。
所述n的取值范围为:100~200。
所述步骤(5)中预先设定频差范围为[-0.5,0.5]。
所述步骤(6)中预先设定的筛选条件具体为:
所选择的两组数据中第一模态对第二模态的耦合量D1i、第二模态的输出零偏值E1i、第二模态对第一模态的耦合量D2i和第一模态的输出零偏值E2i符号均相反,且两组数据中第一模态对第二模态的耦合量D1i和第二模态对第一模态的耦合量D2i绝对值的和最小。
所述步骤(7)中利用步骤(6)中得到的两组数据输入端的输入值计算1步长筛选第一输入端输入值取值区间[bb1,bb2]和第二输入端输入值取值区间[bb′1,bb′2];具体为:
若|b11-b22|=16,则bb1=min{b11,b22},bb2=max{b11,b22};否则,bb1=(b11+b22)/2-8,bb2=(b11+b22)/2+8;
若|b′11-b′22|=16,则bb′1=min{b′11,b′22},bb′2=max{b′11,b′22};否则,bb′1=(b′11+b′22)/2-8,bb′2=(b′11+b′22)/2+8。
本发明与现有技术相比的有益效果是:
(1)本发明先以m步长搜索缩小潜在解范围,然后以1步长搜索确定最后的电压数值,传统搜索算法潜在解有N=65536种可能,本发明中的算法有N′=65536/m/m+289种可能,效率提高了将近m2倍;
(2)本发明采用两个模态频率相近,且两模态间的耦合值小以及零偏输出值小来进行最优数据的选取,这是从陀螺结构非理想性产生的影响得出的判断依据,选取的依据更加合理,选取结果更接近于真实的结果;
(3)本发明对测试环境以及测试仪器的依赖性较小,包括数据的采集、模态的切换等操作均是通过设计相关软件实现的,对硬件设备的要求低,操作效率高。已有方法要用到网络动态信号分析仪,成本较高。
附图说明
图1为本发明中陀螺表头的结构图;
图2为MEMS陀螺谐振环振动方式示意图,其中(a)为谐振环第一模态振动方式示意图,(b)为谐振环第二模态的振动方式;
图3为本发明的系统框图;
图4为本发明的流程图。
具体实施方式
如图1所示为本发明中陀螺表头的结构图,从图1可知,MEMS陀螺振动环的外侧均匀分布八个电极,顺时针依次为电极1~电极8,其中电极2和电极6为第一模态的驱动电极,电极4和电极8为第一模态的检测电极,电极3和电极7为第二模态的驱动电极,电极1和电极5为第二模态的检测电极;图2为MEMS陀螺谐振环振动方式示意图,其中(a)为谐振环第一模态振动方式示意图,(b)为谐振环第二模态的振动方式,第二模态在偏离第一模态45°角方向上振动。
振动环的内圈均匀分布有16个电极实现静电平衡调整,顺时针依次为电极9~电极24,所述16个电极分为4组,第一组包括电极9、电极13、电极17和电极21,电压值记住V1,第二组包括电极10、电极14、电极18和电极22,电压值记住V2,第三组包括电极11、电极15、电极19和电极23,电压值记住V3,第四组包括电极12、电极16、电极20和电极24,电压值记住V4。
所述每一组中4个电极所加电压大小相等,电压在0~28V范围内取值,所述第一组电极和第三组电极的电压由第一平衡调整电路控制,对应第一输入端b,第二和第四组由第二平衡调整电路控制,对应第二输入端b′,第一输入端输入值b和第二输入端输入值b′均为8位二进制数据,其取值范围为0~255;
第一平衡调整电路和第二平衡调整电路中的两个输出端电压只有一组电压可调,由第一平衡调整电路和第二平衡调整电路输入端电压值的高位控制,若第一输入端输入值b的最高位为1,V1为28V(等于陀螺中振动环的偏置电压),V3可调,若第一输入端输入值b的最高位为0,V3为28V,V1可调;若第二输入端输入值b′的最高位为1,V2为28V,V4可调,若第二输入端输入值b′的最高位为0,V4为28V,V2可调,调整电压的大小由低七位决定,取值为0~127,其中127对应28V,0对应0V,转换公式为:电压值=低七位的十进制数/127×28。
图3为本发明的系统框图,第一模态的检测电极4,8的检测信号经ASIC电路处理以后得到信号A和B,第二模态的检测电极(1,5)的检测信号经ASIC电路处理后得到信号D和E。微处理器MCU接收到ASIC发送来的数据,若实现第一模态的恒幅谐振,A反映第一模态振幅大小的信息,B反映第一模态的频率偏离谐振频率程度的信息,MCU对信号A进行AGC自动增益环路数字控制器的处理,得到驱动信号A′,对B进行锁相环路数字控制器的处理,得到锁相信息B′,然后将A′和B′发送给ASIC,经ASIC处理后作用于第一模态的驱动电极2,6实现第一模态的恒幅谐振,为了减小第一模态恒幅谐振的过程以及随机噪声对测量结果的影响,本发明采用100次数据求平均得到第一模态的信息(L1,D1,E1),其中第一模态的频率L1通过上位机程序控制频率计测量得到,D反映第一模态对第二模态的耦合值的大小D1,E反映第二模态的零偏值的大小E1。通过软件切换实现第二模态的恒幅谐振,D反映第二模态振幅大小的信息,E反映第二模态的频率偏离谐振频率程度的信息,MCU对信号D进行AGC自动增益环路数字控制器的处理,得到驱动信号D′,对E进行锁相环路数字控制器的处理,得到锁相信息E′,然后将D′和E′发送给ASIC,经ASIC处理后作用于第二模态的驱动电极(3,7),为了减小第二模态恒幅谐振的过程以及随机噪声对测量结果的影响,本发明采用100次数据求平均得到第二模态的信息(L2,D2,E2),其中第二模态的频率L2通过上位机程序控制频率计测量得到,A反映第二模态对第一模态的耦合值的大小D2,B反映第一模态的零偏值的大小E2。如图4所示为本发明的流程图,从图4可知,本发明提出的一种MEMS陀螺静电平衡电压值的确定方法,所述确定方法由位于MEMS陀螺振动环外侧的平衡调整电路实现,所述平衡调整电路包括第一平衡调整电路和第二平衡调整电路,第一平衡调整电路包括第一输入端、第一输出端和第二输出端,第二平衡调整电路包括第二输入端、第三输出端和第四输出端,具体步骤如下:
(1)根据第一输入端输入值b取值范围和第二输入端输入值b′取值范围确定平衡调整电路的潜在解集S(b,b′),一组b和b′的值为一个潜在解,所述b和b′均为8位二进制数,b的取值范围为0~255,共256种取值,b′的取值范围为0~255,共256种取值,这样传统穷举算法潜在解共有N=256×256=65536种取值;
(2)第一输入端输入值b和第二输入端输入值b′均以m个A/D最小分辨率,即m个LSB为步长,在S(b,b′)中进行取值,获得m步长筛选后的潜在解集为S1(b,b′),S1(b,b′)中第一输入端输入值b和第二输入端输入值b′的取值均为m的整数倍;所述m的取值范围为:8~16;
(3)将步骤(2)中确定的潜在解集S1(b,b′)中第一个值b1和b′1分别作用于平衡调整电路的第一输入端和第二输入端,首先将代表驱动力大小的驱动电压作用于第一模态的驱动电极(2,6),使第一模态实现恒幅谐振,其驱动的幅度和频率由ASIC芯片和MCU闭环控制,记录第一模态的固有谐振频率值L11、第一模态对第二模态的耦合量D11和第二模态的输出零偏值E11,然后在相同平衡调整电路输入值的情况下,将代表驱动力大小的驱动电压作用于第二模态的驱动电极(3,7),使第二模态实现恒幅谐振,记录第二模态固有谐振频率值L21、第二模态对第一模态的耦合量D21和第一模态的输出零偏值E21;获得一组数据P1=(b1,b′1,L11,D11,E11,L21,D21,E21);所述L11,D11,E11,L21,D21,E21均通过n次试验求平均获得,n的取值范围为:100~200;
(4)按顺序取出潜在解集S1(b,b′)中的第一输入端输入值b和第二输入端输入值b′输入值,重复步骤(3),取得j组数据,其中第i组数据记作Pi=(bi,b′i,L1i,D1i,E1i,L2i,D2i,E2i);所述j为S1(b,b′)中潜在解的个数,i∈[1,j];
(5)计算步骤(4)中确定的j组数据中每组数据的频差L1i-L2i,利用每组数据的频差对该组数据进行筛选,获得频差筛选后的数据集合;具体为:若该组数据的频差处于预先设定的频差范围内,则保留该组数据,否则删除该组数据;所述预先设定频差范围为[-0.5,0.5];
(6)在频差筛选后的数据集合中,按照预先设定的筛选条件选择两组数据,筛选获得的第一组数据中第一输入端的输入值为b11,第二输入端的输入值为b′11,第二组数据中第一输入端的输入值为b22,第二输入端的输入值为b′22,所述|b11-b22|<=16;|b′11-b′22|<=16;
所述预先设定的筛选条件具体为:
所选择的两组数据中第一模态对第二模态的耦合量D1i、第二模态的输出零偏值E1i、第二模态对第一模态的耦合量D2i和第一模态的输出零偏值E2i符号均相反,且两组数据中第一模态对第二模态的耦合量D1i和第二模态对第一模态的耦合量D2i绝对值的和最小;
(7)利用步骤(6)中计算得到的两组数据输入端的输入值确定1步长筛选第一输入端输入值取值区间[bb1,bb2]和第二输入端输入值取值区间[bb′1,bb′2];具体为:
若|b11-b22|=16,则bb1=min{b11,b22},bb2=max{b11,b22};否则,bb1=(b11+b22)/2-8,bb2=(b11+b22)/2+8;
若|b′11-b′22|=16,则bb′1=min{b′11,b′22},bb′2=max{b′11,b′22};否则,bb′1=(b′11+b′22)/2-8,bb′2=(b′11+b′22)/2+8;
(8)第一输入端输入值b和第二输入端输入值b′分别以1个A/D最小分辨率为步长在[bb1,bb2]和[[bb′1,bb′2]中取值,获得1步长筛选后的潜在解集为S2(b,b′);
(9)用步骤(8)获得的1步长筛选后的潜在解集S2(b,b′)更新步骤(2)中的S1(b,b′),重复步骤(3)~步骤(4),选择第一模态对第二模态的耦合量D1i、第二模态的输出零偏值E1i、第二模态对第一模态的耦合量D2i和第一模态的输出零偏值E2i绝对值的和最小的一组数据为最终数据,所述最终数据为:Pbest=(bbest,b′best,L1best,D1best,E1best,L2best,D2best,E2best);
(10)利用步骤(9)确定的bbest和b′best,计算出四组电极加的电压值大小,具体为:
若bbest的最高位为1,则V1为28V,V3可调,若bbest的最高位为0,则V3为28V,V1可调;所述调整电压的大小由bbest的低七位决定,转换公式为:电压值=低七位的十进制数/127×28;
若b′best的最高位为1,则V2为28V,V4可调,若b′best的最高位为0,V4为28V,V2可调,所述调整电压的大小由b′best的低七位决定,转换公式为:电压值=低七位的十进制数/127×28。
实施例
以某陀螺电路为例,其实现步骤如下:
(1)确定平衡调整电路的潜在解集S(b,b′),其取值范围为((0,0),(0,1),(0,2)…(0,255),(1,0),(1,1)…(1,255),(2,0),(2,1),(2,2)…(255,255)),共有N种可能取值,N=256×256=65536。
(2)确定8步长静电平衡调整电路的潜在解集S1(b,b′),其中b和b′的取值均为8的倍数,其取值范围为((0,0),(0,8),(0,16)…(0,248),(8,0),(8,8),(8,16)…(8,248)…(248,0),(248,8)…(248,248)),共有N1种可能取值,N1=32×32=1024。
(3)依次将S1中的值作用于平衡调整电路,得到数据集合P如表1所示:
表1
b | b′ | L1 | D1 | E1 | L2 | D2 | E2 |
… | … | … | … | … | … | … | … |
192 | 168 | 3474.11 | 10161 | 936 | 3474.06 | 921 | 617 |
192 | 176 | 3473.955 | 10760 | -967 | 3474.844 | -1307 | 354 |
… | … | … | … | … | … | … | … |
200 | 160 | 3474.091 | -19333 | -9439 | 3473.769 | -678 | -196 |
200 | 168 | 3474.128 | 2190 | -55 | 3474.433 | -243 | 434 |
… | … | … | … | … | … | … | … |
(4)计算每组数据的频差值L1i-L2i,按照频差排序,挑选出频差在[-0.5,0.5]范围内的数据集合,如表2所示:
表2
b | b′ | L1 | D1 | E1 | L2 | D2 | E2 | L1-L2 |
… | … | … | … | … | … | … | … | … |
192 | 168 | 3474.11 | 10161 | 936 | 3474.06 | 921 | 617 | 0.0500 |
208 | 160 | 3474.173 | -16708 | -2190 | 3473.863 | -1169 | -468 | 0.3100 |
200 | 160 | 3474.091 | -19333 | -9439 | 3473.769 | -678 | -196 | 0.3220 |
192 | 160 | 3474.028 | 3008 | -7407 | 3473.562 | -193 | -3 | 0.4660 |
(5)在频差筛选后的数据集合中选择两组数据,P1=(192,168,3474.11,10161,936,3474.06,921,617),P2=(200,160,3474.091,-19333,-9439,3473.769,-678,-196)。其中10161和-19333相反,936和-9439符号相反,921和-678符号相反,617和-196符号相反。
(6)确定1步长搜索的范围。其中b的取值区间由P1中的192和P2中的200确定,bb1=(192+200)/2-8=188,bb2=(192+200)/2+8=204,即b的范围为188~204。b′的取值区间由P1中的168和P2中的160确定,其中bb′1=(168+160)/2-8=156,bb′2=(168+160)/2+8=172,即b′的范围为156~172。则一步长的取值范围为S2((188,156)…(188,172),(189,156)…(189,172)…(204,156)…(204,172)),共有N2种可能取值,N2=17×17=289。
(7)将S2中的值作用于平衡调整电路,得到数据集合如表3所示:
表3
b | b′ | L1 | D1 | E1 | L2 | D2 | E2 |
… | … | … | … | … | … | … | … |
190 | 160 | 3474.478 | 2276 | -494 | 3473.961 | -123 | 48 |
190 | 161 | 3474.49 | 458 | 27 | 3474.04 | -2 | 59 |
190 | 162 | 3474.503 | -1721 | 789 | 3474.11 | 131 | 103 |
… | … | … | … | … | … | … | … |
从表3中挑选出耦合值和零偏输出值小的一组数据Pbest=(190,161,3474.49,458,27,3474.04,-2,59,),bbest=190=0b10111110,b′best=161=0b10100001。根据bbest和b′best即可计算出四组电极加的电压值的大小。bbest的最高位为1,则V1为28V,V3=13.6693。b′best的最高位为1,则V2为28V,V4=7.2756。
传统穷举算法潜在解有N=65536种可能,改进穷举算法有N′=N1+N2=1313种可能,在m=8时效率提高了将近50倍。
Claims (7)
1.一种MEMS陀螺静电平衡电压值的确定方法,所述确定方法由位于MEMS陀螺振动环外侧的平衡调整电路实现,所述平衡调整电路包括第一平衡调整电路和第二平衡调整电路,第一平衡调整电路包括第一输入端、第一输出端和第二输出端,第二平衡调整电路包括第二输入端、第三输出端和第四输出端,其特征在于步骤如下:
(1)根据第一输入端输入值b取值范围和第二输入端输入值b′取值范围确定平衡调整电路的潜在解集S(b,b′),一组b和b′的值为一个潜在解,所述b和b′均为8位二进制数;
(2)第一输入端输入值b和第二输入端输入值b′均以m个A/D最小分辨率,即m个LSB为步长,在S(b,b′)中进行取值,获得m步长筛选后的潜在解集为S1(b,b′),S1(b,b′)中第一输入端输入值b和第二输入端输入值b′的取值均为m的整数倍;
(3)将步骤(2)中确定的潜在解集S1(b,b′)中第一个值b1和b′1分别作用于平衡调整电路的第一输入端和第二输入端,首先使第一模态实现恒幅谐振,记录第一模态的固有谐振频率值L11、第一模态对第二模态的耦合量D11和第二模态的输出零偏值E11,然后使第二模态实现恒幅谐振,记录第二模态固有谐振频率值L21、第二模态对第一模态的耦合量D21和第一模态的输出零偏值E21;获得一组数据P1=(b1,b′1,L11,D11,E11,L21,D21,E21);
(4)按顺序取出潜在解集S1(b,b′)中的第一输入端输入值b和第二输入端输入值b′输入值,重复步骤(3),取得j组数据,其中第i组数据记作Pi=(bi,b′i,L1i,D1i,E1i,L2i,D2i,E2i);所述j为S1(b,b′)中潜在解的个数,i∈[1,j];
(5)计算步骤(4)中确定的j组数据中每组数据的频差L1i-L2i,利用每组数据的频差对该组数据进行筛选,获得频差筛选后的数据集合;具体为:若该组数据的频差处于预先设定的频差范围内,则保留该组数据,否则删除该组数据;
(6)在频差筛选后的数据集合中,按照预先设定的筛选条件选择两组数据,筛选获得的第一组数据中第一输入端的输入值为b11,第二输入端的输入值为b′11,第二组数据中第一输入端的输入值为b22,第二输入端的输入值为b′22,所述b11与所述b22之差的绝对值|b11-b22|<=16;所述b′11与所述b′22之差的绝对值|b′11-b′22|<=16;
(7)利用步骤(6)中得到的两组数据输入端的输入值计算1步长筛选第一输入端输入值取值区间[bb1,bb2]和第二输入端输入值取值区间[bb′1,bb′2];
(8)第一输入端输入值b和第二输入端输入值b′分别以1个A/D最小分辨率为步长在[bb1,bb2]和[bb′1,bb′2]中取值,获得1步长筛选后的潜在解集为S2(b,b′);
(9)用步骤(8)获得的1步长筛选后的潜在解集S2(b,b′)更新步骤(2)中的S1(b,b′),重复步骤(3)~步骤(4),选择第一模态对第二模态的耦合量D1i、第二模态的输出零偏值E1i、第二模态对第一模态的耦合量D2i和第一模态的输出零偏值E2i绝对值的和最小的一组数据为最终数据,所述最终数据为:Pbest=(bbest,b′best,L1best,D1best,E1best,L2best,D2best,E2best);
(10)利用步骤(9)确定的bbest和b′best,计算出四组电极加的电压值大小,具体为:
若bbest的最高位为1,则V1为28V,V3可调,若bbest的最高位为0,则V3为28V,V1可调;所述调整电压的大小由bbest的低七位决定,转换公式为:电压值=低七位的十进制数/127×28;
若b′best的最高位为1,则V2为28V,V4可调,若b′best的最高位为0,V4为28V,V2可调,所述调整电压的大小由b′best的低七位决定,转换公式为:电压值=低七位的十进制数/127×28。
2.根据权利要求1所述的一种MEMS陀螺静电平衡电压值的确定方法,其特征在于:所述步骤(2)中m的取值范围为:8~16。
3.根据权利要求1所述的一种MEMS陀螺静电平衡电压值的确定方法,其特征在于:所述步骤(3)中L11,D11,E11,L21,D21,E21均通过n次试验求平均获得。
4.根据权利要求3所述的一种MEMS陀螺静电平衡电压值的确定方法,其特征在于:所述n的取值范围为:100~200。
5.根据权利要求3所述的一种MEMS陀螺静电平衡电压值的确定方法,其特征在于:所述步骤(5)中预先设定频差范围为[-0.5,0.5]。
6.根据权利要求1所述的一种MEMS陀螺静电平衡电压值的确定方法,其特征在于:所述步骤(6)中预先设定的筛选条件具体为:
所选择的两组数据中第一模态对第二模态的耦合量D1i、第二模态的输出零偏值E1i、第二模态对第一模态的耦合量D2i和第一模态的输出零偏值E2i符号均相反,且两组数据中第一模态对第二模态的耦合量D1i和第二模态对第一模态的耦合量D2i绝对值的和最小。
7.根据权利要求1所述的一种MEMS陀螺静电平衡电压值的确定方法,其特征在于:所述步骤(7)中利用步骤(6)中得到的两组数据输入端的输入值计算1步长筛选第一输入端输入值取值区间[bb1,bb2]和第二输入端输入值取值区间[bb′1,bb′2];具体为:
若|b11-b22|=16,则bb1=min{b11,b22},bb2=max{b11,b22};否则,bb1=(b11+b22)/2-8,bb2=(b11+b22)/2+8;
若|b′11-b′22|=16,则bb′1=min{b′11,b′22},bb′2=max{b′11,b′22};否则,bb′1=(b′11+b′22)/2-8,bb′2=(b′11+b′22)/2+8。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410833945.1A CN104729531B (zh) | 2014-12-26 | 2014-12-26 | 一种mems陀螺静电平衡电压值的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410833945.1A CN104729531B (zh) | 2014-12-26 | 2014-12-26 | 一种mems陀螺静电平衡电压值的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104729531A CN104729531A (zh) | 2015-06-24 |
CN104729531B true CN104729531B (zh) | 2017-10-31 |
Family
ID=53453653
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410833945.1A Active CN104729531B (zh) | 2014-12-26 | 2014-12-26 | 一种mems陀螺静电平衡电压值的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104729531B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110567450B (zh) * | 2019-08-26 | 2021-02-05 | 北京自动化控制设备研究所 | Mems哥氏振动陀螺仪与asic电路的参数匹配方法 |
CN110849388B (zh) * | 2019-09-27 | 2021-12-28 | 北京时代民芯科技有限公司 | 一种基于遗传算法求解mems环陀螺静电平衡调整电压的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3965753A (en) * | 1970-06-01 | 1976-06-29 | Browning Jr Alva Laroy | Electrostatic accelerometer and/or gyroscope radioisotope field support device |
CN101876547A (zh) * | 2009-12-08 | 2010-11-03 | 北京大学 | 一种采用静电平衡梳齿驱动器的水平轴微机械音叉陀螺 |
CN103344227A (zh) * | 2013-06-20 | 2013-10-09 | 上海交通大学 | 静电驱动压电检测体声波谐振三轴微陀螺及其制备方法 |
CN103759722A (zh) * | 2013-12-30 | 2014-04-30 | 北京时代民芯科技有限公司 | 一种环形陀螺仪的静电调整方法及静电调整系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8794068B2 (en) * | 2010-12-01 | 2014-08-05 | Analog Devices, Inc. | Non-degenerate mode MEMS gyroscope |
-
2014
- 2014-12-26 CN CN201410833945.1A patent/CN104729531B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3965753A (en) * | 1970-06-01 | 1976-06-29 | Browning Jr Alva Laroy | Electrostatic accelerometer and/or gyroscope radioisotope field support device |
CN101876547A (zh) * | 2009-12-08 | 2010-11-03 | 北京大学 | 一种采用静电平衡梳齿驱动器的水平轴微机械音叉陀螺 |
CN103344227A (zh) * | 2013-06-20 | 2013-10-09 | 上海交通大学 | 静电驱动压电检测体声波谐振三轴微陀螺及其制备方法 |
CN103759722A (zh) * | 2013-12-30 | 2014-04-30 | 北京时代民芯科技有限公司 | 一种环形陀螺仪的静电调整方法及静电调整系统 |
Non-Patent Citations (2)
Title |
---|
微机械陀螺仪反馈控制器鲁棒性分析及闭环接口检测电路设计;卢海曦等;《东南大学学报(自然科学版)》;20110531;第41卷(第3期);第526-532页 * |
微机械陀螺原理与关键技术;高钟毓;《仪器仪表学报》;19951130;第16卷(第4期);第40-44页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104729531A (zh) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108120451A (zh) | 基于改进pso优化神经网络的硅微加速度计温度补偿方法、系统 | |
CN103335814B (zh) | 风洞中实验模型倾角测量误差数据的修正方法 | |
CN108489512B (zh) | 一种半球谐振陀螺标度因数的补偿标定方法及装置 | |
CN104729531B (zh) | 一种mems陀螺静电平衡电压值的确定方法 | |
CN106323451B (zh) | 通过加速信号获得位移信号的方法及其装置 | |
CN106886628B (zh) | 一种基于索张力不确定性的平面薄膜天线薄膜形状确定方法 | |
CN111578966B (zh) | 一种基于lms算法的半球谐振子特征参数辨识方法 | |
CN113111547A (zh) | 基于缩减基的频域有限元模型修正方法 | |
CN109471061A (zh) | 一种鲁棒处理模型参数误差的接收信号强度差定位方法 | |
CN105258689A (zh) | 一种数字陀螺仪信号控制处理系统 | |
CN114046802B (zh) | 一种光纤陀螺分步温度补偿方法 | |
CN107657073A (zh) | 一种三明治复合材料结构的参数识别方法 | |
CN106441359B (zh) | 一种高精度的微机械陀螺现场模态匹配方法 | |
CN115077564A (zh) | 一种半球谐振陀螺检测通道不平衡误差辨识及补偿方法 | |
CN108228977B (zh) | 一种基于飞行状态参数的直升机振动特征折算方法 | |
CN113865619A (zh) | 一种改善高精度光纤陀螺全温零偏稳定性的方法 | |
CN109388858A (zh) | 基于头脑风暴优化算法的非线性传感器校正方法 | |
Ma et al. | Analysis and compensation of mass imperfection effects on 3-D sensitive structure of bell-shaped vibratory gyro | |
CN106053884A (zh) | 基于坐标旋转法的双轴谐振式加速度计数字控制电路装置 | |
CN106771351B (zh) | 基于数字相频检测法的谐振式加速度计频率锁定电路 | |
CN109029498A (zh) | 一种用于振环硅微陀螺仪模态互易的数字测控系统 | |
CN103759722A (zh) | 一种环形陀螺仪的静电调整方法及静电调整系统 | |
CN101741743B (zh) | 基于网络地址排序的双向流组流方法 | |
CN106571826A (zh) | 一种提高单片机模数转换器动态范围的系统和方法 | |
CN108489602B (zh) | 一种基于模态测试技术的惯性参数识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |