We introduce two families of random tridiagonal block matrices for which the joint eigenvalue distributions can be computed explicitly. These distributions are novel within random matrix theory, and exhibit interactions among eigenvalue coordinates beyond the typical mean-field log-gas type. Leveraging the matrix models, we go on to describe the point process limits at the edges of the spectrum in two ways: through certain random differential operators, and also in terms of coupled systems of diffusions. Along the way we establish several algebraic identities involving sums of Vandermonde determinant products.