1、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:1:1000 tim
2、e（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:2 dw1（i,j）
3、=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 ynfigure（2）;pl
4、ot（time,y-yn,errorfigure（3）;plot（time,dyu）;dyu运行结果为：（2）采用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）=（u（k）-0.9
5、*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*（c_1-c_2）
6、; % 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,仿真结果为：（3）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;xiteo=0.35;
7、 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
