To derive the Yule–Walker equations, we multiply the AR model by $x[n-m]$ and take expectations: (m represents the lag (or time delay)) $$ E[x[n]x[n-m]] = \sum_{k=1}^{p} \phi_k E[x[n-k]x[n-m]] + E[e[n]x[n-m]] $$ Since $e[n]$ is white noise and uncorrelated with past samples: $$ E[e[n]x[n-m]] = 0 \quad \text{for } m \ge 1 $$ Using the definition of autocorrelation, we obtain: $$ R[m] = \sum_{k=1}^{p} \phi_k R[m-k], \quad m = 1,2,\dots,p $$ These equations can be written in matrix form as: $$ \underbrace{ \begin{bmatrix} R[0] & R[1] & \cdots & R[p-1] \\ R[1] & R[0] & \cdots & R[p-2] \\ \vdots & \vdots & \ddots & \vdots \\ R[p-1] & R[p-2] & \cdots & R[0] \end{bmatrix} }_{\text{Known autocorrelation matrix}} \; \underbrace{ \begin{bmatrix} \phi_1\\ \phi_2\\ \vdots\\ \phi_p \end{bmatrix} }_{\text{Unknown AR coefficients}} = \underbrace{ \begin{bmatrix} R[1]\\ R[2]\\ \vdots\\ R[p] \end{bmatrix} }_{\text{Kno...