usercontrail.m

来自「RAIM 完好性分析」· M 代码 · 共 18 行

function [long,lat,h]=usercontrail(t)

%*************************************************** 
%*** 用经纬高表示的用户轨迹
%*** 用户以经度116,纬度40为初始坐标飞机轨迹为平行赤道轨迹，速度300m/s
%*************************************************** 

v = 0;
Re=6378137;
f=1/298.257; 
trans=2*pi/360;
L=40;
RN=Re./sqrt(((1-f).^2).*(sin(trans*L))^2+(cos(trans*L)).^2);
R=RN*cos(trans*L);
w=v/R;
long=(w*t)/trans+116;
lat=40;
h=500;