untitled2.m

来自「就是一件小小的源码」· M 代码 · 共 22 行

% (1)-----选择拟合多项式拟合的阶数 
n=input('输入你所输入的时间序列的个数，n=')
x=1:n;
y=input('请输入需拟合的时间序列，y=');
% 计算多项式的各项系数和拟合值 
m=input('输入多项式拟合的阶数 m = '); 
p=nihe(x,y,m); 
y1=polyval(p,x);
yh=y-y1;
x1=x;
% (3)-----绘制观测数据离散点图和多项式曲线 
plot(x,y,'o',x1,y1,'r') 
title('多项式拟和时间序列')
% (4)-----拟合效果和精度检验 
Q=sum((yh).^2); 
zdwucha=max(yh);
SGM = sqrt(Q / (n - 1)); 
disp('p=');
fprintf (1,'      %3.6f\n ',p) 
fprintf (1,'      剩余平方和      Q = %3.6f \n',Q) 
fprintf (1,'      最大误差        zdwucha = %3.6f \n',zdwucha) 
fprintf (1,'      标准误差        Sigma = %3.6f \n',SGM)