求大神解答。 哪里错了
原步骤如下
functionGM1_1(x0 ) )。
格式龙；
x0=[ 132802，133450，134091，134735，135404，136072，136782，137462，138271 ]；
[m，n]=size(x0；
lamda=x0(1:n-1 )./x0 ) 2:n )；
x1=cumsum(x0 )； %累积
X2=[]；
for i=1:8
x2(I， )=x1 ) I ) x1 ) I1；
结束
B=-0.5.X2；
t=ones (n-1，1 )；
B=[B，t]； 求% b矩阵
yn=x0 ) 2:结束；
p_t=yn/x1(1: ) length(x0(-1 )； 对%原始数据序列X0进行标准平滑性检查，求出对序列X0平滑比p(t )=x0 ) t )/x1 ) t-1 )
a=inv(b.‘b ) B.'YN.)；
a=a(1；
u=a(2；
c=u/a；
b=x0(1)-c；
x=[num2str(b )、‘exp’、’)、num2str(-a )、‘k’、’)、num2str]；
strcatx(k1 )=X；
syms k；
fort=1:length(x0 ) )
k(1，t )=t-1；
结束
k；
y_k_1=bexp(-AK ) c；
forj=1:length(k )-1
y(1，j )=y_k_1) j1 )-Y_k_1(j ) j；
结束
xy=[y_k_1(1)，y ) %预测值
CA=XY-X0； %残差数列
西塔=ca； %残差检验绝对误差序列
XD_Theta=CA./X0； %残差检验相对误差
av=mean(ca )； %残差数列平均值
TEMP0=(ca-av ).^2；
TEMP1=sum(TEMP0)/Length ) ca；
S2=sqrt(temp1)； %绝对误差序列的标准差
av_0=mean(x0 )； %原始序列的平均值
temp_0=(x0-av_0).^2；
temp_1=sum(temp_0)/length ) ca；
S1=sqrt(temp_1)； %原始序列的标准偏差
TempC=S2/S1100； %方差比
c=strcat(num2str ) tempc )，’%’； %事后检查差检查
SS=0.675S1；
ELTA=ABS(ca-av )；
tempn=find (增量=ss )；
N1=Length(tempn )；
N2=Length(ca )；
TempP=N1/N2*100；
p=strcat(num2str ) tempp )，’%’； %后验差检查计算小误差概率