In this article, the authors consider the optimal portfolio on tracking the expected wealth process with liquidity constraints. The constrained optimal portfolio is first formulated as minimizing the cumulate variance between the wealth process and the expected wealth process. Then, the dynamic programming methodology is applied to reduce the whole problem to solving the Hamilton-Jacobi-Bellman equation coupled with the liquidity constraint, and the method of Lagrange multiplier is applied to handle the constraint. Finally, a numerical method is proposed to solve the constrained HJB equation and the constrained optimal strategy. Especially, the explicit solution to this optimal problem is derived when there is no liquidity constraint.