1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
|
import numpy as np import matplotlib.pyplot as plt X=np.random.randint(0,20,(20,10)) b=np.random.randint(0,5,(10,1)) z=np.random.normal(size=[20,1]) y=np.dot(X,b)+z print(y) b_=np.linalg.lstsq(X, y, rcond=None)[0] print(b_) s_b=b.T s_b_=b_.T x = list(range(1,11)) plt.scatter(x, s_b, c='b', marker='o', label='true parameters')
plt.scatter(x, s_b_, c='r', marker='x', label='estimated parameters')
plt.legend()
plt.xlabel('index') plt.ylabel('value') plt.savefig('2.png') plt.show()
|
近期评论