作業(yè)4空間后方交會-文庫吧在線文庫
【正文】 s(Omega)*sin(Kappa)。i++) { double X=R[0][0]*(Data[i][2]Xs)+R[1][0]*(Data[i][3]Ys)+ R[2][0]*(Data[i][4]Zs)。 A[2*i][0]=(R[0][0]*f+R[0][2]*(xxx))/Z。 A[2*i+1][2]=(R[2][1]*f+R[2][2]*(yyy))/Z。j++) { *(AT+j*2*Num+i)=A[i][j]。 Zs+=Xg[2]。A[0][0],2*Num,Xg,1,6,V)。for (i=0。delete []ATA。b){double temp=a。int i,j。 pt++。i++) { for (j=k。i++) { swap(*(m+k*row+i),*(m+is[k]*row+i))。 } } for (i=0。 } }}int r。i++) { swap(*(m+i*row+r),*(m+i*row+is[r]))。 ptemp++?？傮w來說，空間后方交會的知識，在學(xué)習(xí)時就有點(diǎn)懵，后來自己下來又看書加深了理解，把它弄懂了。對自己可以說是又一次提高。搜集資料時，在網(wǎng)上搜索到了一些矩陣求逆的代碼，但是有些也不能實(shí)現(xiàn)，找到合適的后，也要考慮把它運(yùn)用到主程序的問題。delete []m。irow。jrow。j++) { if(j!=k) { *(m+i*row+j)=*(m+i*row+k)**(m+k*row+j)。 } } *(m+k*row+k)=1/(*(m+k*row+k))。 js[k]=j。krow。i++){ for (j=0。,doubleamp。delete []V。R[0][0])。i++){ V[i]=L[i]。} while(fabs(Xg[0])=PRECISION ||fabs(Xg[1])=PRECISION || fabs(Xg[2])=PRECISION ||fabs(Xg[3])=PRECISION || fabs(Xg[4])=PRECISION || (Xg[5])=PRECISION)。 MatrixMul(AT,6,L,1,2*Num,ATL)。 } //求矩陣A的轉(zhuǎn)置矩陣AT for (i=0。 A[2*i][4]=f*sin(Kappa)((xxx)/f)*((xxx)* sin(Kappa)+(yyy)*cos(Kappa))。 xxx=f*X/Z。 R[2][1]=sin(Phi)*sin(Kappa)+cos(Phi)*sin(Omega)*cos(Kappa)。double *Xg=new double[6]。double Phi(0),Omega(0),Kappa(0)。int i,j。fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。j6。j++) { fprintf(fp_out,%10lf\t,*(R+i*3+j))。fprintf(fp_out,\tOmega=%10lf\n,Omega)。 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])。inum。fprintf(fp_out,%\n,temp)。fprintf(fp_out,已知數(shù)據(jù)：\n\n已知點(diǎn)數(shù)：)。Omega, const double amp。 }} }/***********************************************函數(shù)名：OutPut *函數(shù)介紹：*參數(shù)：Q協(xié)因數(shù)陣，m精度，m0單位權(quán)中誤差，6個外*方位元素，旋轉(zhuǎn)矩陣*返回值：int，0成功，1失敗**********************************************/int OutPut(const double *amp。irow。 } m=scale/(Num1)。Data[i][3],amp。double scale=0。}fscanf(fp_input,%lf%lf%lf,amp。m,double amp。 } else break。==order || 39。 } return 1。int InverseMatrix(double *matrix,const int amp。Data,double amp。4．6 組成誤差方程： 逐點(diǎn)計算誤差方程式的系數(shù)和常數(shù)項(xiàng)。即利用攝影測量空間后方交會的方法，獲取相片的6個外方位元素。見表1。typedef double DOUBLE[5]。Data,const double amp。if(InputData(Num,Data,m,f)){ if (Data!=NULL) { delete []Data。 //刷新輸入流 char order=getchar()。==order || 39。}/***********************************************函數(shù)名：InputData
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號-1