使用S_TIDE预报任意时刻水位

首先去https://www.researchgate.net/project/A-non-stationary-tidal-analysis-toolbox-S-TIDE下载工具包，然后使用stide工具包里的s_construct2函数来预报，程序如下：

clear;clc;

%使用kushiro站八大主要分潮调和常数预报瞬时潮位

% this example was first added in S_TIDE v1.20 update

tides={'M2';'S2';'K1';'O1';'N2';'Q1';'P1';'K2'};

dt=1;length=720;

amp=[288.23,133.23,250.80,193.39,35.1736,35.3202,82.4656,36.9610];

pha=[177.14,217.38,24.42,356.23,159.69,345.72,21.14,211.69];

[sumconsti2,consti2]=s_construct2(tides,amp,pha,dt,length,[1994,01,01],42.5);

load kushiro.mat  %the observed water levels at Kushiro start from 1993/01/01

plot(sumconsti2);hold on

plot(kushiro(8761:8761+720)-1807.48,'r')%1807.48 is mean water level

xlim([1 350])

legend('预报','实测');

需要注意的是s_construct2使用的是中间时刻的交点因子和订正角，在1年或者更短的时间内，交点因子和订正角的可以看作是常数，所以使用中间时刻是合理的。但是要预报长期水位（比如1年以上），s_construct2是不合适的，需要使用s_construct3函数，两个函数的使用方法是一致的，但是需要指出的是s_construct3计算了每个时刻的交点因子和订正角，计算量很大，计算速度很慢。