作業(yè)4空間后方交會-在線瀏覽
2025-05-11 07:03本頁面
【正文】 LE[Num]。double scale=0。for (i=0。i++){ //讀取數(shù)據(jù)，忽略編號 if(fscanf(fp_input,%*d%lf%lf%lf%lf%lf, amp。Data[i][1],amp。Data[i][3],amp。 } //單位換算成m Data[i][0]/=。}//如果m未知則歸算其值if (0==m){ for (i=0。i++) { temp[i]=(Data[i][2]Data[i+1][2])/(Data[i][0]Data[i+1][0])+ (Data[i][3]Data[i+1][3])/(Data[i][1]Data[i+1][1])。 } m=scale/(Num1)。delete []temp。}/***********************************************函數(shù)名：MatrixMul *函數(shù)介紹：求兩個矩陣的積,*參數(shù)：Jz1(第一個矩陣),row(第一個矩陣行數(shù)),*Jz2(第二個矩陣),row(第二個矩陣列數(shù)),(第一個*矩陣列數(shù)),(out)JgJz(存放結(jié)果矩陣)*返回值：void**********************************************/void MatrixMul(double *Jz1,const int amp。line,const int amp。irow。jline。 for (int k=0。k++) { temp+=*(Jz1+i*+k)*(*(Jz2+k*line+j))。 }} }/***********************************************函數(shù)名：OutPut *函數(shù)介紹：*參數(shù)：Q協(xié)因數(shù)陣，m精度，m0單位權(quán)中誤差，6個外*方位元素，旋轉(zhuǎn)矩陣*返回值：int，0成功，1失敗**********************************************/int OutPut(const double *amp。m,const double amp。Xs,const double amp。Zs, const double amp。Omega, const double amp。if (!(fp_out=fopen(,w))){ return 1。if (!(fp_input=fopen(,r))){ return 1。fprintf(fp_out,\n空間后方交會程序(C\\C++)\n測繪一班\n 學(xué)號：20080729\n姓名：陳聞亞\n\n)。fprintf(fp_out,已知數(shù)據(jù)：\n\n已知點數(shù)：)。double temp,x,y。num,amp。fprintf(fp_out,%d\n,num)。fprintf(fp_out,%\n,temp)。fscanf(fp_input,%lf%lf%lf,amp。x,amp。fprintf(fp_out,%10lf\t%10lf\t%10lf\n,temp,x,y)。inum。 fscanf(fp_input,%*d%lf%lf%lf%lf%lf, amp。temp[1],amp。temp[3],amp。 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])。fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。fprintf(fp_out,\tXs=%10lf\n,Xs)。fprintf(fp_out,\tZs=%10lf\n,Zs)。fprintf(fp_out,\tOmega=%10lf\n,Omega)。fprintf(fp_out,旋轉(zhuǎn)矩陣：\n)。i3。 for (int j=0。j++) { fprintf(fp_out,%10lf\t,*(R+i*3+j))。}fprintf(fp_out,\n單位權(quán)中誤差：%10lf\n\n,m0)。for (i=0。i++){ fprintf(fp_out,\t)。j6。 } fprintf(fp_out,\n)。for (i=0。i++){ fprintf(fp_out,%10lf\t,m[i])。fprintf(fp_out,************************************** ************************************** ************************************** *********************************\n)。return 0。Num,const DOUBLE *amp。m, const double amp。int i,j。iNum。 Ys+=Data[i][3]。Ys/=Num。double Phi(0),Omega(0),Kappa(0)。double *L=new double[2*Num]。Double6 *A=new Double6[2*Num]。double *ATA=new double[6*6]。double *Xg=new double[6]。 R[0][1]=cos(Phi
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號-1