jingweiduzhuanhua.h

这是大地主题解算时求解大地距离的大椭圆距离的程序实现,用到的头文件我将同时上传

#ifndef JINGWEIDUZHUANHUA
#define JINGWEIDUZHUANHUA
#include<math.h>
#include"jingweidu.h"
#include"tuoqiu.h"
//首先定义平方函数 
long double square(long double s)
{
	return s*s;
}

//大地经纬度转化为地心经纬度
dixinjingweidu daditodixin(dadijingweidu da,tuoqiu tuo)
{
	dixinjingweidu di;
	di.psi=atan(square(1-tuo.bianlv)*tan(da.B));
	di.L=da.L;
	return di;
}
//椭球面坐标到辅助球面坐标的对应关系
qiumianjingweidu fuzhuqiumianzuobiao1(dadijingweidu da)//包括两种对应关系 ，此为一种
{
	qiumianjingweidu qiu;
	qiu.psi=da.B;
	qiu.lambda=da.L;
	return qiu;
}
qiumianjingweidu fuzhuqiumianzuobiao2(dadijingweidu da,tuoqiu tuo)//另一种
{
	long double a=tuo.changzhou;
	long double b=tuo.changzhou-tuo.changzhou*tuo.bianlv;
	qiumianjingweidu qiu;
	qiu.psi=atan(tan(da.B)*b/a);
	qiu.lambda=da.L;
	return qiu;
}
//符号函数
int SGN(long double B)
{
	if(B>0)
		return 1;
	else if(B<0)
		return -1;
	else 
		return 0;
}
#endif