主要内容gydF4y2Ba

重力球谐gydF4y2Ba

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

折叠所有页面gydF4y2Ba

语法gydF4y2Ba

[gx gy gz] =重力球谐(planet_coordinates)gydF4y2Ba
[gx gy gz] =重力球谐(planet_coordinates,degree)gydF4y2Ba
[gx gy gz] =重力球谐(planet_coordinates,model)gydF4y2Ba
[gx gy gz]=重力球坐标（行星坐标、模型、度数）gydF4y2Ba
[gx gy gz]=重力球坐标（行星坐标、模型、度、作用）gydF4y2Ba
[gx gy gz] = gravitysphericalharmonic(planet_coordinates，'Custom'，degree，{datafile dfreader}，action)gydF4y2Ba

描述gydF4y2Ba

默认的行星模型gydF4y2Ba

例子gydF4y2Ba

［gydF4y2BagxgydF4y2BagygydF4y2Ba广州gydF4y2Ba) = gravitysphericalharmonic (gydF4y2Baplanet_coordinatesgydF4y2Ba）gydF4y2Ba实现了基于行星引力势的球谐行星引力的数学表示。此函数用于计算gydF4y2BaNgydF4y2Ba重力值gydF4y2BaxgydF4y2Ba-轴心国，gydF4y2BaygydF4y2Ba-轴线，以及gydF4y2BazgydF4y2Ba-以行星为中心的行星的轴-固定的行星坐标。函数使用gydF4y2Baplanet_coordinatesgydF4y2Ba,一个gydF4y2Ba米gydF4y2Ba-行星中心行星固定坐标的by-3阵列。gydF4y2Ba

［gydF4y2BagxgydF4y2BagygydF4y2Ba广州gydF4y2Ba) = gravitysphericalharmonic (gydF4y2Baplanet_coordinatesgydF4y2Ba，gydF4y2Ba学位gydF4y2Ba）gydF4y2Ba使用gydF4y2Ba学位gydF4y2Ba指定。gydF4y2Ba

指定行星模型gydF4y2Ba

例子gydF4y2Ba

［gydF4y2BagxgydF4y2BagygydF4y2Ba广州gydF4y2Ba) = gravitysphericalharmonic (gydF4y2Baplanet_coordinatesgydF4y2Ba，gydF4y2Ba模型gydF4y2Ba）gydF4y2Ba实现行星模型的数学表示，gydF4y2Ba模型gydF4y2Ba．gydF4y2Ba

［gydF4y2BagxgydF4y2BagygydF4y2Ba广州gydF4y2Ba) = gravitysphericalharmonic (gydF4y2Baplanet_coordinatesgydF4y2Ba，gydF4y2Ba模型gydF4y2Ba，gydF4y2Ba学位gydF4y2Ba）gydF4y2Ba使用gydF4y2Ba学位gydF4y2Ba指定。gydF4y2Ba模型gydF4y2Ba指定行星模型。gydF4y2Ba

［gydF4y2BagxgydF4y2BagygydF4y2Ba广州gydF4y2Ba) = gravitysphericalharmonic (gydF4y2Baplanet_coordinatesgydF4y2Ba，gydF4y2Ba模型gydF4y2Ba，gydF4y2Ba学位gydF4y2Ba，gydF4y2Ba行动gydF4y2Ba）gydF4y2Ba使用指定的gydF4y2Ba行动gydF4y2Ba当输入超出范围时。gydF4y2Ba

自定义行星模型gydF4y2Ba

例子gydF4y2Ba

［gydF4y2BagxgydF4y2BagygydF4y2Ba广州gydF4y2Ba) = gravitysphericalharmonic (gydF4y2Baplanet_coordinatesgydF4y2Ba，gydF4y2Ba“习俗”gydF4y2Ba，gydF4y2Ba学位gydF4y2Ba，{gydF4y2Ba数据文件读取器gydF4y2Ba},gydF4y2Ba行动gydF4y2Ba）gydF4y2Ba实现自定义模型行星的数学表示。gydF4y2Ba数据文件gydF4y2Ba定义行星模型。gydF4y2BadfreadergydF4y2Ba指定的读取器gydF4y2Ba数据文件gydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

全部崩溃gydF4y2Ba

计算重力场中的重力gydF4y2BaxgydF4y2Ba-地球表面赤道上的轴线。此示例使用默认120度EGM2008模型，并执行默认警告操作。gydF4y2Ba

gx=重力球谐（[-6378.137e3 0]）gydF4y2Ba
gx=9.8143gydF4y2Ba

计算地球南极上空25000米的重力。这个例子使用了带有错误动作的EGM96的70度模型。gydF4y2Ba

[gx，gy，gz]=重力球谐（[0 0-6381.751e3]，gydF4y2Ba‘EGM96’gydF4y2Ba，gydF4y2Ba“错误”gydF4y2Ba）gydF4y2Ba
Gx = 0 Gx = 0 Gx = 0gydF4y2Ba

计算赤道上方15000米和北极上方11000米处的重力。此示例使用30阶GMM2B火星模型，并进行警告操作。gydF4y2Ba

P = [2412.648e3 -2412.648e3 0;0 0 3397.2 e3];[gx, gy, gz] =重力球谐(p，gydF4y2Ba“GMM2B”gydF4y2Ba30岁的gydF4y2Ba“警告”gydF4y2Ba）gydF4y2Ba
gx=-2.6085 0 gy=2.6073 0 gz=0.0000-3.6895gydF4y2Ba

使用带有警告动作的120阶EIGEN-GL04C地球模型计算地球南极上空25000米处的重力。gydF4y2Ba

p=[0 0-6381.751e3]；[gx，gy，gz]=重力球谐（p，gydF4y2Ba'特征GL04C'gydF4y2Ba，gydF4y2Ba．..gydF4y2Ba120，gydF4y2Ba“警告”gydF4y2Ba）gydF4y2Ba
Gx = 0 Gx = 0 Gx = 0gydF4y2Ba

计算赤道上空15000米和北极上空11000米的重力。这个例子使用了一个60度的自定义行星模型，没有动作。gydF4y2Ba

P = [2412.648e3 -2412.648e3 0;0 0 3397年e3];[gx, gy, gz] =重力球谐(p，gydF4y2Ba“习俗”gydF4y2Ba,60,gydF4y2Ba．..gydF4y2Ba｛gydF4y2Ba“GMM2BC80_SHA.txt”gydF4y2Ba@astReadSHAFile},gydF4y2Ba“没有”gydF4y2Ba）gydF4y2Ba
gx=-2.60790Gy=2.60670Gz=0.0002-3.6902gydF4y2Ba

输入参数gydF4y2Ba

全部崩溃gydF4y2Ba

行星坐标，指定为gydF4y2Ba米gydF4y2Ba-by-3行星中心行星固定坐标阵列米。的gydF4y2BazgydF4y2Ba-地轴正对着北极。如果gydF4y2Ba模型gydF4y2Ba是gydF4y2Ba“EGM2008”gydF4y2Ba或gydF4y2Ba‘EGM96’gydF4y2Ba(地球)，行星坐标为ECEF坐标。gydF4y2Ba

数据类型：gydF4y2Ba双重的gydF4y2Ba

行星模型，指定为这些值之一。gydF4y2Ba

行星模型gydF4y2Ba 行星gydF4y2Ba
“EGM2008”gydF4y2Ba

地球引力模型2008。行星坐标为ECEF（WGS84）。gydF4y2Ba
‘EGM96’gydF4y2Ba

地球引力模型1996。行星坐标为ECEF (WGS84)。gydF4y2Ba
“10万卢比”gydF4y2Ba

第100度月球模型。gydF4y2Ba
“LP165P”gydF4y2Ba

165度月球模型。gydF4y2Ba
“GMM2B”gydF4y2Ba

戈达德火星2B型。gydF4y2Ba
“习俗”gydF4y2Ba

自定义行星模型，你定义gydF4y2Ba数据文件gydF4y2Ba．gydF4y2Ba

请注意gydF4y2Ba

要部署自定义行星模型，显式地将自定义数据和读取器文件包含到gydF4y2BaMATLABgydF4y2Ba®gydF4y2Ba编译器™gydF4y2Ba（gydF4y2BamccgydF4y2Ba)编译时的命令。例如：gydF4y2Ba

mcc-m MyCustomSphereCalGravityFunction。。。gydF4y2Ba
-a customDataFile -a customReaderFilegydF4y2Ba

对于其他行星模型，请使用gydF4y2BaMATLAB编译器gydF4y2Ba像往常一样。gydF4y2Ba

有关详细信息，请参阅gydF4y2Ba“习俗”gydF4y2Ba．gydF4y2Ba
'特征GL04C'gydF4y2Ba

EIGEN-GL04C复合地球重力场模型。gydF4y2Ba

当输入一个大的PCPF数组和一个高阶值时，您可能会收到内存不足的错误。有关在MATLAB环境中避免内存不足错误的更多信息，请参见gydF4y2Ba解决“内存不足”错误gydF4y2Ba．gydF4y2Ba

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

数据类型：gydF4y2Ba烧焦gydF4y2Ba|gydF4y2Ba一串gydF4y2Ba

调和重力的度和阶，指定为标量。gydF4y2Ba

行星模型gydF4y2Ba 学位和秩序gydF4y2Ba

“EGM2008”gydF4y2Ba

最大度和阶数为2159。gydF4y2Ba

默认度和顺序为120。gydF4y2Ba

‘EGM96’gydF4y2Ba

最大度数和顺序为360。gydF4y2Ba

默认度和顺序为70。gydF4y2Ba

“10万卢比”gydF4y2Ba

最大程度和秩序是100。gydF4y2Ba

默认度和顺序为60。gydF4y2Ba

“LP165P”gydF4y2Ba

最大度数和顺序为165。gydF4y2Ba

默认度和顺序为60。gydF4y2Ba

“GMM2B”gydF4y2Ba

最大度数和顺序为80。gydF4y2Ba

默认度和顺序为60。gydF4y2Ba

“习俗”gydF4y2Ba

最大度数是默认度数和顺序。有关更多信息，请参阅gydF4y2Ba“习俗”gydF4y2Ba．gydF4y2Ba

'特征GL04C'gydF4y2Ba

最大度数和顺序为360。gydF4y2Ba

默认度和顺序为70。gydF4y2Ba

当输入一个大的PCPF数组和一个高阶值时，您可能会收到内存不足的错误。有关在MATLAB环境中避免内存不足错误的更多信息，请参见gydF4y2Ba性能和内存gydF4y2Ba．gydF4y2Ba

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

数据类型：gydF4y2Ba烧焦gydF4y2Ba|gydF4y2Ba一串gydF4y2Ba

自定义行星模型定义，指定为gydF4y2Ba“习俗”gydF4y2Ba．用定义数据文件和相应的阅读器指定行星模型定义。有关详细信息，请参阅gydF4y2Ba数据文件读取器gydF4y2Ba．gydF4y2Ba

数据类型：gydF4y2Ba烧焦gydF4y2Ba|gydF4y2Ba一串gydF4y2Ba

自定义行星模型定义文件和读取器，指定为矢量。gydF4y2Ba数据文件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+ 1)————(gydF4y2Ba学位gydF4y2Ba+1)含归一化球面调和系数矩阵，gydF4y2Ba年代gydF4y2Ba

