rbf.rar_RBF_RBF拟合_rbf拟合函数_函数拟合

function rbf(x,d,test,dtest,unit) %Radial Basis Function for fitting %x stands for train sample %d stands for x's teacher %test stands for test sample %dtest stands for dtest's real output %unit stands for the number of system's hidden units tic dmax=max(x)-min(x); delta=dmax/(sqrt(unit)); tk=linspace(min(x),max(x),unit); len=length(x); t=x'*ones(1,unit)-ones(len,1)*tk; phi=exp(-t.^2/(delta^2)); w=pinv(phi)*d'; lent=length(test); testt=test'*ones(1,unit)-ones(lent,1)*tk; phit=exp(-testt.^2/(delta^2)); ytest=phit*w; toc figure,plot(test,ytest,'.',test,dtest,'r'); title( ['number of rbf: ', num2str(unit),' ;time cost: ',num2str(toc),'sec'] ) end