main.cpp

用来计算航天器再入质点弹道

#include"Model.h"
#include"数值积分.h"
#include"插值.h"
#include "AirModel.h"
//************************************************************
int main()
{
    int				N			=		6;              
	double			STRAT_time	=		0.0;
	double          STOP_time	=		10.0;
	double			Step		=		0.001;
	double          Atitude		=		0.0;
	double			R0			=		6378135.0;
	const  double	PI			=		3.14159265;
	double*			X			=		new double [N]; 
    double*			F			=		new double [N];
    FILE*			fp			=		fopen("result.txt","w");
					
					X[0]		=		(120000.0+R0)/R0;
					X[1]		=		77*PI/180;
					X[2]		=		-15*PI/180;
					X[3]		=		7910/sqrt(9.81*6378135);
					X[4]		=		-1.5*PI/180;
					X[5]		=		53.0*PI/180;
    int				T			=		0;
    while(X[3]*sqrt(9.81*6378135)>800.0)
	{
		RK4(N,Step,X,F);
		Atitude=(X[0]-1.0)*R0;
		T++;
	    fprintf(fp,"%d   %lf   %lf   %lf   %lf\n",T,Atitude,X[3]*sqrt(9.81*6378135),X[0],X[3]);
	}

	delete [] X;
    delete [] F;

	return 0;
}
//**************************************************************