高阶次勒让德计算（跨阶次递推公式）

跨阶次递推勒让德公式

今天卫星大地测量课的一个问题
问题：已知某一点的大地坐标BLH，求该点的高程异常。

（1）从ICGEM网站获取模型数据，并提取l,m,c,s数据。
（2）在c位系数中减去正常重力位系数
（3）将给定点的大地坐标转化为地心坐标
（4）计算给定点处的平均正常重力
（5）Legendre递推公式的计算(标准前向列递推公式)（跨阶次递推公式）
（6）计算给定点处高程异常

matlab代码：

clc;clear;
% B = input('输入大地纬度B: ');
% L = input('输入大地经度L：');
% H = input('输入大地高： ');
H = 0;
GM = 0.3986004415E+15;
R = 0.6378136300E+07;
Degree = 360;
e2=0.00669437999013;
ge=9.7803253359;
k1=0.00193185265246;
Alpha=1/298.257223563;
Beta=0.00344978650684;
Pi=3.1415926535897932384626433;
% B = B*Pi/180;
% L = L*Pi/180;

%读取LMCS数据（读取前删除数据的前20行介绍）
path = 'C:\Users\Administrator\Desktop\EGM2008.gfc';
Data = importdata(pa