clc,clear,close all
y=[20.96;21.40;21.96;21.52;22.39;22.76;23.48;23.66;24.10;24.01;24.54;24.30;25.00;25.64;26.36;26.98;27.52;27.78;28.24;28.78];
x=[127.3;130.0;132.7;129.4;135.0;137.1;141.2;142.8;145.5;145.3;148.3;146.4;150.2;153.1;157.3;160.7;164.2;165.6;168.7;171.7];
rstool(x,y)
a=ones(20,1);
b=0.176x-1.455a;
c=y-b % 残差
c1=c(2:20,1);
c2=c(1:19,1);
plot(c1,c2,‘*’)
xlabel(‘e(i-1)’),
ylabel(‘e(i)’)
hold on
d=0;
d1=-0.15:0.001:0.25;
plot(d,d1)
hold on
plot(d1,d)
% DW的程序:
num=c(1:19)'c(2:20)
den=sum(c(2:20).^2)
rou = num/den
DW=2(1-rou)
c