#!/usr/local/bin/python
import numpy as np
import mympltools as util
import matpyplot
if __name__ == "__main__":
ax = matpyplot.draw('ar1_ex_em.mat')
lines = ax.get_lines()
# put better data limits
lh = lines[0]
x, y = lh.get_xdata(), lh.get_ydata()
alpha = np.diag([0.05, 0.05, 0.1, 0.1])
T = np.array([[1, -1, 0, 0], [-1, 1, 0, 0], [0, 0, 1, -1], [0, 0, -1,
1]]).dot(alpha)
ax.axis((T + np.eye(4)).dot([x.min(), x.max(), y.min(), y.max()]))
# add vertical lines
lbmax = lines[-1].get_xdata()
lines[-1].set_visible(False)
texts = (r'$\theta_{k}$', r'$\theta_{k+1}$', r'$\theta_{k+2}$')
for k in range(len(lbmax)):
ax.axvline(lbmax[k], lw=0.4, color='black', alpha=0.4)
ax.annotate(texts[k],
xy=(lbmax[k], 1),
xytext=(-5, 3),
xycoords=('data', 'axes fraction'),
textcoords='offset points',
size='small')
ax.get_figure().savefig('ar1_ex_em.pdf')
#!/usr/local/bin/python
from mpltools import style
from scipy.io import loadmat
import numpy as np
import matplotlib.pyplot as plt
import mympltools as util
import matpyplot
if __name__ == "__main__":
ax = matpyplot.draw('ballistic_trajectory.mat')
lim = ax.axis("image")
ax.set_xticks(np.arange(0, lim[1] + 1, 10))
ax.set_yticks(np.arange(0, lim[3] + 1, 10))
util.padaxis(ax, 0.03, l=lim)
#ax.grid()
ax.get_figure().savefig('ballistic_trajectory.pdf')
#!/usr/local/bin/python
from mpltools import style
from scipy.io import loadmat
import numpy as np
import matplotlib.pyplot as plt
import mympltools as util
import matpyplot
if __name__ == "__main__":
ax = matpyplot.draw('ar1_ex_b.mat')
x,y = ax.get_lines()[0].get_xdata(),ax.get_lines()[0].get_ydata()
mi = y.argmax()
ax.hold(True)
#ax.plot(x[mi], y[mi],'k*')
#xt = ax.get_xticks()
#xt.append(x[mi])
ax.axvline(x[mi],lw=0.4,color='black',alpha=1.0)
ax.annotate(r'MAP',xy=(x[mi],1),xytext=(-8,2),
xycoords=('data','axes fraction'),
textcoords='offset points',size='x-small')
# ax.annotate('ML/MAP',
# xy=(x[mi], y[mi]), xycoords='data',
# xytext=(-45,-25), textcoords='axes points',
# size='small', alpha=1.0, color='0.2',
# arrowprops=dict(arrowstyle="->",
# color="0.0",
# relpos=(0,0.5),
# shrinkA=5, shrinkB=6,
# patchA=None,
# patchB=None,
# connectionstyle="angle3,angleA=110,angleB=20",
#!/usr/local/bin/python
from mpltools import style
from scipy.io import loadmat
import numpy as np
import matplotlib.pyplot as plt
import mympltools as util
import matpyplot
if __name__ == "__main__":
ax = matpyplot.draw('ballistic_trajectory.mat')
lim = ax.axis("image")
ax.set_xticks(np.arange(0,lim[1]+1,10))
ax.set_yticks(np.arange(0,lim[3]+1,10))
util.padaxis(ax,0.03,l=lim)
#ax.grid()
ax.get_figure().savefig('ballistic_trajectory.pdf')
#!/usr/local/bin/python
import numpy as np
import mympltools as util
import matpyplot
if __name__ == "__main__":
ax = matpyplot.draw('ar1_ex_em.mat')
lines = ax.get_lines();
# put better data limits
lh = lines[0]
x,y = lh.get_xdata(),lh.get_ydata()
alpha = np.diag([0.05, 0.05, 0.1, 0.1])
T = np.array([[1,-1,0,0],[-1,1,0,0],[0,0,1,-1],[0,0,-1,1]]).dot(alpha)
ax.axis((T+np.eye(4)).dot([x.min(),x.max(),y.min(),y.max()]))
# add vertical lines
lbmax = lines[-1].get_xdata()
lines[-1].set_visible(False)
texts = (r'$\theta_{k}$',r'$\theta_{k+1}$',r'$\theta_{k+2}$')
for k in range(len(lbmax)):
ax.axvline(lbmax[k],lw=0.4,color='black',alpha=0.4)
ax.annotate(texts[k],xy=(lbmax[k],1),xytext=(-5,3),
xycoords=('data','axes fraction'),
textcoords='offset points',size='small')
ax.get_figure().savefig('ar1_ex_em.pdf')
