#include
#include
#include
#include
#include
#include
#include
#define?maxsize?5
using?namespace?std;
double?power（double?x?double?n）
{
	double?val?=?1.0;
	while?（n--）
		val?*=?x;
	return?val;
}
const?double?PI?=?3.1415926;
const?double?we?=?7.2921151467?/?power（10?5）;
const?double?GM?=?3.986005*power（10?14）;
class?KK
{
public:
	double?dn?Mo?e?w?Cuc?Cus?Crc?Crs?Cis?Cic?Wo?W?n0;
	double?n?Ek?Ek0?Mk?Vk?Ok?du?di?dr?uk?rk?ik?xk?yk?zk?Lk?Xk?Yk?Zk?Qk?a?UT?io?wo?I?JD?tk?toe?Idot?t;
	int?w_n?Y?M?D?H?min;
	double?sec?a1?a2?a0?aode?L2?GPSn_week?L2P?wxjd?MSB?TGD?IODC?sendtime?wm?ek1?ek2?cosVk?sinVk;
};
void?main（）
{
	int?j?=?1?i?=?0;
	KK?hhd[maxsize];
	ifstream?in;
	in.open（“E:\\fifth\\123.txt“）;
	ofstream?out（“E:\\fifth\\5.txt“）;
	
	in?>>?hhd[j].w_n>>?hhd[j].Y?>>?hhd[j].M?>>?hhd[j].D?>>?hhd[j].H?>>?hhd[j].min?>>?hhd[j].sec
		>>?hhd[j].a0?>>?hhd[j].a1?>>?hhd[j].a2?>>?hhd[j].aode?>>?hhd[j].Crs?>>?hhd[j].dn?>>?hhd[j].Mo?>>?hhd[j].Cuc?>>?hhd[j].e?>>?hhd[j].Cus?>>?hhd[j].a?>>
		hhd[j].toe?>>?hhd[j].Cic?>>?hhd[j].Wo?>>?hhd[j].Cis?>>?hhd[j].I?>>?hhd[j].Crc?>>?hhd[j].w?>>?hhd[j].W?>>?hhd[j].Idot?>>?hhd[j].L2?>>?hhd[j].GPSn_week?>>
		hhd[j].L2P?>>?hhd[j].wxjd?>>?hhd[j].MSB?>>?hhd[j].TGD?>>?hhd[j].IODC?>>?hhd[j].sendtime;
	cout?<	hhd[j].n0?=?sqrt（GM）?/?pow（hhd[j].a?3）;
	hhd[j].n?=?hhd[j].n0?+?hhd[j].dn;

	for?（int?k?=?0;?k	{
		hhd[j].tk?=?k*??60+28800;
		hhd[j].Mk?=?hhd[j].Mo?+?hhd[j].n*hhd[j].tk;
		hhd[j].ek1?=?hhd[j].Mk;

		do{
			
			hhd[j].Ek0?=?hhd[j].ek1;
			hhd[j].ek1?=?hhd[j].Mk?+?hhd[j].e*sin（hhd[j].Ek0）;
		}?while?（fabs（hhd[j].Ek0?-?hhd[j].ek1）>1e-15）;

		hhd[j].Ek?=?hhd[j].ek1;
		//計算觀測時刻的偏近點角Ek
		hhd[j].cosVk?hhd[j].sinVk;
		hhd[j].cosVk?=?（cos（hhd[j].Ek）?-?hhd[j].e）?/?（1?-?hhd[j].e*cos（hhd[j].Ek））;
		hhd[j].sinVk?=?sqrt（1?-?hhd[j].e*hhd[j].e）*sin（hhd[j].Ek）?/?（1?-?hhd[j].e*cos（hhd[j].Ek））;
		if?（hhd[j].cosVk?>?0?&&?hhd[j].sinVk?>?0）
			hhd[j].Vk?=?atan（hhd[j].sinVk?/?hhd[j].cosVk）;
		else?if?（hhd[j].sinVk?>?0?&&?hhd[j].cosVk?			hhd[j].Vk?=?acos（hhd[j].cosVk）;
		else?if?（hhd[j].sinVk?			hhd[j].Vk?=?PI?-?asin（hhd[j].sinVk）;
		else
			hhd[j].Vk?=?2?*?PI?-?acos（hhd[j].cosVk）;

		/*Vk=atan（（sqrt（1-e*e）*sin（Ek））/（cos（Ek）-e））;*///計算真近點角Vk
		hhd[j].Qk?=?hhd[j].Vk?+?hhd[j].w;??????????????????????????????????????????????//計算升交距角Qk
		hhd[j].du?=?hhd[j].Cuc*cos（2?*?hhd[j].Qk）?+?hhd[j].Cus*sin（2?*?hhd[j].Qk）;
		hhd[j].dr?=?hhd[j].Crc*cos（2?*?hhd[j].Qk）?+?hhd[j].Crs*sin（2?*?hhd[j].Qk）;
		hhd[j].di?=?hhd[j].Cic*cos（2?*?hhd[j].Qk）?+?hhd[j].Cis*sin（2?*?hhd[j].Qk）;//計算攝動改正項rurrri
		hhd[j].uk?=?hhd[j].Qk?+?hhd[j].du;
		hhd[j].rk?=?（1?-?hhd[j].e*cos（hhd[j].Ek））*hhd[j].a*hhd[j].a?+?hhd[j].dr;
		hhd[j].ik?=?hhd[j].I?+?hhd[j].di?+?hhd[j].Idot*hhd[j].tk;????????//計算攝動改正的

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-10-25?20:36??作業一程序--C++\Debug\
?????文件??????103424??2015-11-03?11:24??作業一程序--C++\Debug\WGS.exe
?????文件??????677068??2015-11-03?11:24??作業一程序--C++\Debug\WGS.ilk
?????文件?????1305600??2015-11-03?11:24??作業一程序--C++\Debug\WGS.pdb
?????文件?????7471104??2015-11-03?11:26??作業一程序--C++\GNSS作業--C++.sdf
?????文件?????????955??2015-10-24?18:34??作業一程序--C++\GNSS作業--C++.sln
?????文件???????20480??2015-11-03?11:26??作業一程序--C++\GNSS作業--C++.v12.suo
?????目錄???????????0??2015-11-03?11:24??作業一程序--C++\WGS\
?????目錄???????????0??2015-11-03?11:24??作業一程序--C++\WGS\Debug\
?????文件??????388096??2015-11-03?11:24??作業一程序--C++\WGS\Debug\vc120.idb
?????文件??????430080??2015-11-03?11:24??作業一程序--C++\WGS\Debug\vc120.pdb
?????文件????????1466??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.log
?????目錄???????????0??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.tlog\
?????文件????????1242??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.tlog\cl.command.1.tlog
?????文件???????26956??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.tlog\CL.read.1.tlog
?????文件?????????404??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.tlog\CL.write.1.tlog
?????文件????????2162??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.tlog\link.command.1.tlog
?????文件????????2536??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.tlog\link.read.1.tlog
?????文件?????????374??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.tlog\link.write.1.tlog
?????文件?????????184??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS.tlog\WGS.lastbuildstate
?????文件??????279000??2015-11-03?11:24??作業一程序--C++\WGS\Debug\WGS--C++.obj
?????文件????????4084??2015-10-24?19:50??作業一程序--C++\WGS\WGS.vcxproj
?????文件?????????949??2015-10-24?19:50??作業一程序--C++\WGS\WGS.vcxproj.filters
?????文件????????3741??2015-11-03?11:24??作業一程序--C++\WGS\WGS--C++.cpp