ex-6: add EMD vs RL rounds plots

This commit is contained in:
Michele Guerini Rocco 2020-05-19 16:56:26 +00:00
parent 62f03170af
commit 7b8c61e3e5
5 changed files with 11 additions and 8 deletions

View File

@ -3,24 +3,26 @@ import matplotlib.pyplot as plt
def plot(table, title='', log=False): def plot(table, title='', log=False):
plt.figure() plt.figure(figsize=(6.3, 6))
plt.suptitle(title) plt.suptitle(title)
plt.subplot(211) plt.subplot(211)
if log: if log:
plt.xscale('log') plt.xscale('log')
plt.title('EMD distance' + ' '*10, loc='right') plt.title('EMD' + ' '*10, loc='right')
plt.plot(table[0], table[1], color='#92182b') plt.plot(table[0], table[1], color='#92182b')
plt.tick_params(axis='y', labelcolor='#92182b') plt.tick_params(axis='y', labelcolor='#92182b')
plt.ylabel('average', color='#92182b') plt.ylabel('average', color='#92182b')
plt.ticklabel_format(style='sci', axis='y', scilimits=(0, 0)) plt.ticklabel_format(style='sci', axis='y',
scilimits=(0, 0), useMathText=True)
twin = plt.twinx() twin = plt.twinx()
twin.plot(table[0], table[2], color='gray') twin.plot(table[0], table[2], color='gray')
twin.tick_params(axis='y', labelcolor='gray') twin.tick_params(axis='y', labelcolor='gray')
twin.set_ylabel('standard deviation', color='gray') twin.set_ylabel('standard deviation', color='gray')
twin.ticklabel_format(style='sci', axis='y', scilimits=(0, 0)) twin.ticklabel_format(style='sci', axis='y',
scilimits=(0, 0), useMathText=True)
plt.subplot(212) plt.subplot(212)
if log: if log:
@ -29,12 +31,13 @@ def plot(table, title='', log=False):
plt.xlabel('RL rounds') plt.xlabel('RL rounds')
plt.plot(table[0], table[3], color='xkcd:gray') plt.plot(table[0], table[3], color='xkcd:gray')
plt.tight_layout()
plt.rcParams['font.size'] = 12
table = np.loadtxt('ex-6/plots/emd-round.txt') table = np.loadtxt('ex-6/plots/emd-round.txt')
plot(table[:27].T, title='noiseless', log=True) plot(table[:27].T, title='noiseless', log=True)
plot(table[27:47].T, title=r'$\sigma=0.005$') plot(table[27:47].T, title=r'noise at $\sigma=0.005$')
plot(table[47:67].T, title=r'$\sigma=0.01$') plot(table[47:67].T, title=r'noise at $\sigma=0.01$')
plot(table[67:].T, title=r'$\sigma=0.05$') plot(table[67:].T, title=r'noise at $\sigma=0.05$')
plt.show() plt.show()

Binary file not shown.

Binary file not shown.

Binary file not shown.

Binary file not shown.