我见过的最脑残也是最好懂的人工神经网络算法教程（五）

数学中国YY主管

4.5 神经网络的编码（Encoding the Network ）

　　 在本书的开始几章中，你已经看到过怎样用各种各样的方法为遗传算法编码。但当时我并没有向你介绍过一个用实数编码的具体例子，因为我知道我要留在这里向你介绍。我曾经讲到，为了设计一个前馈型神经网络，编码是很容易的。我们从左到右读每一层神经细胞的权重，读完第一个隐藏层，再向上读它的下一层，把所读到的数据依次保存到一个向量中，这样就实现了网络的编码。因此，如果我们有图14所示的网络，则它的权重编码向量将为:

' x* x+ x8 f( Y1 `6 H, c- N4 w, ~

     0.3, -O.8， -O.2， 0.6， O.1， -0.l， 0.4， 0.5

　　在这一网络中，为了简单，我没有把偏移值的权重包括进去。但在实际实现编码时，你必须包含偏移值这个权重，否则你肯定无法获得你所需要的结果。

$ I6 S# u  M' n5 j. t# k& `

/ Z' W. X( J( \% j/ K

图14 为权重编码。

在此之前讲的事情你都懂了吗？好极了，那下面就让我们转来考虑，怎样用遗传算法来操纵已编码的基因吧。

0 Q$ w! g! p* k: J% S4 j, Y

4.6  遗传算法（The Genetic Algorithm ）

        到此，所有的权重已经象二进制编码的基因组那样，形成了一个串，我们就可以象本书早先讨论过的那样来应用遗传算法了。遗传算法（GA）是在扫雷机已被允许按照用户指定的帧数（为了某种缘故, 我下面更喜欢将帧数称作滴答数，英文是ticks）运转后执行的。你可以在ini文件中找到这个滴答数（iNumTicks）的设置。下面是基因组结构体的代码。这些对于你应该是十分面熟的东西了。

0 U) y5 r* o8 B: M$ ^% d5 @" }

Struct SGenome

