谐波分析matlab程序，可直接应用，也调整后应用，比如添加真值约束等(Inputs:
　ts　　=　ts(i)= i,
　　　　　where i=1, ..., ni
　HiLo　= 2-character string indicating rejection of high or low outliers
　low　 = valid range minimum
　high　= valid range maximum (values outside the valid range are rejeced
　　　　　right away)
　fet　 = fit error tolerance (points deviating more than fet from curve
　　　　　fit are rejected)
　dod　 = degree of overdeterminedness (iteration stops if number of
　　　　　points reaches the minimum required for curve fitting, plus
　　　　　dod). This is a safety measure
　delta = small positive number (e.g. 0.1) to suppress high amplitudes

Outputs:

amp　 = returned array of amplitudes, first element is the average of
　　　　the curve
phi　 = returned array of phases, first element is zero
yr　　= array holding reconstructed time series)