bsep.m

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function [out, w, b] = bsep( mix, ep, p, h, m);

if nargin == 4
        m = 0;
end

% make permuted timeindex
i = randperm( cols( mix)); 

% prewhiten data (ensures zero-mean and secord-order independence)
[wmix, v, mn] = prewhite( mix(:,i));

% make initial weights and bias
w = diag(ones( 1, rows( mix)), 0);
b = zeros( 1, rows( mix));

% find unmixing matrix
[w, b] = sepcore( wmix, w, b, ep, p, h, m);

% separate
out = w * mix;