利用matlab实现语音盲分离

本文利用了独立变量分析的算法，用 matlab实现了语音信号的盲分离。这在语音识别，以及未来机器人智能化上起着至关重要的作用

function [W Wefica Wwasobi ISRwa ISRef metoda]= combi（XAR_orderortho）
%
% combines EFICA WASOBI
%

if nargin<3
    ortho=false;
end

%COMMON PREPROCESSING
X = X-mean（X2）*ones（1size（X2））;
C = cov（X‘）;
CC = C^（-1/2）;
x = CC*X;


[Wefica ISRef]=efica（x）;
[Wwasobi WsobiISRwa]= ewasobi（x100.9）;
ISR1=sum（ISRef2）;
ISR2=sum（ISRwa2）;
imin=min（[ISR1 ISR2]）;
if imin（1）    hotove=find（ISR1    ostatni=find（ISR1>=imin（2））;
    W=Wefica（hotove:）;
    if ortho
        Wost=null（Wefica（hotove:））‘;
    else
        Wost=Wefica（ostatni:）;
    end
    y=Wost*x;
    metoda=ones（size（hotove‘））;
else
    hotove=find（ISR2    ostatni=find（ISR2>=imin（1））;
    W=Wwasobi（hotove:）;
    if ortho
        Wost=null（Wwasobi（hotove:））‘;
    else
        Wost=Wwasobi（ostatni:）;
    end
    y=Wost*x;
    metoda=2*ones（size（hotove‘））;
end

while length（ostatni）>1
    [Wefi ISRef0]=efica（y）;
    [WwaAOL_initISRwa0]= ewasobi（y100.9）;
    ISR1=sum（ISRef02）;
    ISR2=sum（ISRwa02）;
    imin=min（[ISR1 ISR2]）;
    if imin（1）        hotove=find（ISR1        ostatni=find（ISR1>=imin（2））;
        W=[W; Wefi（hotove:）*Wost];
        if ortho
            Wost=null（Wefi（hotove:））‘*Wost;
        else
            Wost=Wefi（ostatni:）*Wost;
        end
        y=Wost*x;
        metoda=[metoda ones（size（hotove‘））];
    else
        hotove=find（ISR2        ostatni=find（ISR2>=imin（1））;
        W=[W; Wwa（hotove:）*Wost];
        if ortho
            Wost=null（Wwa（hotove:））‘*Wost;
        else
            Wost=Wwa（ostatni:）*Wost;
        end
        y=Wost*x;
        metoda=[metoda 2*ones（size（hotove‘））];
    end
end
if ~isempty（ostatni）
    W=[W; Wost];
    metoda=[metoda 0];
end

W=W*CC;
Wefica=Wefica*CC;
Wwasobi=Wwasobi*CC;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      10933  2007-03-29 02:55  icalab\acrsobibpf.p

     文件       7603  2007-03-29 02:55  icalab\acsobiro.p

     文件       7166  2007-03-29 02:55  icalab\adv_opt_matrixH.p

     文件      37194  2007-03-29 03:05  icalab\adv_param.p

     文件       2014  2007-03-29 02:55  icalab\als_reg.p

     文件       3344  2007-03-29 02:55  icalab\amuse.p

     文件       1327  2007-03-29 02:55  icalab\amusepermuta.p

     文件     480192  2002-08-15 19:09  icalab\benchmarks\10halo.mat

     文件    1200184  2002-02-18 19:38  icalab\benchmarks\25speakersHALO.mat

     文件    1200184  2002-02-18 17:14  icalab\benchmarks\25speakersNOSYNC.mat

     文件    1792184  2002-02-18 19:38  icalab\benchmarks\64sounds_std.mat

     文件     179799  2007-03-28 11:40  icalab\benchmarks\ABio5.mat

     文件     219141  2007-03-28 11:40  icalab\benchmarks\ABio6.mat

     文件     280200  2002-08-15 19:09  icalab\benchmarks\ABio7.mat

     文件      80208  2002-08-15 19:09  icalab\benchmarks\AC10-7sparse.mat

     文件     160240  2002-11-01 13:33  icalab\benchmarks\ACPos24sparse10.mat

     文件     160248  2002-11-01 13:27  icalab\benchmarks\ACposin10.mat

     文件      80184  2002-11-01 13:31  icalab\benchmarks\ACposvsparse.mat

     文件      80184  2003-03-19 16:02  icalab\benchmarks\ACpwcon10.mat

     文件      80184  2002-08-15 19:09  icalab\benchmarks\ACsin10d.mat

     文件      32216  2002-08-15 19:09  icalab\benchmarks\ACsin4d.mat

     文件      80184  2002-11-01 13:14  icalab\benchmarks\ACsincpos10.mat

     文件      80184  2002-08-15 19:09  icalab\benchmarks\ACsparse10.mat

     文件     448184  2002-08-15 19:09  icalab\benchmarks\acspeech16.mat

     文件      80184  2003-02-14 17:14  icalab\benchmarks\acv10_sin.mat

     文件      80184  2002-08-15 19:09  icalab\benchmarks\ACvsparse10.mat

     文件      17196  2005-12-30 11:40  icalab\benchmarks\AC_2noi_2smooth.mat

     文件      32184  2003-02-14 17:14  icalab\benchmarks\bbb4.mat

     文件     181003  2007-03-28 22:55  icalab\benchmarks\c5signals.mat

     文件      32184  2003-02-14 17:14  icalab\benchmarks\ccc4.mat

............此处省略237个文件信息