重力资料整理MATLAB实验.docVIP

PAGE 重力资料整理上机实验 班级_1250242___ 学号__125024202_ 姓名_侯继伟_ 1、基点网平差实验 下图为三环路基点网，各边段差值如图所示，顺时针方向为正。 假设各边段时间差均为1，试采用线性方程组法求解平差系数，要求用MATLAB编程求解并输出各环路的平差系数。 原理： MATLAB程序： clc; V1=[-2.03 -4.322 1.526 4.836]; V2=[-1.526 -1.486 -0.83 3.828]; V3=[0.83 -1.91 0.918]; dV1=sum(V1); dV2=sum(V2); dV3=sum(V3); a=[4 -1 0;-1 4 -1;0 -1 3] K=inv(a)*[dV1;dV2;dV3] 计算结果： a = 4 -1 0 -1 4 -1 0 -1 3 K = -0.0023 -0.0192 -0.0604 2、布格重力异常值计算 下表为某10个重力测点的地理坐标和相对基点的重力值，试用MATLAB编程求各测点的布格重力异常值。已知基点纬度为33°11′33″，坐标x=0，y=0，海拔高程为14.8米，地表岩石平均密度为2.2g/cm3。（不用作地形校正） 测点号 东向横坐标Y /米 北向纵坐标X /米