A=load（‘F:\Dcor\001.txt‘）;
B=A（:8）*pi/180;
C=A（:9）*pi/180;
D=A（:7）;
T=A（:2）;
for?i=1:5991
D（i）=cos（B（i））*cos（C（i））;
end
E=-2*D-A（:7）;
F=A（:6）+E;
%?
%?scatter（A（:2）A（:6）‘m*‘）
%?xlabel（‘橫????軸：????時(shí)?????????間???????????軸‘）
%?ylabel（‘縱軸：GPS天線相位中心的大地高‘）
%?title（‘時(shí)間和GPS大地高的坐標(biāo)‘）

scatter（A（:2）A（:6）‘c*‘）
hold?on
scatter（A（:2）F‘k.‘）
xlabel（‘橫?????軸：????時(shí)?????????????????間??‘）
ylabel（‘縱軸：GPS大地高（（青）上‘*’），海面高（（黑）下）‘.’‘）
title（‘GPS天線大地高與海面高的關(guān)系‘）

%?TOU=F（1）*ones（10001）;
%?CTOU=F（478）*ones（20001）;
%?TF=[TOU;F（1:1:477）;CTOU;F（478:1:5991）];
%?[NWn]=buttord（1/57601/2455）;
%?[ba]=butter（NWn）;
%?y=filter（baTF）;
%?yy=[y（1001:1:1477）;y（3478:1:8991）];
%?scatter（A（:2）F‘c.‘）
%?hold?on
%?scatter（A（:2）yy‘k*‘）
%?xlabel（‘橫????????軸?：??????