matplotlib作业

import numpy as np
import matplotlib.pyplot as plt
import math


x = np.linspace(0, 2, 1000)
y = (math.e**(-x*x))*np.sin(x-2)*np.sin(x-2)
plt.plot(x, y)
plt.axis([0, 2, -0.5, 1.5])
plt.title('Plotting a function')
plt.xlabel('x')
plt.ylabel('f(x)')
plt.show()

import numpy as np
import matplotlib.pyplot as plt
import math

X = np.random.randint(0, 9, (20, 10))
z = np.random.randint(1, 10, (20, 1))
b0 = np.random.random_sample((10, 1)) * 3 - 1.5
y = np .dot(X, b0) + z
b1 = np.linalg.lstsq(X, y, rcond=None)[0]

x = np.linspace(0, 9, 10)
p1 = plt.scatter(x, b0, marker = 'x', c = 'r')
p2 = plt.scatter(x, b1, marker = 'o', c = 'b')
plt.axis([0, 9, -3, 3])
plt.title('Parameter plot')
plt.xlabel('index')
plt.ylabel('value')
plt.show()

import numpy as np
import matplotlib.pyplot as plt
import math
import scipy.stats as cstats


z = sorted(np.random.normal(1, 1, size = (10000)))
plt.hist(z, bins = 25, density = True)
kernel = cstats.gaussian_kde(z)
plt.plot(z, kernel.pdf(z))
plt.show()