%該程序用加權一階局域法對完好率數據進行進行一步預測
%skyhawk
clear all;
data=load('bk.txt');%讀取數據
whl=data(:,4);
[whsl,lllll]=size(whl);
N=80;%預測點數
m=6;
fch=0.;
for i=whsl-N+1:whsl         %預測后N個點
    for j=1:(i-1)            
        whlsj(j)=whl(j);
    end       
    [y(i),y1(i),y2(i)]=jiaquanyijie(m,whlsj,i-1);%預測第i+1個點 
    fch=fch+(y(i)-whl(i))*(y(i)-whl(i));%計算標準差　
    clear whlsj;
    i;
end 
fch=sqrt(fch)/N

yyy=[whl,y'];
save('bkyc.txt','yyy','-ASCII');

kk=1:whsl;
plot(kk,whl,'b',kk,y,'r')