demix.m

卷积盲源分离

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y = demix(Wt,x)
% this is very slow, better would be something 
% like overlap-save, see Haykin or Oppenheim/Schaefer
% it is assumed that size(Wt) = [ds dx T]
% size(x) = [?? ds]

[ds dx T] = size(Wt);
[rx cx] = size(x);

% check for consistency
if ds ~= cx
error('Wt and x must fit together')
end

% allocate memory for y
y = zeros(rx,ds);

% do the convolution using matlab's built-in filter
for ii = 1:ds
for jj = 1:dx
y(:,ii) = y(:,ii) + filter(Wt(ii,jj,:),1,x(:,jj));
end
end