主要内容gydF4y2Ba

球面调和重力模型gydF4y2Ba

实现行星重力的球面调和表示gydF4y2Ba

  • 库:gydF4y2Ba
  • 航空航天模块/环境/重力gydF4y2Ba

  • 球面调和重力模型块gydF4y2Ba

描述gydF4y2Ba

的gydF4y2Ba球面调和重力模型gydF4y2BaBlock基于行星引力势实现了球面调和行星引力的数学表示。它提供了一种方便的方法来描述行星表面外的球面调和膨胀重力场。gydF4y2Ba

你可以使用球谐来修改球重力的大小和方向(-GM/r)gydF4y2Ba2gydF4y2Ba).最重要或最大的球谐项是第二次纬向谐项J2,它解释了行星的扁率。gydF4y2Ba

如果你想要比球形重力模型更精确的重力值,请使用此块。例如,非大气飞行应用可能需要更高的精度。gydF4y2Ba

限制gydF4y2Ba

  • 该块排除了行星旋转的离心效应,以及进动参考系的效应。gydF4y2Ba

  • 球面调和重力模型适用于径向位置大于行星赤道半径的情况。在行星表面附近或附近的径向位置可能会发生较小的误差。球面调和重力模型不适用于径向位置小于行星表面的情况。gydF4y2Ba

港口gydF4y2Ba

输入gydF4y2Ba

全部展开gydF4y2Ba

从行星中心出发的固定坐标系坐标,指定为一个n × 3矩阵,以选定的单位。矩阵的每一行都是要计算的单独位置。的gydF4y2BazgydF4y2Ba轴是正向北极的。如果gydF4y2Ba中央的身体模型gydF4y2Ba值为gydF4y2BaEGM2008gydF4y2Ba或gydF4y2BaEGM96gydF4y2Ba,这个矩阵包含以地球为中心的地球固定(ECEF)坐标。gydF4y2Ba

当输入一个大的固定框架矩阵和一个高阶值时,您可能会收到内存不足的错误。有关在MATLAB中避免内存不足错误的更多信息gydF4y2Ba®gydF4y2Ba环境,看到gydF4y2Ba解决“内存不足”错误gydF4y2Ba.gydF4y2Ba

当输入一个大型固定框架矩阵时,您可能会收到最大矩阵大小限制。要确定可以在MATLAB环境中为平台创建的最大矩阵或数组,请参见gydF4y2Ba性能和内存gydF4y2Ba.gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

输出gydF4y2Ba

全部展开gydF4y2Ba

的重力值数组gydF4y2BaxgydF4y2Ba设在,gydF4y2BaygydF4y2Ba设在,gydF4y2BazgydF4y2Ba固定框架坐标的轴,以选定的长度单位每秒的平方。矩阵的每一行返回输入矩阵中对应行的计算重力向量。gydF4y2Ba

数据类型:gydF4y2Ba双gydF4y2Ba

参数gydF4y2Ba

全部展开gydF4y2Ba

输入输出单元,具体为:gydF4y2Ba

单位gydF4y2Ba

输入gydF4y2Ba

输出gydF4y2Ba

度量(MKS)gydF4y2Ba

米(m)gydF4y2Ba

米/秒gydF4y2Ba2gydF4y2Ba(米/秒gydF4y2Ba2gydF4y2Ba)gydF4y2Ba

英语gydF4y2Ba

英尺(英尺)gydF4y2Ba

英尺/秒gydF4y2Ba2gydF4y2Ba(英尺/秒gydF4y2Ba2gydF4y2Ba)gydF4y2Ba

编程使用gydF4y2Ba

块参数gydF4y2Ba:gydF4y2Ba单位gydF4y2Ba
类型gydF4y2Ba:特征向量gydF4y2Ba
值gydF4y2Ba:gydF4y2Ba“度量(MKS)”gydF4y2Ba|gydF4y2Ba“英语”gydF4y2Ba
默认的gydF4y2Ba:gydF4y2Ba“度量(MKS)”gydF4y2Ba

超出范围的输入行为,指定为:gydF4y2Ba

价值gydF4y2Ba 描述gydF4y2Ba

没有一个gydF4y2Ba

没有行动。gydF4y2Ba

警告gydF4y2Ba

诊断查看器中的警告,模型模拟继续。gydF4y2Ba

错误gydF4y2Ba

MATLAB返回一个异常,模型仿真停止。gydF4y2Ba

编程使用gydF4y2Ba

块参数gydF4y2Ba:gydF4y2Ba行动gydF4y2Ba
类型gydF4y2Ba:特征向量gydF4y2Ba
值gydF4y2Ba:gydF4y2Ba“没有”gydF4y2Ba|gydF4y2Ba“警告”gydF4y2Ba|gydF4y2Ba“错误”gydF4y2Ba
默认的gydF4y2Ba:gydF4y2Ba“警告”gydF4y2Ba

行星模型,具体为:gydF4y2Ba

中央的身体模型gydF4y2Ba 笔记gydF4y2Ba
EGM2008gydF4y2Ba

