Mercurial > forge
view main/comm/inst/biterr.m @ 9666:67d4cfc5eeb3 octave-forge
comm: update license to GPLv3+
author | carandraug |
---|---|
date | Tue, 13 Mar 2012 04:31:21 +0000 |
parents | b81e11a89ae2 |
children | fe6da9839797 |
line wrap: on
line source
## Copyright (C) 2003 David Bateman ## ## 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/>. ## -*- texinfo -*- ## @deftypefn {Function File} {[@var{num}, @var{rate}] = } biterr (@var{a},@var{b}) ## @deftypefnx {Function File} {[@var{num}, @var{rate}] = } biterr (@var{...},@var{k}) ## @deftypefnx {Function File} {[@var{num}, @var{rate}] = } biterr (@var{...},@var{flag}) ## @deftypefnx {Function File} {[@var{num}, @var{rate} @var{ind}] = } biterr (@var{...}) ## ## Compares two matrices and returns the number of bit errors and the bit ## error rate. The binary representations of the variables @var{a} and ## @var{b} are treated and @var{a} and @var{b} can be either: ## ## @table @asis ## @item Both matrices ## In this case both matrices must be the same size and then by default the ## the return values @var{num} and @var{rate} are the overall number of bit ## errors and the overall bit error rate. ## @item One column vector ## In this case the column vector is used for bit error comparision column-wise ## with the matrix. The returned values @var{num} and @var{rate} are then ## row vectors containing the num of bit errors and the bit error rate for ## each of the column-wise comparisons. The number of rows in the matrix ## must be the same as the length of the column vector ## @item One row vector ## In this case the row vector is used for bit error comparision row-wise ## with the matrix. The returned values @var{num} and @var{rate} are then ## column vectors containing the num of bit errors and the bit error rate for ## each of the row-wise comparisons. The number of columns in the matrix ## must be the same as the length of the row vector ## @end table ## ## This behaviour can be overridden with the variable @var{flag}. @var{flag} ## can take the value 'column-wise', 'row-wise' or 'overall'. A column-wise ## comparision is not possible with a row vector and visa-versa. ## ## By default the number of bits in each symbol is assumed to be give by the ## number required to represent the maximum value of @var{a} and @var{b}. ## The number of bits to represent a symbol can be overridden by the variable ## @var{k}. ## @end deftypefn ## 2003 FEB 13 ## initial release function [num, rate, ind] = biterr (a, b, varargin) if ((nargin < 2) || (nargin > 4)) usage ("[num rate ind] = biterr (a, b [,k [,flag]])"); endif if (ndims (a) > 2 || ndims (b) > 2) error ("biterr: a and b must have at most two dimensions"); endif if (any(any(isinf(a))) || any(any(isnan(a))) || any(any(isinf(b))) || ... any(any(isnan(b))) || !isreal(a) || !isreal(b) || ... any(any((floor(a)) != a)) || any(any((floor(b)) != b)) || ... any(any(a < 0)) || any(any(b < 0))) error ("biterr: a and b must contain only non-negative integers"); endif [ar,ac] = size(a); [br,bc] = size(b); k = max([max(a(:)),max(b(:))]); m = 1; while (k > (2^m-1)) m = m + 1; end if ((ar == br) && (ac == bc)) type = "matrix"; flag = "overall"; c = 1; elseif (any([ar,br] == 1)) type = "row"; flag = "row"; if (ac != bc) error ("biterr: row-wise comparison must have the same number of columns in inputs"); endif if (ar == 1) a = ones(br,1) * a; else b = ones(ar,1) * b; endif elseif (any([ac,bc] == 1)) type = "column"; flag = "column"; if (ar != br) error ("biterr: column-wise comparison must have the same number of rows in inputs"); endif if (ac == 1) a = a * ones(1,bc); else b = b * ones(1,ac); endif else error ("biterr: matrix sizes must match"); endif k = 0; for i =1:length(varargin) arg = varargin{i}; if (ischar(arg)) if (strcmp(arg,"row-wise")) if (strcmp(type,"column")) error ("biterr: row-wise comparison not possible with column inputs"); endif flag = "row"; elseif (strcmp(arg,"column-wise")) if (strcmp(type,"row")) error ("biterr: column-wise comparison not possible with row inputs"); endif flag = "column"; elseif (strcmp(arg,"overall")) flag = "overall"; else error ("biterr: unrecognized string argument"); endif else k = arg; if (k < m) error ("biterr: the symbol size is too small for largest element"); endif endif end if (k == 0) k = m; endif ## Call the core error function to count the bit errors ind = __errcore__(a,b); switch (flag) case 'row', if (strcmp(type,"matrix") && (ac == 1)) num = ind; else num = sum(ind')'; endif rate = num / k / max(ac,bc); case 'column', if (strcmp(type,"matrix") && (ar == 1)) num = ind; else num = sum(ind); endif rate = num / k / max(ar,br); case 'overall', num = sum(sum(ind)); rate = num / k / max(ar,br) / max(ac,bc); otherwise error("impossible"); endswitch endfunction