radicurv.m

来自「惯性导航系统传递对准系统模型程序」· M 代码 · 共 15 行

function [Rm,Rp] = radicurv(lat)

if nargin<1,error('insufficient number of input arguments'),end

EQUA_RADIUS = 6378137.0;    % equatorial radius of the earth; WGS-84
ECCENTRICITY = 0.0818191908426;  % eccentricity of the earth ellipsoid

e2 = ECCENTRICITY^2;
den = 1-e2*(sin(lat))^2;
Rm = (EQUA_RADIUS*(1-e2))/( (den)^(3/2) );%地球短半轴半径。
Rp = EQUA_RADIUS/( sqrt(den) );%地球长半轴半径。
%该函数计算地球椭球模型的长半轴和短半轴半径。该函数的输入量是地球半径EQUA_RADIUS、地球扁心率ECCENTRICITY，计算出地球长半轴半径
%和短半轴半径