網(wǎng)上似乎沒(méi)有準(zhǔn)調(diào)和分析的公開(kāi)程序，而現(xiàn)實(shí)中我們又需要對(duì)短期觀測(cè)（比如25小時(shí)）的數(shù)據(jù)做調(diào)和分析，由于數(shù)據(jù)太短，ttide等經(jīng)典工具包無(wú)法把S2和M2以及K1和O1分開(kāi)，這里我介紹一下stide新添加的s_quasi_HA函數(shù)，可以實(shí)現(xiàn)準(zhǔn)調(diào)和分析，對(duì)于很短的數(shù)據(jù)也可以把四大主要分潮分開(kāi)，這里核心是用了基于衛(wèi)星觀測(cè)開(kāi)發(fā)的潮汐模型EOT20來(lái)提供全球任意點(diǎn)的M2/S2以及K1/O1之間的振幅和遲角關(guān)系，用于調(diào)和分析（即差比關(guān)系）。

實(shí)例如下：

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')