% Andrew Ng機(jī)器學(xué)習(xí)課上使用MATLAB ICA 實(shí)現(xiàn)的雞尾酒會(huì)語音分離
%Cock tail party problem algorithm:
%[W,s,v] = svd((repmat(sum(x.*x,1),size(x,1),1).*x)*x‘);
%audio examples
%http://research.ics.aalto.fi/ica/cocktail/cocktail_en.cgi
[x1, Fs1] = audioread(‘first.wav’);
[x2, Fs2] = audioread(‘second.wav’);

len=min(length(x1),length(x2));
xx = [x1(1:len), x2(1:len)]’;
yy = sqrtm(inv(cov(xx’)))*(xx-repmat(mean(xx,2),1,size(xx,2)));
[W,s,v] = svd((repmat(sum(yy.*yy,1),size(yy,1),1).*yy)*yy’);

a = W*xx; %W is unmixing matrix
subplot(2,2,1); plot(x1); title(‘mixed audio - mic 1’);
subplot(2,2,2); plot(x2); title(‘mixed audio - mic 2’);
subplot(2,2,3); plot(a(1,:), ‘g’); title(‘unmixed wave 1’);
subplot(2,2,4); plot(a(2,:),’r’); title(‘unmixed wave 2’);

audiowrite(‘unmixed1.wav’, a(1,:), Fs1);
audiowrite(‘unmixed2.wav’, a(2,:), Fs1);

總結(jié)

以上是生活随笔為你收集整理的MATLAB实现 ICA 鸡尾酒会语音分离的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。