使用S_TIDE分离潮汐数据里不同分潮示例

首先去（https://www.researchgate.net/project/A-non-stationary-tidal-analysis-toolbox-S-TIDE）下载S_TIDE工具包，并且安装到MATLAB当前目录下。然后运行如下程序，kushiro是日本一个验潮站的逐时水位数据，s_tide输入输出参数的含义见s_tide函数说明以及工具包中文教程，其中输出参数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分潮叠加水位也是如此。

这里只分辨了特定的四大主要分潮，如果想按照瑞利准则自动分辨能分辨的所有分潮，可以用如下命令：

[St,Ht,Gt,coef,xout,ju,Stint,Htint,Gtint]=s_tide(kushiro(1:8767),1,1,'autoselected','autoselected',1,'spline','ols');