在布尔控制网络中，一个节点的状态可以用一个二元变量来描述，即“0”或“1”。不同节点之间的相互影响可以由布尔逻辑函数来描述。所以，布尔控制网络的动态方程 [ 11 ] 为：

其中 x a ( t ) ∈ P , a = 1 , 2 , ⋯ , m 是状态变量， u b ( t ) ∈ P , b = 1 , 2 , ⋯ , n 为控制变量， f a 为 P m + n → P 的逻辑函数， t = 0 , 1 , 2 , ⋯ 是时间序列。

1) 一个较复杂的N节点布尔控制网络的一般形式不方便储存和传输，可以用布尔矩阵对其进行描述(每行只有一个元素为“1”并且其余元素均为“0”的 2 N × 2 N -阶矩阵被称作布尔矩阵)。我们利用布尔网络线性化矩阵 [ 12 ] 的概念对布尔控制网络进行编码及线性化，得到相应的线性化矩阵，就称作该布尔控制网络对应的布尔矩阵。式(1.1)中的布尔控制网络对应一个 ( 2 m + n × 2 m + n ) -阶布尔矩阵B：

B = [ 1 2 3 ⋯ 2 m + n − 1 2 m + n i 1 i 2 i 3 ⋯ i 2 m + n − 1 i 2 m + n ] k i k (1.2)

其中，i_k代表B的第k个行向量的第i_k个元素为1。矩阵中的元素(k,i_k)表示布尔控制网络的第k个节点下一时刻的状态为i_k。

2) 根据代数学的基本知识，任一个 2 m + n -阶的布尔矩阵B存在一个非退化的矩阵W，使得 B ˜ = W B W − 1 ，其中 B ˜ 是如下形式 [ 13 ] ：

B ˜ = ( A 0 A 1 A 2 ⋯ A n O B i 1 0 ⋯ 0 O O B i 2 ⋯ O ⋮ ⋮ ⋮ ⋱ ⋮ O O O ⋯ B i n ) ， B i j = ( C j O O ⋯ O O C j O ⋯ O O O C j ⋯ O ⋮ ⋮ ⋮ ⋱ O O O O O C j )

C j = ( 0 1 0 ⋯ 0 0 0 1 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 ⋯ 1 1 0 0 ⋯ 0 ) (1.3)

其中，A₀是 s 0 × s 0 阶矩阵，A_j是 s 0 × i j s i j 阶矩阵 ( j = 1 , 2 , ⋯ , r ) ，C_j是 i j × i j 阶矩阵 ( j = 1 , 2 , ⋯ , r ) ， { i 1 , i 2 , ⋯ , i r } ⊆ { 1 , 2 , ⋯ , 2 n } ， i 1 < i 2 < ⋯ < i r ， s i j ( j = 1 , 2 , ⋯ , r ) 是正整数，并且有：

在布尔矩阵B中，可以求得B的特征多项式：

Φ ( λ ) = λ s 0 ( λ − 1 ) s 1 ( λ 2 − 1 ) s 2 ⋯ ( λ i r − 1 ) s i r (1.5)

其中 s 0 ≠ 2 m + n , s i ( i = 1 , 2 , ⋯ , 2 m + n ) 是非负整数。

3) 布尔控制网络的状态转移过程可以用有向图的拓扑结构 [ 13 ] 来描述，从而得到其环结构和暂态等动力学性质。一般地，周期环分为两类：一是吸引子，二是伊甸园孤岛。因此，我们可以利用线性化矩阵的特征多项式与动力学性质的关系，对布尔控制网络进行分析和研究。

定理 [ 14 ] 对任意一个布尔网络，其对应的布尔矩阵为B，B所对应的特征多项式为 Φ ( λ ) = λ s 0 ( λ − 1 ) s 1 ( λ 2 − 1 ) s 2 ⋯ ( λ i r − 1 ) s i r ，那么有：

1) 如果 s 0 = 0 ，那么存在 ( s i 1 + s i 2 + ⋯ + s i r ) 个周期环，并且他们都是伊甸园孤岛，它们的周期分别是 i j ，并且满足 ，其中 s i j 是 i j -周期环的数量；

2) 如果 s 0 ≠ 0 ，那么存在 s 0 个暂态和 ( s i 1 + s i 2 + ⋯ + s i r ) 个周期环，并且这些周期环的周期分别是 i j ，进一步地满足 s 0 + i 1 s i 1 + i 2 s i 2 + ⋯ + i r s i r = 2 m + n ，其中 s i j 是 i j -周期环的数量；另外，如果存在正整数 q ( 1 ≤ q ≤ r ) ，使得矩阵中的A_q是一个零矩阵，那么块矩阵 B i j 所对应的 i j -周期环是伊甸园孤岛。

