ica.m

来自「这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的」· M 代码 · 共 17 行

function PI=ica_demo(s)
%205.11.28
tic
A=[1.5 1;1 1.4];
x=A*s;
[N,t]=size(x);
for i=1:1000
%Q=((ones(2,1)*sum(abs(x).^2)).*x)*x'-x*x'*conj(x*x')/t-x*x'*x*x'/t;
Q=x(:,1:t-i)*x(:,i+1:t)'+x(:,i+1:t)*x(:,1:t-i)';
%Q=eye(N);
%Q=x*x';
[W,D]=eig(x*x',Q);
y=W'*x;
r=xsxs(s,y);
PI(i)=xnzb(r);
end
toc