newneg.m

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

function y=newNeg()
%
%
%
%
s1=wavread('ica1');
s2=wavread('ica2');
s3=wavread('ica3');
s4=wavread('ica4');
s1=s1';
s2=s2';
s3=s3';
s4=s4';
s=[s1;s2;s3;s4];
A=[0.7486 0.5624 0.3829 0.4798;0.3741 0.3723 0.2528 0.3683;0.4541 0.7928 0.3429 0.7646;0.0386 0.7952 0.9678 0.3771];
x=A*s;
x=baihua(x);
x1=x(1,:);
x2=x(2,:);
x3=x(3,:);
y=zeros(size(x));
u=0.025;
w=rand(4);
for i=1:length(x)
    y(:,i)=w*x(:,i);
    w=w-u*x(:,i)*tanh(y(:,i)');
    w=w/norm(w);
    p=w*A;
    a=0;
    c=0;
    for k=1:length(p)
        for j=1:length(p)
            a=a+abs(p(k,j))/max(abs(p(k,:)));
            c=c+abs(p(j,k))/max(abs(p(:,k)));
        end
    end
    err(i)=a+c;
end
plot(err);