#include?
#include?
#include?

#define?GM???????????????3.986005e14

void?main（）
{
long?double?toe?=??0.720000000000e+04;				//星歷表參考歷元
long?double?a	=??0.2656013544884677e+08;			//軌道長半徑
long?double?e???=??0.678421219345e-02;				//軌道偏心率
long?double?i0??=??0.958512160302e+00;				//按參考歷元toe計算的軌道傾角
long?double?w???=??-0.258419417299e+01;				//近地點角距
long?double?o0??=?-0.137835982556e+01;				//按參考歷元toe計算的升交點赤徑
long?double?M0??=?-0.290282040486e+00;				//按參考歷元toe計算的平近點角
long?double?delta_n?=??0.451411660250e-08;			//由精密星歷計算得到的衛星平均角速度與按給定參數計算所得的平均角速度之差
long?double?o???=?-0.819426989566e-08;				//升交點赤徑變化率
long?double?I???=?-0.253939149013e-09;				//軌道傾角變化率
long?double?cus?=??0.912137329578e-05;				//緯度俯角的正弦調和項改正的振幅
long?double?cuc?=??0.189989805222e-06;				//緯度俯角的余弦調和項改正的振幅
long?double?cis?=??0.949949026108e-07;				//軌道傾角的正弦調和項改正的振幅
long?double?cic?=??0.130385160446e-07;				//軌道傾角的余弦調和項改正的振幅
long?double?crs?=??0.406250000000e+01;				//軌道半徑的正弦調和項改正的振幅
long?double?crc?=??0.201875000000e+03;				//軌道半徑的余弦調和項改正的振幅
long?double?we?=??7.29211567e-5;					//地球自轉的速率
long?double?t???=?0.720000000000e+04;				//觀測時刻
long?double?n0ntk;
long?double?Mk?Ek??Vk?fk;
long?double?u?r?i?uk?rk?ik?ok?xk?yk;
long?double?Xk?Yk?Zk;
int?diedai;
?????????
n0?=?sqrt（GM/（a*a*a））;????????????????
n?=?n0?+?delta_n;
tk?=?t?-?toe;
Mk?=?M0?+?n*tk;
Ek?=?Mk;
for（diedai=0;?diedai<3;?diedai++）?Ek?=?Mk?+?e*sin（Ek）;??
Vk?=?atan（?sqrt（1.0-e*e）*sin（Ek）/（cos（Ek）-e））;?
fk?=?Vk?+?w;
u?=?cuc*cos（2.0*fk）?+?cus*sin（2.0*fk）;
r?=?crc*cos（2.0*fk）?+?crs*sin（2.0*fk）;
i?=?cic*cos（2.0*fk）?+?cis*sin（2.0*fk）;
uk?=?fk?+?u;
rk?=?a*（1.0-e*cos（Ek））?+?r;
ik?=?i0?+?i?+?I*tk;
xk?=?rk*cos（uk）;
yk?=?rk*sin（uk）;
ok?=?o0?+?（o-we）*tk?-?we*toe;
Xk?=?xk*cos（ok）?-?yk*sin（ok）*cos（ik）;
Yk?=?xk*sin（ok）?+?yk*cos（ok）*cos（ik）;
Zk?=?yk*sin（ik）;
?
cout?<<“衛星運行平均角速度n0：“<cout?<<“校正過的衛星運行平均角速度n：“<cout?<<“歸化時間tk：“<cout?<<“衛星平近點角Mk：“<cout?<<“衛星偏近點角Ek：“<cout?<<“衛星真近點角Vk：“<cout?<<“衛星升交距角fk：“<cout?<<“攝動改正項uri：“<cout?<<“經過攝動改正的升交距角uk：“<cout?<<“經過攝動改正的衛星矢徑rk：“<cout?<<“經過攝動改正的軌道傾角ik：“<cout?<<“衛星在軌道平面坐標系的坐標xk，yk：“<cout?<<“觀測時刻升交點經度ok：“<cout?<<“衛星在地心固定坐標系中的直角坐標為XkYkZk：“<}

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-12-01?22:48??z\
?????目錄???????????0??2015-12-01?22:51??z\GPS衛星位置計算\
?????目錄???????????0??2015-12-01?22:50??z\GPS衛星位置計算\Debug\
?????文件??????229462??2015-12-01?22:50??z\GPS衛星位置計算\Debug\GPS衛星位置計算.exe
?????文件??????254596??2015-12-01?22:50??z\GPS衛星位置計算\Debug\GPS衛星位置計算.ilk
?????文件??????277092??2015-12-01?22:50??z\GPS衛星位置計算\Debug\GPS衛星位置計算.pch
?????文件??????418816??2015-12-01?22:50??z\GPS衛星位置計算\Debug\GPS衛星位置計算.pdb
?????文件???????41984??2015-12-01?22:51??z\GPS衛星位置計算\Debug\vc60.idb
?????文件???????61440??2015-12-01?22:50??z\GPS衛星位置計算\Debug\vc60.pdb
?????文件???????10412??2015-12-01?22:50??z\GPS衛星位置計算\Debug\z.obj
?????文件????????4378??2015-12-01?22:50??z\GPS衛星位置計算\GPS衛星位置計算.dsp
?????文件?????????538??2015-12-01?22:48??z\GPS衛星位置計算\GPS衛星位置計算.dsw
?????文件???????33792??2015-12-01?22:51??z\GPS衛星位置計算\GPS衛星位置計算.ncb
?????文件???????48640??2015-12-01?22:51??z\GPS衛星位置計算\GPS衛星位置計算.opt
?????文件????????1334??2015-12-01?22:50??z\GPS衛星位置計算\GPS衛星位置計算.plg
?????文件????????2784??2015-12-01?22:50??z\GPS衛星位置計算\z.cpp