On the received coefficients of the analytical forms for deuteron wave function in coordinate representation in form r^(l+1)*exp(-A*r) for the modern realistic nucleon-nucleon potentials NijmI, NijmII, Nijm93, Reid93 and Argonne v18 are calculated charge GC(p), quadrupole GQ(p) and magnetic GM(p) deuteron form factors. An original dipole fit for proton and neutron isoscalar electric and magnetic form factors is used for calculations. Theoretical calculations of values of the deuteron form factor are compared with their experimental data of world collaborations (Bates, BLAST, Bonn, JLab, Mainz, Naval Research Lab, NIKHEF, Orsay, Saclay, SLAC, Stanford, VEPP3 and VEPP4) and reviews (Abbott, Boden, Garcon and Karpius). The change of the sign of form factors is in the pulse regions at 4.7-4.9 fm-1 for GC; at 12.8-14.7 fm-1 for GQ; at 6.3-8.1 and 11.4-12.2 fm-1 for GM. The theoretical values for the form factors GQ(0) and GM(0) were calculated at the boundary condition for the momentum at p2=0. Calculated positions of the zero for deuteron form factors are compared with values for other potential models. The peculiarities of parametrizations of deuteron form factors for theoretical approaches and fits of theoretical calculations to experimental values, which are described in the cited literature, are analyzed. Formulas for spherical and quadrupole form factors are writed, which are expressed in terms of the coefficients for the indicated deuteron wave function. At large values of momentum the asymptotics of the deuteron form factors are determined by the coefficients of the analytical forms for deuteron wave function, nucleon isoscalar form factors and the order of momentum in the denominator. The further use of deuteron form factors for obtaining polarization observables in processes involving deuteron as a projectile is discussed.