import numpy
import matplotlib.pyplot as plt
x = numpy.linspace(0, 2, 20)
y = numpy.power(numpy.sin(x-2), 2) * numpy.exp(-1*(x**2))
plt.plot(x, y, 'r-')
plt.xlabel("x label")
plt.ylabel("y label")
plt.title("EX11.1")
plt.show()
import numpy
import matplotlib.pyplot as plt
import scipy.optimize
X = numpy.random.rand(20, 10)
b = numpy.random.rand(10, 1)
z = numpy.random.normal(loc=0, scale=1.0, size=(20, 1))
y = numpy.dot(X, b) + z
B = numpy.dot(numpy.dot(numpy.linalg.inv(numpy.dot(X.T, X)), X.T), y)
plt.title('EX11.2')
plt.xlabel('index')
plt.ylabel('value')
x = numpy.arange(0,10)
plt.scatter(x, B, c='b', marker='o', label='$Estimated coefficients$')
plt.scatter(x, b, c='r', marker='x', label='$true coefficients$')
plt.legend()
plt.show ()
import numpy
import matplotlib.pyplot as plt
import scipy.stats
z = numpy.random.normal(loc=0, scale=5.0, size=10000)
plt.hist(z , bins=25, density = True, color='b')
kernel = scipy.stats.gaussian_kde(z)
x = numpy.linspace(-20, 20, 10000)
plt.title('EX11.3')
plt.plot(x, kernel(x), 'k')
plt.show()