Based on the Longuet-Higgins wave model theory, the previews studies have shown that freak waves can be generated in finite space and time successfully. However, as to generating high nonlinear freak waves, the simulation results will be unrealistic. Therefore, a modified phase modulation method for simulating high nonlinear freak waves was developed. The surface elevations of some wave components at certain time and place are positive by modulating the corresponding random initial phases, then the total surface elevation at the focused point is enhanced and furthermore a freak wave event is generated. The new method can not only make the freak wave occur at certain time and place, but also make the simulated wave surface time series satisfy statistical properties of the realistic sea state and keep identical with the target wave spectrum. This numerical approach is of good precision and high efficiency by the comparisons of the simulated freak waves and the recorded freak waves. Based on the Longuet-Higgins wave model theory, the previews studies have shown that freak waves can be generated in finite space and time successfully. However, as to generate high nonlinear freak waves, the simulation results will be unrealistic. Therefore, a modified phase modulation method for simulating high nonlinear freak waves was developed. The surface elevations of some wave components at certain time and place are positive by modulating the corresponding random initial phases, then the total surface elevation at the focused point is enhanced and furthermore a freak wave event is generated. The new method can not only make the freak wave occur at certain time and place, but also make the simulated wave surface time series satisfaction statistical properties of the realistic sea state and keep the same with the target wave spectrum. of good precision and high efficiency by the comparisons of the simulated freak waves and the recorded freak waves.