水准网间接平差程序设计（C++).docVIP

//////////////////////////////////////////////////// // // // 山东科技大学 测绘学院 // // 闫循鹏 // // QQ 1187774395 // //////////////////////////////////////////////////// //////////////////////////////////////////////////// // visual C++6.0 编译通过 // //////////////////////////////////////////////////// /////////////////////////////////////////////////// // 参考资料 // // 部分网络资料 // // 宋力杰《测量平差程序设计》 // //姚连壁《基于matlab的控制网平差程序设计》 // /////////////////////////////////////////////////// #includeiostream #includefstream #include stdlib.h #includemath.h #include iomanip using namespace std; //////////////////////////////////////////////////////////////////////////class class SZWPC { private: int gcz_zs; //高差总数 int szd_zs; //总点数 int yz_szd_zs; //已知点数 double m_pvv; //[pvv] int *qsd_dh; //高差起点号 int *zd_dh; //高差终点号 char **dm; //点名地址数组 double *gcz; //观测值数组 double *szd_gc; //高程值数组 double *P; //观测值的权 double *ATPA,*ATPL; //法方程系数矩阵与自由项 double *dX; //高程改正数、平差值 double *V; //残差 double m_mu; //单位权中误差 public: SZWPC(); ~SZWPC(); int ij(int i,int j);//对称矩阵下标计算函数 bool inverse(double a[],int n);//对称正定矩阵求逆(仅存下三角元素)(参考他人） void inputdata(char *datafile);//输入原始数据函数 int dm_dh(char *name); //点名转点号 void ca_H0(); //近似高程计算函数 void ca_ATPA(); //法方程组成函数 void ca_dX(); //高程平差值计算函数 void printresult(char *resultfile); //精度估计与平差值输出函数 double ca_V(); //残差计算函数 void zxecpc(char *resultfile);//最小二乘平差函数 }; ////////////////////////////////////////////////////////////////////// // 构造函数 SZWPC::SZWPC() { gcz_zs=0; szd_zs=0; yz_szd_zs=0; } ////////////////////////////////////////////////////////////////////// // 析构函数 SZWPC::~SZWPC() { if(gcz_zs0) { delete []qsd_dh; delete []zd_dh; delete []gcz; delete []P; delete