matlab的lomb函数

function [fPprob] = lomb（thofachifac）
% LOMB（THOFACHIFAC） computes the Lomb normalized periodogram （spectral
% power as a function of frequency） of a sequence of N data points H
% sampled at times T which are not necessarily evenly spaced. T and H must
% be vectors of equal size. The routine will calculate the spectral power
% for an increasing sequence of frequencies （in reciprocal units of the
% time array T） up to HIFAC times the average Nyquist frequency with an
% oversampling factor of OFAC （typically >= 4）.
% 
% The returned values are arrays of frequencies considered （f） the
% associated spectral power （P） and estimated significance of the power
% values （prob）.  Note: the significance returned is the false alarm
% probability of the null hypothesis i.e. that the data is composed of
% independent gaussian random variables.  Low probability values indicate a
% high degree of significance in the associated periodic signal.
% 
% Although this implementation is based on that described in Press
% Teukolsky et al. Numerical Recipes  In C section 13.8 rather than using
% trigonometric rercurrences this takes advantage o