python作业——Matplotlib

Exercise11.1

description:

code:

import numpy as np
import matplotlib.pyplot as plt

x = np. linspace (0, 2, 10000)
y = (np.sin(x - 2) ** 2) * np.exp(-1 * (x**2))

plt.xlim(0, 2)
plt.ylim(0, 1)
plt.xlabel('''pp's x''')
plt.ylabel('''pp's y''')
plt.title('Exercise 11.1')

plt.plot(x, y)
plt.show()

result:

Exercise11.2

description:

code:

import numpy as np
import matplotlib.pyplot as plt
import math
X = np.random.randint(10, 20, (20, 10))
b = np.random.random(10)
z = np.random.random(20)

y = np.dot(X, b) + z

b_hat = np.linalg.lstsq(X, y, rcond=None)[0]

x = list(range(1, 11))
true = plt.scatter(x, b, c='black', marker='o', label='true')
estimated = plt.scatter(x, b_hat, c='r', marker='x', label='estimated')
plt.legend()
plt.show()

result:

Exercise11.3

description:

code:

import numpy as np
import matplotlib.pyplot as plt
import scipy.stats as s

data = np.random.randn(1000)
data.sort()
kde = s.gaussian_kde(data)
plt.hist(data, bins=25, density=True, color='b')
plt.plot(data, kde.pdf(data))
plt.show()

result: