第十二周exercise

import numpy as np

import matplotlib.pyplot as plt



x = np.linspace(0,2,200)
y = np.power(np.sin(x-2),2)*np.exp(-np.power(x,2))

plt.plot(x,y)
plt.xlabel('x')
plt.ylabel('y')
plt.title('f(x)=(sin(x-2)^2)*(e^(-x^2))')
plt.show()

import numpy as np

import matplotlib.pyplot as plt

from scipy import linalg

X = np.random.randint(1,10,size=(10,20))
b =np.random.normal(1,2,20).reshape([20,1])
z =np.random.normal(1,2,10).reshape([10,1])
y = np.dot(X,b)+z

b_ = np.dot(np.dot(linalg.inv(np.dot(X.T,X)),X.T),y)

x = np.linspace(-10,10,20)
plt.scatter(x,b)

plt.scatter(x,b_,c='r')
plt.show()

import numpy as np

import matplotlib.pyplot as plt

from scipy import stats

x = np.arange(-5,5,0.001)
y = stats.norm.pdf(x,0,1)
data = stats.norm.rvs(loc=0,scale=1,size=1000)

plt.plot(x,y,color='red')
plt.hist(data,bins=20,normed=1,color='blue',alpha=0.7)

plt.show()