cal_parameter_back.m

加了动量项和域值

function [WMI_Xin,WIJ_Xin,WJP_Xin,D_WMI,D_WIJ,D_WJP]=Cal_Parameter_Back(WMI_Yuan,WIJ_Yuan,WJP_Yuan,Chuan_Can,VI,VJ,y,X,T,D_WMI,D_WIJ,D_WJP)
M=Chuan_Can(1);%输入点个数
I=Chuan_Can(2);%第I曾神经元个数
J=Chuan_Can(3);%第J层神经元个数
P=Chuan_Can(4);%输出点个数
S_Rate_In=Chuan_Can(5);
S_Rate_Yin=Chuan_Can(6);
S_Rate_Out=Chuan_Can(7);
aa=Chuan_Can(8);
%T=[0.3771 0.60405 0.6600 0.4609 0.1336];
%X=[-0.7 -0.6 -0.5 -0.4 -0.3];
%修正WJP
for p=1:P
    DP(p)=y(p)*(1-y(p))*(T(p)-y(p));
    for j=1:J
        WJP_Xin(j,p)=WJP_Yuan(j,p)+aa*D_WJP(j,p)+S_Rate_Out*DP(p)*VJ(j);
    end
end
D_WJP=WJP_Xin-WJP_Yuan;

%修正WIJ
for j=1:J
    DJ(j)=0;
    for p=1:P
        DJ(j)=DJ(j)+VJ(j)*(1-VJ(j))*DP(p)*WJP_Yuan(j,p);    
    end
end

for j=1:J
    for i=1:I
        WIJ_Xin(i,j)=WIJ_Yuan(i,j)+aa*D_WIJ(i,j)+S_Rate_Yin*DJ(j)*VI(i);
    end
end
D_WIJ=WIJ_Xin-WIJ_Yuan;

%修正WMI
for i=1:I
    DI(i)=0;
    for j=1:J
        DI(i)=DI(i)+VI(i)*(1-VI(i))*DJ(j)*WIJ_Yuan(i,j);    
    end
end

for i=1:I
    for m=1:M
        WMI_Xin(m,i)=WMI_Yuan(m,i)+aa*D_WMI(m,i)+S_Rate_In*DI(i)*X(m);
    end
end
D_WMI=WMI_Xin-WMI_Yuan;