getgpssatellite.m

来自「根据卫星星历及测者位置」· M 代码 · 共 34 行

function GPSsat = getGPSSatellite( n )
% 24 个GPS卫星
GPSsatellites = [ [55  17 26562000 0.01  89 240]; %1
    [55  17 26562000 0.01  89 330] ;%2
    [55  17 26562000 0.01  89  60] ;%3
    [55  17 26562000 0.01  89 150] ;%4
    [55  77 26562000 0.01 179 270] ;%5
    [55  77 26562000 0.01 179   0] ;%6
    [55  77 26562000 0.01 179  90] ;%7
    [55  77 26562000 0.01 179 180] ;%8
    [55 137 26562000 0.01  89 330] ;%9
    [55 137 26562000 0.01  89  30] ;%10
    [55 137 26562000 0.01  89 120] ;%11
    [55 137 26562000 0.01  89 210] ;%12
    [55 197 26562000 0.01 179 330] ;%13
    [55 197 26562000 0.01 179  60] ;%14
    [55 197 26562000 0.01 179 150] ;%15
    [55 197 26562000 0.01 179 240] ;%16
    [55 257 26562000 0.01  89   0] ;%17
    [55 257 26562000 0.01  89  90] ;%18
    [55 257 26562000 0.01  89 180] ;%19
    [55 257 26562000 0.01  89 270] ;%20
    [55 317 26562000 0.01 179  30] ;%21
    [55 317 26562000 0.01 179 120] ;%22
    [55 317 26562000 0.01 179 210] ;%23
    [55 317 26562000 0.01 179 300]; ];%24
while n<24
    n=n+24;
end
if n>24
    n=rem(n,24);%求余数
end
GPSsat = GPSsatellites( n, : );