% DW模型优化后的程序:
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];
y1=y(2:20,1);
y2=y(1:19,1);
y3=y1-0.71133y2;
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];
x1=x(2:20,1);
x2=x(1:19,1);
x3=x1-0.71133x2;
a=ones(19,1);
y4=y3-0.173x3+0.263a % 残差
y5=y4(2:19,1);
y6=y4(1:18,1);
y7=y5-y6;
y8=sum(y7.^2)
y9=sum(y5.^2)
DW= y8/y9
p=1-DW/2
y4 =
-0.0709
0.0412
0.1061
-0.0858
-0.0089
-0.0029
-0.1073
-0.0656
-0.1017
-0.0513
0.0296
0.0091
0.1171
0.0121
0.0486
-0.0395
0.0249
-0.0641
0.0112
y8 =
0.1325
y9 =
0.0728
DW =
1.8207
p =
0.0897