使用S_TIDE分離潮汐數(shù)據(jù)里不同分潮示例
首先去（https://www.researchgate.net/project/A-non-stationary-tidal-analysis-toolbox-S-TIDE）下載S_TIDE工具包，并且安裝到MATLAB當(dāng)前目錄下。然后運(yùn)行如下程序，kushiro是日本一個驗(yàn)潮站的逐時(shí)水位數(shù)據(jù)，s_tide輸入輸出參數(shù)的含義見s_tide函數(shù)說明以及工具包中文教程，其中輸出參數(shù)consti代表各分潮水位
clear;load kushiro.mat %data start from 1993/01/01
[St,Ht,Gt,coef,xout,ju,Stint,Htint,Gtint,consti]=s_tide(kushiro(1:8767),1,1,{'M2';'S2';'K1';'O1'},4,1,'spline','ols');
figure()
plot(consti(1,1:720)+consti(2,1:720))
xlabel('Time(hour)')
title('M2 tide + S2 tide')?
figure()
plot(consti(3,1:720)+consti(4,1:720))
xlabel('Time(hour)')
title('K1 tide + O1 tide')?
可以考到m2和S2分潮的疊加水位有明顯的半月波動，K1和O1分潮疊加水位也是如此。
這里只分辨了特定的四大主要分潮，如果想按照瑞利準(zhǔn)則自動分辨能分辨的所有分潮，可以用如下命令：
[St,Ht,Gt,coef,xout,ju,Stint,Htint,Gtint]=s_tide(kushiro(1:8767),1,1,'autoselected','autoselected',1,'spline','ols');
