verhulst代码

function [Xcerror1error2]=verr（X1k）
% 建立函数[Xcerror1error2]=example9_3_2_3（X1k）
% 其中X1为输入序列，k为预测长度，
% X为预测输出序列，c为后验差检验数，error1为残差，error2为相对误差
format long;
n=length（X1）;
X0=[];
X0（1）=X1（1）;
for i=2:n
    X0（i）=X1（i）-X1（i-1）;             %计算累减生成序列
end
for i=1:n-1
B（i1）=-0.5*（X1（i）+X1（i+1））;   %计算B，Yn
    B（i2）=0.25*（X1（i）+X1（i+1））^2;
    Y（i）=X0（i+1）;
end
alpha=（B‘*B）^（-1）*B‘*Y‘;            %做最小二乘估计
a=alpha（11）;
b=alpha（21）;
c=a*X1（1）;