化工过程智能建模及优化 作业

《化工过程智能建模及优化 作业》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、化工过程智能建模与优化样本训练MATLAB程序如下：（以点（X1，X2）为（-1，1）及相应的y值1为例，阈值与权值控制在0.5-3之间任意取）>>x0=1;y0=1;w110(1)=2;w111(1)=2;w112(1)=2;w120(1)=2;w121(1)=2;w122(1)=2;w210(1)=1.5;w211(1)=1.5;w212(1)=1.5;fort=1:150x1=-1;x2=1;y10=1/(1+exp(x0));y11=1/(1+exp(x1));y12=1/(1+exp(x2));z21=w110(t)*y10+w111(t)*

2、y11+w112(t)*y12;z22=w120(t)*y10+w121(t)*y11+w122(t)*y12;z20=y0;y20=1/(1+exp(y0));y21=1/(1+exp(z21));y22=1/(1+exp(z22));z31=w210(t)*y20+w211(t)*y21+w212(t)*y22;y31(t)=1/(1+exp(z31));Y=1;e(t)=Y-y31(t);dew=-e(t)*1/(exp(z31)*(1/exp(z31)+1)^2);dew210=-e(t)*1/(exp(z31)*(1/exp(z31)+1)^

3、2)*y20;dew211=-e(t)*1/(exp(z31)*(1/exp(z31)+1)^2)*y21;dew212=-e(t)*1/(exp(z31)*(1/exp(z31)+1)^2)*y22;w210(t+1)=w210(t)+5*dew210;w211(t+1)=w211(t)+5*dew211;w212(t+1)=w212(t)+5*dew212;dew110=1/(exp(z21)*(1/exp(z21)+1)^2)*y10*e(t)*1/(exp(z31)*(1/exp(z31)+1)^2)*w211(t);dew111=1/(exp

4、(z21)*(1/exp(z21)+1)^2)*y11*e(t)*1/(exp(z31)*(1/exp(z31)+1)^2)*w211(t);dew112=1/(exp(z21)*(1/exp(z21)+1)^2)*y12*e(t)*1/(exp(z31)*(1/exp(z31)+1)^2)*w211(t);dew120=1/(exp(z22)*(1/exp(z21)+1)^2)*y10*e(t)*1/(exp(z31)*(1/exp(z31)+1)^2)*w212(t);dew121=1/(exp(z22)*(1/exp(z21)+1)^2)*y11

5、*e(t)*1/(exp(z31)*(1/exp(z31)+1)^2)*w212(t);dew122=1/(exp(z22)*(1/exp(z21)+1)^2)*y12*e(t)*1/(exp(z31)*(1/exp(z31)+1)^2)*w212(t);w110(t+1)=w110(t)+5*dew110;w111(t+1)=w111(t)+5*dew111;w112(t+1)=w112(t)+5*dew112;w120(t+1)=w120(t)+5*dew120;w121(t+1)=w121(t)+5*dew121;w122(t+1)=w122(t

6、)+5*dew122;>>end此外，再把（1，1）、（-1，-1）、（1，-1）及相应的y值代入，并对四个计算得到的y值相比较。经过150次重复计算，能得到权值、阈值、相应y值的表格，经过分析后取第103个值，则权值和阈值如下表所示：W110W111W112W120W121W122W210W211W2122.2886972.7847612.2886972.2886972.7847612.288697-6.461640.1217450.121745得到的y值分别为0.849266、0.002256、0.001048、0.849266，与1、-1、-1、

7、1相近。预测程序：任意取（X1，X2），如（2，2），把阈值和权值代入程序，程序如下：>>x0=1;y0=1>>w110=2.288697;w111=2.784761;w112=2.288697;w120=2.288697;w121=2.784761;w122=2.288697;w210=-6.46164;w211=0.121745;w212=0.121745;>>x1=2;x2=2;>>y10=1/(1+exp(x0));>>y11=1/(1+exp(x1));y12=1/(1+exp(x2));>>z21=w110*y10+w111*y11+w11

8、2*y12;>>z22=w120*y10+w121*y11+w122*y12;>>z20=y0