系统稳定性是指系统在某些子系统行为失调时恢复到原来状态的能力或性质。布尔网络的动力学性质模型不依赖于参数，因此可以定性描述布尔控制网络的系统稳定性。

通向吸引子的状态集称为吸引子盆 [ 15 ] ，并且吸引子盆内的各种状态最终都通向吸引子的状态集。吸引子盆能使系统经过适当干扰后仍重复原来的状态循环，所以吸引子可以保持系统的稳定性。在布尔控制网络中，变量经过局部的突变重复原来的状态循环，就称该状态稳定。如果一个布尔控制网络受到干扰后恢复到原来状态的能力越强，系统稳定性就越高，反之系统稳定性越低。另外，吸引盆区域越大说明吸引子的影响力也越大，从而系统的稳定性就越强。

考虑一个布尔控制网络：哺乳动物细胞周期调节的核心网络 [ 16 ] 。该网络包含一个输入CycD和九个状态(蛋白质Rb，E2，CycE，CycA，p27，Cdc20，Cdh1，UbcH10和CycB)，依次取u，x₁，x₂，x₃，x₄，x₅，x₆，x₇，x₈，和x₉，分别表示这十个变量激活(activation)与抑制(inactivation)两个状态，得到相应的布尔逻辑方程：

{ x 1 ( t + 1 ) = [ u ( t ) ∩ x ¯ 3 ( t ) ∩ x ¯ 4 ( t ) ∩ x ¯ 9 ( t ) ] ∪ [ x ¯ 5 ( t ) ∩ u ¯ ( t ) ∩ x ¯ 9 ( t ) ] x 2 ( t + 1 ) = [ x ¯ 1 ( t ) ∩ x ¯ 4 ( t ) ∩ x ¯ 9 ( t ) ] ∪ [ x ¯ 5 ( t ) ∩ x 1 ( t ) ∩ x ¯ 9 ( t ) ] x 3 ( t + 1 ) = x 2 ( t ) ∩ x ¯ 1 ( t ) x 4 ( t + 1 ) = [ ( x 2 ( t ) ∩ x ¯ 1 ( t ) ) ∩ x ¯ 6 ( t ) ∩ [ x 7 ( t ) ∩ x 8 ( t ) ] ¯ ]                                   ∪ [ ( x 4 ( t ) ∩ x ¯ 1 ( t ) ) ∩ x ¯ 6 ( t ) ∩ ( x 7 ( t ) ∩ x 8 ( t ) ) ¯ ] x 5 ( t + 1 ) = [ u ( t ) ∩ x ¯ 3 ( t ) ∩ x ¯ 4 ( t ) ∩ x ¯ 9 ( t ) ]                                   ∪ [ x 5 ( t ) ∩ ( x 3 ( t ) ∩ x 4 ( t ) ) ¯ ∩ u ¯ ( t ) ∩ x ¯ 9 ( t ) ] x 6 ( t + 1 ) = x 9 ( t ) x 7 ( t + 1 ) = [ x ¯ 4 ( t ) ∩ x ¯ 9 ( t ) ] ∪ x 6 ( t ) ∪ [ x 5 ( t ) ∩ x ¯ 9 ( t ) ] x 8 ( t + 1 ) = x ¯ 7 ( t ) ∪ [ x 7 ( t ) ∩ x 8 ( t ) ∩ ( x 6 ( t ) ∩ x 4 ( t ) ∩ x 9 ( t ) ) ] x 9 ( t + 1 ) = x ¯ 6 ( t ) ∩ x ¯ 7 ( t ) (1.6)

该网络对应的布尔矩阵(部分矩阵)如下：

[ 1 2 3 4 ⋯ 1021 1022 1023 1024 408 136 12 12 ⋯ 263 7 15 15 ] k i k (1.7)

所对应布尔矩阵的特征多项式为：

Φ ( λ ) = λ 1020 ( λ 4 − 1 ) (1.8)

可以得到，该布尔控制网络共有1024个状态，只有1个吸引盆，其中1020个暂态，一个4-周期环。该细胞周期调节网络有1020个状态都收敛到44，9，405，199共四个吸引子上。该吸引盆覆盖网络的全部区域，表明该网络稳定程度很高。该网络对应的有向拓扑图(部分)如图1所示。

在该网络中，Rb，E2F，CycE，CycA，p27，Cdc20，Cdh1，UbcH10和CycB九个蛋白质可以组成512个状态。由于输入CycD的控制作用，当输入CycD为0或1时，这512个状态分别对应两种状态转移结果，两种结果可能相同，也可能不同。我们发现当输入CycD从0切换到1时，该调控网络277，279，280，405，407，408，437，440，469，471，501，504这12种状态都降低了272，分别对应5，7，8，133，135，136，165，168，197，199，229，232；这12种状态总共包含128个状态节点。输入CycD从0切换到1时，该调控网络261，263，26，389，391，421，424，453，455，485，488这11种状态都降低了256，分别对应5，7，8，133，135，165，168，197，199，229，232；这11种状态共包含32个状态节点。另外，有352个状态不受输入CycD的影响。可以得到，输入CycD能有规律地控制部分节点的状态转移方式，从而对该网络进行控制。这种情况对于哺乳动物细胞周期调节网络的稳定性来说是有利的，可以保证细胞周期调节功能的稳定性。另外，如果布尔控制逻辑方程中只有某一节点变化，会导致吸引子和周期环会发生较大变化，该网络将会出现新的动力学关系，吸引喷也会呈现新的控制区域。例如，将方程(1.3)中 x 1 ( t + 1 ) 进行如下变化(将 u ( t ) 与 x ¯ 3 ( t ) 之间的关系由逻辑“并”变为逻辑“或”)：

x 1 ( t + 1 ) = [ u ( t ) | x ¯ 3 ( t ) ∩ x ¯ 4 ( t ) ∩ x ¯ 9 ( t ) ] ∪ [ x ¯ 5 ( t ) ∩ u ¯ ( t ) ∩ x ¯ 9 ( t ) ] (1.9)

所对应布尔矩阵的特征多项式为 Φ ( λ ) = λ 1015 ( λ 9 − 1 ) ，形成1015个暂态和一个9-周期环。这说明细

胞周期调节网络某一部分发生变异很容易引起细胞周期调解网络的动力学关系。网络内部结构的改变容易引起内部秩序发生变化，从而引发变异。因此，该细胞周期调节网络是稳定的，但是内部结构的细微变化很可能引起整个网络结构的变化，这不利于该哺乳动物细胞的稳定性。不过，在某种程度上，它容易扩大局部变异的影响，再通过自然选择而发生进化。

图1. 哺乳动物细胞周期调节网络状态转移图

在真实世界里，网络内部中子系统之间关系的突变、环境的瞬间改变容易使网络的演化过程在某些时刻发生状态突变。本文只考虑单个节点的状态发生扰动的情形，并利用另一个布尔控制网络对扰动进行模拟，从而分析新生成网络的动力学性质及稳定性。单个节点发生扰动的布尔逻辑方程：

{ x 1 ( t + 1 ) = g 1 ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) , u 1 ( t ) , u 2 ( t ) , ⋯ , u n ( t ) ) x 2 ( t + 1 ) = g 2 ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) , u 1 ( t ) , u 2 ( t ) , ⋯ , u n ( t ) )                                                                                 ⋮ x m ( t + 1 ) = g m ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) , u 1 ( t ) , u 2 ( t ) , ⋯ , u n ( t ) ) (2.1)