阅读gydF4y2Ba数据文件gydF4y2Ba中指定一个MATLAB函数gydF4y2BadfreadergydF4y2Ba参数。指定的读取器文件取决于文件类型gydF4y2Ba数据文件gydF4y2Ba．gydF4y2Ba

数据文件类型gydF4y2Ba 描述gydF4y2Ba

MATLAB文件gydF4y2Ba

指定MATLABgydF4y2Ba负载gydF4y2Ba例如,函数gydF4y2Ba@loadgydF4y2Ba．gydF4y2Ba

其他文件类型gydF4y2Ba

指定一个自定义的MATLAB阅读器函数。有关自定义阅读器函数的示例，请参见gydF4y2BaastReadSHAFile.mgydF4y2Ba和gydF4y2BaastReadEGMFile.mgydF4y2Ba. 注意这些文件中的输出变量顺序。gydF4y2Ba

例子：gydF4y2Ba{'GMM2BC80_SHA.txt'@astReadSHAFile}gydF4y2Ba

数据类型：gydF4y2Ba双重的gydF4y2Ba

超出量程输入的动作，指定为:gydF4y2Ba

  • “错误”gydF4y2Ba-显示警告并指示输入超出范围。gydF4y2Ba

  • “警告”gydF4y2Ba-显示错误并指示输入超出范围。gydF4y2Ba

  • “没有”gydF4y2Ba—不显示警告或错误信息。gydF4y2Ba

