题目三
代码:
from async_generator import yield_
import numpy as np
from numpy import *
x0=ones(10)
x1=np.array([64.3,99.6,145.45,63.75,135.46,92.85,86.97,144.76,59.3,116.03])
x2=np.array([2, 3, 4, 2, 3, 4, 2, 4, 1, 3])
y=np.array([62.55,82.42,132.62,73.31,131.05,86.57,85.49,127.44,55.25,104.84])
print("X=")
X=stack((x0,x1,x2),axis=1)
print(X)
print("Y=")
Y=vstack(y)
print(Y)
X=mat(X)
Y=mat(Y)
print("W=")
Xt=transpose(X)
mul_x=dot(Xt,X)
MT=linalg.inv(mul_x)
W=dot(MT,Xt)
W=dot(W,y)
W=transpose(W)
print(W)
print(W.shape)
结果: