北航惯性导航综合实验一实验报告.doc

实验一陀螺仪关键参数测试与分析实验加速度计关键参数测试与分析实验二零一三年五月十二日实验一 陀螺仪关键参数测试与分析实验一、 实验目的通过在速率转台上的测试实验，增强动手能力和对惯性测试设备的感性认识；通过对陀螺仪测试数据的分析，对陀螺漂移等参数的物理意义有清晰的认识，同时为在实际工程中应用陀螺仪和对陀螺仪进行误差建模与补偿奠定基础。二、 实验内容 利用单轴速率转台，进行陀螺仪标度因数测试、零偏测试、零偏重复性测试、零漂测试实验和陀螺仪标度因数与零偏建模、误差补偿实验。三、 实验系统组成 单轴速率转台、MEMS 陀螺仪（或光纤陀螺仪）、稳压电源、数据采集系统与分析系统。四、 实验原理 1. 陀螺仪原理陀螺仪是角速率传感器，用来测量载体相对惯性空间的角速度，通常输出与角速率对应的电压信号。也有的陀螺输出频率信号（如激光陀螺）和数字信号（把模拟电压数字化）。以电压表示的陀螺输出信号可表示为： (1-1)式中是与比力有关的陀螺输出误差项，反映了陀螺输出受比力的影响，本实验不考虑此项误差。因此，式（1-1）简化为 (1-2)由（1-2）式得陀螺输出值所对应的角速度测量值： (1-3)对于数字输出的陀螺仪，传感器内部已经利用标度因数对陀螺仪模拟输出进行了量化，直接输出角速度值，即： (1-4)是是陀螺仪的零偏，物理意义是输入角速度为零时，陀螺仪输出值所对应的角速度。且 (1-5)精度受陀螺仪标度因数、随机漂移、陀螺输出信号的检测精度和的影响。通常和表现为有规律性，可通过建模与补偿方法消除，表现为随机特性，可通过信号滤波方法抵制。因此，准确标定和是实现角速度准确测量的基础。五、 陀螺仪测试实验步骤1) 标度因数和零偏测试实验a. 接通电源，预热一定时间；b. 陀螺工作稳定后，测量静止情况下陀螺输出并保存数据；c. 转台正转，测试陀螺仪输出，停转；转台反转，测试陀螺仪输出，停转。在正转和反转时测试陀螺仪输出量，并分别保存数据；d. 改变转台输入角速率重复步骤c，正负角速率的速率档分别不少于5 个（按军标要求是11 个）；e. 转速结束后，当转台静止时，采集陀螺仪输出数据，并保存。f.根据最小二乘法公式 (1-6) (1-7)计算陀螺标度因数和零偏。2) 零漂测试（零偏稳定性）在静止下采集陀螺仪数据，并由测试数计算陀螺仪零偏稳定性。军标中通常的测试时间是1 小时，并对所采集的数据进行1 秒、10 秒及100秒等不同时间的平滑。本实验中可采集数据10 分钟左右，并分别进行1 秒、10 秒及100 秒平滑。按如下公式 (1-8)计算陀螺仪零偏稳定性，并进行比较。3) 零偏重复性测试a. 令转台某角速度200下进行正转，转速平稳后，采集陀螺输出数据，并保存。b. 令转台某角速度-200下进行反转，转速平稳后，采集陀螺输出数据，并保存。c. 按计算陀螺零偏；d. 关掉陀螺电源，并重新启动，重复步骤a、b；e. 重复步骤d 进行3-5 次，共得到陀螺零偏5-7 个；f. 对5-7个陀螺零偏按下式(1-9) (1-9)求均方差，得零偏重复性指标。六、 实验结果1. 数据处理将原始数据剔除后绘图如下2. 计算陀螺标度因数和零偏根据陀螺在10/s，20/s，40/s，60/s ，80/s角速率下正反转的输出，分别求得正转下陀螺的标度因数和零偏，及反转下陀螺的标度因数和零偏，然后求的均值。 = 0.9901 = 0.03583. 零偏稳定性对所采集的数据进行1 秒、10 秒及100秒等不同时间的平滑，如下图。零漂计算结果（1000s平滑）：Bs= 0.0144 4. 零偏重复性以角速度40/s正反转，共采集5组数据组号陀螺零偏10.1878920.1784430.1696540.1845950.19401零偏重复性：0.009309七，实验小结由零漂平滑后的结果可知，对采集的数据平滑时间长可以提高零偏的稳定性。八，源程序%加载数据%Gyro_0end=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_0end.txt);Gyro_0start=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_0start.txt);Gyro_10n=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_10n.txt);Gyro_10p=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_10p.txt);Gyro_20n=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_20n.txt);Gyro_20p=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_20p.txt);Gyro_40n=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_40n.txt);Gyro_40p=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_40p.txt);Gyro_60n=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_60n.txt);Gyro_60p=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_60p.txt);Gyro_80n=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_80n.txt);Gyro_80p=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Gyro_data1 标度因数和零偏测试Gyro_80p.txt);%剔除不合格数据%Gyro_10p=Gyro_10p(find(Gyro_10p9)&(Gyro_10p15)&(Gyro_20p-50)&(Gyro_40n35)&(Gyro_40p50)&(Gyro_60p70)&(Gyro_80p-0.8)&(Gyro_result-50)&(Gyro135)&(Gyro3-2)&(A0-1)&(A600)&(A900)&(A1200)&(A1800)&(A270-1)&(A3001.5); plot (A120, DisplayName, A120, YDataSource, A120); hold all; plot (A150, DisplayName, A150, YDataSource, A150); plot (A180, DisplayName, A180, YDataSource, A180); plot (A210, DisplayName, A210, YDataSource, A210); plot (A240, DisplayName, A240, YDataSource, A240); plot (A270, DisplayName, A270, YDataSource, A270); plot (A30, DisplayName, A30, YDataSource, A30); plot (A300, DisplayName, A300, YDataSource, A300); plot (A330, DisplayName, A330, YDataSource, A330); plot (A60, DisplayName, A60, YDataSource, A60); plot (A90, DisplayName, A90, YDataSource, A90); hold off; figure(gcf)title(剔除数据后,fontsize,12);%5555计算标度因数%5f0=mean(A0);f30=mean(A30);f60=mean(A60);f90=mean(A90);f120=mean(A120);f150=mean(A150);f180=mean(A180);f210=mean(A210);f240=mean(A240);f270=mean(A270);f300=mean(A300);f330=mean(A330);%求正转标度因数%F= f0 f30 f60 f90 f120 f150 f180 f210 f240 f270 f300 f330; %以指向地心为初始零度A= 1 1*cos(pi/6) 1*cos(pi/3) 1*cos(pi/2) 1*cos(pi*2/3) 1*cos(pi*5/6) 1*cos(pi) 1*cos(pi*7/6) 1*cos(pi*4/3) 1*cos(pi*3/2) 1*cos(5/3*pi) 1*cos(22/12*pi);J=1*f0 1*cos(pi/6)*f30 1*cos(pi/3)*f60 1*cos(pi/2)*f90 1*cos(pi*2/3)*f120 1*cos(pi*5/6)*f150 1*cos(pi)*f180 1*cos(pi*7/6)*f210 1*cos(pi*4/3)*f240 1*cos(pi*3/2)*f270 1*cos(5/3*pi)*f300 1*cos(22/12*pi)*f330 ;KG0=(sum(J)-(sum(F)*sum(A)/12)/(sum(A.2)-(sum(A)*sum(F)/12); %KG0=-1.0159%求零偏%55F01=sum(F)/12-KG0*sum(A)/12;%F01=-0.036595%加载静止时的数据%S1=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data2 偏置稳定性测试偏置测试0度.txt);S2=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data2 偏置稳定性测试偏置测试180度.txt);ff0=(mean(S1)+mean(S2)/2; %ff0=-0.011212g%零偏重复性测试D10=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data3 加计偏置重复测试1_0.txt);D11=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data3 加计偏置重复测试1_180.txt);D20=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data3 加计偏置重复测试2_0.txt);D21=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data3 加计偏置重复测试2_180.txt);D30=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data3 加计偏置重复测试3_0.txt);D31=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data3 加计偏置重复测试3_180.txt);D50=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data3 加计偏置重复测试5_0.txt);D51=load(E:惯性器件综合实验惯性导航试验数据1惯导实验1实验数据Acc_data3 加计偏置重复测试5_180.txt);d1=(mean(D10)+mean(D11)/2;d2=(mean(D20)+mean(D21)/2;d3=(mean(D30)+mean(D31)/2;d5=(mean(D50)+mean(D51)/2;d=d1 d2 d3 d5;h0=mean(d);Br=sum(d-h0).2)/3 Br1=sum(Br.0.5);