jdiv_selector.m

implement basic routines in Matlab and Octave (currently using Matlab R13) for simulating adaptive b

function [n, myval] = jdiv_selector(b, r, s, Nmax)

if (size(b, 2) > size(r, 2))
    error('Bit-Vector and Received Vector Size Mismatch')
end

N = size(b, 1);
R = (1/size(r,2))*(r*r');
mu = 0;

for i = 1:N
    mu = mu + ( b(i)*(sum(sort(b) == b(i))/size(b, 1)) )/N;
end


w_av = calculate_av_mvdr(Nmax, s, R);
for d = 1:Nmax
    W_MVDR_AV = w_av(:, d);
    numerator = 0; denominator = 0;
    for i = 1:N
        numerator = numerator + (1/N)*b(i)*real(W_MVDR_AV'*r(:,i));
        denominator = denominator + (1/N)*(b(i)*real(W_MVDR_AV'*r(:,i)) - mu)^2;
    end
    myval(d) = 4*numerator^2/denominator;
end

[x n] = max(myval);