用MATLAB进行大地测量学上的子午线弧长计算

clear;
clc;
a=6378140;
b=6356755.288157528;
e=sqrt(a^2-b^2)/a;
% 确定系数关系
m0=a*(1-e^2);
m2=1.5*e^2*m0;
m4=5.0/4.0*e^2*m2;
m6=7.0/6.0*e^2*m4;
m8=9.0/8.0*e^2*m6;
a0=m0+m2/2+3/8*m4+5/16*m6+35/128*m8;
a2=m2/2+m4/2+15/32*m6+7/16*m8;
a4=m4/8+3/16*m6+7/16*m8;
a6=m6/32+m8/16;
a8=m8/128;
%求赤道到某一纬度的弧长
lat=60; % 北极
lat=lat*pi/180;  % 化为弧度
X=a0*lat-a2/2*sin(2*lat)+a4/4*sin(4*lat)-a6/6*sin(6*lat)+a8/8*sin(8*lat);

该代码是根据纬度的度数来求得赤道到某一纬度的弧长。其中输入的更改可以对lat=60进行更改;

下面是lat=60时的运行结果：

用MATLAB进行大地测量学上的子午线弧长计算相关推荐

赤道子午线弧长反演大地纬度
matlab赤道子午线弧长反演大地纬度话不多说,直接看公式太懒了,直接拍课本的照片 clear clc S=input('请输入由赤道起算的子午线弧长S:'); a=input('请输入长半轴半径 ...
子午线弧长与纬度相互间的转换
本程序是根据武汉大学出版社出版的<大地测量学基础>里的相关知识,对子午线弧长和纬度间相互转换的程序编写.程序内有详细的代码注释.程序是利用python编写,用pycharm编译器编译! i ...
Matlab 如何回退上一步、批量注释以及取消注释
写在这里的初衷,一是备忘,二是希望得到高人指点,三是希望能遇到志同道合的朋友. (今天看到有位博主写的话,感觉很认同) 学习matlab的时候,很多时候想回退到上一步,又找不到对应选项. 想批量注释和 ...
matlab生成exe-在没有安装matlab的电脑上运行matlab生成的exe
在没有安装matlab的电脑上运行matlab生成的exe 本文介绍如何将matlab生成exe,并将该exe在没有安装matlab的电脑上运行,与以前的博文相似地,不说废话,直接干货. 1.编写函数 ...
matlab在常微分方程的应用,MATLAB在常微分方程上应用.ppt
您所在位置:网站首页 > 海量文档 &nbsp>&nbsp计算机&nbsp>&nbspmatlab MATLAB在常微分方程上应用.ppt11页本文 ...
matlab和气象,matlab在气象上的应用
matlab在气象上的应用 ,,在大气科学中的应用,李若楠,一.Matlab界面简介以及一些重要指令和注意事项二.数据处理,,plot相关命令绘图,contour(f)命令绘图,不同投影方式绘图,其 ...
matlab开方分布上分位点,概率密度分布函数和上分位点的数值计算
概率密度分布函数和上分位点的数值计算 1 大连民族学院数学实验报告课程: 数理统计实验题目: 概率密度.分布函数和上分位点的数值计算系别: 理学院专业: 信息与计算科学姓名: 历 ...
Matlab 在图片上画虚线矩形框
matlab在图片上rectangle详解 Rectangle 属性矩形的外观和行为矩形属性控制矩形对象的外观和行为.通过更改属性值,您可以修改矩形的特定方面.使用圆点表示法查询和设置属性. h ...
MATLAB上的GPU加速计算
概述怎样在MATLAB上做GPU计算呢?主要分为三个步骤:数据的初始化.对GPU数据进行操作.把GPU上的数据回传给CPU 一.数据的初始化首先要进行数据的初始化.有两种方法可以进行初始化:一是先 ...

用MATLAB进行大地测量学上的子午线弧长计算

用MATLAB进行大地测量学上的子午线弧长计算相关推荐

最新文章

热门文章