盲源分离程序同样适用于非平稳信号

function [yWtJ] = convbss（xTQKNverboseiterprinteveryetads）;
%CONVBSS demixes convolutively mixed signals.
%
% This implements an algorithm for blind source separation of
% convolutive mixtures that has been published in 
%
%    Lucas Parra Clay Spence 
%    “Convolutive blind source separation of non-stationary sources“ 
%    IEEE Trans. on Speech and Audio Processing pp.320-327 May 2000.
%
% and it has been patented （US patent US6167417）.
%
% The input arguments:
%   x       - the mixed signals （along the columns）
%   T       - the length of the FFT （default==512）
%   Q       - the length of the filter in time domain （default==128）
%   K       - the number of matrices to diagonalize （default==5）
%   N       - the number of intervals used to estimate each
%             cross-power-matrix （default==7）
%   verbose - if >0 the progress is printed out （default==1）
%   iter    - the maximum number of iterations （default==1000）
%   printevery - how often should be something printed （default==10）
%   eta     - the learning rate （default==1.0）
%   ds      - the number of sources to be estimated
%             （default==number_of_columns_of_x）
%
%Copyright notice:
%
%  author: Stefan Harmeling
%          GMD FIRST Berlin
%              Germany
%
%THIS IS UNPUBLISHED PROPRIETARY SOURCE CODE of GMD FIRST Berlin.
%  The purpose of this software is the dissemination of
%  scientific work for scientific use. The commercial
%  distribution or use of this source code is prohibited 
%  as is any modification of its content.
%The copyright notice does not evidence any actual or intended 
%publication of this code.
%
%Term and Termination:
%
%This License shall continue for as long as you use the Software. 
%However it will terminate if you fail to comply with any of its 
%terms and conditions. You agree upon termination to discontinue
%using and to destroy all copies of the Software. 
%The Limitations of Warranties and Liability set out below shall 
%continue in force even after any termination.
%
%Limitation of Warranties and Liability:
%
%THIS SOFTWARE IS PROVIDED BY THE AUTHOR “AS IS“ AND ANY EXPRESS OR
%IMPLIED WARRANTIES INCLUDING BUT NOT LIMITED TO THE IMPLIED
%WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
%HEREBY DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
%ANY DIRECT INDIRECT INCIDENTAL SPECIAL EXEMPLARY OR CONSEQUENTIAL
%DAMAGES （INCLUDING BUT NOT LIMITED TO PROCUREMENT OF SUBSTITUTE 
GOODS
%OR SERVICES; LOSS OF USE DATA OR PROFITS; OR BUSINESS INTERRUPTION）
%HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY WHETHER IN CONTRACT
%STRICT LIABILITY OR TORT （INCLUDING NEGLIGENCE OR OTHERWISE） ARISING 
IN
%ANY WAY OUT OF THE USE OF THIS SOFTWARE EVEN IF ADVISED OF THE
%POSSIBILITY OF SUCH DAMAGE.
%
%----------- 2001 GMD FIRST Berlin - Stefan Harmeling -----------
%



% the size is needed very early
[rxcx] = size（x）;

%