資源描述：
《北航慣性導航大作業(yè).doc》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在行業(yè)資料-天天文庫。
1、慣性導航基礎課程大作業(yè)報告（一）光纖陀螺誤差建模與分析班級：111514：學號2014年5月26日一.系統(tǒng)誤差原理圖二.系統(tǒng)誤差的分析（一）漂移引起的系統(tǒng)誤差1.εx，εy，εz對東向速度誤差δVx的影響clc;clearall;t=1:0.01:25;g=9.8;L=pi/180*39;Ws=2*pi/84.4*60;Wie=2*pi/24;R=g/(Ws)^2;e=0.1*180/pi;mcVx1=e*g*sin(L)/(Ws^2-Wie^2)*(sin(Wie*t)-Wie*sin(Ws*t)/Ws);mcVx2=e*((Ws^2-(Wie^2)*((cos(L)
2、)^2))/(Ws^2-Wie^2)*cos(Ws*t)-(Ws^2)*((sin(L))^2)*cos(Wie*t)/(Ws^2-Wie^2)-(cos(L))^2);mcVx3=(sin(L))*(cos(L))*R*e*((Ws^2)*cos(Wie*t)/(Ws^2-Wie^2)-(Wie^2)*cos(Ws*t)/(Ws^2-Wie^2)-1);plot(t,[mcVx1',mcVx2',mcVx3']);title('Ex,Ey,Ez對Vx的影響');xlabel('時間t');ylabel('Vx(t)');legend('Ex-mcVx1','Ey-mc
3、Vx2','Ez-mcVx3');grid;axissquare;分析：εx，εy，εz對東向速度誤差δVx均有地球自轉周期的影響，εx，εy還會有舒勒周期分量的影響，其中，εy對δVx的影響較大。2.εx，εy，εz對東向速度誤差δVy的影響clc;clearall;t=1:0.01:25;g=9.8;L=pi/180*39;Ws=2*pi/84.4*60;Wie=2*pi/24;R=g/(Ws)^2;e=0.1*180/pi;mcVy1=e*g*(cos(Wie*t)-cos(Ws*t))/(Ws^2-Wie^2);mcVy2=g*sin(L)*e/(Ws^2-Wi
4、e^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t));mcVy3=g*cos(L)*e/(Ws^2-Wie^2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t));plot(t,[mcVy1',mcVy2',mcVy3']);title('Ex,Ey,Ez對Vy的影響');xlabel('時間t');ylabel('Vy(t)');legend('Ex-mcVy1','Ey-mcVy2','Ez-mcVy3');grid;axissquare;分析：εx，εy，εz對北向速度誤差δVy均有地球自轉周期，舒勒周期分量的影響。其中，εx對δVy的影
5、響較大;εy，εz產(chǎn)生的影響幾乎相近。3.εx，εy，εz對東向速度誤差δL的影響clc;clearall;t=1:0.01:25;g=9.8;L=pi/180*39;Ws=2*pi/84.4*60;Wie=2*pi/24;R=g/(Ws)^2;a=Ws^2-Wie^2;b=sin(L);c=cos(L);d=sin(Wie*t);e=cos(Wie*t);f=sin(Ws*t);h=cos(Ws*t);mcL1=Ws^2*0.1/a*(d/Wie-f/Ws);mcL2=(Ws^2*Wie*b/a*(h/Ws^2-e/Wie^2)+b/Wie)*0.1;mcL3=(Ws
6、^2*c*e/Wie/a-Wie*c*h/a-c/Wie)*0.1;plot(t,[mcL1',mcL2',mcL3']);title('Ex,Ey,Ez對mcL的影響');xlabel('時間t');ylabel('mcL(t)');legend('Ex-mcL1','Ey-mcL2','Ez-mcL3');grid;分析：εx，εy，εz對緯度誤差δL均有地球自轉周期的影響，εx還會有舒勒周期分量的影響。4.εx，εy，εz對東向速度誤差δλ的影響clc;clearall;t=1:0.01:50;pi=3.14;g=9.8;L=pi/180*39;Ws=2*pi/8
7、4.4*60;Wie=2*pi/24;R=g/(Ws)^2;a=Ws^2-Wie^2;b=sin(L);c=cos(L);d=sin(Wie*t);e=cos(Wie*t);f=sin(Ws*t);h=cos(Ws*t);mcLONG1=(tan(L)/Wie*(1-e)-Wie*c*h/a)*0.1;mcLONG2=(sec(L)*(Ws^2-Wie^2*c^2)*f/Ws/a-Ws^2*tan(L)*b*d/Wie/a-t*c)*0.1;mcLONG3=(Ws^2*b*d/Wie/a-Wie^2*b*f/a/Ws-t*b)*0.1;p
