//?Code?difference.cpp?:?Defines?the?entry?point?for?the?console?application.
//
#include?“stdafx.h“
#include?
#include?
#include?
#include?
#include?
#include?
#include?
#include???
using?namespace?Eigen;??
using?namespace?std;
const?double?c=299792458;?//light?speed
const?double?GM=3.986005e14;
const?double?we=7.292115e-5;?//地球自轉速度
//-------------------------------------------------------------------
typedef?vector>?Mat;

//民用時轉化為GPS時
void?gpstime（int?yearint?monthint?dayint?hour?int?minute?double?second?int?&week?double?&secofweek）
{
	int?ry[12]={312931303130313130313031}fry[12]={312831303130313130313031};
	if（year?>?80）?
		year?=?1900?+?year;
	else?
		year?=?2000?+?year;
	int?z=0;
	for（int?i1=1980;?i1?	{
		if（i1%4?==?0?&&?i1%400?!=?0）
			z++;????????????????????//記閏年數
	}	
	if（year%4?==?0?&&?year%100?!=?0）
	{
		for（int?i2?=?0;?i2?			day?+=?ry[i2];

	}
	else
	{
		for（int?i3=0;?i3?			day?+=?fry[i3];

	}
	day?+=?（year-1980）*365?+?z?-?5;
	week?=?day/7;
	secofweek=?day%7*24*3600?+?hour*3600?+?minute*60?+?second?;
}

//觀測文件類
class?ReadObsFile?//觀測文件類
{
public:
	double?X0Y0Z0;?//估計坐標
	int?epochNum;
	vectorsvNum;?//衛星個數
	vector?>c1p1p2l1l2d1d2;
	vector?>prn;//衛星的prn號二維向量
	int?week;
	vectorsecofweek;
	void?ReadObsHeader（string?ofile）
	{
		//用文件流讀文件
		ifstream?in_o（ofile.c_str（））;
		if（in_o）
			cout<<“打開O文件成功“<		else
			cerr<<“未找到O文件“<		string?str;
		//定位到初始坐標
		while（str.find（“APPROX?POSITION?XYZ“）?==?string::npos）
		{
			getline（in_ostr）;
		}
		stringstream?ss（str）;	//用字符串流提取出來初始坐標
		ss>>X0>>Y0>>Z0;	//初始坐標提取

		//--------------------------------------------------------------------	
		//定位到數據頭
		while（str.find（“END?OF?HEADER“）?==?string::npos）
		{
			getline（in_ostr）;
		}
		//計算歷元數
		epochNum=0;
		while（in_o）
		{
			getline（in_ostr）;
			for（int?i?=?0;?（size_t）i?				if（str[i]?==?‘G‘）
				{
					epochNum++;
					break;		
				}
		}
	}?//end?of?function?ReadObsHeader

	void?ReadObsData（string?ofile）
	{
		//定位到文件頭
		ifstream?in_o1（ofile.c_str（））;
		string?str;
		//定位到數據頭
		while（str.find（“END?OF?HEADER“）?==?string::npos）
		{
			getline（in_o1str）;
		}
		//用歷元數定義向量的列數
		c1.resize（epochNum）;
		p1.resize（epochNum）;
		p2.resize（epochNum）;
		l1.resize（epochNum）;
		l2.resize（epochNum）;
		d1.resize（epochNum）;
		d2.resize（epochNum）;
		prn.resize（epochNum）;
		svNum.resize（epochNum）;
		secofweek.resize（epochNum）;

		vector?second（epochNum）;
		int?year?=?0month?=?0day?=?0hour?=?0minute?=?0mark?=?0;
		week=0;
		char?g;

		for（int?i?=?0;?i?		{
			//week和secofweek提取
			in_o1>>year>>month>>day>>hour>>minute>>second[i]>>mark;
			//cout<