第十二周作业(Matplotlib专题)

 Exercise 11.1: Plotting a function

代码

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

x=np.linspace(0,2,1000)
y=np.power(np.sin(x-2),2)*np.exp((-1)*np.power(x,2))
fig, ax=plt.subplots()
ax.plot(x,y)
ax.set_xlabel('x')
ax.set_ylabel('y')
ax.set_xlim((0,2))
ax.set_ylim((0,1))
ax.set_title('Exercise 11.1:')
plt.show()

运行结果

 Exercise 11.2: Data

代码

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

X=np.random.rand(20,10)
b=np.random.rand(10,1)
z=np.random.normal(0,1,(20,1))
y=np.dot(X,b)+z

D=np.dot(np.transpose(X),X)
d=np.dot(np.transpose(X),y)
b_=np.dot(np.linalg.inv(D),d)

a=np.linspace(0,9,10)
fig, ax=plt.subplots()
ax.grid(linestyle='-')
True_C,=ax.plot(a,b,'bo')
Est_C,=ax.plot(a,b_,'r+')
plt.legend([True_C,Est_C],['True cofficients','Estimated cofficients'])
ax.set_xlabel('index')
ax.set_ylabel('value')
ax.set_xlim((0,10))
ax.set_ylim((-2,2))
ax.set_title('Exercise 11.2:')
plt.show()

运行结果

 Exercise 11.3: Histogram and density estimation

代码

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

z=np.random.randn(10000)
figure,ax=plt.subplots()
c,d,e=ax.hist(z, bins=25,color='b',density=True)
ax.set_title('Exercise 11.3:')
plt.setp(e, edgecolor='k')
kernel=stats.gaussian_kde(z)
x=np.linspace(-4,4,1000)
y=kernel.pdf(x)
ax.plot(x,y,'k-')
plt.show()

运行结果