The higher-order nonlinear response of random two-phase composites is investigated within the spectral representation theory. Both components are assumed to be weakly nonlinear with a displacement ( D )—electric field ( E ) relation of the form D =ϵ i E +χ i| E | 2 E , where ϵ i is the linear dielectric constant and χ i is the third-order nonlinear response of the ith component ( i=1,2). The effective field-dependent dielectric response of the composites is defined as ϵ ̃ e=ϵ e+χ e| E 0| 2+η e| E 0| 4 , where ϵ e , χ e and η e are the effective linear dielectric constant, the effective third-order and fifth-order nonlinear susceptibilities, respectively. By taking into account the correction to the nonlinear local fields, general expressions are derived for χ e and η e to the second order of χ i . We perform numerical calculations on Au/SiO 2 composites with three typical microstructures, described by the Maxwell–Garnett approximation, Bruggeman effective medium approximation and differential effective medium one. Our formulae are then generalized to systems of arbitrary nonlinearity. In the dilute limit, the formulae reduce to those in previous works, as expected. For large volume fractions, our theoretical results are in good agreement with simulation data in nonlinear random resistor networks.