海洋测绘——GPS在航潮位提取

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（‘横    轴：    时         间           轴‘）
% ylabel（‘纵轴：GPS天线相位中心的大地高‘）
% title（‘时间和GPS大地高的坐标‘）

scatter（A（:2）A（:6）‘c*‘）
hold on
scatter（A（:2）F‘k.‘）
xlabel（‘横     轴：    时                 间  ‘）
ylabel（‘纵轴：GPS大地高（（青）上‘*’），海面高（（黑）下）‘.’‘）
title（‘GPS天线大地高与海面高的关系‘）

% 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（‘横        轴 ：