Multifractal properties of the two-dimensional Ising model are studied by means of the f (␣) spectra of singularities of probability measures supported by energy spectra. These measures are determined through calculations of energy degeneracies for finite-size systems up to 12ϫ12 spins in cases of square and triangular lattices and up to 22ϫ22 spins in the case of the hexagonal lattice. The calculations are performed in an exact manner using the transfer-matrix method. It is argued that, in the thermodynamic limit, the scaling exponent ␣ max associated with the most probable energy of the system takes at the critical temperature a minimum value and, consequently, it is argued that a given system reveals a phase transition at a finite temperature if ␣ max possesses a minimum at the finite temperature. It is also shown that, in the thermodynamic limit, the spectrum f (␣) shrinks when the critical temperature is approached.