GRACE球谐系数阶方差计算+维纳滤波的matlab代码

HDR filename=GIA_n100_mass_0km.txt HDR Institution=JPL / GRACE-TELLUS HDR Variable=Glacial Isostatic Adjustment as Equivalent H2O Thickness Variation Rate up to Degree 100 with no Gaussian Smoothing [Column_3] HDR Units=mm/yr HDR _FillValue=32767.0 HDR Column_1=Longitude HDR Column_2=Latitude HDR Platform=GRACE HDR Sensor=GRACE HDR Date_created=10-Dec-2012 HDR Data_author="F. W. Landerer" <landerer@jpl.nasa.gov> HDR Processing=GIA rates for compressible Earth, computed (by Geruo A) using ICE5G ice history and VM2 viscosity profile HDR Reference=A, G., J. Wahr, and S. Zhong, 2012: Computations of the viscoelastic response of a 3-D compressible Earth to surface loading: an application to Glacial Isostatic Adjustment in Antarctica and Canada, Geophys. J. Int. (in press) HDR LON1_NLONS_DLON=0.5. 360 1. HDR LAT1_NLATS_DLAT=-89.5 180 1. HDR Filter=Gaussian HDR Filter_Width_KM=0 HDR Filter_Max_Degree=100 0.5 -89.5 13.198 0.5 -88.5 14.323 0.5 -87.5 15.158 0.5 -86.5 15.312 0.5 -85.5 14.881 0.5 -84.5 14.257 0.5 -83.5 13.619 0.5 -82.5 12.765 0.5 -81.5 11.377 0.5 -80.5 9.378 0.5 -79.5 7.081 0.5 -78.5 5.043 0.5 -77.5 3.691 0.5 -76.5 3.026 0.5 -75.5 2.718 0.5 -74.5 2.484 0.5 -73.5 2.210 0.5 -72.5 1.761 0.5 -71.5 0.939 0.5 -70.5 -0.225 0.5 -69.5 -1.358 0.5 -68.5 -2.062 0.5 -67.5 -2.252 0.5 -66.5 -2.103 0.5 -65.5 -1.812 0.5 -64.5 -1.511 0.5 -63.5 -1.274 0.5 -62.5 -1.104 0.5 -61.5 -0.964 0.5 -60.5 -0.842 0.5 -59.5 -0.741 0.5 -58.5 -0.646 0.5 -57.5 -0.538 0.5 -56.5 -0.438 0.5 -55.5 -0.360 0.5 -54.5 -0.284 0.5 -53.5 -0.199 0.5 -52.5 -0.129 0.5 -51.5 -0.078 0.5 -50.5 -0.029 0.5 -49.5 0.022 0.5 -48.5 0.058 0.5 -47.5 0.083 0.5 -46.5 0.109 0.5 -45.5 0.129 0.5 -44.5 0.137 0.5 -43.5 0.142 0.5 -42.5 0.147 0.5 -41.5 0.141 0.5 -40.5 0.128 0.5 -39.5 0.118 0.5 -38.5 0.104 0.5 -37.5 0.078 0.5 -36.5 0.049 0.5 -35.5 0.025 0.5 -34.5 -0.006 0.5 -33.5 -0.045 0.5 -32.5 -0.083 0.5 -31.5 -0.119 0.5 -30.5 -0.164 0.5 -29.5 -0.212 0.5 -28.5 -0.256 0.5 -27.5 -0.302 0.5 -26.5 -0.357 0.5 -25.5 -0.411 0.5 -24.5 -0.459 0.5 -23.5 -0.513 0.5 -22.5 -0.575 0.5 -21.5 -0.632 0.5 -20.5 -0.684 0.5 -19.5 -0.744 0.5 -18.5 -0.809 0.5 -17.5 -0.867 0.5 -16.5 -0.925 0.5 -15.5 -0.990 0.5 -14.5 -1.055 0.5 -13.5 -1.114 0.5 -12.5 -1.174 0.5 -11.5 -1.239 0.5 -10.5 -1.299 0.5 -9.5 -1.355 0.5 -8.5 -1.411 0.5 -7.5 -1.464 0.5 -6.5 -1.507 0.5 -5.5 -1.546 0.5 -4.5 -1.584 0.5 -3.5 -1.615 0.5 -2.5 -1.646 0.5 -1.5 -1.697 0.5 -0.5 -1.777 0.5 0.5 -1.883 0.5 1.5 -1.998 0.5 2.5 -2.073 0.5 3.5 -2.010 0.5 4.5 -1.744 0.5 5.5 -1.316 0.5 6.5 -0.871 0.5 7.5 -0.552 0.5 8.5 -0.416 0.5 9.5 -0.426 0.5 10.5 -0.500 0.5 11.5 -0.571 0.5 12.5 -0.609 0.5 13.5 -0.616 0.5 14.5 -0.604 0.5 15.5 -0.586 0.5 16.5 -0.573 0.5 17.5 -0.563 0.5 18.5 -0.554 0.5 19.5 -0.550 0.5 20.5 -0.551 0.5 21.5 -0.552 0.5 22.5 -0.553 0.5 23.5 -0.559 0.5 24.5 -0.570 0.5 25.5 -0.584 0.5 26.5 -0.602 0.5 27.5 -0.621 0.5 28.5 -0.630 0.5 29.5 -0.618 0.5 30.5 -0.584 0.5 31.5 -0.539 0.5 32.5 -0.519 0.5 33.5 -0.604 0.5 34.5 -0.863 0.5 35.5 -1.258 0.5 36.5 -1.635 0.5 37.5 -1.830 0.5 38.5 -1.779 0.5 39.5 -1.525 0.5 40.5 -1.184 0.5 41.5 -0.905 0.5 42.5 -0.793 0.5 43.5 -0.824 0.5 44.5 -0.869 0.5 45.5 -0.837 0.5 46.5 -0.754 0.5 47.5 -0.693 0.5 48.5 -0.680 0.5 49.5 -0.701 0.5 50.5 -0.779 0.5 51.5 -0.957 0.5 52.5 -1.210 0.5 53.5 -1.396 0.5 54.5 -1.331 0.5 55.5 -0.924 0.5 56.5 -0.240 0.5 57.5 0.516 0.5 58.5 1.046 0.5 59.5 1.004 0.5 60.5 0.137 0.5 61.5 -1.496 0.5 62.5 -3.444 0.5 63.5 -5.086 0.5 64.5 -6.021 0.5 65.5 -6.242 0.5 66.5 -5.999 0.5 67.5 -5.591 0.5 68.5 -5.262 0.5 69.5 -5.158 0.5 70.5 -5.297 0.5 71.5 -5.620 0.5 72.5 -6.069 0.5 73.5 -6.572 0.5 74.5 -6.934 0.5 75.5 -6.863 0.5 76.5 -6.212 0.5 77.5 -5.185 0.5 78.5 -4.252 0.5 79.5 -3.826 0.5 80.5 -4.020 0.5 81.5 -4.668 0.5 82.5 -5.521 0.5 83.5 -6.384 0.5 84.5 -7.082 0.5 85.5 -7.430 0.5 86.5 -7.333 0.5 87.5 -6.888 0.5 88.5 -6.307 0.5 89.5 -5.754 1.5 -89.5 13.210 1.5 -88.5 14.351 1.5 -87.5 15.187 1.5 -86.5 15.324 1.5 -85.5 14.875 1.5 -84.5 14.242 1.5 -83.5 13.599 1.5 -82.5 12.717 1.5 -81.5 11.263 1.5 -80.5 9.176 1.5 -79.5 6.807 1.5 -78.5 4.734 1.5 -77.5 3.377 1.5 -76.5 2.714 1.5 -75.5 2.408 1.5 -74.5 2.184 1.5 -73.5 1.938 1.5 -72.5 1.538 1.5 -71.5 0.775 1.5 -70.5 -0.332 1.5 -69.5 -1.423 1.5 -68.5 -2.104 1.5 -67.5 -2.287 1.5 -66.5 -2.134 1.5 -65.5 -1.838 1.5 -64.5 -1.534 1.5 -63.5 -1.299 1.5 -62.5 -1.132 1.5 -61.5 -0.992 1.5 -60.5 -0.868 1.5 -59.5 -0.768 1.5 -58.5 -0.675 1.5 -57.5 -0.569 1.5 -56.5 -0.469 1.5 -55.5 -0.390 1.5 -54.5 -0.314 1.5 -53.5 -0.232 1.5 -52.5 -0.163 1.5 -51.5 -0.112 1.5 -50.5 -0.061 1.5 -49.5 -0.010 1.5 -48.5 0.023 1.5 -47.5 0.047 1.5 -46.5 0.074 1.5 -45.5 0.096 1.5 -44.5 0.102 1.5 -43.5 0.104 1.5 -42.5 0.111 1.5 -41.5 0.109 1.5 -40.5 0.094 1.5 -39.5 0.080 1.5 -38.5 0.067 1.5 -37.5 0.045 1.5 -36.5 0.016 1.5 -35.5 -0.011 1.5 -34.5 -0.041 1.5 -33.5 -0.077 1.5 -32.5 -0.114 1.5 -31.5 -0.152 1.5 -30.5 -0.196 1.5 -29.5 -0.242 1.5 -28.5 -0.286 1.5 -27.5 -0.332 1.5 -26.5 -0.384 1.5 -25.5 -0.436 1.5 -24.5 -0.485 1.5 -23.5 -0.537 1.5 -22.5 -0.596 1.5 -21.5 -0.653 1.5 -20.5 -0.706 1.5 -19.5 -0.764 1.5 -18.5 -0.827 1.5 -17.5 -0.887 1.5 -16.5 -0.946 1.5 -15.5 -1.007 1.5 -14.5 -1.070 1.5 -13.5 -1.130 1.5 -12.5 -1.191 1.5 -11.5 -1.252 1.5 -10.5 -1.310 1.5 -9.5 -1.366 1.5 -8.5 -1.421 1.5 -7.5 -1.471 1.5 -6.5 -1.514 1.5 -5.5 -1.554 1.5 -4.5 -1.590 1.5 -3.5 -1.621 1.5 -2.5 -1.654 1.5 -1.5 -1.704 1.5 -0.5 -1.783 1.5 0.5 -1.893 1.5 1.5 -2.020 1.5 2.5 -2.114 1.5 3.5 -2.086 1.5 4.5 -1.863 1.5 5.5 -1.460 1.5 6.5 -1.000 1.5 7.5 -0.639 1.5 8.5 -0.459 1.5 9.5 -0.436 1.5 10.5 -0.492 1.5 11.5 -0.557 1.5 12.5 -0.599 1.5 13.5