In this chapter, we first use the Gray code representation of the genetic code C = 00, U = 10, G = 11, and A = 01 (C pairs with G, A pairs with U) to generate a sequence of genetic code-based matrices. In connection with these code-based matrices, we use the Hamming distance to generate a sequence of numerical matrices. We then further investigate the properties of the numerical matrices and show that they are doubly stochastic and symmetric. We determine the frequency distributions of the Hamming distances, building blocks of the matrices, decomposition and iterations of matrices. We present an explicit decomposition formula for the genetic code-based matrix in terms of permutation matrices. Furthermore, we establish a relation between the genetic code and a stochastic matrix based on hydrogen bonds of DNA. Using fundamental properties of the stochastic matrices, we determine explicitly the decomposition formula of genetic code-based biperiodic table. By iterating the stochastic matrix, we demonstrate the symmetrical relations between the entries of the matrix and DNA molar concentration accumulation. The evolution matrices based on genetic code were derived by using hydrogen bondsbased symmetric stochastic (2x2)-matrices as primary building blocks. The fractal structure of the genetic code and stochastic matrices were illustrated in the process of matrix decomposition, iteration and expansion in corresponding to the fractal structure of the biperiodic table introduced by Petoukhov (2001a, 2001b, 2005).