S_TIDE是基于T_TIDE開發(fā)的新一代潮汐調(diào)和分析工具包，其可以求解潮汐振幅和遲角的時間變化，同時也對T_TIDE存在的各種問題進(jìn)行了改進(jìn)，比如T_TIDE無法處理長時間數(shù)據(jù)（18.61年），無法處理不均勻采樣數(shù)據(jù)。

這里提供一個S_TIDE與T_TIDE對比的實例，詳細(xì)說明T_TIDE分析時可能出現(xiàn)的問題。程序使用了美國Astoria驗潮站的水位，該數(shù)據(jù)文件已經(jīng)包含在了s_tide工具包里。結(jié)果如下圖所示，可以看到S2分潮振幅在第59年有一個突然的下降，K1分潮振幅有一個突然的上升。這種異常情況是因為當(dāng)年水位觀測的時間計算錯誤，這種錯誤沒有被校正，而是保留在了原始數(shù)據(jù)中。

astoria=ncread('h572a.nc','sea_level');

%1925-01-26到2016-12-31

amp=zeros(91,67,'double');amp(:,:)=NaN;

for i=1:91

? ? ?aa=1+8766*(i-1);bb=8766+aa;

?if sum(isnan(astoria(aa:bb)))<8767*0.25

[NAME,FREQ,TIDECON,XOUT,z0]=t_tide(astoria(aa:bb),'interval',1);

?amp(i,:)=TIDECON(:,1);?

?mwl(i)=z0;

?end

end

plot(amp(:,20),'r.-');hold on;plot(amp(:,41),'k*-');

legend('K1','S2')

xlabel('Year');ylabel('Amplitude(mm)')

下面展示的S_TIDE處理結(jié)果，使用了魯棒擬合robustfit來消除時間錯亂的影響，結(jié)果如下圖所示，異常的突降和突升幾乎沒有了。

amp=zeros(91,67,'double');amp(:,:)=NaN;

for i=1:91

? ? ?aa=1+8766*(i-1);bb=8766+aa;

?if sum(isnan(astoria(aa:bb)))<8767*0.25

[St,Ht,Gt,coef,xout,ju,Stint,Htint,Gtint]=s_tide(astoria(aa:bb),1,1,'autoselected','autoselected',1,'spline','robustfit');

?amp(i,:)=Ht(:,1);?

?mwl(i)=z0;

?end

end

plot(amp(:,20),'r.-');hold on;plot(amp(:,41),'k*-');

legend('K1','S2')

xlabel('Year');ylabel('Amplitude(mm)')