nvs.m

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

% number of visible satellites
sec_day = 24*60; % 24小时内的分钟数
t_span = 1;        % 数据更新间隔分钟数
t=0:t_span:sec_day; %　每分钟更新一次数据，记录24小时内数据
nvs_gps = 6 + 0.2 * rand(size(t));
nvs_gal = 7 + 0.2 * rand(size(t));
plot(t,nvs_gps,'r-')
axis([0,sec_day,0,10])
hold on
plot(t,nvs_gal,'b*')
plot(t,2+sin((t/100)*pi),'g')
xlabel('t');    ylabel('n');
hold off