北航慣性導(dǎo)航大作業(yè)

《北航慣性導(dǎo)航大作業(yè)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、慣性導(dǎo)航基礎(chǔ)課程大作業(yè)報(bào)告（一）光纖陀螺誤差建模與分析班級(jí)：111514姓名:學(xué)號(hào)2014年5月26日?系統(tǒng)誤差原理圖R?代.nf—I

2、)/Ws)；mcVx2=e*((Ws^2?(Wie^2)*((cos(L))x2))/(Ws^2-Wiex2)*cos(Ws*t)?(WsA2)*((sin(L))^2)*cos(Wie*t)/(Wsx2-Wiex2)-(cos(L))A2)；mcVx3=(sin(L))*(cos(L))*R*e*((Wsx2)*cos(Wie*t)/(Wsx2-Wie^2)-(WieA2)*cos(Ws*t)/(WsA2-WieA2)-1)；plot(t,[mcVxl1,mcVx21zmcVx3'])；title(1Ex,EyzEz對(duì)

3、Vx的影響1);xlabel.「時(shí)間t*)；ylabel('Vx(t)1)；legend('Ex-mcVxl1,'Ey-mcVx2'JEz-mcVx3')；grid；axissquare；￡z對(duì)東向速度誤差5Vy的影響分析：ex,ey,￡z對(duì)東向速度誤差§Vx均有地球自轉(zhuǎn)周期的影響，ex,￡y還會(huì)有舒勒周期分量的影響，其中，對(duì)§Vx的影響較大。2.ex,ey,clc;clearall;t二1:0.01:25;滬9.8;L=pi/180*39;Ws=2*pi/84.4*60;Wie=2*pi/24;R=g/(Ws)八2;e

4、=0.1*180/pi;mcVyl=e*g*(cos(Wie*t)-cos(Ws*t))/(WsA2-WieA2);mcVy2=g*sin(L)*e/(WsA2-WieA2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t));mcVy3=g*cos(L)*e/(WsA2-WieA2)*(sin(Wie*t)-Wie/Ws*sin(Ws*t));plot(t/[mcVyl,/mcVy2,/mcVy3,]);title('ExzEy,Ez對(duì)Vy的影響')；xlabelf時(shí)間f);ylabel('Vy(t)');le

5、gend(,Ex-mcVyl7Ey-mcVy2'/,Ez-mcVy3');grid;axissquare;分析：ex,ey,￡z對(duì)北向速度誤差§Vy均有地球自轉(zhuǎn)周期，舒勒周期分量的影響。其中，￡x對(duì)§Vy的影響較大；￡y,￡z產(chǎn)生的影響幾乎相近。3.ex,ey,￡z對(duì)東向速度誤差§L的影響clc;clearall;t=l:0.01:25;g=98;L=pi/180*39;Ws=2*pi/84.4*60;Wie=2*pi/24;R=g/(Ws)A2;a=WsA2-WieA2;b=sin(L);c=cos(L);d=sin(

6、Wie*t);e=cos(Wie*t);f=sin(Ws*t);h=cos(Ws*t);mcLl=WsA2*0.1/a*(d/Wie-f/Ws)；mcL2=(WsA2*Wie*b/a*(h/WsA2-e/WieA2)+b/Wie)*0.1;mcL3=(WsA2*c*e/Wie/a-Wie*c*h/a-c/Wie)*0.1;plot(t/[mcLl'zmcL2'/mcL3']);title('Ex,Ey,Ez對(duì)mcL的影響')；xlabel(咐間t');ylabel('mcL(t)');legend('Ex-mcL17Ey

7、-mcL27Ez-mcL3,);grid;分析：ex,ey,￡z對(duì)緯度誤差6L均有地球自轉(zhuǎn)周期的影響，ex述會(huì)有舒勒周期分量的影響。4.8x,Sy,對(duì)東向速度誤差§入的影響clc;clearall;t=l:0.01:50;pi=3.14;滬9.8;L=pi/180*39;Ws=2*pi/84.4*60;Wie=2*pi/24;R=g/(Ws)A2;a=WsA2-WieA2;b=sin(L);c=cos(L);d=sin(Wie*t);e=cos(Wie*t);f=sin(Ws*t);h=cos(Ws*t);mcLONGl

8、=(tan(L)/Wie*(l-e)-Wie*c*h/a)*0.1;mcLONG2=(sec(L)*(WsA2-WieA2*cA2)*f/Ws/a-WsA2*tan(L)*b*d/Wie/a-t*c)*0.1;mcLONG3=(WsA2*b*d/Wie/a-WieA2*b*f/a/Ws-t*b)*0.1;plot(t