Abstract
The synchronous Boolean network (SBN) is a simple and powerful model for describing, analyzing, and simulating cellular biological networks. This paper seeks a complete understanding of the dynamics of such a model by employing a matrix method that relies on relating the network transition matrix to its function matrix via a self-inverse state matrix. A recursive ordering of the underlying basis vector leads to a simple recursive expression of this state matrix. Hence, the transition matrix is computed via multiplication of binary matrices over the simplest finite (Galois) field, namely the binary field GF(2), i.e., conventional matrix multiplication involving modulo-2 addition, or XOR addition. We demonstrate the conceptual simplicity and practical utility of our approach via an illustrative example, in which the transition matrix is readily obtained, and subsequently utilized (via its powers, characteristic equation, minimal equation, 1-eigenvectors, and 0-eigenvectors) to correctly predict both the transient behavior and the cyclic behavior of the network. Our matrix approach for computing the transition matrix is superior to the approach of scalar equations, which demands cumbersome manipulations and might fail to predict the exact network behavior. Our approach produces result that exactly replicate those obtained by methods employing the semi-tensor product (STP) of matrices, but achieves that without sophisticated ambiguity or unwarranted redundancy.