This article considers the tracking control of unknown nonlinear nonaffine repetitive discrete-time multi-input multi-output systems. Two data-driven iterative learning control (ILC) schemes are designed based on two equivalent dynamic linearization data models of an unknown ideal learning controller, which exists theoretically in the iteration domain. The two control schemes provide ways of selecting learning controllers based on the complexity of the controlled nonlinear systems. The learning control gain matrixes of the two learning controllers are optimized through the steepest descent method using only the measured input-output data of the nonlinear systems. The proposed ILC approaches are pure data-driven since no model information of the controlled systems is involved. The stability and convergence of the proposed ILC approaches are rigorously analyzed under reasonable conditions. Numerical simulation and an experiment based on a Gantry-type linear motor drive system are conducted to verify the effectiveness of the proposed data-driven ILC approaches.