指北方位平臺(tái)慣導(dǎo)系統(tǒng)導(dǎo)航解算報(bào)告(編輯修改稿)
2025-08-28 19:24 本頁面
【文章內(nèi)容簡(jiǎn)介】 為解算指令角速度圖8 jlfw文件中的數(shù)值五 源程序1 matlab源程序clear allclcclose allload(39。F:\39。)format long。%%%%%%%%%%%%%%%%%%%%%%%%變量的定義及預(yù)處理%%%%%%%%%%%%%%%%%%%%%%%%%%%%%J=zeros(1,60001)。%定義經(jīng)度數(shù)組變量L=zeros(1,60001)。%定義緯度數(shù)組變量a=zeros(1,2)。%定義加速度數(shù)組變量用于存放東向和北向的加速度a2=zeros(1,2)。%用于存放有害加速度v=zeros(1,2)。%用于存放東向和北向的速度w2=zeros(3,60001)。%用于存放結(jié)算出的指令角速度t=0::600。%飛機(jī)飛行時(shí)間erro=zeros(3,60001)。%用于存放東北天方向的指令角速度誤差J(1)=。%初始經(jīng)度L(1)=。%初始緯度Re=6378245。%地球長(zhǎng)軸半徑w1=。%地球自轉(zhuǎn)角速度e=1/%地球的橢球度J(1)=J(1)/(180/pi)。L(1)=L(1)/(180/pi)。%將初始經(jīng)緯度轉(zhuǎn)換為弧度%%%%%%%%%%%%%%%%%%%%%%%%%%進(jìn)行解算%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%for i=1:60001 Rx=Re/(1e*(sin(L(i))^2))。%與子午面垂直平面上的主曲率半徑 Ry=Re/(1+2*e3*e*(sin(L(i))^2))。%當(dāng)?shù)刈游缑嬷髑拾霃? a2(1,1)=2*w1*sin(L(i))*v(2)+(v(1)*v(2)/Rx)*tan(L(i))。%東向有害加速度 a2(2,1)=2*w1*sin(L(i))*v(1)v(1)^2/Rx*tan(L(i))。%北向有害加速度 a(1)=f(1,i)+a2(1,1)。%計(jì)算東向加速度 a(2)=f(2,i)+a2(2,1)。%計(jì)算北向加速度 v(1)=a(1)*+v(1)。 v(2)=a(2)*+v(2)。%速度解算 L(i+1)=L(i)+v(1,2)/Ry*。%緯度計(jì)算J(i+1)=
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號(hào)-1