C++高斯正反算.doc

测绘程序设计基础实验七 实验报告 班 级：测绘0903班 姓 名： 学 号： 指导老师： 实验 7 常用测量程序设计 实验目的 1、巩固类的创建与使用 2、掌握数组参数的传递 3、掌握常用测绘程序设计的技巧 二、实验内容 要清楚计算步骤和公式，还有各个参数计算的前后顺序。编写高斯正反算公式。 界面设计 实验步骤 combox。 我为了实现椭球参数的选择使用了combox控件。 添加每行的显示信息 添加combox函数，定义一个全局变量type获得选择的参数 添加一个C++类，编写关于高斯正反算的各个函数 在头文件中定义 在CPP中写函数 定义参数 角度转换 改变参数 调用函数 实验结果 正算： 反算： #include stdafx.h #include Gauss.h #include math.h const double PI=3.1415926535897932384626433832795; double a=6378245; double b=6356863.0187730473; double e2=sqr(a*a-b*b)/a; double e2=sqr(a*a-b*b)/b; double dms2deg(double dms)//度分秒转换度 { double deg; int idegree imin; double dsec; idegree=int(dms); imin=int((dms-idegree)*100); dsec=((dms-idegree)*100-imin)*100; deg=ideree+double(imin/60)+double(dsec/3600)； return deg; } double deg2dms(double deg)//度转换度分秒 { double dms; int idegree imin; double dsec; double dtmp; idegree=int(dms); dtmp=(deg=idegree)*60; immin=int(dtmp); dsec=(dtmp-imin)*60; return dms; } //椭球参数 void ellipd(int type) { if (type==0) { a=6378245; b=6356863.0187730473; e2=0.006693421622965949; e2=0.00673852541468397; }; if (type==1) { a=6378140; b=6356755.2881575287; e2=0.006694384999587952; e2=0.006739501819472927; }; } //正算 double ComptX(double B) { double a0=1+3/4*e2+45/64*pow(e2,2)+350/512*pow(e2,3)+11025/16384*pow(e2,4); double a2=-1/2*(3/4*e2+60/64*pow(e2,2)+525/512*pow(e2,3)+17640/16384*pow(e2,4)); double a4=1/4*(15/64*pow(e2,2)+210/512*pow(e2,3)+8820/16384*pow(e2,4)); double a6=-1/6*(35/512*pow(e2,3)+2520/16384*pow(e2,4)); double a8=1/8*(315/16384*pow(e2,4)); double X=a*(1-e2)*(a0*B+a2*sin(2*B)+a4*sin(4*B)+a6*sin(6*B)+a8*sin(8*B)); return X; } double ComptForward(double B,double L,double cmd,int type,double x,double y) { L=dms2deg(L); cmd=dms2deg(cmd); double l=(L-cmd)*3600/206265; double b=dms2deg(B)*PI/180; ellipd(type) double N=a/sqr(1-e2*sin(b)*s