果然程序還是得多看，多編，然后融會貫通。

目前主要有CSR、GFZ、JPL三大機構(gòu)提供GRACE球諧系數(shù)，ICGEM網(wǎng)站提供的都是以gfc為后綴的文件。

GSM-2_2002095-2002120_GRAC_UTCSR_BA01_0600.gfc

下面的代碼可以讀取gfc文件，并轉(zhuǎn)成mat文件。

點擊生成的mat文件，比如Cnm系數(shù)

其中的C00為1，使用時由于全球質(zhì)量守恒，因此需要將其置0。此外，C10為0，由于GRACE衛(wèi)星對地心運動不敏感，需要補充一階項。另外C20（表中的第三行、第一個）不準確，也需要替換。下面是經(jīng)過處理的數(shù)據(jù)。

后期還需要進行濾波處理。希望對大家有所幫助?。?/p>