主要内容孔侑dF4y2Ba

gravitysphericalharmonic孔侑dF4y2Ba

实现行星重力的球面调和表示孔侑dF4y2Ba

折叠所有页面孔侑dF4y2Ba

语法孔侑dF4y2Ba

[gx gy gz] =重力球谐(planet_coordinates)孔侑dF4y2Ba
[gx gy gz] =重力球谐(planet_coordinates,degree)孔侑dF4y2Ba
[gx gy gz] =重力球谐(planet_coordinates,model)孔侑dF4y2Ba
[gx gy gz] =重力球谐(planet_座标，模型，度)孔侑dF4y2Ba
[gx gy gz] =重力球谐(planet_坐标，模型，度，作用)孔侑dF4y2Ba
[gx gy gz] = gravitysphericalharmonic(planet_coordinates，'Custom'，degree，{datafile dfreader}，action)孔侑dF4y2Ba

描述孔侑dF4y2Ba

默认的行星模型孔侑dF4y2Ba

例子孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba) = gravitysphericalharmonic (孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba）孔侑dF4y2Ba实现了基于行星引力势的球面调和行星引力的数学表示。这个函数计算的数组孔侑dF4y2BaN孔侑dF4y2Ba重力值孔侑dF4y2Bax孔侑dF4y2Ba设在,孔侑dF4y2Bay孔侑dF4y2Ba设在,孔侑dF4y2Baz孔侑dF4y2Ba-以行星为中心的行星的轴-固定的行星坐标。函数使用孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba,一个孔侑dF4y2Ba米孔侑dF4y2Ba-by-3行星中心行星固定坐标数组。孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba) = gravitysphericalharmonic (孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba学位孔侑dF4y2Ba）孔侑dF4y2Ba使用程度和顺序孔侑dF4y2Ba学位孔侑dF4y2Ba指定。孔侑dF4y2Ba

指定的行星模型孔侑dF4y2Ba

例子孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba) = gravitysphericalharmonic (孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba模型孔侑dF4y2Ba）孔侑dF4y2Ba实现了行星模型的数学表示，孔侑dF4y2Ba模型孔侑dF4y2Ba．孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba) = gravitysphericalharmonic (孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba模型孔侑dF4y2Ba，孔侑dF4y2Ba学位孔侑dF4y2Ba）孔侑dF4y2Ba使用程度和顺序孔侑dF4y2Ba学位孔侑dF4y2Ba指定。孔侑dF4y2Ba模型孔侑dF4y2Ba指定行星模型。孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba) = gravitysphericalharmonic (孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba模型孔侑dF4y2Ba，孔侑dF4y2Ba学位孔侑dF4y2Ba，孔侑dF4y2Ba行动孔侑dF4y2Ba）孔侑dF4y2Ba使用指定的孔侑dF4y2Ba行动孔侑dF4y2Ba当输入超出范围时。孔侑dF4y2Ba

自定义行星模型孔侑dF4y2Ba

例子孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba) = gravitysphericalharmonic (孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba“自定义”孔侑dF4y2Ba，孔侑dF4y2Ba学位孔侑dF4y2Ba，{孔侑dF4y2Ba丢失dfreader孔侑dF4y2Ba}，孔侑dF4y2Ba行动孔侑dF4y2Ba）孔侑dF4y2Ba实现自定义模型行星的数学表示。孔侑dF4y2Ba数据文件孔侑dF4y2Ba定义行星模型。孔侑dF4y2Badfreader孔侑dF4y2Ba指定读取器孔侑dF4y2Ba数据文件孔侑dF4y2Ba．孔侑dF4y2Ba

例子孔侑dF4y2Ba

全部折叠孔侑dF4y2Ba

计算重力孔侑dF4y2Bax孔侑dF4y2Ba在地球表面的赤道上。本例使用EGM2008默认的120度模型，并带有默认的警告动作。孔侑dF4y2Ba

Gx =重力球谐([-6378.137e3 0 0])孔侑dF4y2Ba
gx = 9.8143孔侑dF4y2Ba

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

[gx, gy, gz] =重力球谐([0 0 -6381.751e3]，孔侑dF4y2Ba“EGM96”孔侑dF4y2Ba，孔侑dF4y2Ba“错误”孔侑dF4y2Ba）孔侑dF4y2Ba
Gx = 0 Gx = 0 Gx = 0孔侑dF4y2Ba

计算赤道上空15000米和北极上空11000米的重力。本例使用一个带有警告动作的30阶GMM2B Mars模型。孔侑dF4y2Ba

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

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

P = [0 0 -6381.751e3];[gx, gy, gz] =重力球谐(p，孔侑dF4y2Ba“EIGENGL04C”孔侑dF4y2Ba，孔侑dF4y2Ba．..孔侑dF4y2Ba120年,孔侑dF4y2Ba“警告”孔侑dF4y2Ba）孔侑dF4y2Ba
Gx = 0 Gx = 0 Gx = 0孔侑dF4y2Ba

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

P = [2412.648e3 -2412.648e3 0;0 0 3397年e3];[gx, gy, gz] =重力球谐(p，孔侑dF4y2Ba“自定义”孔侑dF4y2Ba现年60岁的孔侑dF4y2Ba．..孔侑dF4y2Ba｛孔侑dF4y2Ba“GMM2BC80_SHA.txt”孔侑dF4y2Ba@astReadSHAFile},孔侑dF4y2Ba“没有”孔侑dF4y2Ba）孔侑dF4y2Ba
Gx = -2.6079 0 gy = - 2.6067 0 gz = 0.0002 -3.6902孔侑dF4y2Ba

输入参数孔侑dF4y2Ba

全部折叠孔侑dF4y2Ba

行星坐标，指定为孔侑dF4y2Ba米孔侑dF4y2Ba-by-3行星中心行星固定坐标阵列米。的孔侑dF4y2Baz孔侑dF4y2Ba轴是正向北极的。如果孔侑dF4y2Ba模型孔侑dF4y2Ba是孔侑dF4y2Ba“EGM2008”孔侑dF4y2Ba或孔侑dF4y2Ba“EGM96”孔侑dF4y2Ba(地球)，行星坐标为ECEF坐标。孔侑dF4y2Ba

数据类型:孔侑dF4y2Ba双孔侑dF4y2Ba

行星模型，指定为这些值之一。孔侑dF4y2Ba

行星模型孔侑dF4y2Ba 地球孔侑dF4y2Ba
“EGM2008”孔侑dF4y2Ba

地球引力模型2008。行星坐标为ECEF (WGS84)。孔侑dF4y2Ba
“EGM96”孔侑dF4y2Ba

地球引力模型1996。行星坐标为ECEF (WGS84)。孔侑dF4y2Ba
“LP100K”孔侑dF4y2Ba

第100度月球模型。孔侑dF4y2Ba
“LP165P”孔侑dF4y2Ba

第165度月亮模型。孔侑dF4y2Ba
“GMM2B”孔侑dF4y2Ba

戈达德火星2B模型。孔侑dF4y2Ba
“自定义”孔侑dF4y2Ba

自定义行星模型，你定义孔侑dF4y2Ba数据文件孔侑dF4y2Ba．孔侑dF4y2Ba

请注意孔侑dF4y2Ba

要部署自定义行星模型，显式地将自定义数据和读取器文件包含到孔侑dF4y2BaMATLAB孔侑dF4y2Ba®孔侑dF4y2Ba编译器™孔侑dF4y2Ba（孔侑dF4y2Ba世纪挑战集团孔侑dF4y2Ba)命令。例如:孔侑dF4y2Ba

mcc - m mycustomsphericalgravityfunction……孔侑dF4y2Ba
-a customDataFile -a customReaderFile孔侑dF4y2Ba

对于其他行星模型，使用孔侑dF4y2BaMATLAB编译器孔侑dF4y2Ba像往常一样。孔侑dF4y2Ba

有关更多信息，请参见孔侑dF4y2Ba“自定义”孔侑dF4y2Ba．孔侑dF4y2Ba
“EIGENGL04C”孔侑dF4y2Ba

EIGEN-GL04C复合地球重力场模型。孔侑dF4y2Ba

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

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

数据类型:孔侑dF4y2Ba字符孔侑dF4y2Ba|孔侑dF4y2Ba字符串孔侑dF4y2Ba

调和重力的度和阶，指定为标量。孔侑dF4y2Ba

行星模型孔侑dF4y2Ba 学位和秩序孔侑dF4y2Ba

“EGM2008”孔侑dF4y2Ba

最大度和阶数为2159。孔侑dF4y2Ba

默认度和顺序是120。孔侑dF4y2Ba

“EGM96”孔侑dF4y2Ba

最大程度和顺序是360。孔侑dF4y2Ba

默认度和顺序为70。孔侑dF4y2Ba

“LP100K”孔侑dF4y2Ba

最大程度和秩序是100。孔侑dF4y2Ba

默认级别和顺序为60。孔侑dF4y2Ba

“LP165P”孔侑dF4y2Ba

最大度和阶数为165。孔侑dF4y2Ba

默认级别和顺序为60。孔侑dF4y2Ba

“GMM2B”孔侑dF4y2Ba

最大程度和顺序是80。孔侑dF4y2Ba

默认级别和顺序为60。孔侑dF4y2Ba

“自定义”孔侑dF4y2Ba

最大程度是默认程度和顺序。有关更多信息，请参见孔侑dF4y2Ba“自定义”孔侑dF4y2Ba．孔侑dF4y2Ba

“EIGENGL04C”孔侑dF4y2Ba

最大程度和顺序是360。孔侑dF4y2Ba

默认度和顺序为70。孔侑dF4y2Ba

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

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

数据类型:孔侑dF4y2Ba字符孔侑dF4y2Ba|孔侑dF4y2Ba字符串孔侑dF4y2Ba

自定义行星模型定义，指定为孔侑dF4y2Ba“自定义”孔侑dF4y2Ba．用定义数据文件和相应的阅读器指定行星模型定义。有关更多信息，请参见孔侑dF4y2Ba丢失dfreader孔侑dF4y2Ba．孔侑dF4y2Ba

数据类型:孔侑dF4y2Ba字符孔侑dF4y2Ba|孔侑dF4y2Ba字符串孔侑dF4y2Ba

自定义行星模型定义文件和阅读器，指定为矢量。孔侑dF4y2Ba数据文件孔侑dF4y2Ba必须包含这些变量。孔侑dF4y2Ba

变量孔侑dF4y2Ba 描述孔侑dF4y2Ba
再保险孔侑dF4y2Ba

行星赤道半径标量(米)孔侑dF4y2Ba
通用汽车孔侑dF4y2Ba

行星引力参数的标量，单位为米的立方每秒的平方(m孔侑dF4y2Ba3.孔侑dF4y2Ba/秒孔侑dF4y2Ba2孔侑dF4y2Ba）孔侑dF4y2Ba
学位孔侑dF4y2Ba

最大次数标量孔侑dF4y2Ba
C孔侑dF4y2Ba

（孔侑dF4y2Ba学位孔侑dF4y2Ba+ 1)————(孔侑dF4y2Ba学位孔侑dF4y2Ba+1)含归一化球面调和系数矩阵，孔侑dF4y2BaC孔侑dF4y2Ba
年代孔侑dF4y2Ba

（孔侑dF4y2Ba学位孔侑dF4y2Ba+ 1)————(孔侑dF4y2Ba学位孔侑dF4y2Ba+1)含归一化球面调和系数矩阵，孔侑dF4y2Ba年代孔侑dF4y2Ba

阅读孔侑dF4y2Ba数据文件孔侑dF4y2Ba中指定一个MATLAB函数孔侑dF4y2Badfreader孔侑dF4y2Ba参数。指定的读取器文件取决于文件类型孔侑dF4y2Ba数据文件孔侑dF4y2Ba．孔侑dF4y2Ba

数据文件类型孔侑dF4y2Ba 描述孔侑dF4y2Ba

MATLAB文件孔侑dF4y2Ba

指定MATLAB孔侑dF4y2Ba负载孔侑dF4y2Ba例如,函数孔侑dF4y2Ba@load孔侑dF4y2Ba．孔侑dF4y2Ba

其他文件类型孔侑dF4y2Ba

指定一个自定义的MATLAB阅读器函数。有关自定义阅读器函数的示例，请参见孔侑dF4y2BaastReadSHAFile.m孔侑dF4y2Ba和孔侑dF4y2BaastReadEGMFile.m孔侑dF4y2Ba．请注意这些文件中的输出变量顺序。孔侑dF4y2Ba

例子:孔侑dF4y2Ba{“GMM2BC80_SHA.txt”@astReadSHAFile}孔侑dF4y2Ba

数据类型:孔侑dF4y2Ba双孔侑dF4y2Ba

超出量程输入的动作，指定为:孔侑dF4y2Ba

  • “错误”孔侑dF4y2Ba—显示警告信息，表示输入超出范围。孔侑dF4y2Ba

  • “警告”孔侑dF4y2Ba—显示错误，表示输入超出范围。孔侑dF4y2Ba

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

数据类型:孔侑dF4y2Ba字符孔侑dF4y2Ba|孔侑dF4y2Ba字符串孔侑dF4y2Ba

输出参数孔侑dF4y2Ba

全部折叠孔侑dF4y2Ba

重力值孔侑dF4y2Bax孔侑dF4y2Ba-以行星为中心的行星固定坐标的轴，返回为数组孔侑dF4y2Ba米孔侑dF4y2Ba重力值的单位是米每秒平方(m/s)孔侑dF4y2Ba2孔侑dF4y2Ba）.孔侑dF4y2Ba

重力值孔侑dF4y2Bay孔侑dF4y2Ba-以行星为中心的行星固定坐标的轴，返回为数组孔侑dF4y2Ba米孔侑dF4y2Ba重力值的单位是米每秒平方(m/s)孔侑dF4y2Ba2孔侑dF4y2Ba）.孔侑dF4y2Ba

重力值孔侑dF4y2Baz孔侑dF4y2Ba-以行星为中心的行星固定坐标的轴，返回为数组孔侑dF4y2Ba米孔侑dF4y2Ba重力值的单位是米每秒平方(m/s)孔侑dF4y2Ba2孔侑dF4y2Ba）.孔侑dF4y2Ba

限制孔侑dF4y2Ba

  • 该函数排除了行星旋转的离心效应和进动参考系的效应。孔侑dF4y2Ba

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

提示孔侑dF4y2Ba

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

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

参考文献孔侑dF4y2Ba

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

David A. Vallado孔侑dF4y2Ba天体动力学基础与应用孔侑dF4y2Ba．纽约:麦格劳-希尔，1997年。孔侑dF4y2Ba

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

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

Lemoine, F. G.， D. E. Smith, D. D. Rowlands, M. T. Zuber, G. A. Neumann, D. S. Chinn。“火星全球探测器的火星重力场改进解(GMM-2B)”。孔侑dF4y2Ba地球物理研究杂志孔侑dF4y2Ba106年,没有。E10(2001年10月25日):23359-23376。孔侑dF4y2Ba

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

帕夫利斯，n.k.， s.a. Holmes, s.c. Kenyon和j.k. Factor。“2160度地球引力模型:EGM2008”。2008年4月13-18日在奥地利维也纳欧洲地球科学联盟大会上发表的论文。孔侑dF4y2Ba

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

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

[10] Hill, k.a.。“天平动点轨道的自主导航”博士羞辱。科罗拉多大学博尔德分校，2007年。孔侑dF4y2Ba

奥斯卡·L·科伦坡孔侑dF4y2Ba球面调和分析的数值方法孔侑dF4y2Ba．报告没有。310.哥伦布:俄亥俄州立大学大地科学系，1981年。孔侑dF4y2Ba

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

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

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

另请参阅孔侑dF4y2Ba

|孔侑dF4y2Ba|孔侑dF4y2Ba

介绍了R2010a孔侑dF4y2Ba

航空航天工具箱文档孔侑dF4y2Ba

万博1manbetx

尝试MATLAB, Si万博1manbetxmulink和其他产品s manbetx 845孔侑dF4y2Ba

得到审判现在孔侑dF4y2Ba