yichuanfahua.m

来自「此包含有遗传算法、粒子群算法、BP算法优化对角递归神经网络的MATLAB程序」· M 代码 · 共 24 行

function [e,y]=yichuanfahua(p,o,endz_indiv)
t=8.*ones(1,8);
 t=cumsum(t);
 t=[0 t];
 xj=zeros(8,1);
 xj_1=xj;
 e=zeros(1,86);
for k=1:86
    Ini = p(k,:)';
    yo=0; %放在这里，为了下面命令yo不会把所有样本的输出全部叠加
     for j=1:8 %隐层节点数
        sj(j)=endz_indiv(1,t(j)+1:t(j+1))*Ini+endz_indiv(1,(64+j))*xj_1(j);
        xj(j)=(1-exp(-sj(j)))/(1+exp(-sj(j)));
        yo=yo+xj(j)*endz_indiv(1,72+j); 
        xj_1(j)=xj(j);
     end  
     y(k)=yo;
     e(k) =  o(k) - y(k); 
end



%[y1]=postmnmx(y,2330.7,2636.8);[o1]=postmnmx(o,2330.7,2636.8);e1=sum(e(1,:).*e(1,:))/86;
%plot([1:86],o1,[1:86],y1,'-r'),axis([1 86 2300 2700]);xlabel('测试样本'); ylabel('丙烯/ppm');legend('实际曲线','遗传算法');