jdiv_blind_selector.m

来自「implement basic routines in Matlab and O」· M 代码 · 共 19 行

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

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


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

[x n] = max(myval);