{

  vector <double> 　 vecWeights;

  double 　　　　　 dFitness;

! a& J+ {! i+ I7 ~0 V

  SGenome():dFitness(0) {}

  SGenome(vector <double> w, double f):vecWeights(w),dFitness(f){}

  //重载'<'的排序方法

      friend bool  operator<(const SGenome& lhs, const SGenome& rhs)

         {

     　      return (lhs.dFitness < rhs.dFitness);

         }

};

3 x# Q4 \  ^7 Y( \" Z

　　从上面的代码你可看出，这一SGenome结构和我们在本书所有其他地方见到的SGenome结构几乎完全一致，唯一的差别就是这里的染色体是一个双精度向量std::vector。因此，可以和通常一样来应用杂交操作和选择操作。但突变操作则稍微有些不同，这里的权重值是用一个最大值为dMaxPerturbation的随机数来搔扰的。这一参数dMaxPerturbation在ini文件中已作了声明。另外，作为浮点数遗传算法，突变率也被设定得更高些。在本工程中，它被设成为0.1。 下面就是扫雷机工程遗传算法类中所见到的突变函数的形式:

3 x" B6 o- D6 i* Q9 P- H/ R

void CGenAlg::Mutate(vector<double> &chromo)

{

  // 遍历权重向量，按突变率将每一个权重进行突变

  for (int i=0; i<chromo.size(); ++i)

  {

     // 我们要骚扰这个权重吗？

     if (RandFloat() < m_dMutationRate)

     {

       // 为权重增加或减小一个小的数量

       chromo += (RandomClamped() * CParams::dMaxPerturbatlon);

     }

  }

}

　　如同以前的工程那样，我已为v1.0版本的Smart Minesweepers工程保留了一个非常简单的遗传算法。这样就能给你留下许多余地，可让你利用以前学到的技术来改进它。就象大多数别的工程一样，v1.O版只用轮盘赌方式选精英，并采用单点式杂交。

# c1 y2 c- H# l5 z  n- T

注意:

当程序运行时，权重可以被演化成为任意的大小，它们不受任何形式的限制。

1 w/ I9 P7 K: O7 a; p
0 y# x$ u& H5 h2 G0 C1 \

4.7  扫雷机类（The CMinesweeper Class ）

　　 这一个类用来定义一个扫雷机。就象上一章描述的登月艇类一样，扫雷机类中有一个包含了扫雷机位置、速度、以及如何转换方向等数据的纪录。类中还包含扫雷机的视线向量（look-at vector）；它的2个分量被用来作为神经网络的2个输入。这是一个规范化的向量，它是在每一帧中根据扫雷机本身的转动角度计算出来的，它指示了扫雷机当前是朝着哪一个方向，如图11所示。 下面就是CMinesweeper扫雷机类的声明:

class CMinesweeper

{

private:

     // 扫雷机的神经网络

      CNeuralNet        m_ItsBrain;

# \+ ~# F6 \) X# H8 |

     // 它在世界坐标里的位置

     SVector2D         m_vPosition;

     // 扫雷机面对的方向

     SVector2D          m_vLookAt;

" P( r- Z: u& P  i

     // 它的旋转(surprise surprise)

     double      m_dRotation;

     double      m_dSpeed;

     // 根据ANN保存输出

     double         m_lTrack,

                    m_rTrack;

0 d! n1 C0 m" p, i+ D5 T9 l2 ^2 G
- K/ e, H! p% c- q( R% _4 l( d

5 }6 H2 s2 w2 g2 _1 W

     m_lTrack和m_rTrack根据网络保存当前帧的输出。

     这些就是用来决定扫雷机的移动速率和转动角度的数值。

+ ?$ B" z) \) g) c  l, k

     // 用于度量扫雷机适应性的分数

     double   　　　 m_dFitness;

0 b: @" y  u) H. K! {* G
: [5 D4 D# U/ q/ Y/ |
9 u% W: c8 K' P

     每当扫雷机找到一个地雷，它的适应性分数就要增加。

% {- m* L: k& H0 v% A! p  V


) l& q; j( n* {, z# C) c

    // 扫雷机画出来时的大小比例

     double      m_dScale;

8 Z5 D' `+ q$ z  V# X; w9 a

    // 扫雷机最邻近地雷的下标位置

     int       m_iClosestMine;

+ K  |3 H9 M1 ^0 V1 B  G9 h, u" o
$ v1 ]/ M: a3 i

7 v$ w, a5 O, c' b# G. P! i2 D

    在控制器类CControl1er中，有一个属于所有地雷的成员向量std::vector。而m_iClosestMine就是代表最靠近扫雷机的那个地雷在该向量中的位置的下标。

- u3 \' }4 p8 P# O
& a/ d/ I% @: P4 g
3 L; w& @, V' t$ _, y, C# O) R% p7 s
1 Q8 e5 O1 I4 J; f6 M! u. w. c

public:

2 S+ v% C& T5 C

    CMinesweeper();

    // 利用从扫雷机环境得到的信息来更新人工神经网

    bool Update(vector<SVector2D> &mines);

    // 用来对扫雷机各个顶点进行变换，以便接着可以画它出来

    void  WorldTransform(vector<SPoint> &sweeper);

: S4 N+ b0 J4 d

    // 返回一个向量到最邻近的地雷

    5Vector2D GetClosestMine(vector<SVector2D> &objects);

( f+ q. T# k. H, x* l+ z

    // 检查扫雷机看它是否已经发现地雷

    int        CheckForMine(vector<SVector2D> &mines, double size);

    void       Reset();

    // ----------------- 定义各种供访问用的函数

    SVector2D  Position()const { return m_vPosition; }

    void       IncrementFitness(double val) { m_dFitness += val; }

    double     Fitness()const { return m_dFitness; }

    void       PutWeights(vector<double> &w) { m_ItsBrain.PutWeights(w); }

    int        GetNumberOfWeights()const

                              { return m_ItsBrain.GetNumberOfWeights(); }

};

: }5 O3 i, ~% i+ _9 K# j5 I

4.7.1 The CMinesweeper::Update Function （扫雷机更新函数）

' B. ~* K& [$ E( m$ }( j

　　需要更详细地向你说明的CMinesweeper类的方法只有一个，这就是Update更新函数。该函数在每一帧中都要被调用，以更新扫雷机神经网络。让我们考察这函数的肚子里有些什么货色:

bool CMinesweeper::Update(vector<SVector2D> &mines)

{

     //这一向量用来存放神经网络所有的输入

     vector<double> inputs;

     //计算从扫雷机到与其最接近的地雷（2个点）之间的向量

     SVector2D vClosestMine = GetClosestMine(mines);

     //将该向量规范化

     Vec2DNormalize(vClosestMine);

4 b" p6 K! ?+ y% i

　　首先，该函数计算了扫雷机到与其最靠近的地雷之间的向量，然后使它规范化。（记住，向量规范化后它的长度等于1。）但扫雷机的视线向量（look-at vector）这时不需要再作规范化，因为它的长度已经等于1了。由于两个向量都有效地化成了同样的大小范围，我们就可以认为输入已经是标准化了，这我前面已讲过了。

( K% `& F- y; Q/ P4 W' z- j& d

     //加入扫雷机->最近地雷之间的向量

     Inputs.push_back(vClosestMine.x);

     Inputs.push_back(vCIosestMine.y);

) H& j. s( L' H) ^- }* @

     //加入扫雷机的视线向量

     Inputs.push_back(m_vLookAt.x);

     Inputs.push_back(m_vLookAt.y);

     //更新大脑，并从网络得到输出

     vector<double> output = m_ItsBrain.Update(inputs);

2 k+ K( N0 ^! G% y1 Y0 k

然后把视线向量，以及扫雷机与它最接近的地雷之间的向量，都输入到神经网络。函数CNeuralNet::Update利用这些信息来更新扫雷机网络，并返回一个std::vector向量作为输出。

, l( |$ D; L4 }% [- ~5 D7 J; l+ n

     //保证在输出的计算中没有发生错误

     if (output.size() < CParams::iNumOutputs)

      {

         return false;

      }

" \5 Y$ G) j6 v

     // 把输出赋值到扫雷机的左、右轮轨

     m_lTrack = output[0];

     m_rTrack = output[1];

+ R; N' z* m$ ?# ]  H- y
4 `; z. r2 I# o1 e8 n5 {

　　 在更新神经网络时，当检测到确实没有错误时，程序把输出赋给m_lTrack和m_rTrack。 这些值代表施加到扫雷机左、右履带轮轨上的力。

     // 计算驾驶的力

     double RotForce = m_lTrack - m_rTrack;

     // 进行左转或右转

     Clamp(RotForce, -CParams::dMaxTurnRate, CParams::dMaxTurnRate);

     m_dSpeed = (m_lTrack + m_rTrack);

9 A: N% R* m9 U* m& Y- o6 ?. p4 I0 U/ [

　　 扫雷机车的转动力是利用施加到它左、右轮轨上的力之差来计算的。并规定，施加到左轨道上的力减去施加到右轨道上的力，就得到扫雷机车辆的转动力。然后就把此力施加给扫雷机车，使它实行不超过ini文件所规定的最大转动率的转动。而扫雷机车的行进速度不过就是它的左侧轮轨速度与它的右侧轮轨速度的和。既然我们知道了扫雷机的转动力和速度，它的位置和偏转角度也就都能更新了。

     //更新扫雷机左右转向的角度

     m_dRotation += RotForce;

     // 更新视线角度

     m_vLookAt.x = -sin(m_dRotation);

     m_vLookAt.y = cos(m_dRotation);

+ [' P) Z1 q/ N: q

     // 更新它的位置

     m_vPosition += (m_vLookAt* m_dSpeed);

  l$ }1 J& a% Q

     // 如果扫雷机到达窗体四周，则让它实行环绕，使它不至于离开窗体而消失

     If (m_vPosition.x > CParams::WindowWidth) m_vPosition.x = 0;

     If (m_vPosition.x < 0) m_vPosition.x = CParams::WindowWidth;

     If (m_vPosition.y > CParams::WindowHeight) m_vPosition.y = 0;

     If (m_vPosition.y < D) m_vPosition.y = CParams::WindowHeight;

　　为了使事情尽可能简单，我已让扫雷机在碰到窗体边框时就环绕折回(wrap)。采用这种方法程序就不再需

要做任何碰撞-响应方面的工作。环绕一块空地打转对我们人来说是一桩非常不可思议的动作，但对扫雷机，这

就像池塘中的鸭子。

& q% q8 H0 L4 L. o: B

      Returen true;

}

4.8  CController Class （控制器类）

CController类是和一切都有联系的类。图15指出了其他的各个类和CController类的关系。

下面就是这个类的定义:

class CController

{

private:

     // 基因组群体的动态存储器（一个向量）

      vector<SGenome>   m_vecThePopulation;

+ u4 o8 J( d, m: D4 f+ r4 f; [& z! S
3 C$ O( z2 p' l. a4 {

      图15  minesweeper工程的程序流程图

4 Q0 n9 b# Q/ h1 S# Z. K
& X/ _# n2 ~# F$ T' V

    // 保存扫雷机的向量

    vector<CMinesweeper>  m_vecSweepers;

    // 保存地雷的向量

    vector<SVector2D>   m_vecMines;

3 y; E1 R) Z9 B5 S) S- @6 x% g

    // 指向遗传算法对象的指针

    CGenAIg*          m_pGA;

    int               m_NumSweepers;

    int               m_NumMines;

5 `" K5 Q5 y) a# y% w0 G

    // 神经网络中使用的权重值的总数

    int               m_NumWeightsInNN;

    // 存放扫雷机形状各顶点的缓冲区

    vector<SPoint>    m_SweeperVB;

" H  ^/ i& A6 O8 h* H

    // 存放地雷形状各顶点的缓冲区

    vector<SPoint>    m_MineVB;

    // 存放每一代的平均适应性分数，供绘图用

    vector<double>    m_vecAvFitness;

    // 存放每一代的最高适应性分

    vector<double>    m_vecBestFitness;

    // 我们使用的各种不同类型的画笔

    HPEN              m_RedPen;

    HPEN              m_BluePen;

    HPEN              m_GreenPen;

    HPEN              m_OldPen;

    // 应用程序窗口的句柄

    HWND              m_hwndMain;

   // 切换扫雷机程序运行的速度

    bool              m_bFastRender;

    // 每一代的帧数（滴答数）

    int               m_iTicks;

    // 代的计数

    int               m_iGenerations;

' k, T4 E! Q; j1 c8 Z- U9 C

    // 窗体客户区的大小

    int               cxClient，cyClient;

    // 本函数在运行过程中画出具有平均-，和最优适应性值的图

    void              PlotStats(HDC surface);

public:

    CController(HWND hwndMain);

  o7 U) r* {4 T& j0 h4 G

    ~CController();

    void   Render(HDC surface);

    void   WorldTransform(vector<SPoint> &VBuffer,

           SVector2D        vPos);

1 b+ Z' Q# `8 _8 v

    bool   Update();

4 K" E) Q- Y+ }6 b. C+ i

    // 几个公用的访问方法

    bool   FastRender() { return m_bFastRender; }

    void   FastRender(bool arg){ m_bFastRender = arg; }

    void   FastRenderToggle() { m_bFastRender = !m_bFastRender; }

};

* c! g9 y; N) v1 d
: h9 s3 n& s% r- o" m$ r; D1 v

当创建CController类的某个实例时，会有一系列的事情发生：

*创建CMinesweeper对象。

*统计神经网络中所使用的权重的总数，然后此数字即被利用来初始化遗传算法类的一个实例。

*从遗传算法对象中随机提取染色体(权重)并(利用细心的脑外科手术)插入到扫雷机的经网络中。

*创建了大量的地雷并被随机地散播到各地。

*为绘图函数创建了所有需要用到的GDI画笔。

*为扫雷机和地雷的形状创建了顶点缓冲区。

所有的一切现都已完成初始化，由此Update方法就能在每一帧中被调用来对扫雷机进行演化。

4 J! S8 X2 Z7 @4 T$ [  M" e
% ?! n! ?  Q0 k

4.8.1  CController::Update Method （控制器的更新方法）

1 V4 c( s2 z3 ], s  X

　　控制器更新方法CController::Update方法（或函数）在每一帧中都要被调用。当调用update函数时，函数的前一半通过对所有扫雷机进行循环，如发现某一扫雷机找到了地雷，就update该扫雷机的适应性分数。由于m_vecThePopulation包含了所有基因组的拷贝，相关的适应性分数也要在这时进行调整。如果为完成一个代（generation）所需要的帧数均已通过，本方法就执行一个遗传算法时代(epoch)来产生新一代的权重。这些权重被用来代替扫雷机神经网络中原有的旧的权重，使扫雷机的每一个参数被重新设置，从而为进入新一generation做好准备。

  bool CController::Update()

{

     // 扫雷机运行总数为CParams::iNumTicks次的循环。在此循环周期中，扫雷机的神经网络

     // 不断利用它周围特有的环境信息进行更新。而从神经网络得到的输出，使扫雷机实现所需的

     // 动作。如果扫雷机遇见了一个地雷，则它的适应性将相应地被更新，且同样地更新了它对应

     // 基因组的适应性。

     if (m_iTicks++ < CParams::iNumTicks)

       {

        for (int i=O; i<m_NumSweepers; ++i)

        {

        //更新神经网络和位置

          if (!m_vecSweepers.Update(m_vecMines))

           {

            //处理神经网络时出现了错误，显示错误后退出

             MessageBox(m_hwndMain, 'Wrong amount of NN inputs!",

               "Error", MB_OK);

1 L; j! o0 Z$ S% q& L' z

             return false;

           }

# P# ]! E8 f3 J9 t

         // 检查这一扫雷机是否已经发现地雷

         int GrabHit = m_vecSweepers.CheckForMine(m_vecMines，

                                                     CParams::dMineScale);

( t+ L0 w3 g3 O9 Q

         if (GrabHit >= 0)

          {

            // 扫雷机已找到了地雷，所以要增加它的适应性分数

            m_vecSweepers.IncrementFitness();

( Y, g8 a5 q1 Z  ~5 g* @, |

            // 去掉被扫雷机找到的地雷，用在随机位置放置的一个新地雷来代替

            m_vecMines[GrabHit] = SVector2D(RandFloat() * cxClient，

                                          RandFloat() * cyClient);

          }

& C" ^) O/ T1 i( O2 p- o$ v

       // 更新基因组的适应性值

        m-vecThePopulation.dFitness = m_vecSweepers.Fitness();

      }

   }

   // 一个代已被完成了。

   // 进入运行遗传算法并用新的神经网络更新扫雷机的时期

   else

    {

     // 更新用在我们状态窗口中状态

     m_vecAvFitness.push_back(m_pGA->AverageFitness());

     m_vecBestFitness.push_back(m_pGA->BestFitness());

# H- m# o. K" Z( [1 D( e5 t( ~

     // 增加代计数器的值

     ++m_iGenerations;

     // 将帧计数器复位

     m_iTicks = 0;

     // 运行GA创建一个新的群体

     m-vecThePopulation = m_pGA->Epoch(m_vecThePopulation);

     // 在各扫雷机中从新插入新的（有希望）被改进的大脑

     // 并将它们的位置进行复位，等

     for(int i=O; i<m_NumSweepers; ++i)

      {m_vecSweepers.m_ItsBrain.PutWeights(m_vecThePopulation.vecWeights);

. c7 R& Q- k! K" z) _/ x  T' {

       m_vecSweepers.Reset();

      }

    }

  returen true;

}

* u8 z& P* N0 M2 o8 x+ G+ E
0 V9 m  B; A3 o9 A6 S& \

概括起来，程序为每一世代做的工作是:

     l．为所有扫雷机和为iNumTicks个帧组织循环，调用Update函数并根据情况增加扫雷机适应值的得分。

     2．从扫雷机神经网络提取权重向量。

     3．用遗传算法去演化出一个新的网络权重群体。

     4．把新的权重插入到扫雷机神经网络。

　 5．转到第1步进行重复，直到获得理想性能时为止。

最后，表3列出了Smart Sweepers工程 v1.0版所有缺省参数的设置值。

( Z7 G" w8 i8 G) C/ _+ V

4.9  运行此程序 (Running the Program )

　　    当你运行程序时，“F”键用来切换2种不同的显示状态，一种是显示扫雷机怎样学习寻找地雷，一种是示在运行期中产生的最优的与平均的适当性分数的统计图表。 当显示图表时，程序将会加速运行。

) [  W9 ~2 M  A* }- g* x

phcs

支持数学建模，支持数学中国！

phcs

支持数学建模，支持数学中国！