代码如下:
import matplotlib.pyplot as plt
import numpy as np
from scipy import stats
#11.1
x=np.linspace(0.0,2.0,50)
y=np.sin(x-2)*np.sin(x-2)*np.exp(-x**2)
plt.figure(1)
plt.xlabel("x")
plt.ylabel("y")
plt.title(r"$f(x) = sin^2(x-2)e^{-x^2}$")
plt.plot(x,y)
plt.show()
得到结果如下:
代码如下:
#11.2
X=np.random.normal(size=(20,10))
b=np.random.normal(size=(10,1))
z=np.random.normal(size=(20,1))
y=np.dot(X,b)+z
be=np.linalg.lstsq(X,y,None)[0]
plt.figure(2)