作業(yè)4空間后方交會(huì)(編輯修改稿)
2025-04-20 07:03 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 )。fprintf(fp_out,\n空間后方交會(huì)程序(C\\C++)\n測(cè)繪一班\n 學(xué)號(hào)：20080729\n姓名：陳聞亞\n\n)。fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。fprintf(fp_out,已知數(shù)據(jù)：\n\n已知點(diǎn)數(shù)：)。int num。double temp,x,y。fscanf(fp_input,%d%lf,amp。num,amp。temp)。fprintf(fp_out,%d\n,num)。fprintf(fp_out,攝影比例尺（0表示其值位置）：)。fprintf(fp_out,%\n,temp)。fprintf(fp_out,內(nèi)方位元素(f x0 y0)：)。fscanf(fp_input,%lf%lf%lf,amp。temp,amp。x,amp。y)。fprintf(fp_out,%10lf\t%10lf\t%10lf\n,temp,x,y)。for (int i=0。inum。i++){ double temp[5]。 fscanf(fp_input,%*d%lf%lf%lf%lf%lf, amp。temp[0],amp。temp[1],amp。temp[2],amp。temp[3],amp。temp[4])。 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])。}fclose(fp_input)。fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。fprintf(fp_out,計(jì)算結(jié)果如下：\n\n外方位元素：\n)。fprintf(fp_out,\tXs=%10lf\n,Xs)。fprintf(fp_out,\tYs=%10lf\n,Ys)。fprintf(fp_out,\tZs=%10lf\n,Zs)。fprintf(fp_out,\tPhi=%10lf\n,Phi)。fprintf(fp_out,\tOmega=%10lf\n,Omega)。fprintf(fp_out,\tKappa=%10lf\n\n,Kappa)。fprintf(fp_out,旋轉(zhuǎn)矩陣：\n)。for (i=0。i3。i++){ fprintf(fp_out,\t)。 for (int j=0。j3。j++) { fprintf(fp_out,%10lf\t,*(R+i*3+j))。 } fprintf(fp_out,\n)。}fprintf(fp_out,\n單位權(quán)中誤差：%10lf\n\n,m0)。fprintf(fp_out,協(xié)因數(shù)陣：\n)。for (i=0。i6。i++){ fprintf(fp_out,\t)。 for (int j=0。j6。j++) { fprintf(fp_out,%20lf\t,*(Q+i*6+j))。 } fprintf(fp_out,\n)。}fprintf(fp_out,\n外方位元素精度：)。for (i=0。i6。i++){ fprintf(fp_out,%10lf\t,m[i])。}fprintf(fp_out,\n)。fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。fclose(fp_out)。return 0。}/***********************************************函數(shù)名：Resection *函數(shù)介紹：計(jì)算*參數(shù)：Num(點(diǎn)數(shù))，Data(數(shù)據(jù))，m,f(焦距),x0,y0*返回值：int，0成功，其它失敗**********************************************/int Resection(const int amp。Num,const DOUBLE *amp。Data,const double amp。m, const double amp。f){double Xs=0,Ys=0,Zs=0。int i,j。//設(shè)置初始值for (i=0。iNum。i++){ Xs+=Data[i][2]。 Ys+=Data[i][3]。}Xs/=Num。Ys/=Num。Zs=m*f。double Phi(0),Omega(0),Kappa(0)。double R[3][3]={}。double *L=new double[2*Num]。typedef double Double6[6]。Double6 *A=new Double6[2*Num]。double *AT=new double[2*Num*6]。double *ATA=new double[6*6]。d
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1