calculate_lms_mvdr.m

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

% Frost LMS Algorithm - From P.810 of Optimum Array Processing

function w_lms = calculate_lms_mvdr(s, r, mu)

numIterations = size(r, 2);
N = size(r,1);

w_lms = zeros(N, numIterations);

w_lms(:,1) = s/(norm(s)^2);
tempResult = (eye(N)-((s*s')/(norm(s)^2)));

for i = 2:numIterations
    w_lms(:,i) = tempResult*(w_lms(:,i-1) - mu*r(:, i)*r(:, i)'*w_lms(:,i-1)) + s/norm(s);
end