作業(yè)4空間后方交會(存儲版)
【正文】 typedef double Double6[6]。 R[0][2]=sin(Phi)*cos(Omega)。iNum。 L[2*i+1]=Data[i][1](f*Y/Z)。 A[2*i+1][1]=(R[1][1]*f+R[1][2]*(yyy))/Z。j6。 Ys+=Xg[1]。MatrixMul(amp。double *mm=new double[6]。delete []AT。a,double amp。double *ptemp,*pt=m。 ptemp++。irow。irow。j++) { if (j!=k) { *(m+k*row+j)*=*((m+k*row+k))。i++) { if(i!=k) { *(m+i*row+k)*=(*(m+k*row+k))。irow。j++) { *ptemp=*pt。此次作業(yè)，便是又一次的挑戰(zhàn)。雖然有些地方是效仿的，不過也是自己要弄懂以后，再試著自己一步步理出來。不斷調(diào)試，不斷改進，這個過程還是有些難度的，最后還實現(xiàn)了用文本導入數(shù)據(jù)計算及文本輸出結(jié)果。return 0。i++){ for (j=0。j++) { swap(*(m+r*row+j),*(m+js[r]*row+j))。 } } } } for (i=0。 for (j=0。 } } } if (0 == max) { return 1。k++){ double max=0。jrow。)**********************************************/int InverseMatrix(double *matrix,const int amp。return 0。delete []L。 VTV+=V[i]*V[i]。//注：協(xié)因數(shù)陣，旋轉(zhuǎn)矩陣等計算本應該使用最后外方位元素值，//由于變換很小忽略double *Q=ATA。 MatrixMul(ATA,6,ATL,1,6,Xg)。i2*Num。 A[2*i][5]=(yyy)。 yyy=f*Y/Z。 R[2][2]=cos(Phi)*cos(Omega)。//迭代計算do { //旋轉(zhuǎn)矩陣 R[0][0]=cos(Phi)*cos(Kappa)sin(Phi)*sin(Omega)*sin(Kappa)。double R[3][3]={}。//設置初始值for (i=0。fclose(fp_out)。j++) { fprintf(fp_out,%20lf\t,*(Q+i*6+j))。 } fprintf(fp_out,\n)。fprintf(fp_out,\tKappa=%10lf\n\n,Kappa)。}fclose(fp_input)。i++){ double temp[5]。fprintf(fp_out,內(nèi)方位元素(f x0 y0)：)。int num。Kappa,const double *R){FILE *fp_out。Q,const double *amp。i++){ for (int j=0。}fclose(fp_input)。Data[i][4])!=5) { return 3。int i。f,amp。f){double x0,y0。 system(cls)。p39。}if (Data!=NULL){ delete []Data。row)。m,double amp。4．7 組成法方程式： 計算法方程的系數(shù)矩陣ATA與常數(shù)項ATL。2作業(yè)思想利用攝影測量空間后方交會的方法求解。4．2 量測控制點的像點坐標： 本次作業(yè)中為已知。5 源程序include include include const double PRECISION=1e5。Num,const DOUBLE *amp。double f(0),m(0)。 fflush(stdin)。R39。return 0。}fscanf(fp_input,%d%lf,amp。f/=1000。i++){ //讀取數(shù)據(jù)，忽略編號 if(fscanf(fp_input,%*d%lf%lf%lf%lf%lf, amp。}//如果m未知則歸算其值if (0==m){ for (i=0。}/***********************************************函數(shù)名：MatrixMul *函數(shù)介紹：求兩個矩陣的積,*參數(shù)：Jz1(第一個矩陣),row(第一個矩陣行數(shù)),*Jz2(第二個矩陣),row(第二個矩陣列數(shù)),(第一個*矩陣列數(shù)),(out)JgJz(存放結(jié)果矩陣)*返回值：void**********************************************/void MatrixMul(double *Jz1,const int amp。 for (int k=0。Xs,const double amp。if (!(fp_input=fopen(,r))){ return 1。num,amp。x,amp。temp[1],amp。fprintf(fp_out,\tXs=%10lf\n,Xs)。i3。for (i=0。for (i=0。Num,const DOUBLE *amp。 Ys+=Data[i][3]。Double6 *A=new Double6[2*Num]。 R[1][0]=co
教學教案相關推薦
鄂ICP備17016276號-1