潮汐调和分析

来自潘海东，最终知识产权属于他本人，若有转发，请署名来源

function [S] = l_calculate_coefficient（IDPntime）
% computes the linear interpolation weights
% Please citing:
% Guo Z. H. Pan A. Cao and X. Lv （2018） A harmonic analysis method adapted to capturing slow variations of tidal amplitudes and phases 
%                                             Cont. Shelf Res. 164（June） 37–44 doi:10.1016/j.csr.2018.06.005.
r = （ntime-1）/（IDP-1）;   
S = zeros（ntimeIDP‘double‘）; 
for t=1:ntime
    k=floor（（t-1）/r）;
    if k==IDP-1
        k=k-1;
    end
    xi=k*r+1;xi1=（k+1）*r+1;
    k1=k+1;k2=k+2;
    S（tk1）=1-（t-xi）/r;
    S（tk2）=（t-xi）/r;   
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-09-08 11:34  s_tide_v1.1\
     文件    21020431  2017-05-05 14:04  s_tide_v1.1\imf.mat
     文件      211912  2017-05-05 14:00  s_tide_v1.1\kushiro.mat
     文件         615  2018-07-29 19:29  s_tide_v1.1\l_calculate_coefficient.m
     文件        1032  2018-08-21 13:28  s_tide_v1.1\s_calculate_coefficient.m
     文件        1788  2018-09-08 11:55  s_tide_v1.1\s_demo.m
     文件         969  2018-07-29 18:42  s_tide_v1.1\s_nodal_correction.m
     文件        2321  2018-09-08 11:18  s_tide_v1.1\s_readme.m
     文件        5327  2018-09-08 14:27  s_tide_v1.1\s_tide.m
     文件        3385  2018-07-29 19:10  s_tide_v1.1\s_tide_error.m
     文件        4395  2003-02-10 02:06  s_tide_v1.1\t_astron.m
     文件       29224  2003-02-11 01:31  s_tide_v1.1\t_constituents.mat
     文件        5688  2010-04-13 14:56  s_tide_v1.1\t_getconsts.m
     文件        4678  2010-04-13 14:56  s_tide_v1.1\t_vuf.m