1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
function testpeph1
C=299792458.0;
dirs='../data/sp3';
td=caltomjd([2010,7,1]);
time=0:30:86400*2-30;
ephs=textread('testpeph1.out');
ephs(ephs==0)=nan;
sat=sprintf('GPS%02d',ephs(1,1));
ephr=readeph(td,time,{sat},dirs,'igs',24,'interp');
clkr=readclk(td,time,{sat},dirs,'igs',24,'interp');
for i=1:length(time)
clkr(i,1)=clkr(i,1)-2*ephr(i,1:3)*ephr(i,4:6)'/C^2;
end
dpos=ephs(:,3:6)-[ephr(:,1:3),clkr(:,1)*1E9];
figure('color','w'), hold on, box on, grid on
plot(time/3600,dpos)
xlabel('time (hr)');
ylabel('error (m)');
xlim(time([1,end])/3600);
ylim([-0.05,0.05]);
legend({'x','y','z','clk'})
text(0.02,0.95,sprintf('STD: X=%.4f Y=%.4f Z=%.4f CLK=%.4fm',...
std(dpos(~isnan(dpos(:,1)),1:3)),std(dpos(~isnan(dpos(:,4)),4))),...
'units','normalized')
title(sprintf('testpeph: interpolation error %s',sat));
moveax
|
