實驗一中使用 GPS 預估星歷估計時，100 分鐘短弧與 240 分鐘短 弧軌道之 3D RMS 值相差無幾均約為±35 公分；使用 GPS 精密星曆 估計後，兩不同時間長度弧段之 3D RMS 值亦相當一致約為±28 公 分。100 分鐘短弧軌道之 3D RMS 值於使用 GPS 精密星曆後較使用 GPS 預估星曆提升 16%，對 240 分鐘短弧軌道而言則提升 24%。
綜合比較 100 分鐘與 240 分鐘短弧軌道成果之內在精度與外在精 度，無論使用 GPS 預估星曆或是精密星曆，兩者之成果顯示精度均 滿足大氣任務所需之±50 公分等級。
與使用 GPS 預估星曆之 CHAMP 成果比較，使用 GPS 精密星曆 以動態法估計之 CHAMP 軌道成果精度約可提升 20%左右，且成果較 穩定，若軌道需要較高之精度要求時，建議可採用 GPS 精密星曆。
然而 GPS 預估星曆可事先獲得，對於須近即時精密低軌衛星軌道計 算的大氣科學研究任務而言，只要 GPS 預估星曆品質穩定，使用 GPS 預估星曆也已可滿足其任務需求。
實驗中可發現當 CHAMP 衛星之訊號發生斷訊時，動態法無法對 該斷訊時刻的 CHAMP 衛星位置進行解算。這是使用動態法求解衛星 軌道所面臨的最大問題，即使有觀測資料，也有可能因為訊號品質不 良而導致求解出精度不佳的成果。因此若要求解連續且穩定的軌道成 果，結合動力積分計算軌道勢在必行。
欲得到更佳的低軌衛星軌道成果以適合其他科學任務需求，例如 探測高精度地球重立場，加入衛星力學模式並同時解算 GPS 衛星軌 道以獲得較 GPS 預估星曆更完善之 GPS 星曆進行高精度定軌將是必 要的。
實驗二僅探討以參數吸收積分器與 PSO 差異之成果，採用的改 善方式則是以動力參數和經驗參數加以吸收兩者差異，成果顯示，平 均 43 個 100 分鐘短弧軌道之徑向、沿軌道方向與橫向 RMS 值均在 5
公分以內，3D RMS 在±10 公分內。此一成果表示，研究中所使用的 積分器產生的短弧軌經參數修正後可與 PSO 達到約±10 公分的位置 精度，未來引入實際觀測資料後，若觀測資料處理得當，則可獲得此 一量級的軌道精度。
參考文獻
太空計畫室, http://www.nspo.gov.tw/
洪志偉，GPS 衛星之雷射測距資料應用於大地絕對坐標之測定，國立 成功大學測量工程研究所碩士論文，台南，1998。
胡明城 和 魯福，現代大地測量學，測繪出版社，北京，1994。
侯瓔瑞，以 Colombo 經驗模式表示低軌衛星軌道擾動之分析，國立 成功大學測量工程研究所碩士論文，台南，2001。
劉佾博，近即時動態法中華三號衛星星曆之測定：使用 MicroLab-1 資料，國立成功大學測量工程研究所碩士論文，台南，1999。
Balmino, G.., and J.B., Barriot, Numerical integration techniques revisted,
Manuscripta geodetica, vol. 15, 1-10, 1989.
Barlier, F., C. Berger, J.L. Falin, G.. Kockarts, and G.. Thuillier, A thermospheric model based on satellite drag data, Ann. Geophys, vol. 34, 9-24, 1978.
Boucher, C., Z. Altamimi, M. Feissel, and P. Sillard, The 1997 International Terrestrial Reference Frame (ITRF97), IERS Technical Note 27, Observatoire de Paris, Paris, France, 1999.
Byun, S.H., and B.E. Schutz, Improving satellite orbit solution using double-differenced GPS carrier phase in kinematic mode, Journal of
Geodesy, vol. 75, 533-543, 2001.
Byun, S.H., Satellite orbit determination using triple-differenced GPS carrier phase in pure kinematic mode, Journal of Geodesy, vol. 76, 569-585, 2003.
Cappellari, J.O., C.E. Velez, and A.J. Fuchs, Mathematical theory of Goddard trajectory determination system, Goddard Space Flight Center Report X582-76-77, 1979.
Cefola, P.J., Equinoctial orbit elements-application to artificial satellite orbits, AIAA PAPER 72-0937, Washington, D.C. U.S.A., 1972.
CDDIS Website：ftp://128.183.103.97
Colombo, O.L., Ephemeris error of GPS satellites. Bull Geod., vol. 60, No. 1, 64-84, 1986.
Colombo, O.L., The dynamic of global positioning system orbits and the determination of precise ephemeris, J. Geophys Res., vol. 94, 9167-9182, 1989.
Cui, C. and M. Mareyen, Gauss’s equations of motion in terms of Hill variables and first application tonumerical integration of satellite orbits,
Manu. Geod., vol.17, 155-163, 1992.
Escobal, P.R., Methods of orbit determination, Krieger Publishing Company, Malabar, Florida, 1975.
GFZ Website：http://op.gfz-potsdam.de/champ/index_CHAMP.html
Gill, E., Smooth Bi-polynomial Interpolation of Jacchia Atmospheric Densities for Efficient Satellite Drag Computation, DLR-GSOC IB 96-1, German Aerospace Center (DLR), 1996.
Goad, C.C., and L. Goodman, A modified Hopfield tropospheric refraction correction model, Proceedings of the Fall Annual Meeting of
the American Geophysical Union, San Francisco California, December,
12-17, 1974.Goad, C.C., D.A. Grejner-Brzezinska, and M. Yang, Determination of high-precision GPS orbit using triple-differencing technique, J. Geodesy, vol.70, 655-662, 1996.
Hedin, A.E., MSIS-86 thermospheric model, J. Geophy. Res., vol. 92, 4649-4662, 1987.
Hedin, A.E., Extension of the MSIS thermosphere model into the middle and lower atmosphere, J. Geophys. Res., vol. 96, 1159, 1991.
Heiskanen, W.A., and H. Moritz, Physical Geodesy, W.H. Freeman and Co., San Francisco, 1967.
IGS Website：http://igscb.jpl.nasa.gov
Jacchia, L.G., Revised Static Models of the Thermosphere and Exosphere with Empirical Temperature Profiles, SAO Special Report 332, Cambridge, 1971.
Jacchia, L.G., The upper atmosphere, Smithon. Astrophy. Obs., 20, 1977.
Krogh, F.T., Variable Order Integrators for the Numerical Solution of Ordinary Differential Equations, Manual, JPL, Pasadena, California, 1969.
Lee, L.C.(ed), Special issue on COSMIC, Terr. Atm. Ocean, vol. 11, No. 1, 2002.
Lemoine, F.G., The Development of Joint NASA GSFC and National Imagery and Mapping Agency (NIMA) Geopotential Model EGM96, NASA/TP-1998-206861, 1998.
McCarthy, D.D., IERS Standards(1992), IERS Technical Note 13.
Observatoire de Paris, 1992.
McCarthy, D.D., IERS Standards(1996), IERS Technical Note 21.
Observatoire de Paris, 1996.
McCarthy, J.J. et al., GEODYN II System Operation Manual, vol.1, NASA/Goddard Space Flight Center, Greenbelt, 1993.
Milani, A., A.M. Nobili, and P. Farinella, Non-Gravitational Perturbations and Satellite Geodesy, Adam Hilger, Bristol, 1987.
Montenbruck, O., and E. Gill, Satellite Orbit-Models, Methods, and Application, Springer, New York, 2001.
Reigber, Ch., G.. Balmino, P. Schwintzer, R. Biancale, A. Bode, J.-M.
Lemoine, R. Koenig, S. Loyer, H. Neumayer, J.-C., Marty, F. Barthelmes, F. Perosanz, and S.Y. Zhu, A high-quality global gravity field model from CHAMP GPS tracking data and accelerometry (EIGEN-1S), Geophys.
Res. Lett., vol. 29, No.14, 10.1029/2002GL015064, 2002.
Robert, C.E., Densities based upon Jacchia’s 1970 models, Cel. Mech., vol. 4, 368-377, 1971.
Seeber, G., Satellite Geodesy: Foundations, Methods, and Application, Waler de Gruyter, Berlin, New York, 1993.
Springer, T.A., and G. Beutler, Towards an Official IGS Orbit by Combining the Results of All IGS Processing Centers, Proceedings of the
1993 IGS Workshop, 242-249, 1993.
Svehla, D., and M. Rothacher, Kinematic and reduced-dynamic precise orbit determination of low earth orbits, Advances in Geosciences, vol. 1, 47-56, 2003.
The Astronomical Almanac: for the Year 1994, Washington, U.S.
Government Printing Office, 1994.
Velez, C.E., Calculation of precision satellite with equinoctial elements (VOP formulation), in A. Dold and B. Eckmann (eds.), Lecture notes in Mathematics, vol. 362, Springer-Verlag, New York, 184-206, 1972.
Wahr, J.M., The forced nutations of an elliptical, rotating, elastic, and oceanless earth, Geophysical Journal Royal Astronomical Society, vol. 64, 705-727, 1981.
Wu, S.C., T.P. Yunck and C.L. Thornton, Reduced-dynamic technique for precise orbit determination of low Earth satellites, J. Guidance, Control,
and Dynamics, vol.14, No. 1, 24-30, 1991.
Yunck, T.P., S.C. Wu, J.T. Wu and C.L. Thornton, Precise tracking of remote sensing satellites with the global positioning system, IEEE
Transactions of Geoscience and Remote Sensing, vol. 28, No. 1, 108-116,
1990.附 錄 A 研究開發之程式說明
本研究中使用了動態法估計低軌衛星軌道與模擬動力參數估 計，因此程式發展也分為兩個主軸，第一部分為動態法部分，第二部 分則為衛星受力模式與參數估計，所有的程式碼均以 Fortran 語法開 發撰寫，以下將介紹主要開發程式之功能。
I 動態法相關程式：
newstore：讀入低軌衛星與地面追蹤站之 GPS 觀測資料 storeeph：讀入 GPS 廣播星曆
storepre：讀入 GPS 精密星曆或預估星曆
example2：地面追蹤站單點定位以求得時錶誤差改正值 abs_kin：低軌衛星單點定位以求得近似坐標與時錶誤差改正 leo_inerp：將原本 10 秒一筆的近似坐標內插為 30 秒一筆之坐標 leo_td：將一個地面追蹤站與低軌衛星組成三次差分觀測資料 leo_apriori：內插極運動參數
leo_reduce_batche：地面追蹤站之各項誤差改正 leo_chol2_batch：法方程式約化
leo_estim_batche：最小二乘平差估計低軌衛星坐標
champatt：讀入 CHAMP 衛星姿態參數，並轉換至慣性坐標系 ecef2inertial：慣性坐標與地心地固系坐標系統相互轉換
II 模擬動力參數估計相關程式：
dyncol：建立衛星運動方程式、力學模式積分、軌道預估、參數估計 之主程式
generate：設定輸入與成果輸出檔之副程式
integrat：讀入積分參數與起始時刻、起始位置與速度之副程式 perturb：整合各項力學模式之副程式
voe：整合各項變分方程式之副程式 grvaccl：地球非球體擾動模式之副程式 lunaccl：多體擾動模式之副程式
solaccl：太陽輻射壓模式之副程式 kepaccl：中心力模式之副程式 dragaccl：大氣阻力模式之副程式 tidgrav：海潮與固體潮模式之副程式 dvdq：數值積分之副程式
xyz2arc：慣性卡氏坐標系與衛星本體坐標系統相互轉換之副程式 intsubpwvoe8：動力參數與經驗參數估計之副程式