Practical performances of controller design methods strongly depend on relevancy of identified models. Fractional order system models promise advantage of more accurate modeling of real systems. This study presents a discussion on utilization of two fundamental numerical solution methods of fractional calculus in identification problems of One Noninteger Order Plus Time Delay with one pole (NOPTD-I) models. The identification process is carried out by estimating parameters of a NOPTD-I type transfer function template so that the step response of the NOPTD-I model can sufficiently fit experimental step response data. In the study, step responses of NOPTD-I models are numerically calculated according to two fundamental methods, which are Mittag-Leffler (ML) function and Grunwald-Letnikov (GL) definition. Particle swarm optimization (PSO) algorithm is used to perform data fitting. Illustrative examples are presented to evaluate model parameter estimation performances of these two methods for synthetically generated noisy test data. An experimental study is conducted for modeling pitch rotor of TRMS to compare experimental performances.