Numerical studies are reported to support the idea to explain the particle production in late inflationary era based on the parametric resonance with a noise effect. Two nonperturbative renormalization group formulations are used to numerically calculate the time evolution of particle numbers. Firstly, the dynamical renormalization group (DRG) method is applied to sum up the secular contributions. Secondly, we derive an exact evolution equation of the particle number which turned out to be possible surprisingly owing to the presence of the noise effect. Our numerical results show a drastic enhancement of the particle production in agreement with earlier works qualitatively. A comparison is made of numerical results based on two methods. Our work provides nonperturbative and quantitative methods to evaluate the particle production for the inflationary cosmology.