单片空间后方交会编程

《单片空间后方交会编程》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、单片空间后方交会编程实习报告指导老师：袁修孝付迎春PRS032陈钧20033259002481.实习目的用程序设计语言（VisualC++或者C语言）编写一个完整的单片空间后方交会程序，通过对提供的试验数据进行计算，输出像片的外方位元素并评定精度。本实验的目的在于让我们深入理解单片空间后方交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过程。通过上机调试程序加强动手能力的培养，通过对实验结果的分析，增强学生综合运用所学知识解决实际问题的能力。2.实习内容利用一定数量的地面控制点，根据共线条件方程求解像片外方位元素。3.数据准备1.已

2、知航摄仪的内方位元素：fk=153.24mm，x0=y0=0.0mm，摄影比例尺为1:50000；2.4个地面控制点的地面坐标及其对应像点的像片坐标：点号像片坐标(mm)地面点坐标(m)xyXYZ1-86.15-68.9936589.4125273.322195.172-53.4082.2137631.0831324.51728.693-14.78-76.6339100.9724934.982386.50410.4664.4340426.5430319.81757.31程序运行后的输出结果如下：输入控制点的像点坐标:第1个点的像点坐标为：（-86.15，-68.99

3、）第2个点的像点坐标为：（-53.40，82.21）第3个点的像点坐标为：（-14.78，-76.63）第4个点的像点坐标为：（10.46，64.43）输入控制点的地面坐标:第1个点的地面坐标为：（36589.41，25273.32，2195.17）第2个点的地面坐标为：（37631.08，31324.51，728.69）第3个点的地面坐标为：（39100.97，24934.98，2386.50）第4个点的地面坐标为：（40426.54，30319.81，757.31）4.操作步骤画流程图，上机调试程序并打印结果。88程序实现：#include"stdio.h"#i

4、nclude"stdlib.h"#include"iostream.h"#include#defineN4//观测点数为4voidmult(double*m1,double*m2,double*result,inti_1,intj_12,intj_2)//矩阵相乘{inti,j,k;for(i=0;i

5、_matrix(double*m1,intn)//矩阵求逆{int*is,*js;inti,j,k,l,u,v;doubletemp,max_v;is=(int*)malloc(n*sizeof(int));js=(int*)malloc(n*sizeof(int));if(is==NULL

6、

7、js==NULL){printf("outofmemory!");return(0);}for(k=0;kmax

8、_v){max_v=temp;is[k]=i;js[k]=j;}}if(max_v==0.0){free(is);free(js);printf("inversisnotavailble!");return(0);}if(is[k]!=k)for(j=0;j

9、m1[l]=1.0/m1[l];for(j=0;j=0;k--){if(js[k]!=k)for(j=0;j