文件名称:gps平差程序代码 矩阵运算
文件大小:85KB
文件格式:TXT
更新时间:2012-06-07 00:49:54
gps平差 代码 矩阵
int adj::doadj() // 普通最小二乘平差
{
MAT APA,AT;
AT=A.T();
APA=AT*P*A;
N_1=APA.inverse1();
int flag;
if(APA.R()==APA.GetRow())
flag=1;
else flag=0;
if(flag!=1)
{
this->flag=0;
return 0;
}
MAT AX=A.T()*P*l;
X=N_1*AX;
AX=A*X;
for(int i=0;i