第一篇：GPS实验报告
卫星导航定位算法与程序设计
实验报告
实验一 时空基准转换
一、实验目的1、加深对时空系统及其之间转换关系的理解
2、掌握常用时空基准之间的转换模型与软件实现
二、实验内容
1、编程实现GPS起点1980年1月6日0时对应的儒略日
2、编程实现2011年11月27日对应的GPS周数与一周内的秒数
3、在WGS84椭球的条件下，编程实现当中央子午线为117度时，计算高斯坐标
x=3548910.811290287，y=179854.6172135982对应的经纬度?
4、在WGS84椭球的条件下, 表面x=-2408000,y=4698000,z=3566000处地平坐标系坐
标为:e=704.8615,n=114.8683,u=751.9771的点对应的直角坐标?
三、实验过程
1、这是测试 站心坐标系 空间直角坐标系 相互转换高斯正算和高斯反算 的主程
序功能模块
（1）调用了enu2xyz 站心坐标向空间坐标的转换函数
（2）调用了xyz2enu空间坐标向站心坐标转换的函数
（3）调用了高斯反算函数gauss_fansuan2、时间转换主程序功能模块
disp('1980年1月6日0时对应的儒略时:');
jd1 = julday(1980,1,6,0);
disp(jd1);
disp('2011年11月27日对应的GPS周和周内秒:');
jd2 = julday(2011,11,27,0);
[week,sow] = gps_time(jd2);
disp('week:')
disp(week);
disp('sow:');
disp(sow);
四、实验体会
本次实验中我对坐标转换问题有了更多的了解，同时对matlab的一些基本函数更加了解，对时间转换内容的熟悉程度也有所加深。通过本次实验我学会了matlab程序的调试方法，为后继学习打下了基础。
实验二 RINEX文件读写
一、实验目的1、深入了解RINEX文件格式
2、进一步提高MATLAB程序设计能力
3、掌握N文件、O文件、SP3文件的基本读写技巧
二、实验内容
1、任选IGS站，下载N文件、O文件与SP3文件
2、编程实现N文件的读入，并采用中文标注出主要参数的名称及作用
3、编程实现O文件的读入，并采用中文标注出主要参数的名称及作用
4、编程实现SP3文件的读入，并采用中文标注出主要参数的名称及作用
三、实验过程
1、N文件读取主程序模块功能
（1）调用rinexe函数将文件SITE247J.01N中的数据按一定的格式读入到文件
eph.dat中
rinexe('SITE247J.01N','eph.dat');
（2）调用get_eph函数将文件eph.dat中的数据读入到变量Eph中
Eph = get_eph('eph.dat');
（3）输出得到的Eph变量的内容到命令窗口
disp(Eph);
（4）rinexe函数及get_eph函数具体注释见源代码
2、SP3文件、O文件读取主程序模块功能
读取O文件主要过程：
（1）打开文件，获取头文件的信息。利用anheader获取观测值类型、天线高及偏心及找到相对应观测类型的标识码和判断文件是否结束的标识码。
（2）在函数体中利用findstr()及循环语句逐行查找字符串“END OF HEADER”，找到头文件结束的位置。
（3然后逐行在头文件中寻找字符串“ANTENNA: DELTA H/E/N”、”# / TYPES OF OBSERV”等，得到天线高及偏心高矩阵、观测值类型、观测值个数和观测类型的标识码。
（4）获取观测历元的信息。利用fepoch_0获得某一个历元的观测时间和卫星编号矩阵及文件结束标识符。然后可以从卫星的编号矩阵sats中求得卫星数。
（5）获取某一观测类型（如程序中的P2）的信息。利用grabdata.m函数,根据卫星数NoSv和观测类型数NoObs返回某历元的5*5“obs”观测数据矩阵,再利用循环获得所有历元的观测数据矩阵。
（6）利用fobs_typ.m获取观测类型为P2的观测值矩阵中对应列的信息。
涉及的主要函数及功能：
anheader.m：分析RINEX文件的头文件，输出观测类型和天线设置类型，及找到相应观测类型的标识码。
fepoch_0.m:在打开的RINEX文件中利用识别标志fid寻找下一历元。从历元序列中生成历元时间、卫星编号、文件结束符。只有O文件才会被处理。
grabdata.m：在选定的历元读取编号卫星对应的观测值。
fobs_typ.m：返回观测值矩阵中的列i，其中包含观测类型“Type”。
读取SP3文件主要过程：
（1）打开文件，然后读取文件，利用ReadSP3.m以矩阵的形式返回卫星的轨道信息、卫星数目、及文件头中的各种参数值。
（2）读取精密轨道信息。利用for循环，得到卫星的编号矩阵以及其对应的轨道的位置及种差数据。
（3）利用FormSP3.m根据指定的卫星编号以及起止时间读取SP3文件，并返回相关信息,以矩阵的形式输出。
涉及的主要函数及功能：
ReadSP3.m：读取SP3精密卫星轨道数据文件
调用方式：[SP3data,numsat,header]= ReadSP3(filename)
SP3data:轨道信息
Numsat：卫星数目
Header：文件头
FormSP3.m：读取SP3精密卫星轨道数据文件，并按卫星编号存为矩阵形式。
调用方式：[Time,SP3X,SP3Y,SP3Z,Clk,remark,new_prn]= FormSP3(filename,PRN,sep,lep)根据指定的卫星编号以及起止时间读取SP3文件filename，并返回相关信息。
四、实验体会
本实验主要是学习对N文件、SP3文件、O文件的读取，在实验过程中我学习到了很多数据的读取方法，这个实验之前我也遇到过好多数据读取方面的问题，学习了这个程序代码后，发现它读取数据的很多思路非常的巧妙。通过本次实验，各种文件格式数据的结构有
了深入的了解。
实验三卫星轨道计算
主要过程：
整个程序的实现分三个模块实现：时间转换模块；读广播星历文件(N文件)模块；计算卫星位置模块。
1、时间转换模块：
将O文件的某一历元对应的格里高利历时(实验中的时间：2001,9,4,9+40/60)转换为儒略时，然后将儒略时转换为GPS时得到对应的周数(WN)和周内秒数(TOW)。
主要涉及的函数及功能：
1）Julday.m：将格里高利历时转换为儒略时。
2）gps_time.m:将儒略时转换为GPS时的到对应的周数和周内秒数。
2、读星历文件模块：
在求解卫星位置时，第一需要利用O文件中每个历元的历元时刻t。在计算某时刻卫星的位置时，这里的某时刻便是O文件历元时刻t。第二需要利用读取的每个历元不同的卫星PRN号。根据PRN号和历元时刻t在广播星历N文件中搜索相同的卫星PRN号、合适的历元时刻，利用其对应的数据，计算卫星的位置。
先利用rinexe.m函数读取N文件将导航信息中的每颗卫星的21个参数以矩阵的格式存入文件eph.dat中，然后利用get_eph.m函数将文件eph.dat中的星历信息存入矩阵Eph中。
主要涉及的函数及功能：
1）rinexe.m：读取rinexe导航信息文件，将信息以矩阵的格式写入矩阵。矩阵中有21行，每一列存储一颗卫星的信息。
调用方式：rinexe(ephemerisfile, outputfile)，ephemerisfile为文件名(N文件要放在rinexe函数所在的目录下)
outputfile为文本文件名，将读到N文件的数据存放在此文本文件中。
2）get_eph.m:星历文件中包含的星历信息被存入矩阵eph中。矩阵行数为21行，列数为有星历信息的卫星个数。
调用方式: eph = get_eph(ephemeridesfile)
Ephemeridesfile——文件名，将存放在此文本文件中的导航电文数据提取出来，为计算卫星的位置做准备。该函数返回值为eph星历矩阵(为21*7的矩阵)，它为Matlab内部格式。
3、计算卫星位置模块：
主要用到四个函数，分别为读O文件函数、eph星历矩阵的选择函数、修复函数(修复GPS时间超限或者溢出)、计算卫星的位置函数。
根据要寻找的历元所包含的卫星PRN号，在N文件中搜索对应的卫星PRN号、合适的历元时刻，利用其对应的数据，计算出卫星的位置。
利用find_eph.m函数从eph矩阵中找出相应卫星编号的卫星的参数所在的列数，接着利用satpos.m函数来该颗卫星对应的X、Y、Z坐标值，即在地心地固坐标系中的坐标。
主要涉及的函数及功能：
1）find_eph.m：从Eph矩阵中选中用于计算的列数。
调用方式：col_Eph(t)= find_eph(Eph,svs(t),sow)
Eph——星历矩阵，存放计算卫星位置所需要的N文件的数据。
svs(t)——svs中存放O文件中某个历元观测到的所有卫星PRN号。t为循环控制，svs(t)为svs中的某刻卫星的编号。
Sow——选定的历元的周内秒数值。
返回值——col_Eph矩阵中存放着从Eph矩阵中选中用于计算的列数。
2）X=satpos(tx_GPS,Eph(:,k))
tx_GPS——历元的格里高利历时转化为GPS时后的秒数sow（即所谓的归化后的时间。
Eph(:,k)——Eph星历矩阵中的对应的列的数据。
返回值——卫星在地心地固坐标系中的坐标
在本实验中的具体调用方式：
sat(1:3,t)= satpos(sow,Eph(:,col_Eph(t)))
其中Eph(:,col_Eph(t))为上一步中用col_Eph.m函数并找到的相应卫星在Eph矩阵中对应列的数据。
3）tk = check_t(t-toe)
t-toe——为儒略日
返回值为修复后的儒略日。
4）计算卫星在地固坐标系中的三维坐标的在程序中的具体实现如下： satp(1,1)= x1*cos(Omega)-y1*cos(i)*sin(Omega);
satp(2,1)= x1*sin(Omega)+y1*cos(i)*cos(Omega);
satp(3,1)= y1*sin(i);%计算卫星地心固定坐标系中的直角坐标 附：
1、在编计算每颗卫星位置的程序时，每颗卫星都需要存储导航电文中的21个参数，分别为：
svprn,af2,M0,roota,deltan,ecc,omega,cuc,cus,crc,i0,idot,cic,cis,Omega0,Omegadot,toe,af0,af1,toc。详细情况可以参照代码进行调试理解。
2、实验内容中要求的第3条“根据卫星位置计算公式编写主函数，同时调用时间转换模块、星历读取等的子函数来共同完成卫星位置的计算，最后输出结果”，具体的实现情况可以参照已经给出的代码。
实验四 伪距定位及DOP值
一、实验目的1、掌握伪距离定位的基本原理并了解其公式推导基本流程
2、熟练掌握线性化的基本公式及内涵
3、掌握DOP计算的基本步骤并绘制时段内星座DOP的分布
二、实习内容
1、独立调试一段线性化程序，程序的输出为A矩阵，用于最小二乘求解与DOP值的计算；
2、构建A矩阵并提取对应DOP计算所需元素
3、空间直角坐标系转换至大地直角坐标系的R矩阵的编写
4、绘制指定接收机位置及指定时段内的星座DOP分布图
三、实验过程
1、主程序功能模块
（1）调用函数mpgen计算多路径效应误差；
（2）调用函数llh2xyz将用户坐标转换为WGS4坐标系统下的坐标；
（3）调用函数loadgps加载gps数据；
（4）调用函数olspos计算方向余弦矩阵H；
（5）调用函数plot绘制DOP值图像
2、olspos函数功能模块
根据伪距观测值和卫星坐标计算用户坐标和方向余弦矩阵H；
该函数还利用循环进行迭代平差；
while((iter
for N = 1:numvis
pr0 = norm(svxyzmat(N,:)-estusr(1:3));
y(N,1)= prvec(N)estusr(4);
end
H = hmat(svxyzmat,estusr(1:3));
beta = Hy;
estusr=estusr+beta';
iter=iter+1;
end3、hmat函数模块
根据伪距方程线性化后的函数，求出三个该正数的系数，然后组建h矩阵;
(Xj-X)2(Yj-Y)2(Zj-Z)2
(Xj,Yj,Zj)为第j颗卫星在空中的位置；
(X,Y,Z)为接收机空间的位置。
线性化后得到：
Xj-X0Yj-Y0Zj-Z0 00dx0dy0dz
四、实验体会
通过本次实验的学习，我对伪距离定位的基本原理并及其公式推导基本流程有了更清楚的了解。同时也掌握了线性化的基本公式及内涵、精度因子计算的基本步骤，能够绘制了时段内星座DOP的分布图。
第二篇：GPS实验报告
实验一：GPS静态测量实验
实验目的：
1、掌握天宝GPS接收机的操作。
2、掌握GPS静态相对定位数据采集方法。
3、掌握卫星预报软件的使用方法。
4、掌握数据传输与后台处理软件的使用方法。实习任务：对已有控制点进行多时段静态测量 实验步骤：
放置脚架，对中整平，安置好仪器。 量取天线高
打开接收机电源，接收机跟踪大于4颗以上卫星时，卫星指示灯慢闪；打开数据记录灯；此时开始记录数据。（注：一定要保证数据记录灯亮，否则没有记录数据） 认真填写外业记录表
结束测量时，先关闭数据记录灯，再关闭接收机电源。
2.静态数据内业处理
（1）接收机的数据传输
关于外业观测数据的传输，比较特别的是，Trimble 5700接收机的数据传输需要安装Data Transfer数据传输软件才能实现传输。
(2)将trimble接收机的数据文件转成RINEX格式
安装好Convert to RINEX软件后，运行，选择好要转换的trimble数据文件，如图：
点击“编辑”，对相关参数进行设置，选择观测方法为“护圏的中心”，并根据外业观测记录表，填好初始天线高，点击“改正”即可。设置完成后，就是进行格式转换了。
（3）HGO软件，新建项目，选择相应的坐标系统 如图：
(4)处理基线观测残差序列图和基线处理
这一个环节，主要是通过查看基线的残差序列图来初步判断该基线的质量好坏。质量控制只作为了解，是基线解算质量的三个恒量标准，即比率（ratio）、均方根（RMS）。我们主要通过屏蔽某段信号或者某颗卫星的信号来使得ratio值和 RMS值增大，ratio值越大越好，信号好的话，ratio值一般在50-100之间。RMS值越小越好，信号好的话，RMS值一般会在0.005左右。如图:
未屏蔽信号前：ratio值为16，RMS值为0.0062
屏蔽信号，修复完成后，ratio值为99，RMS值为0.0045.这里要特别注意的是，在屏蔽卫星信号时，必须保证同一观测时段必须有不少于4颗卫星。同时，在屏蔽基线时必须保证每个测站至少与两条基线相连。
按照同样的方法处理每天基线，使其合格。对于多次处理仍未合格的基线，可改变其高度截止角和采样间隔的参数使得其合格。
（5）闭合环和重复基线处理
1）同步环处理：
由于2009版的GPS测量规范对同步环的要求太高了，所以在实际工作当中，我们不要求每个同步环都合格。
2）异步环处理：
必须保证每个异步环都合格。对于不合格的异步环，选中该异步环，对构成该异步环的若干基线进行再处理，方法跟上面处理基线的一样，先处理ratio值比较小，而RMS值比较大的基线，重新解算该基线，再看异步环有没有合格。按照这个方法，直到所有的的异步环合格。
（6）平差
1)选取特定的控制点
为保证整网的点位精度均匀，起算点一般应均匀地分布在GPS网的周围,避免所有的起算点分布在网中一侧的情况。故这次中，我们组选取了45、46、50这三个点作为控制点。当在HGO中把这三个点转为控制点后，我们要在其属性里，根据任务书给出的已知点坐标，依次输入坐标，需要注意的是由于我们已知坐标是国家80坐标系下的，所以在HGO的对话框里选择的是“国家坐标”，另外，由于我们的点位都位于同一投影带上，所有东坐标都不需要输带号（不在同一投影带才需要输入换带坐标）。如图：
2)平差：先进行平差报告设置
依次进行各项平差，如自由网平差，二维约束平差。平差完成后，系统对话框就会显示平差的结果，若X2检验不通过，我们可以在平差设置里输入系统提示的松弛因子来使其通过。若Tau检验未通过，我们可以点击生成报告，查看未通过的基线，再处理该基线（处理方法同上面的基线处理方法），最后再平差即可。平差通过后，我们可以导出自己的项目总结报告。
实验二：RTK测量
实验目的：掌握利用GPS-RTK技术进行地形测绘和放样等基本测绘工作的方法，培养和提高利用所学GPS理论知识解决实际问题的能力；加深对测绘事业的理解，树立为测绘事业奋斗的精神；培养积极进取、勇于探索、团结协作、对工作一丝不苟、吃苦耐劳的品德和作风，提高业务素质和竞争能力。实验任务：
1、完成控制点的测量
2、进行坐标系之间的换算
3、进行点放样
实验步骤：
1、把手机卡装到仪器上，开机设置仪器参数模式为【GSM移动台】。2.把手部连接到仪器上，在手部中点击{GPS】，搜出仪器对应的型号，然后连接，接着输入广工的基准站IP地址，手簿显示屏上方会显示单点，伪距，浮动，固定，当显示为固定时表示接收机已接受到服务器的信号，并信号稳定，可以开始测量。
3、新建项目,在菜单中选择项目，然后新建一个项目
新建项目完成后，打开该项目即可。设置好参数就可以开始了。
4.控制点测量以及参数解算
在我们使用RTK时，我们没有启用任何参数所得到的直角坐标往往是不准确的，所以我们还有一个很重要的步骤就是求解参数，参数有几种，工程上常用的一般就是四参数和高程拟合，求解四参数之前，我们必要的条件是：在测区至少要有两个以上的已知点。
那么本次作业中，求解过程如下：首先用移动站在15号和16号两个已知点上采集两个没有参数的坐标，接着退出测量界面，进入参数界面，选择“参数计算”。在完成参数计算后，我们还选择了16号作为检验点，检验结果符合要求后，我们才正式开始下面的其他测量工作，此时得到的坐标就是和已知点相同坐标系的我们需要的直角坐标。如图：
5.放样
首先，在测量界面点击左上角下拉菜单，选择“点放样”，进入放样模式 在放样的时候，我们直接进入点库中选中我们之前已经设计好的放样点进行放样，按照手簿提示的方位和距离来移动我们的移动站，从而找到该点进行放样。其他点的放样一样按照这个方法。
需要注意的是，在放样的时候，需要解的状态为“固定”，并且值要不大于3.0。
四．实验总结：
通过这个RTK实验，自己更加了解了RTK的使用方法和操作流程，我学会了运用所学知识解决处理简单问题的方法与技巧，一个人的知识和能力只有在实践中才能发挥作用,才能得到丰富、完善和发展。
第三篇：静态GPS实验报告
GPS静态测量实习报告 学院：
专业：
学号：
姓名：
组号：土木与水利工程学院测绘工程09-1班王震阳20094176C6
一、实验目的和要求
⑴了解静态GPS测量系统的组成部分。
⑵了解静态GPS相对定位模式的作业方法。
⑶了解GPS观测数据在计算机上的处理过程。
二、实验仪器和工具
1.外业：南方s82静态GPS接收机1套（7台）、已充电锂电池7块、对点器基座7套、铝合金三脚架7个。
2.内业：数据传输电缆1根、数据处理软件（南方
GPS）（含采集器与计算机通讯软件、基线向量处理软件、网平差及坐标转换软件）
三、实验环境
外业: 合肥工业大学南区
内业：南方GPS数据处理软件
四、实验成果
见《C6组GPS成果报告》
五、实验小结
可以说我们本次GPS静态测量实习，一共进行了三次尝试才最终将成果打印出来。前两次的失败有仪器的原因有人为的原因。第一次失败是因为学校的GPS仪器注册码过期导致无法继续工作而中途废弃。回到宿舍后我们给南方公司要了注册码将班级的所有仪器均注册了一下才使仪器可以使用。第二次失败是因为我们在观测的过程中有一台仪器没有电了，我们换了一节电池后，回来发现基线解算不合格。经过多次尝试还是没有合格后，我们组决定重新再来，进行了第三次GPS测量，这一次我们在前两次失败的经验教训之上终于完成了合格的解算工作。
正真正的实习过程中总会出现各种各样的错误，有些可能是我们经验不足导致的，有些可能是外界原因导致的。反正不管什么原因我们都要有能力去探索解决问题的方法。
第四篇：GPS 实验报告
GPS实习报告
学院：资源与环境学院 专业： 测绘工程 班级： 姓名：
学号： 指导老师：
第一部分：GPS控制网外业观测设计
一 项目概况
这次实习是由张老师带领我们在学校及其周围进行GPS测量。本次实习包括GPS静态观测，GPS实习差分动态测量（RTK），手持便携GPS导航的认识与应用。这些实习巩固、扩大和加深我们从课堂上所学理论知识，获得测量工作的初步经验和基本技能，着重培养我们的独立工作能力，进一步熟练掌握测量仪器的操作技能，提高运用理论及计算能力。工程测区：华北水利水电学院，东风渠，黄河测绘局等地区。
测区概况：地势较为平坦，测区内高楼比较多。有一条河渠横穿过测区。
项目开始时间：2013年3月9日 项目结束时间：2013年3月17日
二 技术依据
（1）采用西安80坐标系
（2）《全球定位系统（GPS）测量规范》 三 现有测绘资料
（1）《全球定位系统（GPS）测量规范》
（2）已知点位信息
四 选点情况
本次选点是由张老师选定，共有8个点位，分布在学校及其周围。点位远离大功率无线电发射机，避开了大型建筑物，减少了多路径效应，视野开阔，适合进行GPS观测。具体情况如下：
观测方案及质量控制方法 1，采样间隔15s，截至高度角为15°。
2，观测期间不得在天线附近50m范围内使用电台，不得在10m范围内使用对讲机或手机。
3，用于GPS基准网观测的接收机必须是符合GPS规范要求的双频机，其标称精度为5mm+-1*10
-6。
五 观测方案及质量控制方法
本次实习的精度为E级，由三台GPS同步接收机，组成同步闭合环进行同步观测，每一时段观测50分钟以满足精度要求。每时段观测前后都应量取一次仪器高。
为提高观测质量，同一时段观测过程中不允许出现接收机重新启动，进行自测试，改变卫星高度角，改变数据采样间隔，关闭或删除文件等操作。
五 提交成果资料
（1）点之记
（2）GPS静态测量外业观测手簿（3）GPS实习报告
第二部分：GPS网控制网技术总结
一 外业观测情况
本次测量时采用的是，Trimble GPS接收机三台。采用同步观测的相对定位方法，从而保证了卫星星历误差、卫星钟误差、电离层延迟等误差的强相关性，通过差分的方法来消除这些误差。观测时为了保证测量的精度时段长度规定为50分钟。按照静态定位的测量原理，测量时观测的最少卫星数位四个。外业观测时需要对GPS接收机进行以下设置：
(1)调度安排，确定每台接收机观测的测站，开机时间，搬站情况。
(2)观测组按调度表规定的时间进行作业，保证同步观测同一卫星组。
(3)每时段开机前，量取天线高，并及时记录测站名、年月日、时段号、天线高等信息。关机后再量取一次天线高作校核。
(4)仪器工作过程中，作业人员对照指示灯工作状况说明，判断仪器是否正常工作。
(5)一个时段观测过程中，不得进行以下操作：关闭接收机，又重新开机；进行自测试；改变卫星高度角；改变数据采样间隔；改变天线位置；
(6)观测院在作业期间不得擅自离开测站，并应防止仪器收到震动，防止人或其他物体靠近天线，遮挡卫星信号。
(7)接收机在观测过程，不应在接收机旁使用对讲机；雷雨天气过境应关机停测，卸下天下以防雷击。
(8)应记录雨、晴、阴、云等天气。二 数据处理
本次实习的数据处理是应用南方GPS软件进行外业观测数据处理，将数据导入该软件后先进行基线结算，这将由计算机自动完成，然后再对不合格数据进行调整，然后输入已知点数据进行网平差，最后输出成果。本次成果第一次基线解算时出现里超限数据，需要改动卫星高度角或数据采样间隔来调整以满足精度要求。三 数据处理结果平差后坐标和点位精度
ID X坐标 Y坐标 rms(mm)dx(mm)dy(mm)点 名
00B2 3852636.764 469463.678 1.640 1.171 1.148 00B2 00G2 3852371.932 470192.135 1.304 1.020 0.812 00G2 GA02 3853810.984 469690.367 1.166 0.787 0.860 GA02 GA01 3854657.352 469424.311 0.000 0.000 0.000 GA01 GA03 3854114.425 468956.490 0.000 0.000 0.000 GA03 00B6 3854761.598 468733.367 1.893 1.341 1.336 00B6 GA05 3853798.608 470095.209 0.000 0.000 0.000 GA05 00G4 3855827.674 468324.367 2.382 1.714 1.654 00G4 1 参数拟合高程 28.858824 内符合精度中误差±27.077(mm)拟合后高程残差
点号 正常高(高程)大地高 正常高(拟合)差值 Rms(mm)GA01
91.258
62.411
91.270 0.012 1.377 GA03
92.333
63.445
92.303-0.030 1.363 GA05
90.200
61.362
90.221 0.021 1.452 拟合高程
ID 正常高(高程)大地高 Rms(mm)点 名 00B2 91.275 62.416 0.757 00B2 00G2 92.773 63.914 1.290 00G2 GA02 90.844 61.985 1.291 GA02 00B6 93.791 64.932 1.493 00B6 00G4 93.703 64.844 1.590 00G4 ID 坐标 X 坐标 Y 高 程 x y h 点 名
GA01
3854657.352
469424.311 91.258 * * * GA01 GA03
3854114.425
468956.490 92.333 * * * GA03 GA05
3853798.608
470095.209 90.200 * * * GA05 00B2
3852636.764
469463.678 91.275 00B2 00B6
3854761.598
468733.367 93.791 00B6 00G2
3852371.932
470192.135 92.773 00G2 00G4
3855827.674
468324.367 93.703 00G4 GA02
3853810.984
469690.367 90.844 GA02 7 参数结果
Dx平移(米):-35.106 Dy平移(米): 80.362 Dz平移(米): 64.242 Rx 旋转(秒): 1.144568 Ry 旋转(秒):-0.609468 Rz 旋转(秒): 1.452463 SF尺度(ppm): 4.573765 第三部分：动态测量实习总结 本次实习是在校园内进行，基准站和流动站都是Trimble公司产品，这次实习是在学校内选定五个点位，先在学校内一开矿地处架设基准站，使用GPS观测手簿连接基准站，在连接到基准站后再连接流动站，并且在GPS观测手簿中输入这五个点位的已知坐标，利用流动站先选取三个点位进行校准，在利用另外两个点位进行放样，通过这次实习让我了解了利用RTK进行测量的步骤，锻炼了我的实践能力，为我以后的工作打下了良好的基础。
第四部分：GPS导航测量总结
这次实习利用GPS导航机进行路线导航，距离测量等，这次我利用GPS导航手簿量测了我宿舍到GA05的距离，发现GPS导航手簿的强大功能及给我们生活带来的便利。
第五部分：实习总结
本次GPS实习虽然仅仅只有一个星期的时间，主要包括两大部分——GPS静态测量和动态测量。其中静态测量包括实地勘踏选点、调度方案的确定、正式外业数据采集、最后也就是最重要的内业处理。动态测量包括利用GPS RTK技术进行RTK测量以及利用RTK进行放样操作。实习过程中我们收获还是十分丰富的。对GPS静态测量和动态测量有了深刻的理解，将之前停留在书本知识阶段延伸到了实际操作阶段，实习过程中，不仅对原理有了更进一步的理解，还在操作中熟练了步骤，为我们以后的工作打下了良好的基础。
当然，这次实习也出现了一些问题，实习之后更要认真总结，我相信这也必将成为我们日后工作和学习中的宝贵财富。
第五篇：北航物理实验GPS模拟与声源定位实验报告
声源定位C++程序
#include
cout << “声源坐标为：（”< if(acos(D/sqrt(pow(A,2)+pow(B,2)))+angle2 < pai/2)else angle1 = angle2x1*(pow(x2,2)+pow(y2,2)-pow(c*dt2,2));B = y2*(pow(x1,2)+pow(y1,2)-pow(c*dt1,2))c*dt2*(pow(x1,2)+pow(y1,2)-pow(c*dt1,2));double x0,x1,x2,y0,y1,y2,t0,t1,t2,t3,A,B,D,angle1,angle2,dt1,dt2,r;x0=0,x1=0,x2=0.45;y0=0,y1=0.3,y2=0;int i = 8;while(i--){ cout <<“请输入第<<8-i<<”组时间数y据Y“<< endl;cin >> t0 >> t1 >> t2>> t3;t0 = t0 / 1000000;t1 = t1 / 1000000;t2 = t2 / 1000000;dt1 = t1t0;} // end main; 程序运行: 2、GPS模拟 C++程序： #include void GPS(double x[10] , double y[10] , double t[10] , double x1 , double y1 , int num){ double a1,a2,a3,a4,b1,b2,mult,D=0,E=0;int i;while(1){ a1 = a2 = a3 = a4 = b1 = b2 = 0;for(i = 0;i < 10;i ++){ a1 = a1-(3*pow((x[i]-x1),2)+pow((y[i]/1000-y1),2)-pow(c*t[i]/1000000,2));a2 = a2+2*(x[i]/1000-x1)*(y1-y[i]/1000);a3 = a2;a4 = a4-(3*pow((y[i]-y1),2)+pow((x[i]/1000-x1),2)-pow(c*t[i]/1000000,2)); b1 = b2 = b1+(pow((x[i]/1000-x1),2)+pow((y[i]/1000-y1),2)-pow(c*t[i]/1000000,2))*(x1-x[i]/1000);b2+(pow((x[i]/1000-x1),2)+pow((y[i]/1000-y1),2)-pow(c*t[i]/1000000,2))*(y1-y[i]/1000);} //end for;mult = a3/a1;a4-= a2*mult;b2-= b1*mult;E = b2/a4;D =(b1-a2*E)/a1;x1 += D;y1 += E;if((pow(D,2)+pow(E,2))< min)break;} // end while;x1 = x1*1000;y1 = y1*1000;cout <<”点P“<< num <<”定位于: X = “<< fixed << setprecision(8)<< x1;cout <<”(mm), Y = “<< fixed << setprecision(8)<< y1 <<”(mm)“< GPS(x,y,t1,0,0,1);GPS(x,y,t2,0,0.3,2);GPS(x,y,t3,0.45,0,3);cin>>yy;return 0;int yy;double t1[10] , t2[10] , t3[10] , x[10] , y[10]; cout << ”请输入时间数据Y(us)与坐标数据Y(mm):“ << endl;int i;for(i = 0;i <= 9;i ++){ cout << ”第“ << i+1 << ”组" << endl;cin >> t1[i] >> t2[i] >> t3[i] >> x[i] >> y[i];} // end for; } // end main; 运行结果：
点击咨询 