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

        % File name: Gabor.m 
        % Notes:     應用RBF網絡對 Gabor 進行函數逼近
        %            檢驗反向傳播的逼近效果
        %            最簡單的簡單的梯度下降算法

        % Version: 1.0
        % Author:  bug zhang
        % Date:    2005.1
        %======================================================================     
        
        close all	%關閉所有的窗口
        clear all	%清除所有的變量
        k=input('訓練集點數=k,k=[3000]   ');
        if isempty(k),
           k=3000;
        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;
       
        X1=rand(1,k);
        X2=rand(1,k);
        X=[X1;X2];
        Y=pi/2.*exp(-2.*(X1.^2+X2.^2)).*cos(2*pi.*(X1+X2));
        
        net=newrb(X,Y,0,S,L,5);
        
        X1=rand(1,30);
        X2=rand(1,30);
        X=[X1;X2];
        Y1=pi/2.*exp(-2.*(X1.^2+X2.^2)).*cos(2*pi.*(X1+X2));
        Y2=sim(net,X);
        
        R=linspace(min(Y1),max(Y1),1000);
        T=R;
        
        figure(2)
        plot(R,T,'-k',Y1,Y2,'+r')