工程計算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 個整數a ,a , ,an 1 2 組成的序列, a n i | |£ ,1 £ i £ n。如果對于i £ j ,有
0 = å
=
j
k i
k a ,則稱序列區間i i j a , a , , a +1 為一個零和區間,相應的區間長度為j-i+1。
給定n 個整數a ,a , ,an 1 2 組成的序列, a n i | |£ ,1 £ i £ n。如果對于i £ j ,有
0 = å
=
j
k i
k a ,則稱序列區間i i j a , a , , a +1 為一個零和區間,相應的區間長度為j-i+1。