求教采用RBF网络进行辨识的程序

piggybaby

RT 因为对matlab是初学，所以好多东西都不懂。老师让我做个RBF用于传感器故障诊断的课题。第一步先让我做个RBF辨识模型。这个是实际系统模型  G(s)= 8/[（1+15S）（1+15S）（1+20S）（1+20S）]  
  还给了我个用BP做的辨识程序(付上），让我看看写个RBF的给她。但是还是没有思路，希望高手能够帮忙解答下。我的e-mail:xiaolong841009@163.com

function nni
%produce values of modle through simulation
clear;
sim('bscy');
%initiate weight
y=delaysig(yf',1,2);
u=delaysig(uf',1,2);
[w1,b1,w2,b2]=initff([y;u],18,'tansig',1,'purelin');

%for i=1:4
    sim('bscy');
    y=delaysig(yf',1,2);
    u=delaysig(uf',1,2);
    disp_frep=100;
    max_epoch=30000;
    err_goal=0.01;
    lr=0.1;
    tp=[disp_frep,max_epoch,err_goal,lr];
    [w1,b1,w2,b2,epoc,tr]=trainbpx(w1,b1,'tansig',w2,b2,'purelin',[y;u],yf',tp);
    ploterr(tr,err_goal);
    %end


%simulate
a=simuff([y;u],w1,b1,'tansig',w2,b2,'purelin');
newplot;
plot(a,'r');
hold on
plot(yf','b');
hold on
e=yf'-a;
plot(e,'k');
save wb w1 b1 w2 b2 ;