组合导航卯酉圈

ref:https://baike.baidu.com/item/%E5%8D%AF%E9%85%89%E5%9C%88/3033623?fr=aladdin

卯酉圈(Prime Vertical)就是指地平坐标系中的大圆。即与子午圈相垂直的地平经圈 ,它与地平圈相交于东点和西点。

过椭球面上一点的法线,可作无限个法截面,其中一个与该点子午面相垂直的法截面同椭球面相截形成的闭合的圈称为卯酉圈。在下图1中,PEE′即为过点P的卯酉圈 [1] 

图1卯酉圈示意图

卯酉圈是一个非常拗口的名称,该名称的由来可参考下图2所示的罗盘简化图。卯酉其实就是东西向的意思。卯酉圈就是经过法线、东西走向的平面与椭球面的交线形成的闭合圈。

图2 罗盘简化图

卯酉圈的曲率半径用N表示。

为了推求N的计算公式,如图1所示,过P点作以O'为中心的平行圈PHK的切线PT,该切线位于垂直于子午面的平行圈平面内。因卯酉圈也垂直于子午面,故PT也是卯酉圈在P点处的切线,即PT垂直于Pn。所以PT是平行圈PHK及卯酉圈PEE'在P点处的公切线。

由麦尼尔定理知,假设通过曲面上一点引两条截弧,一为法截弧,一为斜截弧,且在该点上这两条截弧具有公共切线,这时斜截弧在该点处的曲率半径等于法截弧的曲率半径乘以两截弧平面夹角的余弦。

子午圈曲率半径M=N^3*b^2/a^4

close all clear clc disp('***** 基于EKF的位置速度观测组合导航程序 *****'); disp('Step1:加载数据;'); load IMU_data200.mat %惯导原始数据 load Reference_data.mat %GPS测量数据 disp('Step2:初始化参数;'); %% 一些导航参数常数项 WIE = 7.292115e-5; % 地球自转角速度 r0 = 6378137.0; % 地球半径 EE = 0.0818191908426; % 偏心率 d2r = pi/180; % degree to radian r2d = 180/pi; % radian to degree dh2rs = d2r/3600; % deg/h to rad/s %% 导航坐标系下初始化姿态,速度,位置 yaw = (0)*pi/180;%航向角 pitch = 0*pi/180;%俯仰角 roll = 0*pi/180;%滚动角 cbn=eul2dcm(roll,pitch,yaw); cnb=cbn'; q=dcm2quat(cbn)'; Vn=0;%北向速度 Ve=0;%东向速度 Vd=0;%地向速度 V_last=[Vn Ve Vd]'; Lati = 31.4913627505302*pi/180;%纬度 Longi= 120.849577188492*pi/180;%经度 Alti = 6.6356;%高度 sampt0=1/200;%惯导系统更新时间 Rn = r0*(1-EE^2)/(1-EE^2*(sin(Lati))^2)^1.5; %子午圈曲率半径 Re = r0/(1-EE^2*(sin(Lati))^2)^0.5; %卯酉曲率半径 g_u = -9.7803267711905*(1+0.00193185138639*sin(Lati)^2)... /((1-0.00669437999013*sin(Lati)^2)^0.5 *(1.0 + Alti/r0)^2); g = [0 0 -g_u]';%重力 g0=9.80665; %% 卡尔曼滤波P、Q、R设置 % P的设置 std_roll = (5)*d2r; std_pitch = (5)*d2r; std_yaw = (60)*d2r; std_vel = 0.1; std_pos = 5; std_gyro = 3*0.5*dh2rs; % 陀螺随机漂移0.5度/小时 std_acc = 3*0.15e-3*g0; % 加表零偏0.15mg Pfilter = diag([std_roll^2 std_pitch^2 std_yaw^2 std_vel^2 std_vel^2 std_vel^2 (std_pos/3600/30/57.3)^2 (std_pos/3600/30/57.3)^2 std_pos^2 std_gyro^2 std_gyro^2 std_gyro^2 std_acc^2 std_acc^2 std_acc^2]); % Q的设置 std_Wg = 0.15*(2.909*1e-4); % 陀螺漂移噪声,度/根号小时转化成rad/根号秒 std_Wa = 0.21/60/3; % 加表漂移噪声 Qkf = diag([std_Wg^2 std_Wg^2 std_Wg^2 std_Wa^2 std_Wa^2 std_Wa^2]); G = zeros(15, 6); F = zeros(15); F_i=zeros(9,9); F_s=zeros(9,6); H = zeros(6,15); H(1:3,4:6) = eye(3); H(4:6,7:9) = eye(3); % R的设置 R = diag([std_vel^2 std_vel^2 std_vel^2 (std_pos/3600/30/57.3)^2 (std_pos/3600/30/57.3)^2 (std_pos)^2]);
05-29
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

hjwang1

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值