地球重力场计算EGM1996

地球重力场，计算地球各个地方重力，适用于大地测量专业

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;
using System.IO;

namespace EGM1996
{
    public partial class Form1 : Form
    {
        public Form1（）
        {
            InitializeComponent（）;
        }
        double[] C = new double[361 361];
        double[] S = new double[361 361];//用于存放球谐系数
        double[] P = new double[361 361];//用于存放球谐系数
        double[] dP = new double[361 361];//用于存放球谐系数

        double[] L2 = new double[20 2];//用于存放经纬度
        string[] L3 = new string[20];//用于存放经纬度


        double a = 6378137.0000;//椭球长半轴
        double e2 = 0.0066943799013;//椭圆第一偏心率
        double GM = 3.986004415E14;//地球引力常数
        double H = 0;




        int h = 0;//确定阶数

        public double DMS2Rad（double x）
        {
            int z = 0;
            if （x > 0）
            {
                z = 1;
                x = x * z;
            }
            else
            {
                z = -1;
                x = x * z;
            }

            x = x * 10000;

            double m = x % 100;
            x = （x - m） / 100;

            double f = x % 100;
            x = （x - f） / 100;

            double d = x;
            x = （d + f / 60 + m / 3600） * Math.PI / 180;

            return x * z;
        }

        //大地转地心
        public double DDDX（double B double R double N double H double e2）
        {
            double Y = （N * （1 - e2） + H） * Math.Sin（B）;
            double phi = Math.Asin（Y / R）;
            return phi;

        }

        //求勒让德函数
        public void LRD（double[] P double[] dP double phi int N）
        {
            P[0 0] = 1.0; P[1 0] = Math.Sqrt（3） * Math.Sin（phi）; P[1 1] = Math.Sqrt（3） * Math.Cos（phi）;

            dP[0 0] = 0.0; dP[1 0] = Math.Sqrt（3） * Math.Cos（phi）; dP[1 1] = -Math.Sqrt（3） * Math.Sin（phi）;

            for （int i = 2; i < N; i++）
            {
                for （int j = 0; j <= i; j++）
                {
                    if （j == i）
                    {
                        double n1 = 0.0 n2 = 0.0;
                        n1 = （2 * i + 1）; n2 = 2 * i;
                        P[i j] = Math.Sqrt（（n1 / n2）） * Math.Cos（phi） * P[i - 1 j - 1];
                        dP[i j] = Math.Sqrt（（n1 / n2）） * （Math.Cos（phi） * dP[i - 1 j - 1] - Math.Sin（phi） * P[i - 1 j - 1]）;

                    }
                    else if （j == i - 1）
                    {
                        double n3 = 0.0 n4 = 0.0;
                        n3 = 4 * i * i - 1; n4 = i * i - j * j;
                        P[i j] = Math.Sqrt（（n3 / n4）） * Math.Sin（phi） * P[i - 1 j];
                        dP[i j] = Math.Sqrt（（n3 / n4）） * （Math.Sin（phi） * dP[i - 1 j] + Math.Cos（phi） * P[i - 1 j]）;

                    }
                    else if （j != i 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5612  2012-10-21 12:23  EGM1996\EGM1996\Properties\Resources.resx

     文件        249  2012-10-21 12:23  EGM1996\EGM1996\Properties\Settings.settings

     文件       1364  2012-10-21 12:23  EGM1996\EGM1996\Properties\AssemblyInfo.cs

     文件       2864  2012-10-21 12:23  EGM1996\EGM1996\Properties\Resources.Designer.cs

     文件       1092  2012-10-21 12:23  EGM1996\EGM1996\Properties\Settings.Designer.cs

     文件       3760  2012-10-21 12:32  EGM1996\EGM1996\EGM1996.csproj

     文件        490  2007-07-21 01:33  EGM1996\EGM1996\bin\Debug\EGM1996.vshost.exe.manifest

     文件      14328  2013-03-16 21:40  EGM1996\EGM1996\bin\Debug\EGM1996.vshost.exe

     文件      18432  2012-11-25 18:32  EGM1996\EGM1996\bin\Debug\EGM1996.exe

     文件      34304  2012-11-25 18:32  EGM1996\EGM1996\bin\Debug\EGM1996.pdb

     文件        180  2012-10-22 19:42  EGM1996\EGM1996\obj\Debug\EGM1996.Form1.resources

     文件        180  2012-10-21 12:32  EGM1996\EGM1996\obj\Debug\EGM1996.Properties.Resources.resources

     文件      18432  2012-11-25 18:32  EGM1996\EGM1996\obj\Debug\EGM1996.exe

     文件        564  2013-03-16 21:40  EGM1996\EGM1996\obj\Debug\EGM1996.csproj.FileListAbsolute.txt

     文件        847  2012-10-22 19:42  EGM1996\EGM1996\obj\Debug\EGM1996.csproj.GenerateResource.Cache

     文件      34304  2012-11-25 18:32  EGM1996\EGM1996\obj\Debug\EGM1996.pdb

     文件        488  2012-10-21 12:23  EGM1996\EGM1996\Program.cs

     文件       6012  2012-10-22 19:42  EGM1996\EGM1996\Form1.resx

     文件      10915  2012-10-22 19:43  EGM1996\EGM1996\Form1.Designer.cs

     文件      14825  2012-10-23 13:59  EGM1996\EGM1996\Form1.cs

     文件        911  2012-10-21 12:23  EGM1996\EGM1996.sln

    ..A..H.     13824  2013-03-16 21:40  EGM1996\EGM1996.suo

     文件    5619962  2012-08-30 10:57  EGM1996\egm96.gfc

     文件        542  2012-10-14 22:01  EGM1996\jwd.txt

     文件       1045  2012-10-23 13:59  EGM1996\1.txt

     目录          0  2012-10-21 12:23  EGM1996\EGM1996\obj\Debug\TempPE

     目录          0  2012-10-21 13:15  EGM1996\EGM1996\obj\Debug\Refactor

     目录          0  2012-10-21 12:23  EGM1996\EGM1996\bin\Debug

     目录          0  2012-10-21 12:23  EGM1996\EGM1996\obj\Debug

     目录          0  2012-10-21 12:23  EGM1996\EGM1996\Properties

............此处省略7个文件信息