研究潮汐季节循环需要注意的事项

代码所需要的数据kushiro.mat在s_tide工具包里https://www.researchgate.net/project/A-non-stationary-tidal-analysis-toolbox-S-TIDE
clear
load kushiro.mat
for i=1:36

aa=1+720*(i-1);bb=720*i;
[NAME,FREQ,TIDECON,XOUT]=t_tide(kushiro(aa:bb),'interval',1);
K1(i)=TIDECON(6,1);
S2(i)=TIDECON(12,1);
end
plot(K1);hold on;plot(S2,'r') 结果如下图所示，可以看到K1和S2都存在显著的半年变化。需要注意的是逐月的潮汐调和分析是没法把P1和K1区分，以及K2和S2区分。而P1的频率恰好处在K1的半年变化上，而K2的频率恰好处在S2的半年变化上。所以下图实线所示的K1和S2振幅的半年变化是虚假的。

下面的程序使用了差比关系，根