We propose amixed spectral method for heat transfer in unbounded domains, using generalised Hermite functions and Legendre polynomials. Some basic results on the mixed generalised Hermite-Legendre orthogonal approximation are established, which plays important roles in spectral methods for various problems defined on unbounded domains. As an example, the mixed generalised Hermite-Legendre spectral scheme is constructed for anisotropic heat transfer. Its convergence is proven, and some numerical results demonstrate the spectral accuracy of this approach.