欢迎来到天天文库
浏览记录
ID:31946926
大小:55.50 KB
页数:10页
时间:2019-01-29
《水准网间接平差程序设计(c++)》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、可编辑版//////////////////////////////////////////////////////visualC++6.0编译通过///////////////////////////////////////////////////////////////////////////////////////////////////////////参考资料////部分网络资料////宋力杰《测量平差程序设计》////姚连壁《基于matlab的控制网平差程序设计》/////////////////////
2、////////////////////////////////#include#include#include#include#includeusingnamespacestd;//////////////////////////////////////////////////////////////////////////classclassSZWPC{private:intgcz_zs;//高差总数intszd_zs;
3、//总点数intyz_szd_zs;//已知点数doublem_pvv;//[pvv]int*qsd_dh;//高差起点号int*zd_dh;//高差终点号char**dm;//点名地址数组double*gcz;//观测值数组double*szd_gc;//高程值数组double*P;//观测值的权double*ATPA,*ATPL;//法方程系数矩阵与自由项double*dX;//高程改正数、平差值double*V;//残差doublem_mu;//单位权中误差public:SZWPC();~SZWPC();in
4、tij(inti,intj);//对称矩阵下标计算函数boolinverse(doublea[],intn);//对称正定矩阵求逆(仅存下三角元素)(参考他人)Word完美格式可编辑版voidinputdata(char*datafile);//输入原始数据函数intdm_dh(char*name);//点名转点号voidca_H0();//近似高程计算函数voidca_ATPA();//法方程组成函数voidca_dX();//高程平差值计算函数voidprintresult(char*resultfile);/
5、/精度估计与平差值输出函数doubleca_V();//残差计算函数voidzxecpc(char*resultfile);//最小二乘平差函数};////////////////////////////////////////////////////////////////////////构造函数SZWPC::SZWPC(){gcz_zs=0;szd_zs=0;yz_szd_zs=0;}////////////////////////////////////////////////////////////////
6、////////析构函数SZWPC::~SZWPC(){if(gcz_zs>0){delete[]qsd_dh;delete[]zd_dh;delete[]gcz;delete[]P;delete[]V;}if(szd_zs>0){delete[]szd_gc;delete[]ATPA;delete[]ATPL;delete[]dX;Word完美格式可编辑版for(inti=0;i7、////////////////////////////////////////////////////////////////对称矩阵下标计算函数intSZWPC::ij(inti,intj){return(i>=j)?i*(i+1)/2+j:j*(j+1)/2+i;}////////////////////////////////////////////////////////////////////////////对称正定矩阵求逆(仅存下三角元素)(参考他人)boolSZWPC::inverse(double8、a[],intn){double*a0=newdouble[n];for(intk=0;k
7、////////////////////////////////////////////////////////////////对称矩阵下标计算函数intSZWPC::ij(inti,intj){return(i>=j)?i*(i+1)/2+j:j*(j+1)/2+i;}////////////////////////////////////////////////////////////////////////////对称正定矩阵求逆(仅存下三角元素)(参考他人)boolSZWPC::inverse(double
8、a[],intn){double*a0=newdouble[n];for(intk=0;k
此文档下载收益归作者所有