graphkla.m.svn-base

Bayesian Surprise Toolkit for Matlab T. Nathan Mundhenk, Laurent Itti

% //////////////////////////////////////////////////////////////////// %
% The Baysian Surprise Matlab Toolkit - Copyright (C) 2004-2007        %
% by the University of Southern California (USC) and the iLab at USC.  %
% See http://iLab.usc.edu for information about this project.          %
% //////////////////////////////////////////////////////////////////// %
% This file is part of the Baysian Surprise Matlab Toolkit             %
%                                                                      %
% The Baysian Surprise Matlab Toolkit is free software; you can        %
% redistribute it and/or modify it under the terms of the GNU General  %
% Public License as published by the Free Software Foundation; either  %
% version 2 of the License, or (at your option) any later version.     %
%                                                                      %
% The Baysian Surprise Matlab Toolkit is distributed in the hope       %
% that it will be useful, but WITHOUT ANY WARRANTY; without even the   %
% implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR      %
% PURPOSE.  See the GNU General Public License for more details.       %
%                                                                      %
% You should have received a copy of the GNU General Public License    %
% along with the iBaysian Surprise Matlab Toolkit; if not, write       %
% to the Free Software Foundation, Inc., 59 Temple Place, Suite 330,   %
% Boston, MA 02111-1307 USA.                                           %
% //////////////////////////////////////////////////////////////////// %
%
% Primary maintainer for this file: T. Nathan Mundhenk <[e-mail]>
%

function Z = graphkla(a1,a2,interval,grain)

matsize = interval/grain;

Z = zeros(matsize);
fprintf('Computing gamma kl values\n');
x = 1;
for b1 = 0.1: grain: interval
    fprintf('.');
    y = 1;
    for b2 = 0.1: grain: interval
        Z(x,y) = klgamma(a1,a2,b1,b2);
        y = y + 1;
    end
    x = x + 1;
end

fprintf('Done \n');

h = surfl(Z);
figure(gcf); % subplot(2,1,2);
colormap hsv;
shading interp;
set(h,'EdgeColor','k');
%light('Position',[50,50,50]);
%lighting phong;
%material([0.4,0.6,0.5,30]);
%set(h,'FaceColor',[0.7 0.7 0],'BackFaceLighting','lit')
%view([30,25]);
%set(gca,'CameraViewAngleMode','Manual');
axis([0 matsize 0 matsize 0 100])
%set(gca,'ZTickLabel','Negative||Positive');
%set(gca,'PlotBoxAspectRatio',[2.5 2.5 1]);
xlabel('b Axis');
ylabel('b` Axis');
zlabel('Surprise Value');
title('Surprise by beta values');
set(gcf,'PaperPositionMode','auto');
print -dps2;