测量程序设计实验报告实验名称：大地坐标与空间直角坐标的换算实验四 大地坐标与空间直角坐标的换算一、实验目的编写大地坐标与空间直角坐标相互转换的程序，并对格式化文件数据进行计算,验证程序
二、实验内容：1、大地坐标向空间直角坐标换算转换公式： （1）其中：L 为经度，B 为纬度，h 为大地高,为卯酉圈曲率半径，为第一偏心率，a 为旋转椭球长半轴，b 为短半轴
WGS84 椭球参数：长半轴 a = 6378137 扁率 f = 1/298
257223563根据上式创建以 geo2xyz 命名的函数,函数输入输出格式为 [x, y， z］ = geo2xyz （L， B, h)2、空间直角坐标向大地坐标换算根据式（1）推导大地坐标向空间直角坐标转换公式： 注意计算纬度时需要用到迭代，可用作为初始值
创建以 xyz2geo 命名的函数，函数输入输出格式为 ［L， B， h] = xyz2geo (x, y, z) 三、实验步骤1、 大地坐标向空间直角坐标换算主程序：%％大地坐标向空间直角坐标换算%函数的输入输出格式为［x，y，z］=geo2xyz(L，B，h)[filename，pathname］ = uigetfile(’*
txt’,'请选择打开的数据文件’)；file = ［pathname， filename];data = importdata(file);L=data
data（:，1);B=data
data（：,2）;h=data
data（：，3）;［x,y，z]=geo2xyz(L，B，h）；A=［x,y，z]；A=A’;[filename_out,pathname_out］ = uiputfile('＊
txt'，’请选择要输出数据文件'）；fileout = [pathname_out， filename_out］;
