gps网平差程序说明.txt

GPS网平差计算

一. 全局变量说明
      MAX 140                  为定义矩阵的最大维数  //l=3
      PI 3.14159265358979312   圆周率                //l=4     
      rou (180.0*60*60/PI)     一弧度的秒数          //l=5

二. 工具函数说明
    1 矩阵求逆函数

       int inverse(double C[][MAX],double B[][MAX],int n)  //l=149
       n 求逆矩阵的维数
	   返回0说明矩阵不可求逆；
       返回1说明矩阵求逆完成,B为返回逆阵

    2 矩阵乘积函数

       void AXB(double A[][MAX],double B[][MAX],double C[][MAX],int m,int n,int k) //l=42
       计算矩阵乘积A*B=C 
	   m,n 为A阵的行列数
	   n,k 为B阵的行列数  
	   void AXB(double a,double A[][MAX], double aA[][MAX],int m,int n)       //l=65
	   计算数乘矩阵 
       void AXB(double A[][MAX],double B[][1],double C[][1],int m,int n)         //l=53
	   计算矩阵乘列向量

    3 矩阵转置函数
       void AT(double A[][MAX],double AH[][MAX],int m,int n)               //l=73
       void AT(double A[][1],double AH[][MAX],int m)                       //l=80  
       列向量转置为行向量
       void AT(double A[][MAX],double AH[][1],int m)                       //l=86
       行向量转置为列向量
    
    4 矩阵屏幕输出
      void matdis(double A[][MAX],int n,int m)                             //l=132
      void matdis(double A[][1],int n) //l=142
      输出列向量 
    5 矩阵文件输出
      void matout(double A[][1],int n,ofstream out)         //l=410 向文件输出向量
      void matout(double A[][MAX],int n,int m,ofstream out) //l=416 向文件输出矩阵
三. 基本平差工具函数
    1 计算法方程系数阵 
      void ATPA(double A[][MAX],double P[][MAX],double ATPA[][MAX],int m,int n)//l=93
    
	2 计算法方程常数项
	  void ATPL(double A[][MAX],double P[][MAX],double L[][1],double ATPL[][1],int m,int n)
           
	3 计算单位权中误差辅助数据
	  double VPV(double V[][1],double P[][MAX],int m)    //l=119
               
四. 基本平差结构与函数
    1 基本平差结构  
      struct adj //l=205

    2 平差结构数据键盘输入 
      void ksetadj(adj &a) //l=222
    
	3 平差结构数据文件输入（name 为文件名）
	  int fsetadj(adj &aa,char name[20]) //l=251
      返回1说明数据输入成功，返回0说明数据输入失败

    4 平差处理
      int  doadj(adj &a) //l=286
	  返回1说明数据处理成功，返回0说明数据平差失败

	5 稳健估计   
      int rubust(adj &a)  //l=375
      返回1说明数据处理成功，返回0说明数据平差失败

    6 平差结果屏幕输出
      void adjdis(adj &aa) //l=311
	  
	7 平差结果文件输出（name 为文件名）  
      int foutadj(adj &aa, char name[20])  //l=338
	  返回1说明数据输出成功，返回0说明数据输出失败

    8 极大权法平差函数
	  int doadj(adj &a,int knownp)  //l=584
	  knownp 为已知数据的个数 
五. gps网平差结构与函数
    1 gps点结构
      struct Gpspnt   //l=428

    2 基线向量结构 
      struct xiangliang   //l=436                                       
      
    3 gps网结构
      struct Gpsnet     //l=444   

    4 文件输入gps网数据（fname为文件名）
      int finGpsnet(Hnet &a,char *fname)   //l=455
	  包括数据输入、极大权法平差数据准备、并将该数据保存在date.txt文件中 
     
    5 极大权法平差处理及数据输出（outfile为输出文件名）
      int Gpsadj(Gpsnet &a,char *outfile)  //l=610  

    6 高程网平差主函数---用户使用
      int Gdoadj(Gpsnet &a,char *infilename,char* outfilename)  //l=690
      infilename   输入文件文件名
	  outfilename  输出文件名