例1：布尔控制网络I

布尔控制网络逻辑方程：

{ x 1 ( t + 1 ) = x 2 ( t ) x 2 ( t + 1 ) = x 3 ( t ) x 3 ( t + 1 ) = u ( t ) ∪ [ x 2 ( t ) ∩ x 3 ( t ) ] (2.2)

扰动布尔逻辑方程：

{ x 1 ( t + 1 ) = x ¯ 2 ( t ) x 2 ( t + 1 ) = x ¯ 3 ( t ) x 3 ( t + 1 ) = x 1 ( t ) ∩ x 3 ( t ) (2.3)

1) 若不考虑节点扰动的影响，该布尔控制网络对应的布尔矩阵：

[ 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 1 2 3 4 5 6 8 8 1 2 3 4 5 6 8 8 ] k i k (2.4)

其所对应布尔矩阵的特征多项式为：

Φ ( λ ) = λ 9 ( λ − 1 ) 7 (2.5)

可以得到，该布尔控制网络有7个不动点和9个暂态。该网络对应的有向拓扑图如图2所示。

2) 若对16个节点k分别进行扰动，得到结果如表1所示。

在该布尔控制网络中，对于任一节点扰动所对应布尔矩阵的特征值都为1，对应特征多项式的形式为：

Φ ( λ ) = λ 16 − n ( λ − 1 ) n , n = 6 , 7 (2.6)

不过，不同节点扰动对应的布尔矩阵不一定相同。

例2：布尔控制网络II

布尔控制网络逻辑方程：

