本文共 604 字，大约阅读时间需要 2 分钟。

例题：

y=input('请输入数据');n=length(y);yy=ones(n,1);yy(1)=y(1);for i=2:n    yy(i)=yy(i-1)+y(i)endB=ones(n-1,2);for i=1:(n-1)    B(i,1)=-(yy(i)+yy(i+1))/2;    B(i,2)=1;endBT=B';for j=1:(n-1)    YN(j)=y(j+1);endYN=YN';A=inv(BT*B)*BT*YN;a=A(1);u=A(2);t=u/a;t_test=input('输入需要预测的个数');i=1:t_test+n;yys(i+1)=(y(1)-t).*exp(-a.*i)+t;yys(1)=y(1);for j=n+t_test:-1:2    ys(j)=yys(j)-yys(j-1);endx=1:n;xs=2:n+t_test;yn=ys(2:n+t_test);plot(x,y,'^r',xs,yn,'*-b');det=0;for i=2:n    det=det+abs(yn(i)-y(i));enddet=det/(n-1);disp(['百分绝对误差为：',num2str(det),'%']);    disp(['预测值为：',num2str(ys(n+1:n+t_test))]);

转载地址：http://tigwi.baihongyu.com/