The numerical methods of finding the optimal position of optical phase conjugator (OPC) and the optimal fiber dispersions are proposed, which are able to effectively compensate overall channels in $8{ imes}40$ Gbps WDM system with non zero-dispersion shifted fiber (NZ-DSF) as an optical fiber. And BER characteristics in the system with two induced optimal parameters are compared with those in the system with the currently used mid-span spectral inversion (MSSI) in order to confirm the availability of the proposed methods. It is confirmed that two optimal parameters depend on each other, but less related with the searching procedure. The methods proposed in this research will be expected to alternate with the method of making a symmetrical distribution of power and local dispersion in real optical link which is a serious problem but the condition in the case of applying the OPC into multi-channels WDM system.