{ x 1 ( t + 1 ) = [ x 2 ( t ) ∩ x 3 ( t ) ] ↔ x 5 ( t ) x 2 ( t + 1 ) = [ x ¯ 1 ( t ) ] ∪ x 3 ( t ) x 3 ( t + 1 ) ) = x ¯ 4 ( t ) x 4 ( t + 1 ) = [ x 1 ( t ) ∩ x 2 ( t ) ] ∪ x 3 ( t ) x 5 ( t + 1 ) = x 1 ( t ) ↔ [ x 3 ( t ) ∪ x 4 ( t ) ] (2.7)

图2. 布尔控制网络I无节点扰动的状态转移图

表1. 布尔控制网络I不同扰动的吸引子数量

扰动布尔逻辑方程：

{ x 1 ( t + 1 ) = x ¯ 2 ( t ) ∩ x 5 ( t ) x 2 ( t + 1 ) = x ¯ 2 ( t ) ∩ [ x 2 ( t ) ∪ x 3 ( t ) ] x 3 ( t + 1 ) = x ¯ 2 ( t ) ∩ u ( t ) ∩ [ x 1 ( t ) ∪ x 4 ( t ) ] x 4 ( t + 1 ) = x ¯ 2 ( t ) ∩ u ( t ) ∩ x 1 ( t ) x 5 ( t + 1 ) = x ¯ 2 ( t ) ∩ x ¯ 3 ( t ) (2.8)

1) 若不考虑节点扰动的影响，该布尔控制网络对应的布尔矩阵(部分矩阵)：

[ 1 2 3 4 ⋯ 61 62 63 64 30 30 14 14 ⋯ 11 11 27 27 ] k i k (2.9)

该网络对应的有向拓扑图(部分)如图3所示。

其所对应布尔矩阵的特征多项式为：

Φ ( λ ) = λ 62 ( λ 2 − 1 ) (2.10)

可以得到，该布尔控制网络有一个2-周期环和62个暂态。若对64个节点k分别进行扰动，得到结果如表2所示。

每一节点扰动对应布尔矩阵的特征多项式如表3所示。

上述两个例子说明不同节点的扰动所对应的布尔控制网络特征多项式有较大差异，并且改变了原网络的动力学性质，系统的稳定性也随之发生变化。在例2.2中，节点6的扰动使原网络产生了一个2-周期环和一个不动点；对节点15进行扰动，产生了一个2-周期环和一个3-周期环，出现了两个吸引盆；两

图3. 布尔控制网络II无节点扰动的状态转移图

表2. 布尔控制网络II不同扰动的吸引子数量

表3. 布尔控制网络II的不同扰动的特征多项式

注：若特征多项式相同，对应的布尔矩阵不一定相同。

者所对应布尔控制网络的系统稳定性都发生了较大变化。因此，在布尔控制网络中，单个节点发生扰动将对网络的演化方式有较大影响，从而影响对应布尔控制网络的动力学性质，改变系统的稳定性。

若性质全同或相近的两个或多个动力系统从不同的初始条件出发，因系统间的相互作用使动力系统各自的演化状态逐步接近，并且最终演化到全同的状态，那么就称这些动力系统同步。在两个确定布尔网络的同步问题中，产生命令的布尔网络称作主动布尔网络；带有控制输入及通过状态反馈 [ 17 ] 接收命令的布尔网络称作从动布尔网络 [ 18 ] 。

主动布尔网络动态方程：

{ x 1 ( t + 1 ) = f 1 ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) ) x 2 ( t + 1 ) = f 2 ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) )                             ⋮ x m ( t + 1 ) = f m ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) ) (2.11)

其中 x a ( t ) ∈ P ( a = 1 , 2 , ⋯ , m ) 是状态变量， f a : P m → P ( a = 1 , 2 , ⋯ , m ) 为逻辑函数， t = 0 , 1 , 2 , ⋯ 是时间序列。

从动布尔网络的动态方程：

