【正文】 f%lf%lf,amp。Data[i][3],amp。irow。Omega, const double amp。fprintf(fp_out,%\n,temp)。 fprintf(fp_out,%3d\t%10lf\t%10lf\t%10lf\t%10lf\t%10lf\n, i+1,temp[0],temp[1],temp[2],temp[3],temp[4])。j++) { fprintf(fp_out,%10lf\t,*(R+i*3+j))。fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。double Phi(0),Omega(0),Kappa(0)。 R[2][1]=sin(Phi)*sin(Kappa)+cos(Phi)*sin(Omega)*cos(Kappa)。 A[2*i][4]=f*sin(Kappa)((xxx)/f)*((xxx)* sin(Kappa)+(yyy)*cos(Kappa))。 MatrixMul(AT,6,L,1,2*Num,ATL)。i++){ V[i]=L[i]。delete []V。i++){ for (j=0。 js[k]=j。j++) { if(j!=k) { *(m+i*row+j)=*(m+i*row+k)**(m+k*row+j)。irow。搜集資料時，在網(wǎng)上搜索到了一些矩陣求逆的代碼，但是有些也不能實現(xiàn)，找到合適的后，也要考慮把它運用到主程序的問題?？傮w來說，空間后方交會的知識，在學(xué)習(xí)時就有點懵，后來自己下來又看書加深了理解，把它弄懂了。i++) { swap(*(m+i*row+r),*(m+i*row+is[r]))。 } } for (i=0。i++) { for (j=k。int i,j。delete []ATA。A[0][0],2*Num,Xg,1,6,V)。j++) { *(AT+j*2*Num+i)=A[i][j]。 A[2*i][0]=(R[0][0]*f+R[0][2]*(xxx))/Z。 R[1][0]=cos(Omega)*sin(Kappa)。 Ys+=Data[i][3]。for (i=0。i3。temp[1],amp。num,amp。Xs,const double amp。}/***********************************************函數(shù)名：MatrixMul *函數(shù)介紹：求兩個矩陣的積,*參數(shù)：Jz1(第一個矩陣),row(第一個矩陣行數(shù)),*Jz2(第二個矩陣),row(第二個矩陣列數(shù)),(第一個*矩陣列數(shù)),(out)JgJz(存放結(jié)果矩陣)*返回值：void**********************************************/void MatrixMul(double *Jz1,const int amp。i++){ //讀取數(shù)據(jù)，忽略編號 if(fscanf(fp_input,%*d%lf%lf%lf%lf%lf, amp。}fscanf(fp_input,%d%lf,amp。R39。double f(0),m(0)。5 源程序include include include const double PRECISION=1e5。2作業(yè)思想利用攝影測量空間后方交會的方法求解。m,double amp。}if (Data!=NULL){ delete []Data。 system(cls)。f,amp。Data[i][4])!=5) { return 3。i++){ for (int j=0。Kappa,const double *R){FILE *fp_out。fprintf(fp_out,內(nèi)方位元素(f x0 y0)：)。}fclose(fp_input)。 } fprintf(fp_out,\n)。fclose(fp_out)。double R[3][3]={}。 R[2][2]=cos(Phi)*cos(Omega)。 A[2*i][5]=(yyy)。 MatrixMul(ATA,6,ATL,1,6,Xg)。 VTV+=V[i]*V[i]。return 0。jrow。 } } } if (0 == max) { return 1。 } } } } for (i=0。i++){ for (j=0。不斷調(diào)試，不斷改進(jìn)，這個過程還是有些難度的，最后還實現(xiàn)了用文本導(dǎo)入數(shù)據(jù)計算及文本輸出結(jié)果。此次作業(yè)，便是又一次的挑戰(zhàn)。irow。j++) { if (j!=k) { *(m+k*row+j)*=*((m+k*row+k))。irow。double *ptemp,*pt=m。delete []AT。MatrixMul(amp。j6。 L[2*i+1]=Data[i][1](f*Y/Z)。 R[0][2]=sin(Phi)*cos(Omega)。i++){ Xs+=Data[i][2]。}fprintf(fp_out,\n外方位元素精度：)。for (i=0。temp[0],am