地球-是美国国家地理空间情报局(NGA)最新的地球球面调和引力模型。这个块提供了重力模型的WGS-84版本。如果你需要使用地球的旧标准,你可以使用EGM96行星模型。gydF4y2Ba

EGM96gydF4y2Ba 地球gydF4y2Ba
LP100KgydF4y2Ba

月球-根据计算轨道所需的计算时间来确定月球轨道是最好的。该行星模型与LP165P在大约同一年创建,数据类似。gydF4y2Ba

LP165PgydF4y2Ba

月球-是最好的延长月球任务轨道精度。该行星模型与LP100K在大约同一年创建,数据类似。gydF4y2Ba

GMM2BgydF4y2Ba

火星gydF4y2Ba

自定义gydF4y2Ba

允许您指定自己的行星模型。此选项启用gydF4y2Ba中心体MAT-filegydF4y2Ba参数。gydF4y2Ba

EIGENGL04CgydF4y2Ba

地球-支持重力场金宝app模型,EIGEN-GL04C (gydF4y2Bahttp://icgem.gfz-potsdam.de/tom_longtimegydF4y2Ba).这是EIGEN-CG03C的升级版。gydF4y2Ba

有关中心物体固定坐标系的更多信息,请参见gydF4y2Ba算法gydF4y2Ba.gydF4y2Ba

当你定义自己的行星模型时gydF4y2Ba学位gydF4y2Ba参数限制为int16的最大值。当输入较大的度数时,可能会收到内存不足的错误。有关在MATLAB环境中避免内存不足错误的更多信息,请参见gydF4y2Ba解决“内存不足”错误gydF4y2Ba.gydF4y2Ba

依赖关系gydF4y2Ba

将该参数设置为gydF4y2Ba自定义gydF4y2Ba使gydF4y2Ba中心体MAT-filegydF4y2Ba.gydF4y2Ba

编程使用gydF4y2Ba

块参数gydF4y2Ba:gydF4y2BaptypegydF4y2Ba
类型gydF4y2Ba:特征向量gydF4y2Ba
值gydF4y2Ba:gydF4y2Ba“EGM2008”gydF4y2Ba|gydF4y2Ba“EGM96”gydF4y2Ba|gydF4y2Ba“LP100K”gydF4y2Ba|gydF4y2Ba“LP165P”gydF4y2Ba|gydF4y2Ba“GMM2B”gydF4y2Ba|gydF4y2Ba“自定义”gydF4y2Ba|gydF4y2Ba“EIGENGL04C”gydF4y2Ba
默认的gydF4y2Ba:gydF4y2Ba“EGM2008”gydF4y2Ba

谐模型的程度,指定为标量:gydF4y2Ba

中央的身体模型gydF4y2Ba 推荐程度gydF4y2Ba 最大程度上gydF4y2Ba

EGM2008gydF4y2Ba

120gydF4y2Ba

2159gydF4y2Ba

EGM96gydF4y2Ba

70gydF4y2Ba

360gydF4y2Ba

LP100KgydF4y2Ba

60gydF4y2Ba

One hundred.gydF4y2Ba

LP165PgydF4y2Ba

60gydF4y2Ba

165gydF4y2Ba

GMM2BgydF4y2Ba

60gydF4y2Ba

80gydF4y2Ba

EIGENGL04CgydF4y2Ba

70gydF4y2Ba

360gydF4y2Ba

编程使用gydF4y2Ba

块参数gydF4y2Ba:gydF4y2Ba学位gydF4y2Ba
类型gydF4y2Ba:特征向量gydF4y2Ba
值gydF4y2Ba:gydF4y2Ba标量gydF4y2Ba
默认的gydF4y2Ba:gydF4y2Ba“120”gydF4y2Ba

包含自定义行星模型定义的中央天体mat文件。的gydF4y2Baaerogmm2b.matgydF4y2Ba文件在Aerospace Blockset™是自定义行星模型的默认mat -文件。gydF4y2Ba

该文件必须包含:gydF4y2Ba

变量gydF4y2Ba 描述gydF4y2Ba
再保险gydF4y2Ba

行星赤道半径的标量,单位为米(m)。gydF4y2Ba

通用汽车gydF4y2Ba

行星引力参数的标量,单位为米的立方每秒的平方(mgydF4y2Ba3.gydF4y2Ba/秒gydF4y2Ba2gydF4y2Ba)gydF4y2Ba

学位gydF4y2Ba

最大次数的标量。gydF4y2Ba

CgydF4y2Ba

(gydF4y2Ba学位gydF4y2Ba+ 1)————(gydF4y2Ba学位gydF4y2Ba+1)含归一化球面调和系数矩阵,gydF4y2BaCgydF4y2Ba.gydF4y2Ba

年代gydF4y2Ba

(gydF4y2Ba学位gydF4y2Ba+ 1)————(gydF4y2Ba学位gydF4y2Ba+1)含归一化球面调和系数矩阵,gydF4y2Ba年代gydF4y2Ba.gydF4y2Ba

当使用较大的值时gydF4y2Ba学位gydF4y2Ba,您可能会收到内存不足的错误。有关在MATLAB环境中避免内存不足错误的更多信息,请参见gydF4y2Ba解决“内存不足”错误gydF4y2Ba.gydF4y2Ba

依赖关系gydF4y2Ba

要启用该参数,请设置gydF4y2Ba中央的身体模型gydF4y2Ba来gydF4y2Ba自定义gydF4y2Ba.gydF4y2Ba

编程使用gydF4y2Ba

块参数gydF4y2Ba:gydF4y2Ba数据文件gydF4y2Ba
类型gydF4y2Ba:特征向量gydF4y2Ba
值gydF4y2Ba:gydF4y2Ba“aerogmm2b.mat”gydF4y2Ba| MAT-filegydF4y2Ba
默认的gydF4y2Ba:gydF4y2Ba“aerogmm2b.mat”gydF4y2Ba

算法gydF4y2Ba

的gydF4y2Ba球面谐波重力gydF4y2Ba块在中央机构的固定坐标系坐标系统:gydF4y2Ba

  • 地球-固定坐标系是地球中心固定坐标系(ECEF)。gydF4y2Ba

  • 月球-固定坐标系坐标系为主轴系统(PA),由JPL行星星历DE403确定方向。gydF4y2Ba

  • 火星-固定的坐标系是由磁极的旋转方向和原子午线定义的gydF4y2Ba[14]gydF4y2Ba.gydF4y2Ba

参考文献gydF4y2Ba

[1] Gottlieb, Robert G.,“快速重力,重力偏微分,重力归一化,重力梯度扭矩和磁场:推导,代码和数据”。nasa - cr - 188243。休斯顿,德克萨斯州:美国宇航局林登B.约翰逊航天中心,1993年2月。gydF4y2Ba

[2] Vallado,大卫。gydF4y2Ba天体动力学基础与应用gydF4y2Ba.纽约:麦格劳-希尔,1997年。gydF4y2Ba

国防部世界大地测量系统1984,它的定义及其与地方大地测量系统的关系。尼玛TR8350.2。gydF4y2Ba

作者简介:柯诺普里夫,a.s., W. Asmar, E. Carranza, W. l . Sjogren, D.N. Yuan。“作为探月者任务的最新重力模型,”gydF4y2Ba伊卡洛斯gydF4y2Ba, 150,不。1(2001):队。gydF4y2Ba

Lemoine, F. G., D. E. Smith, D.D. Rowlands, M.T. Zuber, G. A. Neumann, d.s. Chinn。“火星全球探测器的火星重力场改进解(GMM-2B)”。gydF4y2Ba地球物理研究杂志gydF4y2Ba106, np E10 (October 25, 2001): pp 23359-23376。gydF4y2Ba

肯扬,J. Factor, N. Pavlis和S. Holmes。"朝向下一个地球引力模型"勘探地球物理学家协会第77届年会,圣安东尼奥,德克萨斯州,2007年9月23-28日。gydF4y2Ba

[7]帕夫利斯,n.k., S.A. Holmes, S.C. Kenyon, J.K. Factor,《2160度地球引力模型:EGM2008》。2008年4月13日至18日在奥地利维也纳举行的欧洲地球科学联盟大会上发表。gydF4y2Ba

[8] Grueber, T.和A. Köhl。“用gps水准和海洋学分析验证EGM2008重力场”。发表于2008年6月23-27日在希腊查尼亚举行的国际重力、大地水准面与地球观测国际研讨会。gydF4y2Ba

[9] Förste, C., Flechtner等,“结合卫星任务和高程/重力测量地面数据的平均全球重力场模型- EIGEN-GL04C”。gydF4y2Ba地球物理研究摘要gydF4y2Ba8、03462、2006。gydF4y2Ba

[10] Hill, k.a.。“天平动点轨道的自主导航”博士论文,科罗拉多大学博尔德分校,2007。gydF4y2Ba

奥斯卡·L·科伦坡球面调和分析的数值方法俄亥俄州立大学地理科学系报告第310号,俄亥俄州哥伦布市。, 1981年3月。gydF4y2Ba

奥斯卡·L·科伦坡"用两颗卫星绘制全球重力地图"荷兰大地测量委员会第7期,第3期,荷兰代尔夫特,1984。,大地测量科学系报告。310号报告。哥伦布:俄亥俄州立大学,1981年3月。gydF4y2Ba

布兰登·琼斯“重力场评估和估计的有效模型”。博士论文,科罗拉多大学,博尔德,2010。gydF4y2Ba

[14]国际测量联合会/国际测量小组地图坐标和旋转要素工作组的报告:1991年。gydF4y2Ba

扩展功能gydF4y2Ba

C / c++代码生成gydF4y2Ba
使用Simulink®Coder™生成C和c++代码。金宝appgydF4y2Ba

介绍了R2010agydF4y2Ba