%計算基于最大李雅譜諾夫方法的預測值
function [x_1,x_2]=pre_by_lya(m,lmd,whlsj,whlsl,idx,min_d)
% x_1 - 第一預測值， x_2 - 第二預測值，
% m －嵌入維數，lmd - 最大李雅譜諾夫值，whlsj - 數據數組，whlsl - 數據個數， idx - 中心點的最近距離點位置， min_d - 中心點與最近距離點的距離

%相空間重構
LAST_POINT = whlsl-m+1;
for j=1:LAST_POINT           
    for k=1:m
        Y(k,j)=whlsj(k+j-1);
    end
end


a_1=0.;
for k=1:(m-1)
    a_1=a_1+(Y(k,idx+1)-Y(k+1,LAST_POINT))*(Y(k,idx+1)-Y(k+1,LAST_POINT));  % 此處Y(k+1,LAST_POINT)實際上就是Y(k,LAST_POINT+1)
end

deta=sqrt(min_d^2*2^(lmd*2)-a_1);
if (isreal(deta)==0) | (deta>Y(m,idx+1)*0.001)
    deta=Y(m,idx+1)*0.001
end
x_1=Y(m,idx+1)+deta;
x_2=Y(m,idx+1)-deta;