用matlab计算卫星的瞬时位置,运用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=')