攝影測量學(xué)單像空間后方交會編程實習(xí)報告(參考版)
2024-08-15 00:19本頁面
【正文】 fprintf(fp,k的精度:%lf\n,mi[5])。 fprintf(fp,q的精度:%lf\n,mi[3])。 fprintf(fp,Ys的精度:%lfm\n,mi[1])。 fprintf(fp,單位權(quán)中誤差的絕對值:%lfm\n,m0)。 fprintf(fp,%.5lf %.5lf %.5lf\n,c1,c2,c3)。 fprintf(fp,%.5lf %.5lf %.5lf\n,a1,a2,a3)。 // fprintf(fp,迭代次數(shù):%d\n,n)。 //保存到文件，結(jié)果文件默認保存在程序目錄 FILE* fp。 printf(w的精度:%lf\n,mi[4])。 printf(Zs的精度:%lfm\n,mi[2])。 printf(Xs的精度:%lfm\n,mi[0])。 printf(\n外方位元素解:\nWs=%.2lf Ys=%.2lf Zs=%.2lf\nq=%.5lf w=%.5lf k=%.5lf\n\n, Xs, Ys, Zs, q, w, k)。 printf(%.5lf %.5lf %.5lf\n,b1,b2,b3)。 printf(旋轉(zhuǎn)矩陣R:\n)。 } //printf(%lf\n,m0)。 i 6。 Q = (()*A).Inverse()。 //單位權(quán)中誤差的絕對值m0以及第i個未知數(shù)的中誤差mi m0 = sqrt(((A*Xx L).Transpose()*(A*Xx L)).ToDouble() / (2 * 4 6))。 n100)。 }while ((abs(Xx[3][0]) || abs(Xx[4][0]) || abs(Xx[5][0]) ) amp。 //printf(%d\n,n)。 k += Xx[5][0]。 q += Xx[3][0]。 Ys += Xx[1][0]。 //求逆 Xx = (ATA_ * AT) * L。 //轉(zhuǎn)置 ATA = AT * A。 L[2 * i + 1][0] =y[i] (f*Y / Z)。 A[2 * i + 1][5] = xx。 A[2 * i + 1][3] = xx*sin(w) (yy*(xx*cos(k) yy*sin(k)) / f f*sin(k))*cos(w)。 A[2 * i + 1][1] = (b2*f + b3*yy) / Z。
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1