工程計(jì)算MATLAB code to calculate the reorthogonalized sine tapers
input: N = the length of the time series data to be tapered
p = the number of tapers requested
I = the gap structure a vector of length N
I(t) = 1 if there is data at time t, t=1, ..., N
I(t) = 0 if there is a gap at time t
output: X = N-by-p vector of the reorthogonalized sine taper
給定n 個(gè)整數(shù)a ,a , ,an 1 2 組成的序列， a n i | |£ ，1 £ i £ n。如果對(duì)于i £ j ，有
0 = å
=
j
k i
k a ，則稱序列區(qū)間i i j a , a , , a +1 為一個(gè)零和區(qū)間，相應(yīng)的區(qū)間長(zhǎng)度為j-i+1。
給定n 個(gè)整數(shù)a ,a , ,an 1 2 組成的序列， a n i | |£ ，1 £ i £ n。如果對(duì)于i £ j ，有
0 = å
=
j
k i
k a ，則稱序列區(qū)間i i j a , a , , a +1 為一個(gè)零和區(qū)間，相應(yīng)的區(qū)間長(zhǎng)度為j-i+1。
