hifi_f16_aerodata.c

一个F-16仿真程序

		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CXq_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CX1420_ALPHA2_205.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CX1420_ALPHA2_205.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CYr_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CY1620_ALPHA2_407.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CY1620_ALPHA2_407.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CYp_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CY1520_ALPHA2_409.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CY1520_ALPHA2_409.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CZq_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CZ1420_ALPHA2_305.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CZ1420_ALPHA2_305.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CLr_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CL1620_ALPHA2_607.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CL1620_ALPHA2_607.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CLp_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CL1520_ALPHA2_609.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CL1520_ALPHA2_609.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CMq_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CM1420_ALPHA2_105.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CM1420_ALPHA2_105.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CNr_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CN1620_ALPHA2_507.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CN1620_ALPHA2_507.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _delta_CNp_lef(double alpha){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 1; 
	double x[1];	
	FILESIZE = 14;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 14;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA2();
		fp = fopen("CN1520_ALPHA2_509.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CN1520_ALPHA2_509.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _Cy_r30(double alpha, double beta){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 2; 
	double x[2];	
	FILESIZE = 380;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 20;	
		ndinfo.nPoints[1] = 19;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA1();
		X[1] = getBETA1();
		fp = fopen("CY0720_ALPHA1_BETA1_405.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CY0720_ALPHA1_BETA1_405.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
	x[1] = beta;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _Cn_r30(double alpha, double beta){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 2; 
	double x[2];	
	FILESIZE = 380;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 20;	
		ndinfo.nPoints[1] = 19;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA1();
		X[1] = getBETA1();
		fp = fopen("CN0720_ALPHA1_BETA1_503.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CN0720_ALPHA1_BETA1_503.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
	x[1] = beta;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */


double _Cl_r30(double alpha, double beta){
	static int flag = 0;
	static double *DATA = (double*) NULL;
	static double **X;
	static ND_INFO ndinfo ;
	
	FILE *fp;
	double data;
	int i,FILESIZE;
	int nDimension = 2; 
	double x[2];	
	FILESIZE = 380;	

	/* Initialise everything when this function is called for the first time */
	if(flag==0){
		flag = 1;	/* Set to FILE_READ_TAG */
		DATA = (double*) malloc(FILESIZE*sizeof(double)); 
		ndinfo.nDimension = nDimension;
		ndinfo.nPoints = intVector(nDimension);
		ndinfo.nPoints[0] = 20;	
		ndinfo.nPoints[1] = 19;	
		X = (double **) malloc(nDimension*sizeof(double*));

		X[0] = getALPHA1();
		X[1] = getBETA1();
		fp = fopen("CL0720_ALPHA1_BETA1_603.dat","r");
		if(fp==(FILE*) NULL)
			mexErrMsgTxt("Cannot find file CL0720_ALPHA1_BETA1_603.dat in current directory");

		for(i=0;i<FILESIZE;i++){
			fscanf(fp,"%lf",&data);
			DATA[i] = data;
			}
		fclose(fp);
		}

	x[0] = alpha;
	x[1] = beta;
    return (interpn(X,DATA,x,ndinfo));
}/* End of function(...) */