We study fermionic one-matrix, two-matrix and D-dimensional gauge invariant matrix models. In all cases we derive loop equations which unambiguously determine the large-N solution. For the one-matrix case the solution is obtained for an arbitrary interaction potential and turns out to be equivalent to the one for the hermitian one-matrix model with a logarithmic potential and, therefore, belongs to the same universality class. The explicit solutions for the fermionic two-matrix and D-dimensional matrix models are obtained at large N (or in the spherical approximation) for the quadratic potential.