GPS卫星-坐标计算

View1.Items.Add(li);Р M[i] = double.Parse(arrs[i]);Р }Р double t0e=460800.00; Р double t;Р double[] XK = new double[21];Р double[] YK = new double[21];Р double[] ZK = new double[21];Р for (int l = 0; l < 21;l++ )Р {Р //计算平均角速度Р double GM = 3986004.418e008;Р double n, n0;Р n0 = Math.Sqrt(GM / (Math.Pow(M[17], 6)));Р n = M[12] + n0;Р //规划时刻Р t=t0e+l*60;Р double tk = t-t0e;Р double Mk = M[13] + n * tk;Р // 迭代计算平近点角的计算Р double Ek, Ek1;Р Ek = Mk;Р Ek1 = Mk + M[15] * Math.Sin(Ek);Р doР {Р Ek = Ek1;Р Ek1 = Mk + M[15] * Math.Sin(Ek);Р }Р while (Math.Abs(Ek1 - Ek) > 1e-15);Р //计算真近点角Р double Vk = Math.Atan(((Math.Sqrt(1 - M[15] * M[15]) * Math.Sin(Ek))) / (Math.Cos(Ek) - M[15]));Р //反正切值的象限处理Р double A1 = Math.Atan(((Math.Sqrt(1 - M[15] * M[15]) * Math.Sin(Ek))));Р double A2 = Math.Cos(Ek) - M[15];