clear;clc;x0=[127627 128453 129227 129988 130756];n=length(x0);lamda=x0(1:n-1)./x0(2:n)range=minmax(lamda)x1=cumsum(x0)for i=2:nz(i)=0.5*(x1(i)+x1(i-1));endB=[-z(2:n)',ones(n-1,1)];Y=x0(2:n)';u=B\Yx=dsolve('Dx+a*x=b','x(0)=x0');x=subs(x,{'a','b','x0'},{u(1),u(2),x1(1)});yuce1=subs(x,'t',[0:n-1]);digits(6),y=vpa(x) %為提高預(yù)測精度，先計(jì)算預(yù)測值，再顯示微分方程的解yuce=[x0(1),diff(yuce1)]epsilon=(x0-yuce)./x0 %計(jì)算殘差delta=abs(epsilon./x0) %計(jì)算相對誤差rho=1-(1-0.5*u(1))/(1+0.5*u(1))*lamda %計(jì)算級比偏差值