We study the nonlinear discrete random Schrödinger equation \( i \frac{\partial }{\partial t} u+(\varepsilon _n \Delta +V) u+\delta |u|^{2 p} u =0 \quad \left( p \in \mathbb {N}^{+}\right) \) and discrete random wave equation \( u_{tt}+(\varepsilon _n \Delta +V) u+\delta u^{2 p+1} =0 \quad (p\in \mathbb {N}^+) \) on \(\mathbb {Z}^{d} \times [0, \infty )\) , where \(0< \delta <\varepsilon \ll 1,\) \(\Delta \) is the discrete Laplacian and V is the random potential, \(\left| \varepsilon _n\right| \le \varepsilon e^{-\varrho |n|}\) with \( \varrho >0\) . We fix the random potential V in a good set and choose the small amplitudes as parameters to prove a KAM theorem for finding linearly stable quasi-periodic solutions of nonlinear random Schrödinger equation and wave equation.