GPS水准高程拟合报告
GPS高程拟合方法及精度分析
GPS高程拟合方法及精度分析GPS（全球定位系统）是一种通过卫星进行定位的导航系统，它通过接收地面上的GPS 接收器收集到的卫星信号来确定接收器的位置。
GPS系统不仅可以提供经度和纬度等位置信息，还可以提供高程信息。
在实际应用中，由于各种误差的存在，GPS高程数据往往需要进行拟合处理，以提高其精度。
GPS高程拟合方法主要有以下几种：1.大地水准面拟合法：该方法假设地球上存在一个水准面，通过高程数据与该水准面的差值来进行拟合。
大地水准面拟合法可以根据地球椭球体模型进行，也可以根据区域地形特征进行。
2.多项式拟合法：该方法通过将GPS高程数据与多项式函数进行拟合，来估算出真实的地理高程。
多项式拟合法常用的模型有一次、二次和三次多项式，其拟合误差随着多项式的阶数增加而减小。
3.高斯滤波法：该方法考虑到GPS高程数据的时序性，通过滤波算法对数据进行平滑处理，以提高高程数据的精度。
高斯滤波法利用高斯函数对数据进行加权平均，同时考虑到观测误差的方差，使得滤波结果更加符合实际情况。
1.接收器误差：GPS接收器的误差包括时钟误差、接收机硬件误差等，这些误差会直接影响到GPS高程数据的精度。
2.卫星误差：卫星的轨道误差、钟差误差等因素也会对GPS高程数据的精度产生影响。
3.大气误差：由于大气对GPS信号的传播会产生延迟和折射等误差，因此对GPS高程数据的精度也会有一定的影响。
4.数据后处理方法：不同的数据后处理方法对GPS高程数据的精度有着较大的影响。
合理选择数据处理方法可以提高GPS高程数据的精度。
为了提高GPS高程数据的精度，在采集数据时需要注意选择合适的接收器和卫星，并进行数据后处理以减小误差。
还可以通过与地面高程标志点对照来校正高程数据，以获得更高的精度。
利用水准点进行GPS高程拟合几种算法及精度分析
利用水准点进行GPS高程拟合几种算法及精度分析摘要:利用已知水准点高程，进行GPS大地高向正常高转换，其精度受似大地水准面、已知点高程和GPS网点的大地高三种误差的影响关键词:GPS高程拟合；计算方法；拟合精度。
Abstract: The use of known standard point elevation for GPS Height Conversion to normal, its accuracy is affected by geoid, elevation and GPS networks known point of the earth-and high-impact errorsKeywords: GPS Elevation Fitting; calculation method; fitting accuracy.一、水准高程与GPS程在水准测量中采用的正常高h是地面点沿铅垂线到似大地水准面的高度，即以不规则的有起伏的重力等位面为基准面，具有严格的物理意义，正常高可由水准测量结合重力测量得出；而GPS所测量的高程是沿法线方向到WGS84椭球面的高度H（大地高），即以简单的数学曲面为基准面，具有明确的几何意义但缺乏物理意义。
这两种基准面是不一致的，它们之间的差值称为高程异常，其关系式为：ξ=H一h式中:ξ—高程异常，表示似大地水准面参考椭球面的距离;H—大地高;h—正常高。
采用静态方法进行GPS测量后，由GPS三维平差可得到施测点的大地高，同时在所施测GPS控制网中联测部分水准点，则这些点的大地高H、正常高h 是已知的，即可求得这些点的高程异常。
在一定范围内可以认为高程异常变化平缓，但在此范围内高程异常不为常数，因此可以使用一些数学函数来进行拟合，求得能反映GPS网控制范围中高程异常变化的函数，然后通过内插求得GPS网点中个点的高程异常，从而得到控制点的拟合水准高程。
研究GPS高程拟合的意义如下：①通过GPS控制点的大地高通过拟合精确求定正常高；②求定高精度的似大地水准面；③求定不同位置的准确高程异常值。
GPS拟合高程测量
GPS 拟合高程测量一、GPS 拟合高程测量，仅适用于平原或丘陵地区的五等及以下等级高程测量。
二、GPS 拟合高程测量宜与GPS平面控制测量一起进行。
三、GPS 拟合高程测量的主要技术要求，应符合下列规定：1 GPS 网应与四等或四等以上的水准点联测。
联测的GPS 点，宜分布在测区的四周和中央。
若测区为带状地形，则联测的GP S 点应分布于测区两端及中部。
2 联测点数，宜大于选用计算模型中未知参数个数的1．5 倍，点间距宜小于10km。
3 地形高差变化较大的地区，应适当增加联测的点数。
4 地形趋势变化明显的大面积测区，宜采取分区拟合的方法。
5 GPS 观测的技术要求，应按本规范3．2 节的有关规定执行；其天线高应在观测前后各量测一次，取其平均值作为最终高度。
四、GPS 拟合高程计算，应符合下列规定：1 充分利用当地的重力大地水准面模型或资料。
2 应对联测的已知高程点进行可靠性检验，并剔除不合格点。
3 对于地形平坦的小测区，可采用平面拟合模型；对于地形起伏较大的大面积测区，宜采用曲面拟合模型。
4 对拟合高程模型应进行优化。
5 GPS 点的高程计算，不宜超出拟合高程模型所覆盖的范围。
五、对GPS 点的拟合高程成果，应进行检验。
检测点数不少于全部高程点的10％且不少于3个点；高差检验，可采用相应等级的水准测量方法或电磁波测距三角高程测量方法进行，其高差较差不应大于30D mm(D 为检查路线的长度，单位为km)。
1)导线的布设导线的布设有闭合导线、附合导线及支导线三种基本形式，如图所示。
3．支导线从一个高级点C和CD边的已知方位角出发，延伸出去的导线C、9、10、11称为支导线。
由于支导线只具有必要的起始数据，缺少对观测数据的检核，因此，只限于在图根导线和地下工程导线中使用。
GPS拟合高程精度分析
随着全球定位系统（ＧＰＳ）技术的成熟和仪器价格的大幅下降，该技术在测绘方面已普及。
众所周知，全球定位系统（ＧＰＳ）测定平面位置能达到很高的精度，但其测定的高程精度，由于受到高程异常的制约，尚有许多问题需要研究，以便充分发挥全球定位系统（ＧＰＳ）技术的优越性。
ＧＰＳ测定的是从地面到ＷＧＳ－８４椭球的大地高Ｈ，而我国采用的是正常高ｈ，两者存在一个差值，即高程异常，而高程异常与勘探区的地形、地层结构等因素有关。
本文依据六个煤田勘探区的实测资料，对拟合的ＧＰＳ高程在平地、丘陵、山地所达到的精度进行分析，并提出应用范围和可能提高精度的浅见，供研究探讨。
１ＧＰＳ高程异常的求解１．１勘探区ＧＰＳ网概况勘探区ＧＰＳ网施测于０４－０６年，采用Ｔｒｉｍｂｌｅ５７００双频ＧＰＳ接收机，依据０１《全球定位系统（ＧＰＳ）测量规范》，按Ｄ、Ｅ级布网，各区概况如表１。
外业解算出的基线，经过同步环和异步环及复测基线检验，其误差，在限差１／３以内达８０％，其余在１／２以内，证明精度是可靠的。
１．２ＧＰＳ高程异常及变化量求解ＧＰＳ高程异常，大地高Ｈ，正常高ｈ，它们之间的关系为ξ＝Ｈ－ｈ（１）Ｈ由ＧＰＳ测量求得，网中每个点均可获得，ｈ由水准或者三角高程测得。
依据公式（１）求出高程异常值。
Δξ＝ξ－ξｉ（２）其中：ξｉ参考点高程异常。
依据（２）式，相对中部一点（参考点）计算的高程异常变化量见表２。
经过数据统计分析，平地第四纪覆盖区，似大地水准面比较平缓，其变化幅度每公里４ｃｍ左右；丘陵区特别是梯田底部与丘陵的上部，似大地水准面有曲率显示，其变化幅度每公里７ｃｍ左右；在山区，可能受地貌相关性影响，似大地水准面变化较大，其变化幅度每公里１０ｃｍ左右。
从六个勘探区高程异常变化趋势看，东西方向大，南北方向小。
测区号参考点高程异常变化量／ｃｍ２３４５６７８９１０１１１２１３１４Ｋ１１１９．１１５．３１２．６８．１４．０１８．７１６．５１３．８１１．６９．３６．８５．４２．９Ｋ２１１５．７１３．２１０．４７．８５．０１８．７１．７１４．２１３．０１０．７８．４４．９２．３Ｋ３１１９．９１５．０１０．６７．２２．０２．６１２．１８．９３．９１．２Ｋ４１２０．４１７．０１１．３７．３２．１１６．５１６．４１２．５５．２２．４Ｋ５１１９．８１２．０５．８０．７１９．１１９．２５．２０．５Ｋ６１３２．３２５．３１８．６１．８３０．１３．４１３．１１．４表２高程异常变量统计表Ｔａｂｌｅ２Ｈｅｉｇｈｔａｎｏｍａｌｙｖａｒｉａｂｌｅｓｓｔａｔｉｓｔｉｃｓ勘探区编号地形类别面积／ｋｍ２相对高差／ｍＧＰＳ点数联测水准点数三等四等Ｋ１平地１３００１０６２１７０Ｋ２平地９５０２５５１４１０Ｋ３丘陵７００２９０４５４７Ｋ４丘陵６００２８０４４５６Ｋ５山地８５０６８０４９３５Ｋ６山地８００１０３０５４４６表１各区概况统计表Ｔａｂｌｅ１Ｓｔａｔｉｓｔｉｃｓｏｆｅａｃｈａｒｅａ＇ｓｇｅｎｅｒａｌｓｉｔｕａｔｉｏｎ作者简介：雒养社（１９６３—），男，陕西户县人，高级工程师，１９８７年毕业于西安矿业学院，现从事煤田地质测绘工作。
GPS高程拟合及其在公路勘察中的应用研究的开题报告
GPS高程拟合及其在公路勘察中的应用研究的开题报告一、课题背景随着现代化道路交通的快速发展，公路建设以及勘察工作的效率和准确性要求越来越高。
GPS(Global Positioning System，全球卫星定位系统)技术的应用为公路勘察带来了新的思路和方法。
GPS高程拟合是指通过GPS测量数据来拟合大地高程模型的方法，可以在较高精度范围内获取大地高程并反映公路线路的高程情况。
因此，探索GPS高程拟合技术在公路勘察中的应用具有重要意义。
二、研究目的本研究旨在探究GPS高程拟合技术在公路勘察中的应用，研究目的主要包括以下三方面：1. 研究GPS高程拟合原理及其适用范围。
通过对GPS高程拟合原理的深入探究，确定其适用范围和技术要求。
2. 建立高程拟合模型，并对其进行验证。
通过采集一组GPS测量数据，建立高程拟合模型，进行实验验证。
3. 探索GPS高程拟合技术在公路勘察中的应用。
将GPS高程拟合技术应用于公路勘察，比较其效果和传统勘察方法的差异，探索其应用前景。
三、研究内容1. GPS高程拟合原理的研究。
归纳总结GPS高程拟合的基本原理，结合GPS定位技术原理，建立高程拟合模型。
2. GPS测量数据采集及处理。
利用GPS接收机采集一组GPS测量数据，进行数据处理和分析。
3. 高程拟合模型建立及验证。
将采集得到的GPS测量数据应用于高程拟合模型，验证该模型的有效性。
4. GPS高程拟合技术在公路勘察中的应用。
将GPS高程拟合技术应用于公路勘察中实际案例，比较其效果与传统勘察方法的差异，并探索其优势和应用前景。
四、研究方法1. 文献调研。
通过查阅相关文献资料，对GPS高程拟合技术及其在公路勘察中的应用进行综合了解。
2. 采集GPS测量数据。
使用GPS接收机采集一组GPS测量数据，并进行数据预处理和分析，确定处理方案。
3. 建立GPS高程拟合模型，并验证模型的可靠性和精度。
4. GPS高程拟合技术在公路勘察中的应用。
GPS水准高程拟合报告-推荐下载
实验分析：
1 本实验中可以选择两种差值公式算法
Δh(dB)散点图
由以上散点图可以看出以下特点： A 采用 dB dL 作为变量值求得的插值函数所插值的结果与真是值差距普遍集中
在-0.4 到-0.6 之间，显然这种拟合结果的精度并不是很高，另外我们不得不 怀疑这种拟合结果存在某种系统误差因为其散点图具有一定的偏向性，同时 作为公共点已知数据的十个点的差值结果亦如想象的那样理想，这证明我们 的差值函数参数的解算没用问题。 B 采用 B L 作为变量值求得的插值函数所插值的结果与真是值差距普遍集中在0.4 到 0.4 之间，显然这种拟合结果的精度也不是很高，但较第一种差值有 了一定的提高，这种插值结果的散点图分布较为随机不太可能存在系统误差， 同时作为公共点已知数据的十个点的差值结果亦如想象的那样理想，这证明 我们的差值函数参数的解算没用问题。 总结以上拟合结果尽管其精度不尽相同但是两种方法的精度显然不能与传统 的水准测量相媲美，完全不能满足实际工程水准的需要。 2 下面采用拉格朗日 n 次插值对原数据在进行处理
公式一
式中 dB=B-B0;dL=L-L0;B0=1/n∑B;L0=1/n∑L,n 为 GPS 观测点的数量。
利用其中一些具有水准观测资料的公共点上的的大地高和正高可以计算出这些 点的大地水准面差距。利用这些公共点的观测资料求得公式一的参数，再利用 求得的公式进行其他点的大地水准面差距内插，和正高的拟合；
GPS 水准高程拟合报告
GPS高程拟合方法及精度分析
GPS高程拟合方法及精度分析引言随着全球定位系统(GPS)的普及和发展，GPS技术在地球科学、工程测量和导航定位等领域得到了广泛的应用。
GPS高程的测量和拟合在地球科学研究和工程测量中扮演着重要的角色。
对GPS高程拟合方法及其精度进行深入的研究和分析具有重要的意义。
一、GPS高程拟合方法GPS高程的测量是通过GPS卫星信号和接收机接收时间的差值来计算得到的。
在GPS测量中，精确的高程测量是非常重要的。
高程拟合是指根据已知的GPS观测数据，通过一定的数学模型和算法，来拟合出地球表面上各点的高程值。
目前常用的GPS高程拟合方法主要包括差分GPS法、动态大地水准面模型法和GNSS/地球重力模型法。
1. 差分GPS法差分GPS法是基于参考站和移动站测量GPS信号的相位和码距的差值来进行高程测量的方法。
该方法可以减小大气层等误差对高程测量的影响，提高高程测量的精度。
差分GPS法广泛应用于工程测量和导航领域，具有较高的精度和实用性。
2. 动态大地水准面模型法动态大地水准面模型法是基于大地水准面模型预测的高程值和GPS观测数据进行拟合的方法。
通过使用大地水准面模型，可以对GPS测量中的大气层延迟和其他误差进行校正，提高高程测量的精度。
该方法适用于地球科学研究领域，可以得到更为精确的高程值。
二、GPS高程拟合精度分析GPS高程拟合的精度是衡量其可靠性和实用性的重要指标。
在GPS高程拟合过程中，需要对其精度进行综合分析和评估。
1. 精度影响因素GPS高程拟合的精度受到多种因素的影响，主要包括大气层延迟、接收机误差、地形和重力效应、卫星轨道误差等。
这些因素会对GPS高程拟合的精度产生影响，需要在实际应用中进行综合考虑和分析。
2. 精度评估方法针对GPS高程拟合的精度进行评估，可以采用单点定位和差分定位、统计分析和误差分析等方法。
通过对GPS观测数据和拟合结果进行综合分析和评估，可以得到GPS高程拟合的精度水平和可靠性。
GPS高程拟合方法及精度分析
GPS高程拟合方法及精度分析
GPS全球卫星定位系统（Global Positioning System）是一种全球性的导航系统，它可以利用卫星进行高精度的位置定位。
然而，GPS定位的高程精度受到多种因素的影响，
包括GPS接收机本身、信号传输路径等，因此需要对GPS高程进行拟合处理以提高其精
度。
GPS高程拟合方法主要包括差值法、插值法和回归分析法三种。
差值法是根据GPS测量到的位置信息和地面标高测量值之差，通过差值运算来得到GPS高程测量值。
差值法具有计算简单、速度快的特点，但局限性较大，不能解决在GPS
定位时所遇到的某些问题，例如多径效应等。
插值法需要用周围已知高程数据进行插值计算，以得出该位置的高程。
插值法的精度
与区域内高程数据的分布稠密程度有关，一般来说，在数据较为密集的情况下，插值法的
精度较高，反之则不佳。
回归分析法将GPS测量到的位置信息与实测标高之间的相关性进行线性拟合，由此推
导出每个位置的GPS高程测量值。
回归分析法的精度受到模型的影响，模型的构建需要考
虑影响因素的相互作用和相关度。
实际应用中，GPS高程拟合方法的选择需要结合实际情况进行决策。
在拟合方法上，
一般建议采用回归分析法，因为它可以分析其他影响因素，并将其纳入模型中，从而提高
精度。
在应用上，需要结合当地的天气、地形和信号传输情况等因素进行多次测量和比对，以提高GPS高程的精度。
总体而言，在选择GPS高程拟合方法时，应考虑实际需求和精度要求，从而选择适合
自己的方法。
此外，对GPS高程的整体监测和维护也是提高其精度的重要措施。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
GPS水准高程拟合报告实验目的：1掌握GPS水准高程拟合的基本原理，了解高精度GPS水准的研究意义；2能够利用Matlab编程实现几何内插法拟合GPS水准高程；实验内容：利用Matlab编程实现几何内插法拟合GPS水准高程，并作内插结果分析实验原理：1大地水准面，参考椭球面，正高，大地高之间的几何关系A 正高的定义是：由地面点沿通过该点的铅垂线到大地水准面的距离。
一般用符号Hg表示。
B 大地高的定义是：由地面点沿通过该点的椭球面法线到椭球面的距离。
也称为椭球高，一般用符号H表示。
大地高是一个纯几何量，不具有物理意义。
同一个点，在不同的基准下，具有不同的大地高。
利用GPS，可以测定地面点的WGS-84中的大地高。
C 大地水准面差距：大地水准面到椭球面的距离，称为大地水准面差距，记为hg （或N）。
如上图可以看出大地高和正高之间的关系：H=Hg+ hg2几何内插法原理几何内插法是通过一些既进行了GPS观测又具有水准资料的点上的大地水准面差距，采用平面或曲面拟合，配置三次样条等内插方法，得到其他点上的大地水准面差距从而反算这些点上的正高。
3二多项式拟合N=a0+a1*dB+a2*dL+a3*dB2+a4*dL2+a5*dL*dB 公式一式中dB=B-B0;dL=L-L0;B0=1/n∑B;L0=1/n∑L,n为GPS观测点的数量。
利用其中一些具有水准观测资料的公共点上的的大地高和正高可以计算出这些点的大地水准面差距。
利用这些公共点的观测资料求得公式一的参数，再利用求得的公式进行其他点的大地水准面差距内插，和正高的拟合；实验步骤：1输入已知点的GPS观测值和相应的正常高构成矩阵B,L,H,h,分别是纬度矩阵，经度矩阵，大地高矩阵，正高矩阵；2计算dB=B-B0;dL=L-L0;B0=1/n∑B;L0=1/n∑L，构成矩阵矩阵dB,dL和大地水准面差距矩阵N=H-h;3将以上计算得到的矩阵代入公式一经过间接平差求得相应的参数a i,这样就能构成一个确定的多项式二；4输入未知待求点的GPS观测值构成矩阵BB,LL,HH，计算相应的dBB,dLL；5将dBB,dLL矩阵代入多项式二，解算出对应点的大地水准面差距NN矩阵；6反算各点的正高h=H-NN;7对计算得到的正高，大地水准面差距做对比分析；实验分析：1本实验中可以选择两种差值公式算法（1）N=a0+a1*dB+a2*dL+a3*dB2+a4*dL2+a5*dL*dB（2）N=a0+a1*B+a2*L+a3*B2+a4*L2+a5*L*B采用公式（1）的插值结果如下：Δh(dB)散点图注：Δh(dB)是插值点的水准资料与插值结果的差值采用公式（2）的插值结果如下：Δh(dB)散点图由以上散点图可以看出以下特点：A 采用dB dL作为变量值求得的插值函数所插值的结果与真是值差距普遍集中在-0.4到-0.6之间，显然这种拟合结果的精度并不是很高，另外我们不得不怀疑这种拟合结果存在某种系统误差因为其散点图具有一定的偏向性，同时作为公共点已知数据的十个点的差值结果亦如想象的那样理想，这证明我们的差值函数参数的解算没用问题。
B采用 B L作为变量值求得的插值函数所插值的结果与真是值差距普遍集中在-0.4到0.4之间，显然这种拟合结果的精度也不是很高，但较第一种差值有了一定的提高，这种插值结果的散点图分布较为随机不太可能存在系统误差，同时作为公共点已知数据的十个点的差值结果亦如想象的那样理想，这证明我们的差值函数参数的解算没用问题。
总结以上拟合结果尽管其精度不尽相同但是两种方法的精度显然不能与传统的水准测量相媲美，完全不能满足实际工程水准的需要。
2 下面采用拉格朗日n次插值对原数据在进行处理拉格朗日三次差值残差分布拉格朗日四次差值残差分布拉格朗日五次差值残差分布拉格朗日六次差值残差分布注：1由于五次和六次插值的效果较差，五次插值仅选取了拟合效果相对较好的前39项作图，六次插值仅选取前25项；2拟合多项式参数求解并未进行平差解算；分析拉格朗日插值的几幅残差图发现在四次插值时其残差分布与Δh(dB)散点图比较接近，其他的几种高次插值效果都比较差，而且随着插值次数的提高，拟合效果也来越差其残差成指数增加，由此可见多项式插值对GPS水准拟合并不是一种好的拟合方式；对此我们可以尝试用空间插值，例如克里金插值法。
实验总结：与常规的水准测量而言GPS水准具有效率高费用低等特点，可以在大范围内进行高程测量，对于大型工程而言传统的大地水准将耗费大量的人力物力财力，且观测周期非常长，如果GPS可行那么将为大范围水准测量开辟崭新的空间时代，同时高精度的GPS水准也将改变传统大地测量在高程基准面确定上的难题，为测绘作业提供高精度的大地水准面。
但是就目前而言GPS水准高程的精度问题是其进行工程应用的难题，其原因主要是大地水准面差距的精度低，目前我们还难以获得较高精度的全球大地水准面，解决GPS水准的主要方法主要是几何插值法，就本次实验可以看出多项式插值的精度并不能满足实际需要，但是综合目前各方面的实际情况，GPS水准可以满足三等水准的要求。
对于一些特殊的工程而言要获取高精度的GPS水准采用的处理方式是，在局部区域内建立若干GPS长期观测基站利用这些基站的观测值和IGS精密星历解算该区域的大地水准面，以获得精度较高的大地水准面差距从而能获得精度较高的高程，这方面的例子主要是港珠澳大桥项目，不过由于成本太高耗时太长一般的项目难以企及，不具有工程实用性。
附Matlab程序代码%原始数据B=[4160726,4154896,4159270,4148735,4157177,4153413,4157949,4151556,41 55573,4157027]'L=[469370,473772,473622,474714,474356,476134,476233,477565,476918,477814]'H=[43.085,47.677,37.002,56.859,43.904,78.045, 41.697,51.519,33.57, 32.715]' h=[15.395,19.758,9.137,28.900,15.993,50.071,13.721,23.492,5.574,4.685]'%定义dB=B-B0dB=zeros(10,1);dL=zeros(10,1);B0=0;L0=0;for i=1:10B0=B0+B(i,1);L0=L0+L(i,1);endB0=B0/10.0;L0=L0/10.0;for i=1:10dB(i,1)=B(i,1)-B0;dL(i,1)=L(i,1)-L0;end%V=XA-N;求相应数据X=zeros(10,6);%初始化X，赋值for i=1:10X(i,1)=1.0;X(i,2)=dB(i,1);X(i,3)=dL(i,1);X(i,4)=dB(i,1)*dB(i,1);X(i,5)=dL(i,1)*dL(i,1);X(i,6)=dL(i,1)*dB(i,1);end%求NN=zeros(10,1);for i=1:10N(i,1)=H(i,1)-h(i,1);end%间接平差求A；X'PV=0 V=XA-N P=E；A=zeros(6,1);s=X'*X;t=inv(s);u=X'*N;A=t*u;%读入要计算的原始数据BB=[4149569,4152803,4155546,4159369,4157612,4147678,4155194,4159612,4 151539,4157232,4150583,4158508,4160739,4156273,4154803,4157419,4163806,4 163116,4146497,4153030,4143589,4148280,4150278,4156103,4146715,4149128,4 154168,4155429,4144712,4148524,4150861,4154013,4146430,4145664,4148161,4143343,4145375,4144379,4142039, 4139017]'LL=[479940,479819,479795,479113,479861,482166,482154,482603,482291,483 485,485202,485355,485636,486046,487445,488532,488447,490930,486936,487566 ,488117,487969,488166,490412,490503,490823,491023,491935,492338,494523,49 4172,493828,495124,498357,497419,500731,500943,504035,504085,504544]' HH=[35.852,35.099,33.292,32.411,32.857,40.737,34.026,33.772,35.689,33.418, 65.076,39.009,32.066,105.537,128.872,32.045,37.03,34.798,49.149,84.576,44.983,58.323,86.432,36.14,59.754,88.67,100.893,76.64,49.964,39.018,52.246,71.3 9,33.804,32.367,32.36,64.272,58.95,31.911,37.269,37.935]'BB0=0;LL0=0;for i=1:40BB0=BB0+BB(i,1);LL0=LL0+LL(i,1);endBB0=BB0/40.0;LL0=LL0/40.0;dBB=zeros(40,1);dLL=zeros(40,1);XX=zeros(40,6);%计算各点大地水准面差距NNfor i=1:40dBB(i,1)=BB(i,1)-BB0;dLL(i,1)=LL(i,1)-LL0;XX(i,1)=1.0;XX(i,2)=dBB(i,1);XX(i,3)=dLL(i,1);XX(i,4)=dBB(i,1)*dBB(i,1);XX(i,5)=dLL(i,1)*dLL(i,1);XX(i,6)=dLL(i,1)*dBB(i,1);end%MM是10个已知点的大地水准面差距的插值矩阵，o是10个已知值的正高矩阵NN=zeros(40,1);NN=XX*A;MM=X*A;hh=HH-NN;o=H-MM;。