clear?all
load?rr_inteval
load?avg_rr_inteval;
load?avg_rscope;
P=3;m=1;
n=10;
N=1;
x=[avg_rscope;avg_rr_inteval;rr_inteval（4）];
d=[0.003;0.5;1.00];
W=rand（Nn）;
WW=rand（nm）;
a=ones（1n）;
for?j=1:n
????b（j）=j*P/n;
end
epoch=1;
epo=2000;
error=0.1;
err=0.01;
delta=1;
lin=0.5;
while?（error>=err&epoch<=epo）
????u=0;
????for?p=1:P
????????for?j=1:n
????????????u=0;
????????????for?k=1:m
????????????????u=u+WW（jk）*x（pk）;
????????????end
????????????net（pj）=u;
????????end
????end
????for?p=1:P
????????for?j=1:n
????????????u=net（pj）;
????????????u=（u-b（j））/a（j）;
????????????phi（pj）=cos（1.75*u）*exp（-u*u/2）;
????????end
????end
????for?p=1:P;
????????for?i=1:N
????????????u=0;
????????????for?j=1:n
????????????????u=u+W（ij）*phi（pj）;
????????????end
????????????y（pi）=delta*abs（u）;
????????end
????end
????u=0;
????for?p=1:P
????????for?i=1:N
????????????u=u+abs（d（pi）*log（y（pi））+（1-d（pi）*log（1-y（pi））））;
????????end
????end
????err