Deep Semi-NMF源代码

function [ Z H dnorm ] = deep_seminmf （ X layers varargin ）

% Process optional arguments
pnames = { ...
    ‘z0‘ ‘h0‘ ‘bUpdateH‘ ‘bUpdateLastH‘ ‘maxiter‘ ‘TolFun‘ ...
    ‘verbose‘ ‘bUpdateZ‘ ‘cache‘ ‘gnd‘ ...
};


X = bsxfun（@rdivideXsqrt（sum（X.^21）））;

num_of_layers = numel（layers）;

Z = cell（1 num_of_layers）;
H = cell（1 num_of_layers）;

dflts  = {0 0 1 1 500 1e-5 1 1 1 0};

[z0 h0 bUpdateH bUpdateLastH maxiter tolfun verbose bUpdateZ cache gnd] = ...
        internal.stats.parseArgs（pnamesdfltsvarargin{:}）;

if  ~iscell（h0）
    for i_layer = 1:length（layers）
        if i_layer == 1
            % For the first layer we go linear from X to Z*H so we use id
            V = X;
        else 
            V = H{i_layer-1};
        end
        
        if verbose
            display（sprintf（‘Initialising layer #%d with k=%d with size（V）=%s...‘ i_layer layers（i_layer） mat2str（size（V））））;
        end
        if ~iscell（z0）
            % For the later layers we use nonlinearities as we go from
            % g（H_{k-1}） to Z*H_k
            [Z{i_layer} H{i_layer} ~] = ...
                 seminmf（V ...
                     layers（i_layer） ...
                     ‘maxiter‘ maxiter ...
                     ‘bUpdateH‘ true ‘bUpdateZ‘ bUpdateZ ‘verbose‘ verbose ‘save‘ cache ‘fast‘ 1）; 
        else
            display（‘Using existing Z‘）;
            [Z{i_layer} H{i_layer} ~] = ...
                 seminmf（V ...
                     layers（i_layer） ...
                     ‘maxiter‘ 1 ...
                     ‘bUpdateH‘ true ‘bUpdateZ‘ 0 ‘z0‘ z0{i_layer} ‘verbose‘ verbose ‘save‘ cache ‘fast‘ 1）; 
        end
    end

else
    Z=z0;
    H=h0;
    
    if verbose
        display（‘Skipping initialization using provided init matrices...‘）;
    end
end


dnorm0 = cost_function（X Z H）;
dnorm = dnorm0 + 1;

if verbose
    display（sprintf（‘#%d error: %f‘ 0 dnorm0））;
end

%% Error Propagation
if verbose
    display（‘Finetuning...‘）;
end
H_err = cell（1 num_of_layers）;


for iter = 1:maxiter  
    H_err{numel（layers）} = H{n