user_alg2.m

来自「Matlab盲信号分离很有用的工具」· M 代码 · 共 32 行

function W = user_alg2( x );
% ICALAB user can replace this function by own algotrithm computing demixing matrix W
%
% This is only dummy function
%
%   input:
%    x - signals vector, each signal is in different row
%
%   output:
%    W - separation matrix
%        if user algorithm use prewhitening, matrix W = W * Q

fprintf( '\nYou can insert here your own algorithm\n\n' );

M = 20;

W = eye( size(x,1) );

N = size(x, 1);

for i = 1:N
   h(i,:) = aryule(x(i,:), M);
end

h = mean(h);

for i = 1:N
  tX1(i,:) = filter(h, 1, x(i,:));
end

W = jadeop(tX1);