振动整流法辨识加速度计非线性误差模型系数
来源:爱问旅游网
第34卷第9期 2013年9月 仪 器 仪 表 学 报 Chinese Journal of Scientific Instrument Vo1.34 No.9 Sep.2013 振动整流法辨识加速度计非线性误差模型系数 球 鲁金瑞,任顺清,王常虹 (哈尔滨工业大学空间控制与惯性技术研究中心哈尔滨150080) 摘要:为了研究精确标定加速度计在线振动与高g的环境下的误差模型系数的方法。首先介绍了双轴离心机的特点。在反 转平台上设计了双轴夹具,然后精确计算了加速度计3个轴的输入比力。经过主轴和反转平台轴正反角位置静态实验,以及主 轴与反转平台轴以等值反方向的匀角速率旋转的动态实验,通过振动整流法,并补偿掉静态均值,精确辨识了加速度计的误差 模型系数。依据仿真结果可得:该实验方法能够抑制离心机误差对标定精度的影响,其中反转平台的静动态失准角对标定影响 最大,应该提高离心机静动态失准角的测试精度,并加以补偿。 关键词:振动整流;加速度计;误差模型系数;双轴离心机;失准角 中图分类号:U666.1 文献标识码:A 国家标准学科分类代码:590.30 Identification method of accelerometer nonlinear error model coefficients with vibration rectification Lu Jinrui,Ren Shunqing,Wang Changhong (Space Control and Inertial Technology Research Center,Harbin Institute of Technology,Harbin 150080,China) Abstract:In ord er to study the methods of calibrating the accelerometer error model coefficients accurately under linear vibration and high g composite environment,the characteristics of the dual—axis centrifuge are described at first.A dual-axis fixture was designed and fixed on the counter—rotating platform.And the input speciifc forces of the three axes of the accelerometer were calculated accurately.Then,the static test was performed,in which the spindle axis rotates to equal-・spaced angular position and the axis of the counter・-rotating platform rotates reversely to corresponding equal position.The average value of the static test was compensated with vibration rectification method.The dynamic test was performed,in which the spindle axis rotates with constant angular rate and the counter-rotating platform axis rotates in the opposite direction with the same angular rate;then the accelerometer error model coefficients were identified accurately.The test simulation results show that the influence of the cen— trifuge error on calibration accuracy is suppressed greatly using this test method;and the influence of static and dynamic misalignment angles of the counter-rotating platform on the calibration accuracy is quite large and can not be ignored.The measurement accuracy of the centrifuge static and dynamic misafignment angles should be improved and compensated. Keywords:vibration rectification;aceelerometer;error model coefifcient;dual—axis centrifuge;misalignment angle 收稿日期:2013-05 Received Date:2013-05 基金项目:十二・五预研项目(51309050202)资助 第9期 鲁金瑞等:振动整流法辨识加速度计非线性误差模型系数 2131 1 引 言 2离心机误差 2.1离心机的坐标系建立 伴随着当代航空航天技术的进步,要求对于惯性器 件在复合的近真实的环境下,进行更为精确的误差模型 系数标定,补偿后可提高导航精度。这为双轴离心机的 应用提供了很好的契机(能够提供线振动与高g(重力加 速度,单位m/s )的复合力学环境)。线振动与高g的环 境条件为飞行器的发射、变轨、大机动等提供复合力学环 境,国内这方面的研究工作还比较少。本文依据现有的 在单轴离心机上惯性器件标定研究基础,展开探索性的 研究。 图1给出了双轴离心机的结构简图以及相应的坐标 系,其中:基准坐标系OoX。Y。 即为地理坐标系,主轴坐标 系0 。Ylz 相对地理坐标系存在位姿变换,方位轴坐标系 02X y2 :相对主轴坐标系存在位姿变化( 为反转平台轴 与主轴的同步跟踪角度差),夹具坐标系O3X Y z 相对于 方位轴坐标系存在平移变化,加速度计坐标系O4X y4 是 夹具绕X 轴转动0角,然后在新形成的坐标系中夹具绕 Y 轴转动叼角。夹具体相当于一个双轴位置台,通过这2 个角位置可以改变加速度计的三轴对于反转平台的指 向。图中叩为双轴台夹具 轴的转动角度,0为双轴台夹 具 轴的转动角度。 (一 +△∞)什 精密离心机产生高于1 g的精密的向心加速度、谐波加 速度等,用来标定加速度计的高阶非线性系数u 、陀螺仪的 静态误差模型系数 、惯性组合的2,4个误差模型系数 等。 文献[3]给出了加速度计误差模型在离心机上的试验设计。 文献[4]分析了离心机误差对加速度计标定精度的影响,使 用的离心机不带反转平台。文献[5]设计了在反转平台离心 机上标定陀螺加速度表的方法,解决了陀螺加速度计测试 时,离心机大臂旋转引起的牵连运动影响的问题。带反转平 台的离心机,也称双轴离心机 ,一般有2个平行轴系,即主 轴轴系和带反转平台的轴系,主轴和反转平台轴都以匀角速 率旋转,2个角速率矢量大小相等方向相反。它有2个特点, 一(局部放大) 1主轴:2反转平台轴;3工作台:4夹具;5加速度计 是在反转平台上角速率相对喷性空间很小,仅有地球自转 角速率分量;二是与反转平台固联的各点比力均为振动形 图1 双轴精密离心机结构简图 Fig.1 The schematic structure diagram of dual—axis centrifuge 式,且与反转平台上的安装位置无关。此类离心机有类似于 线振动台的功能。文献[8]采用谐波分析法在双轴离心机上 标定陀螺仪的与比力相关的误差模型系数的方法。文献 2.2双轴离心机的误差源 主要研究双轴离心机的误差包括:主轴铅垂度 、 [10—11]给出了线振动台标定加速度计的非线性系数的方 +△ (动态失准角)、 ),反转平台 法,线振动台仅在一个方向振动,它的振动位移与加速度的 卢。,轴线平行度(精度不容易保证。文献[6,12一l8]给出了加速度计的误差模 对准误差 ,离心机动态半径△ ,主轴与反转平台轴的 型以及相关的标定方法。本文将利用双轴离心机的线振动 速率对称性△ (∞=一n+△ )。表1给出了以上7个 台功能,考虑离心机误差源,精确计算加速度计的输入比力, 误差源标称姿态给出的相邻坐标系间的齐次变换 用振动整流法来标定加速度计的非线性误差模型系数,并将 分析离心机误差对标定系数的影响关系。 阵 。 表1 离心机误差项、相关坐标系及其齐次变换式 Table 1 The errors,related coordinate frame and its homogeneous transformation of the centrifuge 变换阵为: 由表1可知,夹具坐标系相对于基准坐标系的齐次 = 1。2=[ ] (1) 式中:A为夹具体坐标系相对于地理坐标系的姿态矩阵; 2132 仪器仪表学报 第3 4卷 w为夹具坐标系原点相对地理坐标系的位移矢量。 3加速度计的输入比力计算 3.1加速度的计算公式 依据加速度的合成定理,在地理坐标系中,加速度可 以表示为: 口 。= + ×W+2.o0× + ×( ×w,) (2) 式中: 为地球相对于惯性空间的自转角速率,在基准 坐标系下的表示。 3.2加速度坐标系下,输入比力的计算 再通过姿态矩阵的转化,求得夹具坐标系中的离心 机加速度分量与重力加速度分量(以下计算均略去二阶 小量): 口:,=A 口 o g :Alr[0 0] (3) 如图1所示,改变夹具体的2个角度 、田,输入比力 将重新分配为: : 厂11 0 0 一sin叼] l0 c。s 0 Il (口: +g ) 一sin COS∞ (4) 输入比力的计算结果为: a =一sin叩[(R (O/l+△ 1)+2尺n(【, COS L・ COS,Ot)/g+1]一sin,cos ,7[一sin(Ao)t+ ) sin + COS(Atot+ )((尺+AR) +2R ̄xo sin )COS,Or]/g+ COS ̄)COS 77[一sin(Awt ) COS 一cos(Atot+ )(( +AR) +2尺 幻 sin )sin ]/g (5) a =COS 7/sin 0[(R ( l+△ 1)+2RDxo COS L・ COS-r )/g+1]一(sin(bcos 0一COS cbsin"r/sin )[一 sin(Atot+ )R COS f 一cos(Awt+ )((R十AR) + 2RDxo ・sin L)sin ]/g一(COS c0s 0+sin ̄bsin+/sin )[一sin(Atot+7)尺 sin +cos(Awt+ )((R+ △R) +2 力 ・sin L)COS,Or]/g (6) a =COS叼c0s [(R (OLl+△ 1)+2RO ̄o COS L・ COS,Ot)/g+1]+(COS ̄bsin0一sin ̄bsin Vcos )[一 sin(htot+y)R ̄sin +c0s(Awt+ )((R+AR) + 2R ・sin L)COS ]/g+(sin ̄bsin +COS 6sin叼c0s ) [一sin(Aeot+ )R cos 一c0s(Atot+ )((R+△ ) + 2冗 sin L)sin伫]/g (7) 式(5)~(7)的推导结果中不含 、 、r ,所以夹具 坐标系的加速度与其在反转平台上的安装位置无关,可 以省去调心。 当z轴作为输入轴时(图2中C7安装方式),则输入 轴上比力的精确表达式可以通过式(7)得: ac7z=COS叼c0s 0(R』 (Otl+△ 1)+2ROw COS L・ COS /g+1)一(COS: S d ̄sin C 0一sin 6sin叼c0s 0){ o+( 1+ 0 n △ 1)cos( ̄t+ )+卢1—1●●●●●●J........。........ si..Ln( /2t+y)一[COS(Atot+y)((尺+ △ ) +2R/ ̄ sin C r L)cos.1 S +R,( ̄sin(Awt+T)sin ̄t]/g}一 (sin ̄bsin +COS ̄喝0 n 叩 bsin ycos叶 0){30+卢lcos( ̄t+ )一 0【osin( + )一[COS(△∞£+ )((R+AR)O 1 0 + 2R 幻 sin L)sin 一R sin(Atot+ )cos ̄Qt]/gl 4 加速度计误差模型及振动整流法 加速度计输出的误差模型 可表示为: E=K1( +a +K2a + a 十Kppa:+Ko。a:+ K a a +K a a。+8) (8) 式中:E为加速度计输出值,单位:V;a 、a。、Ⅱ 为加速度计 输入轴、输出轴、摆轴上的输入比力分量,单位:g; 为 零位偏置,单位:g;K 为标度因数,单位:V/g; 、K 、K。 为各轴二阶非线性系数,单位:g/g ; 为三阶非线性系 数,单位:g/g ; 、K 为交叉耦合系数,单位:g/g ;s为 随机误差,单位:g。 振动整流误差分析是指给加速度计输入单一的定频 正弦振动、一系列频率的正弦振动或特定的随机振动,将 加速度计在振动中输 的直流平均值减去静态时输出的 直流平均值再除以振动幅值的交流有效值的平方 ( ),以得到所需的系数 。 在振动实验中,加速度计输出的平均值: E =Kl[ +(口 ) +K2(。 ) +K3(0 )3 + K即(r王 ): +K∞(ao) +K (ai0 ) +K (ai 。) +8] 以图2中的c7安装方式为例说明振动整流法的求 取过程,将加速度计的输入比力代人到误差模型中,并采 用整周积分取平均的方法,可以得到加速度计的整流输 出为(忽略小量): E(7):K1[ +( ( 1+△ 1)/g+1)+ + ÷K ( ) +÷ 。( /g) 十s] (9) 取加速度计输出的整周平均值与静态输出的平均值 的差,即: Ec(n)=E(n)一E静态均值 振动整流后的输出,扣除静态实验中输出值的均值 之后,c7安装方式的最终结果为: Ec(7)=K1[ +尺 (Oll+Aa1)/g+ 1 、…/ ) +÷ 。(R ̄/g) + ] (10) 下面讨论如何通过振动整流实验的方法将加速度计 第9期 鲁金瑞等:振动整流法辨识加速度计非线性误差模型系数 2133 的高阶系数辨识出来,并且建立双轴离心机误差项与加 速度计误差系数辨识精度之间的关系。 5 加速度计在双轴离心机上的试验设计 依据加速度计测试标准 川,采用振动整流法(加 速度计输入比力采用积分方式确定),通过初始安装位置 的设计,确定其误差模型方程的非线性项。具体的安装 方式如图2所示,图中a为向心加速度, 轴表示输入 轴,PA轴为摆轴,OA轴为输出轴。所设计9种初始安装方 式。 +。 f 初始时刻,夹具坐标系与主轴示意图 fg厶 翳 H g IA , I {s 1s 『g .I—— ,■—, 刚I g 图2安装位置示意图 Fig.2 The schematic diagram of the installing positions 由图2可知:初始安装位置通过反转平台的同步跟 踪与旋转双轴夹具可以实现c1~c9的位置,具体转动 过程如表2所示。 表2从初始位置到C1~C9位置,需旋转的角度 Table 2 The angles required to rotate from initial position to C1一C9 positions 利用静态测试法剔除部分误差项,在图2中每种初 始安装位置下,即双轴离心机的2个轴均处于位置状态, 主轴每转动15。、反转轴反向转动15。后停止转动,并记 录一次加速度计输出;转动一周后,求出这一周输出的平 均值。 一 在双轴离心机实验时,对于图2中每种初始安装位 置下,振动整流后求出加速度计的输出,并减去静态试验 时对应的输出均值。这样就相当于消除了本实验中由于 重力加速度、加速度计安装误差对加速度计输入比力的 影响(而加速度计安装误差项的有些误差是不能准确测 量的,这对提高加速度计的标定精度具有重要的意义)。 与Ec(7)同理,可得: Ec(1)=K1{Ko+ ( ( 。+Aa,) )+ + +÷d 十 +÷ 2K +s} (11) Ec(2):K{ + ( ( 。+Aoq)/g)+÷ + 警 +丢dl2 + d12K。一÷ Kip+s} (12) Ec(3):K.{Ko一 ( ( +Aa,)/g)+扣 一 +÷d12 + 。一÷ Kip+s} (13) Ec(4):Ki{ + ( ( +Aa1)/g)+丢d. + dl2 +÷ + +÷dl2 +s} (14) ,1 1 (5)=K1{ 一等( ( +Aa )/g)+÷d。 一 dl2 +÷ +÷ +÷d1 2Ki。+s }(15) Ee(6)= { 一譬( ( 。+Aa1)/g)+丢 一 警 +÷ +÷ 一÷dl2 + } (16) c(8):K {Ko+÷d +÷d1 2Ko。+s}(17) Ec(9):K.{ +÷d +÷d +s}(18) 式中: dl=R ̄Z/g,d2=尺 (Oll+△ 1)/g (19) 式中: 为离心机主轴的转速,R为离心机的两竖直转动 轴之间的距离,O/ 、△ 分别为反转平台的静、动态失准 角. 莺力加速度 2134 仪器仪表学报 第3 4卷 依据式(11)~(12)可以解算出 后,再通过 式(13)一(14)可以结算出 ,式(15)~(16)可以解 算出 ,式(10)、(17)~(18)可以解算出 、 、 。。。 在初始安装方式下,采用整周积分的方式,依据 下面依据图2方式试验依据式(20)~(25)给出仿 真算例,阐述离心机误差与模型系数之间的关系,从而 说明双轴离心机的误差对于加速度计参数标定的 影响。 式(10)~(18)、图2C1~C9,进而可以推算出: =(Ec(7)一Ec(8)一Ec(9)一d2Kl+ (20) 6计算算例 工作台的纬度为45.732 8。,双轴离心机各误差均为 K )/(K d ) K3:((Ec(4)一Ec(3)+Ec(2)一Ec(1))/K.一 42a:)/(3√2d /4) =(21) (22) (23) (24) (25) 10 量级(对于位移误差相当于10~m=10 Ixm,角度误差 相当于10一rad一2 。设某双轴离心机标称半径为1 m,离 心机的转速应为3-rr rad/s=540。/s,此时产生的向心加速 ( c(7)一Ec(8)+Ec(9)一d K1一 K )/(K d, ) K。。=(Ec(7)+Ec(8)一Ec(9)一K d:一 K。))/(K d。 ) K =2(Ec(2)一Ec(1))/(K d ) K 。:2(Ec(5)一Ec(6))/(K d。 ) 度约为9g;分别对加速度计进行9位置振动整流实验,记 录相应的加速度计输出。依据文献[4]给出的加速度误差 模型数据,转换到本实验用的模型,得出对应的系数:K = 0.69 V/g、 :1.449 3×10一 g,K2=一2.26×10一。g/g 、 由式(2o)~(25)可以清晰的看到,离心机误差项 K3=1.48 X 10一 g/g 、K =一1.12×10~g/g2K、= 。。对标定系数的影响。具体计算时,输入轴上的比力精确到 一1.45×10一 g/g 、Ki。=5.68×10 g/g 、Ki。=5.51× 阶小量,摆轴与输出轴上的比力只计算标称值(一阶小 由于采用振动整流方法扣除了静态平均值,旋转角 10 g/ 。首先利用式(5)~(7)计算相应位置的加速度 计输入值,再依据式(20)一(25)可以计算出加速度计的 高次误差系数。表3给出的数据仿真分析中,双轴离心 机各误差单独作用下对加速度计高次项误差系数的影响 大小。 量与误差模型系数的乘积为二阶小量,可以忽略)。 度误差对标定结果的影响可以忽略不计,故本实验对于 旋转角度的精度要求不高,但要保证角位置的稳定性。 表3双轴离心机误差与加速度计模型系数标定精度之间的关系 Table 3 Relationship between dual-axis centrifuge error and model coeficient calfibration accuracy of accelerometer 从表3可以看出,离心机误差项对加速度计高次项误差 反时,反转平台上,任意一点相对惯性空间的角速度约为 0,能够隔离掉由于主轴旋转引入的牵连旋转运动;带反 系数有一定的影响。依据所设计的测试方法,双轴离心机的 误差OL。 、 ,对标定系数的影响均可忽略;静态与动 态半径误差AR对标定系数的影响也可忽略;动态失准角 As 与静态失准角Ot.的测试误差对加速度计误差模型系数 转的离心机2个转动轴的转速精度相对容易控制,能够 给出精确的谐波加速度。提出了利用振动整流法在辨识 加速度计的非线性误差模型系数的应用,从推导的公式 与仿真结果可得如下结论: 的标定影响很大,对个别模型系数的影响甚至比真实值大2 个数量级,应提高其测试精度,并需要对此项进行补偿,以提 高标定的精度。实际计算时,考虑或不考虑地球自转角速率 对非线性误差模型系数标定精度无影响。 1)由于振动整流效应和补偿静态值,对加速度计的 夹具姿态误差要求可以放宽。 2)采用振动整流法,并补偿静态输出后,能够最大限 度的抑制离心机误差对加速度计标定精度的影响。依据 所设计的测试方法,地球自转对标定精度的影响可以忽 略,离心机主要误差项对于标定精度的影响均可忽略,反 转平台的静动态失准角影响不可忽略,对静动态失准角 7 结 论 针对双轴离心机的特点:既能提供线加速度又能够 提供谐波加速度,当主轴与反转轴转速大小相等方向相 第9期 鲁金瑞等:振动整流法辨识加速度计非线性误差模型系数 2135 必须进行精密测试后加以补偿。 3)对于加速度计的误差模型中的非线性系数标定, 振动整流法的优势比较明显,只需要补偿离心机的动态 失准角,消除了加速度计安装误差对标定精度的影响,也 GAO,X ZH,JIANG X X,CHEN W,et a1.Error model identification method of aceelerometer based on semi--par-- ametric regression『J 1.Journal of Chinese Inertial Tech— nology,2012,20(3):352—357. 消除了离心机工作半径误差的影响。 4)线振动台的位置精度难于控制,所以谐波加速度 精度不高,标定的误差模型系数的精度难于保证;而双轴 离心机的谐波加速度计精度主要决定于旋转角速率精 度,易于控制。所以本文提供的双轴离心机上的测试方 法比用线振动台测试时的精度要高。 参考文献 [1] std 836TM-2001,Gyro and accelerometer panel of the IEEE aerospace and electronic system society,IEEE rec— ommended practice for precision centirfuge testing of lin— ear accelerometers[S].New York,The Institute of Elec— trical and Electronics Engineers,2001. [2] std 836 -2009,Gyro and accelerometer panel of the IEEE aerospace and electronic system society IEEE rec— ommended practice for precision centirfuge testing of lin— ear accelerometers[s].The Institute of Electircal and Electronics Engineers,2009. [3] 苏宝库,李丹东.加速度计精密离心机试验的优化设 计[J].中国惯性技术学报,2010,18(5):620—624. SU B K,LI D D.Optimal design of accelerometer test on precision centirfuge[J].Journal of Chinese Inertial Tech— nology,2010,18(5):620-624. [4] 王世明,任顺清.精密离心机误差对石英加速度计误 差标定精度分析[J].宇航学报,2012,33(4),520—526. WANG SH M,REN SH Q.Relationship between calibra— tion accuracy of error model coefifcients of accelerometer and errors of precision centirfuge[J].Journal of Astronau— tics,2012,33(4),520—526. [5] 乔永辉,刘雨,苏宝库,等.陀螺加速度计误差模型系 数离心机测试方法研究[J].宇航学报,2007,28(4): 854—859,931. QIAO Y H,LIU Y,SU B K,et a1.Research on testing method of error model coefifcients of gyro accelerometer on centirfuge[J].Journal of Astronautics,2007,28(4): 854・859,931. [6]ZADFIYA RAJESH L,SWAMY K B M,SEN,SIDDHAR. THA.An RF based centirfuge calibrator for MEMS accel- erometer testing f C 1.ICECE 2010-6th International Con— ference on Electrical and Computer Engineering,2010: 412.414. [7] 高贤志,蒋效雄,陈伟,等.一种基于半参数回归的加 速度计误差模型辨识方法[J].中国惯性技术学报, 2012,20(3):352.357. [8] 王世明,任顺清,王常虹.在反转平台离心机上标定陀 螺仪的误差模型系数的方法[J].中国惯性技术学报, 2012,20(5):625-630. WANG SH M,REN SH Q,WANG CH H.Calibration of gyro error model coefifcients on precision centirfuge with counter—rotating platform[J].Journal of Chinese Inertila Technology,2012,20(5):625—630. [9]张华强,赵剡,陈雨.捷联惯性导航系统整体标定新方 法[J].北京航空航天大学学报,2012,38(4): 459—463. ZHANG H Q,ZHAO Y,CHEN Y.New system calibration method for strapdown inertila navigation system[J].Jour— nal of Bering University of Aeronautics and Astronautics, 2012,38(4):459—463. [10] Std 1293—1998,IEEE standard speciifcation recommended format guide and test procedure for linear,single—axis, nongroseopic aceelerometers[S].The Institute of Electri— cal and Electronics Engineers,1998. [1 1] 中国人民解放军总装备部.单轴摆式伺服线加速度计 试验方法[S].中华人民共和国国家军用标准, FL 6615. General Armaments Department of the People’s Liberation Army.Test methods for single—axis pendulous servo linear accelerometers『S].The National Military Standard of the People’s Republic of China,FL 6615. [12]刘骏跃,陈明.SAW加速度计量程调整及计量标定 方法研究[J].仪器仪表学报,2006,27(9): l052 1056. LIU J Y.CHEN M.On adjusting measure range&cali— bration for the SAW accelerometer『J 1.Chinese Journal of Scientiifc Instrument,2006,27(9):1052—1056. [13]乔仁晓,孟晓风,李皎,等.离心机法校准加速度计不 确定度的统计模拟算法[J].电子测量与仪器学报, 2010,24(5):450—453. QIAO R X,MENG X F,LI J,et a1.Statistical simulation algorithm to analyzing the uncertainty in calibration of ac— celerometers by centirfuge[J].Journal of Electronic Measurement and Instrument,2010,24(5):451—453. [14]张春熹,魏渊,张唏,等.数字闭环加速度计系统模型 分析与校正设计[J].仪器仪表学报,2012,33(12): 2833—2839. ZHANG CH X,WEI Y。ZHANG X,et a1.Model analysis and emendation design of digitla close—loop accelerometer system[J].Chinese Journal of Scientiifc Instrument. 2136 仪器仪表学报 第3 4卷 2012,33(12):2833—2839. 『15 1 KHALIYULIN R.Experimental calibration and va ̄dation of a micromechanical accelerometer model suited for sys— ten—level design[C].Procedia Chemistry,Proceedings of the Eurosensors 23rd Conference.2009. [16]任春华,潘英俊,刘进江,等.基于多温感测的加速度 计组件的系统补偿研究[J].仪器仪表学报,2008, 29(10):2130—2134. REN CH H,PAN Y J,LIU J J,et a1.Study on system compensation of quartzose flexible accelerometer module based on multi—temperature sensing[J].Chinese Journal of Scientiifc Instrument,2008,29(1O):2130-2134. [17] 许高斌,朱华铭,陈兴.不等高梳齿电容式三轴MEMS 加速度传感器[J].电子测量与仪器学报,2011, 25(8):704-710. XU G B.ZHU H M,CHEN X.Tri—aixs capacitive aceeler— ation sensor with different height combs『J].Journal of Electronic Measurement and Insturment,201 1,25(8): 704-710. [18] 王新龙,申功勋,何乃刚.一种有效的加速度计静态模 型辨识方法[J].仪器仪表学报,2003,24(1):57-60. WANG X L,SHEN G X,HE N G.An efficient discrimi— nation method of the accelerometer static model para--me-- _~ters f J 1.Chinese Journal一 of Sci~ent~iif一c I一nstur_ ment,2003, 24(1):57-60. 作者简介 Ren Shunqing received bachelor degree in 1989,master de— gree in 1994 and doctor degree in 5999 all from Harbin Institute of Technology;he is a professor and doctoral supervisor in Space Control and Inertiat Techno ̄gy Research Cen ̄r.HIT.His main research directions are precision measurement technology,ine ̄ial technology and system identification.