%盲信號分離clc;clear;%EASI參數(shù)u=6e-3;%RLS參數(shù)beta=0.99;%源信號生成，數(shù)據(jù)長度4000n=0:0.0001:0.4;S=zeros(5,4001);S(1,:)=hardlims(cos(2*pi*155*n));S(2,:)=sin(2*pi*800*n);S(3,:)=sin(2*pi*300*n+6*cos(2*pi*60*n));S(4,:)=sin(2*pi*90*n);S(5,:)=rand(1,length(n))*2-1;A=rand(5,5);X=A*S;[W_EASI ECT_EASI]=EASI(X,A,u);%EASI[W_RLS ECT_RLS]=RLS(X,A,beta,0);%普通梯度RLS[W_NRLS ECT_NRLS]=RLS(X,A,beta,1);%自然梯度RLS%畫ECT曲線figure;plot(1:4001,ECT_EASI,'r',1:4001,ECT_RLS,'g',1:4001,ECT_NRLS,'b');legend('EASI（LMS）','普通梯度RLS','自然梯度RLS');ylabel('ECT');%恢復(fù)信號Y_EASI=W_EASI*X;Y_RLS=W_RLS*X;Y_NRLS=W_NRLS*X;%分離效果對比figure;Range=1:200;for ind=1:5    subplot(4,5,ind);    plot(n(Range),S(ind,Range));    ylim([-2,2]);    title(strcat('source',num2str(ind)));endfor ind=1:5    subplot(4,5,ind+5);    plot(n(Range),Y_EASI(ind,Range));    ylim([-2,2]);    title(strcat('EASI',num2str(ind)));endfor ind=1:5    subplot(4,5,ind+10);    plot(n(Range),Y_RLS(ind,Range));    ylim([-2,2]);    title(strcat('RLS',num2str(ind)));endfor ind=1:5    subplot(4,5,ind+15);    plot(n(Range),Y_NRLS(ind,Range));    ylim([-2,2]);    title(strcat('RLS_N',num2str(ind)));end