matlab源音频分离的声音应用背景当2个或更多的源正在讨论的麦克风，它变得很难找到的区别的来源.zip

clc clear all; close all; [x1,fs]=wavread('input 1 ex2.wav'); [x2,fs]=wavread('input 2 ex2.wav'); x=[x1 x2]; [rx,cx] = size(x); T = 512; hT = T/2+1; %needed for fft Q = 128; K = 5; N = floor(rx/(T*K)); iter = 2000; eta = 0.95; ds = cx; dx = cx; % do the calculation disp('estimate the cross-power spectrum of x') Rx = zeros(dx,dx,hT,K); XX = zeros(dx,dx,hT); for k=1:K for n=1:N indx = (k-1)*T*N + (n-1)*T + (1:T); X=fft(x(indx,:)); X=X(1:hT,:); for row = 1:dx for col = 1:dx XX(row,col,:) = reshape(conj(X(:,col)).*X(:,row),[1 1 hT]); end end Rx(:,:,:,k) = Rx(:,:,:,k) + XX; end end Rx = Rx/(N*T); Rxk = sum(Rx,4); m = 1./repmat(sum(sum(Rxk.*conj(Rxk),2),1),[ds dx 1]); Wt = zeros(ds,dx,T); Wt(:,:,1) = eye(ds,dx); Wf = calcfft(Wt); while iter > 0 dWf = zeros(ds,dx,hT); hWf = permute(conj(Wf),[2 1 3]); for k = 1:K WfRx = multiply3(Wf,Rx(:,:,:,k)); Eomega = multiply3(WfRx,hWf); dWf = dWf + multiply3(Eomega,WfRx); end dWf = 2*dWf; dWf = m.*dWf; Wf = Wf - eta*dWf; Wt = calcifft(Wf); Wt(:,:,Q+1:end) = 0; Wf = calcfft(Wt); iter = iter - 1; end [ds dx T] = size(Wt); [rx cx] = size(x); y = zeros(rx,ds); for p = 1:ds for q = 1:dx y(:,p) = y(:,p) + filter(Wt(p,q,:),1,x(:,q)); end end fprintf('done\n'); wavwrite(y(:,1),16000,'output 1 ex2.wav'); wavwrite(y(:,2),16000,'output 2 ex2.wav');