Лабы / ml6

X = 20:3:80
Y = [
2.62
2.41
2.24
2.11
1.91
1.81
1.64
1.54
1.44
1.36
1.28
1.21
1.16
1.11
1.07
1.04
1.002
0.982
0.946
0.94
0.92
];
Y = Y'

M = [sum(X.^6) sum(X.^5) sum(X.^4) sum(X.^3);
sum(X.^5) sum(X.^4) sum(X.^3) sum(X.^2);
sum(X.^4) sum(X.^3) sum(X.^2) sum(X);
sum(X.^3) sum(X.^2) sum(X) 21];
A = [sum(Y.*X.^3) sum(X.^5) sum(X.^4) sum(X.^3);
sum(Y.*X.^2) sum(X.^4) sum(X.^3) sum(X.^2);
sum(Y.*X) sum(X.^3) sum(X.^2) sum(X);
sum(Y) sum(X.^2) sum(X) 21];
B= [sum(X.^6) sum(Y.*X.^3) sum(X.^4) sum(X.^3);
sum(X.^5) sum(Y.*X.^2) sum(X.^3) sum(X.^2);
sum(X.^4) sum(Y.*X) sum(X.^2) sum(X);
sum(X.^3) sum(Y) sum(X) 21];
C = [sum(X.^6) sum(X.^5) sum(Y.*X.^3) sum(X.^3);
sum(X.^5) sum(X.^4) sum(Y.*X.^2) sum(X.^2);
sum(X.^4) sum(X.^3) sum(Y.*X) sum(X);
sum(X.^3) sum(X.^2) sum(Y) 21];
D = [sum(X.^6) sum(X.^5) sum(X.^4) sum(Y.*X.^3);
sum(X.^5) sum(X.^4) sum(X.^3) sum(Y.*X.^2);
sum(X.^4) sum(X.^3) sum(X.^2) sum(Y.*X);
sum(X.^3) sum(X.^2) sum(X) sum(Y)]
a = det(A)/det(M)
b = det(B)/det(M)
c = det(C)/det(M)
d = det(D)/det(M)
V = a*X.^3+ b*X.^2+c*X+d-Y
SSE =sum(V.^2)
S_A = [sum(X.^4) sum(X.^3) sum(X.^2);
sum(X.^3) sum(X.^2) sum(X);
sum(X.^2) sum(X) 21];
SB = [sum(X.^6) sum(X.^4) sum(X.^3);
sum(X.^4) sum(X.^2) sum(X);
sum(X.^3) sum(X) 21];
SC = [sum(X.^6) sum(X.^5) sum(X.^3);
sum(X.^5) sum(X.^4) sum(X.^2);
sum(X.^3) sum(X.^2) 21];
SD = [sum(X.^6) sum(X.^5) sum(X.^4);
sum(X.^5) sum(X.^4) sum(X.^3);
sum(X.^4) sum(X.^3) sum(X.^2)];
SE = SSE./19
sa = det(S_A)*SE/det(D)
sb = det(SB)*SE/det(D)
sc = det(SC)*SE/det(D)
sd = det(SD)*SE/det(D)
COV = cov(X,Y)
RO = corrcoef(X,Y)