function adjmatrix = LPHVG(time_series,L)
[num,len,~] = size(time_series);
time_series = abs(fft(time_series))/len;
min_value=min(time_series);
if (min_value<0)
    time_series = time_series+abs(min_value);
end
[num,len,~] = size(time_series);
adjmatrix = zeros(len,len);
for k = 1:len
    Y = zeros(1,L+1);
    for l = k+1:len
        if abs(k-l)<=L+1
            adjmatrix(k,l)=1;adjmatrix(l,k)=1;
            if Y(1,1)<time_series(l)
               Y(1,1)=time_series(l);
               Y=sort(Y);
            end
        elseif time_series(k)>Y(1,1) && time_series(l)>Y(1,1)
            adjmatrix(k,l)=1;adjmatrix(l,k)=1;
            Y(1,1)=time_series(l);
            Y=sort(Y);
        end
    end
end