gmnewxm0.m

来自「这是一个有预测功能的程序」· M 代码 · 共 14 行

function [Xm0]=gmnewXm0(x,T,N)
n=length(x);
m=N+n;
xT_1=gmT(T-1,x);
xT=gmT(T,x);
q=getab(xT_1,xT);
a=q(1);
b=q(2);
for i=1:m
  xm1(i)=(x(1)-b/a)*exp(-a*(i-1))+b/a;
end
  xm0=xm1;
Xm0=gm_T(T,xm0);