工程計算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
Procedure TSPSA:
begin
init-of-T { T為初始溫度}
S={1，……，n} {S為初始值}
termination=false
while termination=false
begin
for i=1 to L do
begin
generate(S′form S) { 從當前回路S產生新回路S′}
Δt:=f(S′))-f(S) {f(S)為路徑總長}
IF(Δt<0) OR (EXP(-Δt/T)>Random-of-[0,1])
S=S′
IF the-halt-condition-is-TRUE THEN
termination=true
End
T_lower
End
End
一臺精密儀器的工作時間為 n 個時間單位。與儀器工作時間同步進行若干儀器維修程序。一旦啟動維修程序，儀器必須進入維修程序。如果只有一個維修程序啟動，則必須進入該維修程序。如果在同一時刻有多個維修程序，可任選進入其中的一個維修程序。維修程序必須從頭開始，不能從中間插入。一個維修程序從第s個時間單位開始，持續 t個時間單位，則該維修程序在第s+t-1個時間單位結束。為了提高儀器使用率，希望安排盡可能少的維修時間。
