主要内容孔侑dF4y2Ba

gravitysphericalharmonic孔侑dF4y2Ba

实现行星重力的球面谐波表示孔侑dF4y2Ba

折叠所有页面孔侑dF4y2Ba

语法孔侑dF4y2Ba

[gx gy gz] =重力球面谐波(planet_coordinates)孔侑dF4y2Ba
[gx gy gz] =重力球面谐波(行星坐标，度)孔侑dF4y2Ba
[gx gy gz] = gravitysphericalharmonic(planet_coordinates,model)孔侑dF4y2Ba
[gx gy gz] =重力球面谐波(行星坐标，模型，度)孔侑dF4y2Ba
[gx gy gz] = gravitysphericalharmonic(行星坐标，模型，程度，动作)孔侑dF4y2Ba
[gx gy gz] = gravitysphericalharmonic(planet_coordinates，'Custom'，degree，{数据文件dfreader}，action)孔侑dF4y2Ba

描述孔侑dF4y2Ba

默认行星模型孔侑dF4y2Ba

例子孔侑dF4y2Ba

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

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba] =地心引力孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba学位孔侑dF4y2Ba）孔侑dF4y2Ba使用的程度和顺序孔侑dF4y2Ba学位孔侑dF4y2Ba指定。孔侑dF4y2Ba

指定行星模型孔侑dF4y2Ba

例子孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba] =地心引力孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba模型孔侑dF4y2Ba）孔侑dF4y2Ba实现了行星模型的数学表示，孔侑dF4y2Ba模型孔侑dF4y2Ba．孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba] =地心引力孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba模型孔侑dF4y2Ba，孔侑dF4y2Ba学位孔侑dF4y2Ba）孔侑dF4y2Ba使用的程度和顺序孔侑dF4y2Ba学位孔侑dF4y2Ba指定。孔侑dF4y2Ba模型孔侑dF4y2Ba指定行星模型。孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba] =地心引力孔侑dF4y2Baplanet_coordinates孔侑dF4y2Ba，孔侑dF4y2Ba模型孔侑dF4y2Ba，孔侑dF4y2Ba学位孔侑dF4y2Ba，孔侑dF4y2Ba行动孔侑dF4y2Ba）孔侑dF4y2Ba使用指定的孔侑dF4y2Ba行动孔侑dF4y2Ba当输入超出范围时。孔侑dF4y2Ba

自定义行星模型孔侑dF4y2Ba

例子孔侑dF4y2Ba

［孔侑dF4y2Bagx孔侑dF4y2Ba孔侑孔侑dF4y2Ba广州孔侑dF4y2Ba] =地心引力孔侑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 = gravitysphericalharmonic([-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 gy = 0 gz = 9.7552孔侑dF4y2Ba

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

P = [2412.648e3 -2412.648e3 0;0 0 3397.2e3];[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 gy = 0 gz = 9.7552孔侑dF4y2Ba

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

P = [2412.648e3 -2412.648e3 0;0 0 3397e3];[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

地球重力场组合模型eigenen - 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

行星赤道半径标量，单位为米(m)孔侑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

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

[2]瓦拉多，大卫A。孔侑dF4y2Ba天体动力学基础及应用“，孔侑dF4y2Ba．纽约:麦格劳-希尔，1997年。孔侑dF4y2Ba

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

[4]科诺普利夫、A. S.、S. W.阿斯马尔、E.卡兰扎、W. L.斯约根、袁东宁。“最近的引力模型作为月球勘探者任务的结果，伊卡洛斯”150,no。1(2001): 1 - 18。孔侑dF4y2Ba

[5]莱莫因、F. G.、D. E.史密斯、D. D.罗兰、M. T.祖伯、G. A.诺伊曼和D. S.秦。“来自火星全球探测者的火星重力场(GMM-2B)的改进解决方案。”孔侑dF4y2Ba地球物理研究杂志孔侑dF4y2Ba106年,没有。E10(2001年10月25日):23359-23376。孔侑dF4y2Ba

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

[7]帕夫利斯，N.K, S. A.霍姆斯，S. C.凯尼恩和J. K. Factor。2160度的地球引力模型:EGM2008在欧洲地球科学联盟大会上发表的论文，维也纳，奥地利，2008年4月13-18日。孔侑dF4y2Ba

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

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

[10]希尔，堪萨斯州。"振动点轨道的自主导航"博士羞辱。科罗拉多大学博尔德分校，2007年。孔侑dF4y2Ba

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

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

[13]琼斯，布兰登A。重力场评价和估计的有效模型博士羞辱。科罗拉多大学博尔德分校，2010年。孔侑dF4y2Ba

[14]孔侑dF4y2BaIAU/IAG关于地图坐标和旋转元件工作组的报告:1991年孔侑dF4y2Ba．孔侑dF4y2Ba

另请参阅孔侑dF4y2Ba

|孔侑dF4y2Ba|孔侑dF4y2Ba

在R2010a中引入孔侑dF4y2Ba

航天工具箱文档孔侑dF4y2Ba

金宝app

试试MATLAB、Sim金宝appulink和其他产品下载188bet金宝搏孔侑dF4y2Ba

现在就去审判吧孔侑dF4y2Ba