基于Matlab模拟地球重力场

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法  神经网络预测 雷达通信  无线传感器

信号处理 图像处理 路径规划 元胞自动机 无人机  电力系统

⛄ 内容介绍

针对地形起伏较大地区纯粹采用数值逼近方法进行GPS高程转换精度不高的事实,提出了一种在"移去-拟合-恢复"思路的基础上,综合利用EGM2008地球重力场模型和数值逼近算法进行GPS高程转换的方法.通过贵州某段高速公路首级测量控制网实例验证了该方法的有效性和实用性,证实该方法在地形起伏较大的重山区能够达到较高的精度.

⛄ 完整代码

load grav67.doc

x=1;

while x<=1173

f(x,1)=grav67(x,2);

l(x,1)=grav67(x,3);

DgF(x,1)=grav67(x,5);

f1(x,1)=f(x,1)*pi/180;

x=x+1;

end

for x=1:1173

    go67(x,1)=(9.78031846*(1+(0.0053024*sin(f1(x,1))^2)-(0.0000059*sin(2*f1(x,1))^2)))*(10^5);

    go80(x,1)=(9.78032677*((1+(0.001931851353*sin(f1(x,1))^2))/(1-(0.0066943800229*sin(f1(x,1))^2))^(1/2)))*(10^5);

    Dg80(x,1)=DgF(x,1)+go67(x,1)-go80(x,1);

end

[x,y]=meshgrid(-120:(1/12):-114,56:(1/12):60);

z=griddata(f,l,Dg80,y,x,'v4');

load grav80.txt

for i=1:1173

    EIGEN(i,1)=Dg80(i,1)-grav80(i,5);

    

end

o=griddata(f,l,EIGEN,y,x,'v4');

contourf(x,y,z),colorbar,title('GRS80(1)')

figure

contourf(x,y,o),colorbar,title('EIGEN(1)')

load grav80rtm.txt

for i=1:1173

    MRT(i,1)=EIGEN(i,1)-grav80rtm(i,5);

end

figure

lo=griddata(f,l,MRT,y,x,'v4');

contourf(x,y,lo),colorbar

⛄ 运行结果

⛄ 参考文献

[1]张兴福, 沈云中. 应用卫卫跟踪数据恢复地球重力场的模拟分析[J]. 同济大学学报：自然科学版, 2005, 33(10):4

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料