一、问题描述
二、算法核心思想分析
三、代码及运行结果
a.py
import xlrd
import numpy as np
# 读取数据
def read_data(k):
x = []
data = xlrd.open_workbook("lab2_data.xlsx")
table = data.sheets()[0]
rows = table.nrows
for i in range(1, rows):
row_value = table.row_values(i)
if row_value[3] == k:
x.append(row_value)
return x
# 计算均值
def get_u(x):
u = np.mean(x)
return u
# 计算方差
def get_sigmal(x):
sigmal = np.cov(x)
return sigmal
def main():
u = [0, 0, 0]
sigmal = [0, 0, 0]
data = read_data(1)
for i in range(3):
xi = [x[i] for x in filter(lambda x: x, data)]
print("x", i + 1, xi)
u[i] = get_u(xi)
print("u", i + 1, u[i])
sigmal[i] = get_sigmal(xi)
print("sigmal", i + 1, sigmal[i])
if __name__ == '__main__':
main()
类ω1中的3个特征xi的均值和方差分别为:
b.py
import xlrd
import numpy as np
# 读取数据
def read_data(k):
x = []
data = xlrd.open_workbook("lab2_data.xlsx")
table = data.sheets()[0]
rows = table.nrows
for i in range(1, rows):
row_value = table.row_values(i)
if row_value[3] == k:
x.append(row_value)
return x
# 计算均值
def get_u(x):
u = np.mean(np.mat(x).T, axis=0) # 求每列的均值
return u
# 计算协方差
def get_sigmal(x):
sigmal = np.cov(x)
# sigmal = np.cov(np.mat(x).T)
return sigmal
def main():
u = [0, 0, 0]
sigmal = [0, 0, 0]
data = read_data(1)
for i in range(3):
xi = [x[i] for x in filter(lambda x: x, data)], [x[(i+1) % 3] for x in filter(lambda x: x, data)]
print("x", i + 1, (i+1) % 3 + 1, "\n", np.mat(xi))
u[i] = get_u(xi)
print("u", i + 1, (i+1) % 3 + 1, "\n", np.mat(u[i]))
sigmal[i] = get_sigmal(xi)
print("sigmal", i + 1, (i+1) % 3 + 1, "\n", np.mat(sigmal[i]))
if __name__ == '__main__':
main()
类ω1中任意两个特征组合的均值和方差为:
c.py
import xlrd
import numpy as np
# 读取数据
def read_data(k):
x = []
data = xlrd.open_workbook("lab2_data.xlsx")
table = data.sheets()[0]
rows = table.nrows
for i in range(1, rows):
row_value = table.row_values(i)
if row_value[3] == k:
x.append(row_value)
return x
# 计算均值
def get_u(x):
u = np.mean(x, axis=0)
return u
# 计算协方差
def get_sigmal(x):
sigmal = np.cov(np.mat(x).T)
return sigmal
def main():
data = read_data(1)
xi = [x[:3] for x in filter(lambda x: x, data)]
print("x", "\n", np.mat(xi))
u = get_u(xi)
print("u", "\n", np.mat(u))
sigmal = get_sigmal(xi)
print("sigmal", "\n", np.mat(sigmal))
if __name__ == '__main__':
main()
类ω1中3个特征组合的均值和方差为:
d.py
import xlrd
import numpy as np
# 读取数据
def read_data(k):
x = []
data = xlrd.open_workbook("lab2_data.xlsx")
table = data.sheets()[0]
rows = table.nrows
for i in range(1, rows):
row_value = table.row_values(i)
if row_value[3] == k:
x.append(row_value)
return x
# 计算均值
def get_u(x):
u = np.mean(x)
return u
# 计算方差
def get_sigmal(x):
sigmal = np.cov(x)
return sigmal
def main():
u = [0, 0, 0]
sigmal = [0, 0, 0]
data = read_data(2)
for i in range(3):
xi = [x[i] for x in filter(lambda x: x, data)]
print("x", i + 1, "\n", xi)
u[i] = get_u(xi)
print("u", i + 1, "\n", u[i])
sigmal[i] = get_sigmal(xi)
print("sigmal", i + 1, "\n", sigmal[i])
if __name__ == '__main__':
main()
类ω2中均值和协方差矩阵中的三个参数分别为:
如有错误请指正