1、clear all;xite=0.5;alfa=0.5;w1=rands(2,6); % value of w1,initially by randomw1_1=w1;w1_2=w1;w2=rands(6,1); % value of w2,initially by randomw2_1=w2;w2_2=w2_1;dw1=0*w1;x=0,0;u_1=0;y_1=0;I=0,0,0,0,0,0 % input of yinhanceng cellIout=0,0,0,0,0,0 % output of yinhanceng cellFI=0,0,0,0,0,0ts=0.001;for k=1:
2、1:1000 time(k)=k*ts; u(k)=0.5*sin(3*2*pi*k*ts); y(k)=(u_1-0.9*y_1)/(1+y_12); for j=1:6 I(j)=x*w1(:,j); Iout(j)=1/(1+exp(-I(j); end yn(k)=w2*Iout; %output of network e(k)=y(k)-yn(k); % error calculation w2=w2_1+(xite*e(k)*Iout+alfa*(w2_1-w2_2); % rectify of w2 FI(j)=exp(-I(j)/(1+exp(-I(j)2); for i=1:
3、2 dw1(i,j)=e(k)*xite*FI(j)*w2(j)*x(i); % dw1 calculation w1=w1_1+dw1+alfa*(w1_1-w1_2); % rectify of w1 % jacobian information yu=0; yu=yu+w2(j)*w1(1,j)*FI(j); dyu(k)=yu; x(1)=u(k); x(2)=y(k); w1_2=w1_1; w2_2=w2_1; u_1=u(k); y_1=y(k);endfigure(1);plot(time,y,r,time,yn,b);xlabel(timesylabel(y and ynfi
4、gure(2);plot(time,y-yn,errorfigure(3);plot(time,dyu);dyu运行结果为:(二)采用RBF网络仿真 网路结构为2-4-1,采样时间1ms,输入信号,权值的初值随机取值,高斯基函数初值,仿真m文件程序如下:%RBF simulationclear allalfa=0.05;b=1.5*ones(4,1);c=0.5*ones(2,4);w=rands(4,1);w_1=w;w_2=w_1;c_1=c;c_2=c_1;b_1=b;b_2=b_1d_w=0*w;d_b=0*b;2000; u(k)=0.5*sin(1*2*pi*k*ts); y(k)
5、=(u(k)-0.9*y_1)/(1+y_12); x(2)=y_1;4 h(j)=exp(-norm(x-c(:,j)2/(2*b(j)*b(j); ym(k)=w*h em(k)=y(k)-ym(k); d_w(j)=xite*em(k)*h(j); d_b(j)=xite*em(k)*w(j)*h(j)*(b(j)-3)*norm(x-c(:,j)2; d_c(i,j)=xite*em(k)*w(j)*h(j)*(x(i)-c(i,j)*(b(j)-2); w=w_1+d_w+alfa*(w_1-w_1); b=b_1+d_b+alfa*(b_1-b_2); c=c_1+d_c+alfa
6、*(c_1-c_2); % Jacobian information yu=yu+w(j)*h(j)*(c(1,j)-x(1)/b(j)2; w_2=w_1; w_1=w; c_2=c_1; c_1=c; b_2=b_1; b_1=b;,time,ym,time(s)y and ymplot(time,y-ym,plot(time,dyu,仿真结果为:(三)DRNN网路仿真的初值随机取值,取%DRNN simulationwd=rands(7,1);wo=rands(7,1);wi=rands(3,7);xj=zeros(7,1);xj_1=xj;xitei=0.35;xited=0.35;x
7、iteo=0.35; u(k)=0.5*sin(4*pi*k*ts); Ini=u_1,y_1,17 sj(j)=Ini*wi(:,j)+wd(j)*xj(j); xj(j)=(1-exp(-sj(j)/(1+exp(-sj(j); Pj=0*xj; Pj=wo(j)*(1+xj(j)*(1-xj(j)*xj_1(j); Qij=0*wi;3 Qij(i,j)=wo(j)*(1+xj(j)*(1-xj(j)*Ini(i); ymk=0; ymk=ymk+xj(j)*wo(j); ym(k)=ymk; e(k)=y(k)-ym(k); wo=wo+xiteo*e(k)*xj; wd=wd+xited*e(k)*Pj; wi=wi+xitei*e(k)*Qij; xj_1=xj;time/s
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1