首先去https://www.researchgate.net/project/A-non-stationary-tidal-analysis-toolbox-S-TIDE下載工具包，然后使用stide工具包里的s_construct2函數(shù)來(lái)預(yù)報(bào)，程序如下：

clear;clc;

%使用kushiro站八大主要分潮調(diào)和常數(shù)預(yù)報(bào)瞬時(shí)潮位

% 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('預(yù)報(bào)','實(shí)測(cè)');

需要注意的是s_construct2使用的是中間時(shí)刻的交點(diǎn)因子和訂正角，在1年或者更短的時(shí)間內(nèi)，交點(diǎn)因子和訂正角的可以看作是常數(shù)，所以使用中間時(shí)刻是合理的。但是要預(yù)報(bào)長(zhǎng)期水位（比如1年以上），s_construct2是不合適的，需要使用s_construct3函數(shù)，兩個(gè)函數(shù)的使用方法是一致的，但是需要指出的是s_construct3計(jì)算了每個(gè)時(shí)刻的交點(diǎn)因子和訂正角，計(jì)算量很大，計(jì)算速度很慢。