X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=octave_packages%2Ftsa-4.2.4%2Fhisto.m;fp=octave_packages%2Ftsa-4.2.4%2Fhisto.m;h=26819da90f928a8bfbbf79872909f684dd1d192a;hb=f5f7a74bd8a4900f0b797da6783be80e11a68d86;hp=0000000000000000000000000000000000000000;hpb=1705066eceaaea976f010f669ce8e972f3734b05;p=CreaPhase.git

diff --git a/octave_packages/tsa-4.2.4/histo.m b/octave_packages/tsa-4.2.4/histo.m
new file mode 100644
index 0000000..26819da
--- /dev/null
+++ b/octave_packages/tsa-4.2.4/histo.m
@@ -0,0 +1,73 @@
+function [H,X]=histo(Y,Mode)
+% HISTO calculates histogram for each column
+% [H,X] = HISTO(Y,Mode)
+% 	   
+%   Mode 
+%	'rows' : frequency of each row
+%	'1x'   : single bin-values 
+%	'nx'   : separate bin-values for each column
+%   X  are the bin-values 
+%   H  is the frequency of occurence of value X 
+%
+% HISTO(Y) with no output arguments:
+%	plots the histogram bar(X,H)
+%
+% more histogram-based results can be obtained by HIST2RES2  
+%
+% see also:  HISTO, HISTO2, HISTO3, HISTO4
+%
+% REFERENCE(S):
+%  C.E. Shannon and W. Weaver "The mathematical theory of communication" University of Illinois Press, Urbana 1949 (reprint 1963).
+
+%	$Id: histo.m 5148 2008-06-27 10:36:37Z schloegl $
+%	Copyright (C) 1996-2002,2008 by Alois Schloegl <a.schloegl@ieee.org>	
+%    	This is part of the TSA-toolbox 
+%	http://hci.tugraz.at/~schloegl/matlab/tsa/
+%
+%    This program 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 3 of the License, or
+%    (at your option) any later version.
+%
+%    This program 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 this program.  If not, see <http://www.gnu.org/licenses/>.
+
+if nargin<2,
+        Mode='1x';
+end;
+Mode=lower(Mode);
+
+if strcmp(Mode,'rows')
+        R = histo4(Y);
+        
+elseif strcmp(Mode,'column')
+        R = histo4(Y');
+        R.X = R.X';
+        
+elseif strcmp(Mode,'1x')
+        R = histo3(Y);
+        
+elseif strcmp(Mode,'nx')
+        R = histo2(Y);
+        
+end;
+
+H = R.H;
+X = R.X;
+if nargout == 0,
+        if any(size(X)==1),
+                if exist('OCTAVE_VERSION')<5,
+                        bar(R.X,R.H,'stacked');
+                else
+                        bar(R.X,R.H);   
+                end
+        else
+                warning('2-dim X-values not supported\n')
+		%bar3(R.X,R.H);
+        end;
+end;