第十二周作业（matplotlib）

import numpy as np 
import matplotlib.pyplot as plt 
from scipy import stats 
#exercise 11.1
def f(x):
	return np.power((np.sin(x - 2)), 2) * np.exp(-(np.power(x, 2)))
x = np.linspace(0, 2 , 100)        #利用数量大的离散点来绘制连续图像
y = f(x)
plt.subplot(221)          #构造第一个子图
plt.plot(x, y,label = '$sin^2xe^-x^2$')
plt.xlabel('x')
plt.ylabel('y')
plt.title('My Function')
plt.legend()          #将label显示出来
#exercise 11.2
X = np.random.random((20,10))
b = np.random.random((10,1))
z = np.random.normal(0, 1, size = (20,1))
y = np.dot(X, b) + z
X = np.mat(X)
b1 = np.dot(X, b) - y
plt.subplot(222)
plt.plot(b, 'rx', label = 'b')
plt.plot(b1, 'bo', label = 'b1')
plt.ylim(-2, 2)
plt.legend()
#exercise 11.3
z = np.random.normal(loc = 0, scale = 1, size = 10000)
y = stats.norm.pdf(z)
plt.subplot(223)
plt.hist(z, bins = 25, density = True, color = 'b')
plt.subplot(224)
plt.hist(y, bins = 25, density = True, color = 'r')
plt.show()

result: