lomb scargle

谱分析方法一种，用于分析非等间隔信号最有效最方便的方法

function indata = inputtolomb
%CREATE TEST DATA FOR lombscargle.m

%Create time points
t=[0:pi/32:8*pi];

%Keep random time points discard the rest
%Note: To increase the number of points retained decrease the multiplier （and vice versa）
discpct = 30; %Percent of data points to discard
ndiscards = floor（（discpct/100）*length（t））;
keeps = randperm（length（t））;
keeps = sort（keeps（ndiscards+1:end））; %Discard 25% of the samples selected randomly
t=t（keeps）;

%USER-SPECIFIED FREQUENCIES AND AMPLITUDES （for testing purposes）
%Frequencies （Hz） = [0.9  3.3  4.0  4.5]
%Amplitudes       = [2.0  1.8  3.0  1.5]
foft = 2.0*sin（2*pi*0.9*t） + 1.8*sin（2*pi*3.3*t） + 3.0*sin（2*pi*4*t） + 1.5*sin（2*pi*4.5*t）;
indata=[t‘ foft‘];
%lombscargle（indata）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         763  2002-08-29 18:37  inputtolomb.m
     文件        1335  2009-05-26 13:51  license.txt
     文件       21682  2002-11-27 20:41  lombscargle.m