function?[oHDEV]=calhdev（tausPeriodreadings）??%相位數據計算hadamard方差，speriod是采樣時間間隔，reading輸入相位數據序列
oHDEV=[];?
for?i=1:length（tau）
????oHDEV（1i）=calculateHDEV（tau（1i）sPeriodreadings）
end
function?[oH]?=?calculateHDEV（tausPeriodreadings）
N?=?nume