By using vertex operator representation of polynomial Lie algebra bâ(n), a theoretical interpretation of coupled BKP hierarchy is given, where corresponding Hirota bilinear equations are derived. Then based upon this, a wave function matrix is introduced, so that we can construct the corresponding dressing operator matrix and obtain the coupled BKP constraints and Sato equation, which allows us to define the Lax operator matrix and obtain the Lax equation. It is found that the coupled BKP hierarchy is a special case of matrix KP theory. Lastly, we present the BĂ€cklundâDarboux transformations for the coupled BKP equations from both fermionic and bosonic pictures.