close all	%關(guān)閉所有的窗口
        clear all	%清除所有的變量
        k=input('訓(xùn)練集點(diǎn)數(shù)=k,k=[100]   ');
        if isempty(k),
           k=100;
        end;
        L=input('最大隱層神經(jīng)元個(gè)數(shù)=L,L=[20]   ');
        if isempty(L),
           L=20;
        end;
        S=input('歸一化參數(shù)=S,S=[0.5]   ');
        if isempty(S),
           S=0.5;
        end;

        X=10*rand(1,k);
        Y=sin(X).*exp(-X);
       
        
        net=newrb(X,Y,0,S,L,5);
        
        
        
        X2=linspace(0,10,1000);
        Y2=sin(X2).*exp(-X2);
        
        X1=10.*rand(1,30);
        Y1=sin(X1).*exp(-X1);
        R=sim(net,X1);
        
        figure(1)
        plot(X2,Y2,'-k',X1,R,'.b',X1,Y1,'+r')
        grid on
        xlabel('x in [0,10]')
        ylabel('y=sin(x).*exp(-x)')
        title('簡(jiǎn)單bp網(wǎng)絡(luò)應(yīng)用于函數(shù)逼近')