#!/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')