坐标换算和GPS时间计算

GPS课第一次作业，懒得编2个EXE，全弄一起了，分享下。。。

界面如下：

依照惯例分享下代码。。。BS下无穷无尽的伸手党。。。

/// <summary>
        /// 计算大地经度L
        /// </summary>
        /// <param name="x"></param>
        /// <param name="y"></param>
        /// <returns></returns>
        public double cal_l(double x, double y)
        {
            return Math.Atan(y / x);
        }

    /// &lt;summary&gt;
    /// 计算原始的B
    /// &lt;/summary&gt;
    /// &lt;param name="b"&gt;&lt;/param&gt;
    /// &lt;returns&gt;&lt;/returns&gt;
    public double cal_orginb(double b)
    {
        return a / Math.Sqrt(1 - e * square(Math.Sin(b)));
    }

    /// &lt;summary&gt;
    /// 计算新的B
    /// &lt;/summary&gt;
    /// &lt;param name="n"&gt;&lt;/param&gt;
    /// &lt;param name="h"&gt;&lt;/param&gt;
    /// &lt;returns&gt;&lt;/returns&gt;
    public double cal_nb(double n, double h)
    {
        double a = z * (n + h);
        double ba = Math.Sqrt(square(x) + square(y)) * (n * (1 - e) + h);
        return Math.Atan(a/ba);
    }

    //计算高程
    public double cal_h(double n, double b)
    {
        //double tmp = Math.Sin(b);
        return (z / Math.Sin(b)) - n * (1 - e);
    }

    /// &lt;summary&gt;
    /// 运算过程
    /// &lt;/summary&gt;
    /// &lt;param name="err"&gt;&lt;/param&gt;
    /// &lt;returns&gt;&lt;/returns&gt;
    public double[] ConvertProcess(double err)
    {
        try
        {

            double[] tmp = new double[3];
            double n = 0;
            double b = originb;
            double h = 0;
            double rb = b;
            tmp[0] = ConvertToDegree(cal_l(x, y));
            if (tmp[0]&lt;0)
            {
                tmp[0] += 180;
            }
            do 
            {
                b = rb;
                n = cal_orginb(b);
                h = cal_h(n, b);
                rb = cal_nb(n, h);

            } while (Math.Abs(rb - b) &gt; err);

            tmp[1] = ConvertToDegree(rb);
            tmp[2] = h;

            return tmp;

        }
        catch (System.Exception ex)
        {
            return null;
        }
    }

算B用迭代法，误差以参数形式给出，可以自定义误差范围。
下面是GPS时间的计算：
    /// &lt;summary&gt;
    /// 儒略日的计算
    /// &lt;/summary&gt;
    /// &lt;param name="year"&gt;&lt;/param&gt;
    /// &lt;param name="month"&gt;&lt;/param&gt;
    /// &lt;param name="day"&gt;&lt;/param&gt;
    /// &lt;returns&gt;&lt;/returns&gt;
    public static double jd(double year, double month, double day)
    {
        int G = 0;
        double n = 0;
        if(year * 372 + month *31 +Math.Floor(day)&gt;=588829)
            G=1; //判断是否为格里高利历日1582*372+10*31+15
        if (G==1)
        {
            n = Math.Floor(year / 100); 
            n = 2 - n + Math.Floor(n / 4); //加百年闰
        }
        if (month == 1 || month == 2)
        {
            month += 12;
            year -= 1;
        }
        return Math.Floor(365.25 * (year + 4716)) + Math.Floor(30.6001 * (month + 1)) + day + n - 1524.5;
        //return Convert.ToInt64(day - 32075 + 1461 * (year + 4800 + (month - 14) / 12) / 4 + 367 * (month - 2 - (month - 14) / 12 * 12) / 12 - 3 * ((year + 4900 + (month - 14) / 12) / 100) / 4);
    }

    public static string weekday(double jd) 
    {
        int no = Convert.ToInt32(Math.Floor((jd + 1.5) % 7));
        switch(no)
        {
        case 0:
            return "星期日";
        case 1:
            return "星期一";
        case 2:
            return "星期二";
        case 3:
            return "星期三";
        case 4:
            return "星期四";
        case 5:
            return "星期五";
        case 6:
            return "星期六";
        default:
            return null;
        }
    }</pre>

 
作业都要编程。。。真费事。。。THU坑爹。。。