固体潮 matlab,基于MATLAB改进设计的固体潮改正程序及其应用

科技的不断进步和发展,新的方法、软件都通算,得到每一个重力测点的固体潮改正值,并将其与过计算机不断的实现,MATLAB　以独特的优势广泛一定区域使用同一个固体潮改正经纬度计算所得到的运用于各个领域,尤其是数据的计算方面具有计算速每一点的固体潮改正值进行比较,结果采用每一点的度快,编程简单,受到广大计算机爱好者的青睐。重经纬度计算该点的固体潮改正值比一定区域使同一个力勘探作为一种物理勘探方法,越来越广泛的运用,潮改经纬度计算得到每一点的固体潮改正值更准确,而且测量和计算的精度要求也越来越高,但是由于固提高了固体潮改正精度。体潮改正值较小,研究人员较少,目前固体潮改正大1　固体潮改正的原理多数在EXCEL中实现　,郜晓亮等用MATLAB实现由于地球、太阳、月亮的相对位置不断变化会了固体潮的改正计算,董良等将计算公式中天文参数引起重力值的变化,把这两者作用在地球上某一点的标准历元为J1900.0改为最新公布的J2000.0系统,和称为重力固体潮理论值。在重力勘探中,在计算测用MATLAB实现了固体潮的改正,改进前后固体潮量点的绝对重力值时,先要进行固体潮改正,消除固之差为±0.001×10-8m/s2。虽然有所改进,但差值很体潮的影响。固体潮改正值主要计算公式为小,这很小的差值可能还是由于数据位数的取舍所造　　　　　　　　　　　(1)成,重力勘探目前所使用的计算固体潮的方式仍是把工区中心点(或将工区分成几块,每一块选择其中心　　　　　　(2)点)的经纬度作为固体潮改正的经纬度,这样计算较为简单,计算量也大大的减少,但是降低了固体潮改正的精度。　　　　　　　　(　3)笔者将每一个重力测点实测得到的经纬度坐标,　　　　　　　　　　(　　4)作为固体潮改正的经纬度,用MATLAB实现编程计

式中:为固体潮改正值;为潮汐因子,取1.16;为地心至月心的平均距离;为月心至地心的距离;为地心至日心的平均距离;为日心至地心的距离;为月亮对测点的地心天顶距;为太阳对测点的地心天顶距;为测点纬度;为测点地心纬度。具体计算参照参考文献[10]。2　程序设计图1　计算数据格式示意图程序设计主要按照上文中固体潮改正原理中公式进行编写,在计算中主要有数据的读取、计算时刻替换为“0”,然后再点乘即得到了当年月计算的天数,的儒略世纪数T的计算、计算结果的输出等三个模块再加上计算之日以前本月的天数,就是所求的总天数。较为复杂,每个模块的设计及使用作为重点介绍,其2.3　计算结果的输出余均按公式进行编写,在计算程序中体现。计算得到的固体潮改正值,需要输出到文件,为2.1　数据的读取了方便,就将计算结果输出到数据输入的原文档的另一数据的显示格式采用长格式,以便对比分析:个表格中,名称为“固体潮改正值”,采用的计算函format　long;由于重力勘探数据量较大,用Microsoft　数为xlswrite('GUTICHAO.xls',gtc,'固体潮改正值','A1')。Excel　工作表整理较为方便,故该程序以Microsoft　3　计算结果的对比分析Excel　工作表为计算数据的存储文件,在计算时,读为验证程序计算的正确性,笔者选用了《重力取Microsoft　Excel　工作表的计算数据,再进行固体潮调查技术规范(1∶5万)》[10]中固体潮改正的示例:改正计算。文件名为GUTICHAO.xls;数据格式如图经度为93°00′,纬度为31°20′,日期为20031所示,A-G列依次为点号、日期、时间、重力值、年5月6日,时间为19∶45,计算的结果与示例结近区地改值、固体潮改正经