We introduce the first random matrix model of a complex β-ensemble. The matrices are tridiagonal and can be thought of as the non-Hermitian analogue of the Hermite β-ensembles discovered by [I. Dumitriu and A. Edelman, Matrix models for beta ensembles, J. Math. Phys. 43 (2002) 5830–5847]. The main feature of the model is that the exponent β of the Vandermonde determinant in the joint probability density function (j.p.d.f.) of the eigenvalues can take any value in ℝ+. However, when β=2, the j.p.d.f. does not reduce to that of the Ginibre ensemble, but it contains an extra factor expressed as a multidimensional integral over the space of the eigenvectors.