bsep.m

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

function [sdb] = bsep(x,y)
% bsep - Evaluates the quality of Blind Signal Separation
%
% usage: sdb = bsep(x,y);
%
% x(:,i,j) - The i^th microphone signal with only the j^th source active
% y(:,i,j) - The i^th separated signal with only the j^th source active
%      sdb - The attenuation achieved by BSS in dB
%
% See also: http://www.ele.tue.nl/ica99/
%      and: "Evaluation of Blind Signal Separation", in proc.
%            int. workshop on Independent Component Analysis and Blind
%            Signal Separation"

% Author: D.W.E. Schobben ([e-mail])
% rev. 1.02 @ 990105

if nargin ~=2 | length(size(x)) ~= 3 | length(size(y)) ~= 3,
  error('Please supply two 3D arrays')
end

J = size(x,2);

if size(x,3) ~= J | size(y,2) ~= J | size(y,3) ~= J,
  error('Incorrect array size')
end

for j=1:J,
  ymj = mean(y(:,j,j).^2);
  ymji = mean( ( sum(y(:,j,:),3) - y(:,j,j) ).^2 );
  sdb(j) = 10*log(ymj/ymji);
end