count.m

MIMO and VBLAST( please see the documentation)

function varargout = count(data_src,data_est,varargin)
%COUNT Error counter.
%   ERR = COUNT(DATA_SRC,DATA_EST) simply counts an error occurrence
%   of data estimations.
%
%   [ERR,SER] = COUNT(DATA_SRC,DATA_EST) same as above, but also gives
%   a symbol error rate in decibels.
%
%   [ERR,SER] = COUNT(DATA_SRC,DATA_EST,'EchoOn') displays a counting
%   error process as full-report in table form.
%
%   See also DETECT, SOURCE.

%   Copyright 2001-2002 Kamil Anis, [e-mail]
%   Dept. of Radioelectronics, 
%   Faculty of Electrical Engineering
%   Czech Technical University in Pragu
%   $Revision: 2.0 $  $Date: 2002/10/23 17:33:28 $
%   --
%   <additional stuff should go here>


[fr_length,space_dim,frames] = size(data_est);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% BODY BEGIN %%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

err = sum(sum(data_src ~= data_est));

if err ~= 0
	ser = err / (frames * fr_length);
  ser = 10 * log10(ser);
  rel = 100 - ser * 100;
else
  ser = -inf;
  rel = 100;
end

varargout = {err,ser};

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% BODY END %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if (isempty(varargin) == 0) & (varargin{end} == 'EchoOn')
	[indent,gap,name] =  iprompt('COUNT:');
	
	str1 = sprintf('%1d',err);
	str2 = sprintf('%1.0f',ser);
	str3 = sprintf('%1.0f',rel);
	
	disp(' ');
	disp([name,gap,'Total errors counted -> ',str1,'.']);
	disp([indent,'Symbol error rate -> ',str2,' [dB].']);
	disp([indent,'System reliability -> ',str3,' %.']);
	disp(' ');
end