p1.cpp

综合利用三次法和blewitt法检测与修复周跳

// P1.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#define MAXNUM		5000
#define OBSERVE  "F:\\GPSSOFT\\zhout\\OCC-2-5"
#include "Preprocess.h"
#include <stdio.h>
#include <iostream>
using namespace std;
main()
{


	FILE *fpObsData;
	OBSDAT	data[ MAXNUM ];
	CALCDAT	cal[ MAXNUM ];
	TOL tol = { 20, 4, 3, -4.0e7, +4.0e7, 1.0e7, 4.0e7, 6.0, 0.01, 4.0, 1.0  };

	int result;


	unsigned int b = 100;

	long int i=0;

	double temp;

	if ( ( fpObsData = fopen( OBSERVE, "r" ) ) == NULL )
		printf( "\nObservational Data Failed !\n\n" );

	while(1)
	{
		if ( fscanf( fpObsData, "%ld \t%lf \t%lf \t%lf \t%lf \t%lf \n", &data[i].epoch, &data[i].phi1, &data[i].phi2, &temp, &data[i].P1, &data[i].P2 ) == EOF )
				break;
		i++;
	}

		fclose( fpObsData );		
		result = Preprocess( data, cal, i, tol, &b );		
	return 0;
}