4-1预测与回归

xyuanshi=[0.18095 0.19382 0.18382 0.14627];
n=length(xyuanshi);
%%%%%%%对数据做均匀化处理%%%%%%
x(1)=(3*xyuanshi(1)+xyuanshi(2))/4;
x(n)=(3*xyuanshi(n)+xyuanshi(n-1))/4;
for i=2:(n-1)
  x(i)=(xyuanshi(i-1)+xyuanshi(i+1)+2*xyuanshi(i))/4;
end

x1(1)=x(1);
for i=2:n
    x1(i)=x1(i-1)+x(i);
end

%生成并使用
for i=1:(n-1)
    x2(i)=(x1(i)+x1(i+1))/2; %生成累加矩阵
end
%计算待定参数的值
Y=x(2:n)';
B(:,1)=-x2';
B(:,2) = 1;
alph=(B'*B)^(-1)*B'*Y;

X1end(1) = x1(1);
for i=2:(n + 1)   %只推测后1个数据，可以从此修改
    X1end(i)=(X1end(1) - (alph(2)/alph(1)) )*exp(-alph(1)*(i-1)) + alph(2)/alph(1);
end

%Xlend：预测结果
%xend：实际相对误差
Xend(1) = X1end(1);   
for i=2:(n+1)   %只推测后1个数据，可以从此修改
    Xend(i)=X1end(i)-X1end(i-1); %得到预测出来的数据
end
Errorxd=(Xend(1:n)-x)./x;