说明：  谐波分析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)