Генералов Кирилл: Новая страница: «import numpy as np import math f = open("order.txt", "r") X = [] Y = [] for line in f: q, x1, x2, x3, x4, x5, x6, y = map(float, line.split()) x = [] …»

2019-06-07T14:16:44Z

Новая страница: «import numpy as np import math f = open("order.txt", "r") X = [] Y = [] for line in f: q, x1, x2, x3, x4, x5, x6, y = map(float, line.split()) x = [] …»

Новая страница

import numpy as np
import math
f = open("order.txt", "r")
X = []
Y = []
for line in f:
q, x1, x2, x3, x4, x5, x6, y = map(float, line.split())
x = []
x.append(x1)
x.append(x2)
x.append(x3)
x.append(x4)
x.append(x5)
x.append(x6)
x.append(x1 * x1)
x.append(x1 * x2)
x.append(x1 * x3)
x.append(x1 * x4)
x.append(x1 * x5)
x.append(x1 * x6)
x.append(x2 * x2)
x.append(x2 * x3)
x.append(x2 * x4)
x.append(x2 * x5)
x.append(x2 * x6)
x.append(x3 * x3)
x.append(x3 * x4)
x.append(x3 * x5)
x.append(x3 * x6)
x.append(x4 * x4)
x.append(x4 * x5)
x.append(x4 * x6)
x.append(x5 * x5)
x.append(x5 * x6)
x.append(x6 * x6)
X.append(x)
Y.append(y)
f.close()
Mx = []
for j in range(27):
Mx1 = np.sum(X[j]) / len(X)
Mx.append(Mx1)
My = np.sum(Y)/len(Y)
print(Mx)
print("Мат.ожидание(Mx) =", Mx)
Dx = []
Dx11 = 0
Dy = 0

for i in range(len(X)):
Dy += (Y[i] - My) ** 2
for i in range(27):
Dxs = 0.0
for j in range(len(X)):
Dxs += (X[j][i] - Mx[i]) ** 2
print(Dxs)

Dx.append(Dxs/len(X))

print("Дисперсия(Dx) = ", Dx)
kx = []
for i in range(27):
kxs = 0.0
for j in range(len(X)):
kxs += ((X[j][i] - Mx[i]) * (Y[j] - My))
kx.append(kxs)
print("Коэффициент корреляции (kx)=", kx)
r = []
rABS = []
for i in range(len(kx)):
r1 = (kx[i] / (len(X)-1)) / (math.sqrt(Dx[i]) * math.sqrt(Dy))
r.append(r1)
rABS.append(math.fabs(r1))
print("r =", r)
print(np.max(rABS))
print("Первый индекс = ", rABS.index(np.max(rABS)) + 1)

Генералов Кирилл N2 - История изменений

Генералов Кирилл: Новая страница: «import numpy as np import math f = open("order.txt", "r") X = [] Y = [] for line in f: q, x1, x2, x3, x4, x5, x6, y = map(float, line.split()) x = [] …»