摄影测量单片空间后方交会程序《C++》

摄影测量单片空间后方交会程序《C++》摄影测量单片空间后方交会程序《C++》摄影测量单片空间后方交会程序《C++》摄影测量单片空间后方交会程序《C++》编制仅供参考审核批准生效日期地址：电话：传真:邮编:输入文件形式如下：C++源程序如下：#include<iostream>#include<fstream>#include<cmath>#include<string>#include<iomanip>usingnamespacestd;constintn=6;voidinverse(doublec[n][n]);template<typenameT1,typenameT2>voidtranspose(T1*mat1,T2*mat2,inta,intb);template<typenameT1,typenameT2>voidmulti(T1*mat1,T2*mat2,T2*result,inta,intb,intc);template<typenameT>voidinput(T*mat,inta,intb);template<typenameT>voidoutput(T*mat,char*s,inta,intb);intmain(){ ofstreamoutFile; (5); doublex0=,y0=;doublefk=;xt"); while(!()) {for(inti=0;i<a;i++) for(intj=0;j<b;j++) inFile>>mat[i][j]; } (); return;}template<typenameT>voidoutput(T*mat,char*s,inta,intb){cout<<setw(15)<<s<<endl; for(inti=0;i<a;i++) {for(intj=0;j<b;j++) cout<<setw(13)<<mat[i][j]; cout<<endl; } return;}voidinverse(doublec[n][n]){inti,j,h,k;doublep;doubleq[n][12];for(i=0;i<n;i++)//构造高斯矩阵for(j=0;j<n;j++)q[i][j]=c[i][j];for(i=0;i<n;i++)for(j=n;j<12;j++) {if(i+6==j)q[i][j]=1;elseq[i][j]=0;}for(h=k=0;k<n-1;k++,h++)//消去对角线以下的数据for(i=k+1;i<n;i++) {if(q[i][h]==0)continue;p=q[k][h]/q[i][h];for(j=0;j<12;j++){q[i][j]*=p;q[i][j]-=q[k][j];} }for(h=k=n-1;k>0;k--,h--)//消去对角线以上的数据for(i=k-1;i>=0;i--){if(q[i][h]==0)continue;p=q[k][h]/q[i][h];for(j=0;j<12;j++){q[i][j]*=p;q[i][j]-=q[k][j];}}for(i=0;i<n;i++)//将对角线上数据化为1{p=q[i][i];for(j=0;j<12;j++)q[i][j]*=p;}for(i=