准调和分析程序：S_TIDE实现短期观测的准调和分析

网上似乎没有准调和分析的公开程序，而现实中我们又需要对短期观测（比如25小时）的数据做调和分析，由于数据太短，ttide等经典工具包无法把S2和M2以及K1和O1分开，这里我介绍一下stide新添加的s_quasi_HA函数，可以实现准调和分析，对于很短的数据也可以把四大主要分潮分开，这里核心是用了基于卫星观测开发的潮汐模型EOT20来提供全球任意点的M2/S2以及K1/O1之间的振幅和迟角关系，用于调和分析（即差比关系）。

实例如下：

load Qingdao.mat

 [S,H,G,coef,xout,ju,coefint,Ra,Rb]=s_quasi_HA(wl2(1:25),'1997/08/01 00:00:00',25,1,35.23,119.33);

 plot(wl2(1:25),'r.');hold on;plot(xout,'k')