作業(yè)4空間后方交會(完整版)
【正文】 *函數(shù)介紹：從文件（）中讀取數(shù)據(jù),*文件格式如下： *點(diǎn)數(shù) m(未知寫作0)* 內(nèi)方位元素(f x0 y0)*編號 x y X Y Z*實(shí)例：4 0 0 01 2 3 4 *參數(shù)：(in/out)Num(點(diǎn)數(shù))，*(in/out)Data(存放數(shù)據(jù)),m,f,x0,y0*返回值：int ，0成功，1文件打開失敗，2控制點(diǎn)個*數(shù)不足，3文件格式錯誤**********************************************/int InputData(int amp。Num,amp。if (m0 || f0){ return 3。Data[i][0],amp。iNum1。row,double *Jz2, const int amp。k。Ys,const double amp。}fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。temp)。y)。temp[2],amp。fprintf(fp_out,\tYs=%10lf\n,Ys)。i++){ fprintf(fp_out,\t)。i6。i6。Data,const double amp。}Xs/=Num。double *AT=new double[2*Num*6]。 R[1][1]=cos(Omega)*cos(Kappa)。 double Y=R[0][1]*(Data[i][2]Xs)+R[1][1]*(Data[i][3]Ys)+ R[2][1]*(Data[i][4]Zs)。 A[2*i][1]=(R[1][0]*f+R[1][2]*(xxx))/Z。 A[2*i+1][3]=(xxx)*sin(Omega)(((yyy)/f)* ((xxx)*cos(Kappa)(yyy)*sin(Kappa)) f*sin(Kappa))*cos(Omega)。 } } //求ATA MatrixMul(AT,6,amp。 Phi+=Xg[3]。double VTV=0。i6。delete []ATL。a=b。ptemp=matrix。 }}int k。jrow。 } } //列交換 if (js[k]!=k) { for (i=0。irow。//恢復(fù)行列for (r=row1。 } }}ptemp=matrix。 pt++。而編程方面所面臨的問題就比較多了，上一學(xué)期就有過一個類似的作業(yè)，當(dāng)時沒有攻破編程這一題，便用的Excel解決，但事后也沒有繼續(xù)研究，實(shí)為遺憾。 17 。其實(shí)主要總結(jié)起來就一個比較難的點(diǎn)：編程實(shí)現(xiàn)矩陣的運(yùn)算尤其是矩陣求逆。delete []js。for (i=0。r){ if (js[r]!=r) { for (j=0。jrow。i++) { swap(*(m+i*row+k),*(m+i*row+js[k]))。 is[k]=i。for (k=0。irow。}/***********************************************函數(shù)名：InverseMatrix *函數(shù)介紹：求矩陣的逆（高斯約當(dāng)法） *輸入?yún)?shù)：(in/out)matrix（矩陣首地址），*（in）row（矩陣階數(shù)）*輸出參數(shù)：matrix（原矩陣的逆矩陣）*返回值：int ，0成功，1失敗*調(diào)用函數(shù)：swap(doubleamp。delete []mm。}OutPut(Q,mm,m0,Xs,Ys,Zs,Phi,Omega,Kappa,amp。i2*Num。 Kappa+=Xg[5]。 if(InverseMatrix(ATA,6)) return 1。 A[2*i+1][5]=(xxx)。 A[2*i][3]=(yyy)*sin(Omega)(((xxx)/f)* ((xxx)*cos(Kappa)(yyy)*sin(Kappa))+ f*cos(Kappa))*cos(Omega)。 double xxx,yyy。 R[2][0]=sin(Phi)*cos(Kappa)+cos(Phi)*sin(Omega)*sin(Kappa)。double *ATL=new double[6]。Zs=m*f。f){double Xs=0,Ys=0,Zs=0。}fprintf(fp_out,\n)。 for (int j=0。j3。fprintf(fp_out,\tPhi=%10lf\n,Phi)。temp[4])。for (int i=0。fprintf(fp_out,攝影比例尺（0表示其值位置）：)。fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。Phi,const double amp。 } *(JgJz+i*line+j)=temp。,double *JgJz){for (in
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號-1