It has been resulted that Lagrange multiplier method with statistical approach was superior to traditional harmonic balance method in identifying the nonlinear loudspeaker parameters when output signals were contaminated with Gaussian random noise. We have known that the displacement-dependent characteristic values of nonlinear parameters identified by traditional harmonic balance method were estimated less than original values by the increase of output noise and the stiffness coefficients were very sensitive to output noise. Also, by the sensitivity analysis we have verified that the harmonic distortions in acoustic radiation was mainly due to nonlinearity of force factor caused by uneven magnetic fields and that reducing the nonlinearity of damping coefficients were very effective for improving second harmonic distrotion of acoustic radiation.