数据类型：gydF4y2Ba烧焦gydF4y2Ba|gydF4y2Ba一串gydF4y2Ba

输出参数gydF4y2Ba

全部崩溃gydF4y2Ba

中的重力值gydF4y2BaxgydF4y2Ba-以行星为中心的行星固定坐标的轴，返回为数组gydF4y2Ba米gydF4y2Ba重力值的单位是米每秒平方(m/s)gydF4y2Ba2gydF4y2Ba).gydF4y2Ba

中的重力值gydF4y2BaygydF4y2Ba-以行星为中心的行星固定坐标的轴，返回为数组gydF4y2Ba米gydF4y2Ba重力值的单位是米每秒平方(m/s)gydF4y2Ba2gydF4y2Ba).gydF4y2Ba

中的重力值gydF4y2BazgydF4y2Ba-以行星为中心的行星固定坐标的轴，返回为数组gydF4y2Ba米gydF4y2Ba重力值的单位是米每秒平方(m/s)gydF4y2Ba2gydF4y2Ba).gydF4y2Ba

局限性gydF4y2Ba

  • 该函数不包括行星旋转的离心效应和进动参考系的效应。gydF4y2Ba

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

提示gydF4y2Ba

  • 当输入一个大的PCPF数组和一个高阶值时，您可能会收到内存不足的错误。有关在MATLAB环境中避免内存不足错误的更多信息，请参见gydF4y2Ba性能和内存gydF4y2Ba．gydF4y2Ba

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