{ s 1 ( t + 1 ) = g 1 ( s 1 ( t ) , s 2 ( t ) , ⋯ , s m ( t ) , u 1 ( t ) , u 2 ( t ) , ⋯ , u n ( t ) ) s 2 ( t + 1 ) = g 2 ( s 1 ( t ) , s 2 ( t ) , ⋯ , s m ( t ) , u 1 ( t ) , u 2 ( t ) , ⋯ , u n ( t ) )                                                                               ⋮ s m ( t + 1 ) = g m ( s 1 ( t ) , s 2 ( t ) , ⋯ , s m ( t ) , u 1 ( t ) , u 2 ( t ) , ⋯ , u n ( t ) ) (2.12)

s a ( t ) ∈ P ( a = 1 , 2 , ⋯ , m ) 是状态变量， u b ( t ) ∈ P ( b = 1 , 2 , ⋯ , n ) 是控制变量， g a : P m + n → P ( a = 1 , 2 , ⋯ , m ) 为逻辑函数， t = 0 , 1 , 2 , ⋯ 是时间序列。

状态反馈控制的动态方程为：

{ u 1 ( t ) = h 1 ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) , s 1 ( t ) , s 2 ( t ) , ⋯ , s m ( t ) ) u 2 ( t ) = h 2 ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) , s 1 ( t ) , s 2 ( t ) , ⋯ , s m ( t ) )                                                                       ⋮ u m ( t ) = h m ( x 1 ( t ) , x 2 ( t ) , ⋯ , x m ( t ) , s 1 ( t ) , s 2 ( t ) , ⋯ , s m ( t ) ) (2.13)

h a : P 2 m → P ( a = 1 , 2 , ⋯ , m ) 为逻辑函数， t = 0 , 1 , 2 , ⋯ 是时间序列。

例1：考虑下面两个布尔网络，其中布尔网络(2.14)称为主动布尔网络I，布尔网络(2.15)称为从动布尔网络I。

{ x 1 ( t + 1 ) = x 2 ( t ) ∩ x 3 ( t ) x 2 ( t + 1 ) = x 1 ( t ) ∩ x ¯ 2 ( t ) (2.14)

{ s 1 ( t + 1 ) = s 1 ( t ) ∪ ¯ s 2 ( t ) ∩ u 1 ( t ) s 2 ( t + 1 ) = s 1 ( t ) ∪ ¯ u 2 ( t ) (2.15)

主动布尔网络对应的布尔矩阵：

[ 1 2 3 4 2 1 2 4 ] k i k (2.16)

对应的有向拓扑图如图4所示。

对应布尔矩阵的特征多项式为：

Φ ( λ ) = λ ( λ 2 − 1 ) ( λ − 1 ) (2.17)

主动布尔网络有一个暂态、一个不动点和一个2-周期环。

从动布尔网络对应的布尔矩阵：

[ 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 2 1 2 1 2 1 4 3 1 2 3 4 1 2 1 2 ] k i k (2.18)

对应的有向拓扑图如图5所示。

所对应布尔矩阵的特征多项式为：

Φ ( λ ) = λ 14 ( λ 2 − 1 ) (2.19)

可以得到，该布尔控制网络有一个2-周期环和14个暂态，主从动布尔网络对应的布尔矩阵中有三个相同的元素：节点1、节点2和节点3。主从布尔网络都包含由节点1和节点2构成的2-周期环。所以该主从动布尔网络实现了部分同步。

例2：再考虑下面两个布尔网络，其中布尔网络(2.20)称为主动布尔网络II，布尔网络(2.21)称为从动布尔网络II。

{ x 1 ( t + 1 ) = x 2 ( t ) ∩ x 3 ( t ) x 2 ( t + 1 ) = x ¯ 1 ( t ) x 3 ( t + 1 ) = x 2 ( t ) ∪ x 3 ( t ) (2.20)

{ s 1 ( t + 1 ) = s 2 ( t ) ∩ s 3 ( t ) s 2 ( t + 1 ) = s ¯ 1 ( t ) ∪ u ( t ) s 3 ( t + 1 ) = s 2 ( t ) ∪ s 3 ( t ) (2.21)

图4. 主动布尔网络I的状态转移图

图5. 从动布尔网络I的状态转移图

图6. 布尔网络II的状态转移图

状态反馈控制： u ( t ) = [ x 1 ( t ) ∪ ¯ s 1 ( t ) ] ∪ [ x 2 ( t ) ∪ ¯ s 2 ( t ) ]

主动布尔网络II对应的布尔矩阵：

[ 1 2 3 4 5 6 7 8 3 4 4 8 1 2 2 6 ] k i k (2.22)

所对应布尔矩阵的特征多项式为：

从动布尔网络II对应的布尔矩阵：

[ 1 2 3 4 5 6 7 8 3 4 4 8 1 2 2 6 ] k i k (2.24)

所对应布尔矩阵的特征多项式为：

Φ ( λ ) = λ 4 ( λ 4 − 1 ) (2.25)

可以得到，主动布尔网络II和从动布尔控制网络II都有四个暂态和一个4-周期环。主从动布尔网络的布尔简化矩阵完全相同，其有向拓扑图也完全相同，如图6所示。即，二者有相同的动力学行为，实现完全同步。