反向传播matlab代码

 %反向传播
        delta_sparse= beta*(-rho_ae./rho_hat+(1-rho_ae)./(1-rho_hat));    
        delta3=(layer3_dataTrans-layer1_data).* sigmoidInv(layer3_data);
        delta2=(Weight2'*delta3+delta_sparse).* sigmoidInv(layer2_data);
       
        Weight2 = Weight2 - alpha * delta3 * layer2_dataTrans' - lambda * Weight2;
        Weight1 = Weight1 -  alpha  * delta2 * layer1_data' - lambda * Weight1;            
        b1=b1- alpha *delta2;
        b2=b2- alpha *delta3;