参考文献gydF4y2Ba

戈特利布，r.g。快速重力，重力偏微分，重力归一化，重力梯度扭矩和磁场:推导，代码和数据。gydF4y2Ba技术报告NASA承包商报告188243gydF4y2Ba. 休斯顿：美国宇航局林登·B·约翰逊航天中心，1993年2月。gydF4y2Ba

[2] 巴尔拉多，大卫A。gydF4y2Ba天体动力学基础与应用gydF4y2Ba．纽约:麦格劳-希尔，1997年。gydF4y2Ba

[3] 国防测绘局。gydF4y2Ba国防部1984年世界大地测量系统，它的定义及其与地方大地测量系统的关系gydF4y2Ba．TR 8350.2, 2版。费尔法克斯，弗吉尼亚州:DMA, 1991年9月1日。gydF4y2Ba

A. S. S. Asmar, E. Carranza, W. L. Sjogen, D. N. Yuan。“月球探勘者任务的最新重力模型，伊卡洛斯”，第150期。1(2001):队。gydF4y2Ba

[5] 莱蒙，F.G.，D.E.史密斯，D.D.罗兰兹，M.T.祖伯，G.A.纽曼和D.S.钦。“火星全球勘测员提供的火星重力场（GMM-2B）的改进解决方案。”gydF4y2Ba地球物理研究杂志gydF4y2Ba第106号E10（2001年10月25日）：23359-23376。gydF4y2Ba

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

[7] 帕夫利斯，N.K.，S.A.霍姆斯，S.C.肯扬和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·科伦坡gydF4y2Ba球面调和分析的数值方法gydF4y2Ba．报告没有。310.哥伦布:俄亥俄州立大学大地科学系，1981年。gydF4y2Ba

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

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

［14］gydF4y2Ba国际测量联盟/国际测量组制图座标和旋转要素工作组的报告:1991年gydF4y2Ba．gydF4y2Ba

另请参阅gydF4y2Ba

|gydF4y2Ba|gydF4y2Ba

R2010a中引入gydF4y2Ba

航空航天工具箱文档gydF4y2Ba

金宝app

仿真软件MAT金宝appLAB,その他の製品をお試しくださいgydF4y2Ba

無料評価版のダウンロードgydF4y2Ba