mexican.m

RBF程序

        %======================================================================
        % Copyright (C) 2005
        % All rights reserved.

        % File name: Mexican.m 
        % Notes:     应用RBF网络对 2D Mexican Hat 进行函数逼近
        %            检验反向传播的逼近效果
        %            最简单的简单的梯度下降算法

        % Version: 1.0
        % Author:  bug zhang
        % Date:    2005.1
        %======================================================================
        
        close all	%关闭所有的窗口
        clear all	%清除所有的变量
        k=input('训练集点数=k,k=[5000]   ');
        if isempty(k),
           k=5000;
        end;
        L=input('最大隐层神经元个数=L,L=[40]   ');
        if isempty(L),
           L=40;
        end;
        S=input('归一化参数=S,S=[0.5]   ');
        if isempty(S),
           S=0.5;
        end;
     

        X=4*pi*rand(1,k)-2*pi;
        Y=sin(abs(X))./(abs(X));
       
        
        net=newrb(X,Y,0,S,L,5);
        
        
        X2=linspace(0,4*pi,1000)-2*pi;
        Y2=sin(abs(X2))./(abs(X2));
        
        X1=4*pi*rand(1,100)-2*pi;
        Y1=sin(abs(X1))./(abs(X1));
        R=sim(net,X1);
        
        figure(2)
        plot(X2,Y2,'-k',X1,R,'.b',X1,Y1,'+r')
        grid on
        xlabel('x in [0,10]')
        ylabel('y=sin(x).*exp(-x)')
        title('简单bp网络应用于函数逼近')