Convergence Analysis of Numerical Schemes for Liquid Crystals Based on the Invariant Energy Quadratization Method