We develop a matrix model to describe bilayered quantum Hall fluids for a series of filling factors. Considering two coupling layers, and starting from a corresponding action, we construct its vacuum configuration at ν = q i K −1 ij q j , where K ij is a 2 × 2 matrix and q i is a vector. Our model allows us to reproduce several well-known wave functions. We show that the wave function Ψ (m,m,n) constructed years ago by Yoshioka, MacDonald and Girvin for the fractional quantum Hall effect at filling factor 2 m+n and in particular Ψ (3,3,1) at filling 1 2 can be obtained from our vacuum configuration. The unpolarized Halperin wave function and especially that for the fractional quantum Hall state at filling factor 2 5 can also be recovered from our approach. Generalization to more than 2 layers is straightforward. *