matlab计算gps卫星的坐标位置

clear
clc
format long
tp=input('tp=');
toc=input('toc=');
a0=input('a0=');
a1=input('a1=');
a2=input('a2=');
toe=input('toe=');
M0=input('M0=');
a=input('长半径a=');
deltan=input('卫星平均角速度之差deltan=');
e=input('e=');
w=input('w=');
Cuc=input('Cuc=');
Cus=input('Cus=');
Cic=input('Cic=');
Cis=input('Cis=');
Crc=input('Crc=');
Crs=input('Crs=');
i0=input('i0=');
I=input('轨道倾角变化率I=');
OM0=input('OM0=');
OM=input('升交点赤径变化率OM=');
tt=a0+a1*(tp-toc)+a2*(tp-toc);
t=tp-tt;
tk=t-toe;
u=3.986005e14;
n0=(sqrt(u))/(a*a*a);
n=n0+deltan;