% pltvbar.m ---------------------------------------------------
% This script demonstrates the use of vbar and ebar to plot
% vertical bar plots and error bar plots respectively.
% Some other things to note about pltvbar are:
% - The first vbar in the argument list plots two functions on a single
% trace (green) with the 1st function (phase1) defining the position
% of the bottom of the vertical bars and the 2nd function (phase2)
% defining the position of the tops of the bars.
% - The second vbar in the list plots 3 functions (called serp, bell1, and
% bell2). The 3 columns of the first vbar argument define the x coordinates
% for those three functions. The next argument (0) indicates that the bottom
% of all the vertical bars is at y=0. The last vbar argument gives the
% y coordinate for each of the 3 functions (one function per column).
% - The next pair (after the 'LineWidth') plots two traces corresponding to the
% two columns of poly23. The 1st column is a 2nd order polynomial and the
% 2nd column is 3rd order.
% - The ebar function creates two error bar traces, the 1st defined by the 1st
% column of each of the 3 arguments and the 2nd defined by the 2nd column.
% - The 'screensize' property is used to make the plot fill a fixed fraction
% of the space available.
% - The 'TraceID' argument is used to assign names for each trace that are
% appropriate for the data being displayed.
% - The 'LineWidth' argument appears in the middle of the plt call to change
% the width of only the traces defined earlier in the argument list.
% - The 'AxisPos' argument is used to widen the TraceID box to make room for
% the longer than usual trace ID names.
% - The Options' argument is used to enable the menu bar and to add the "Print"
% tag to the menu box.
% - The axis FontSize,FontWeight,FontAngle,Xtick, and Ytick properties are
% modified by including a plt argument formed by adding a "+" in front of the
% property name to indicate an axis property.
% ----- Author: ----- Paul Mennen
% ----- Email: ----- paul@mennen.org
x = (0:.25:6)';
serp = 19 * x ./ (x.*x+1);
bell1 = 12 * exp(-(x-2.5).^2);
bell2 = 12 * exp(-(x-4).^2);
t = (0:.01:1)';
phase1 = 1.4 * exp(-2*t) .* sin(20*t) + 20;
phase2 = 3 * t .* cos(5*pi*(1-t).^3) + 16;
p2 = polyval([-0.7 4 8],x);
p3 = polyval([-0.35 3.5 -9 28],x);
poly23 = [p2 p3]; % [quadradic cubic]
errval = [0*x+1 2*cos(7+x/2)];
Id = {'phase 1&2' 'serpentine' 'bell curve 1' 'bell curve 2' ...
'quadradic' 'cubic' 'err-constant', 'err-cosine'}; % trace IDs
set(0,'Units','pixels'); sz = get(0,'screensize'); % get screen size
sfrac = [.7 .6] .* sz(3:4); % use 70% of the screen width and 60% of the height
plt('TraceID',Id,'Position',[10 50 sfrac],'AxisPos',[1 1 1 1 1.4],...
'FigName','pltvbar','FigBKc',[.1 .1 .2],'xlim',[-.1 6.2],...
vbar(6*t,phase1,phase2),...
vbar([x x+.05 x+.1],0,[serp bell1 bell2]),...
'LineWidth',3,'+Xtick',0:2:6,'+Ytick',[12 13.7 14.7 21.2 26.8],...
x,poly23,ebar([x x],poly23,errval),'Options','Menu +Print',...
'+FontSize',11,'+FontWeight','bold','+FontAngle','Italic');
% A figure somewhat similar to what is created by the above plt
% command could also be generated using the standard MatLab plot
% function (not that I would really recommend this) as follows:
% --------------------------------------------------------------
% hold on;
% h1 = plot(vbar(6*t,phase1,phase2),'LineWidth',3);
% h2 = plot(vbar([x x+.05 x+.1],0,[serp bell1 bell2]),'LineWidth',3);
% h3 = plot(x,poly23);
% h4 = plot(ebar([x x],poly23,errval));
% set([h1;h2;h3;h4],{'color'},num2cell([get(gca,'ColorOrder'); .7 0 0],2));
% set(gcf,'Name','pltvbar','Position',[10 50 1000 600]);
% set(gca,'xlim',[-.1 6.2]);
% legend(Id{:},-1);
% --------------------------------------------------------------
