How is the minimum variance unbiased estimator derived for some RV $X\sim \mathcal{N}(0,\sigma^2)$, given $\mathbf{Y}=X\cdot \mathrm{ones}_{N\times 1}+\mathbf{w}$, where $\mathbf{w}$ is $N$-dimensional multivariate-normal with mean $[0, \dots, 0]^T$ and covariance matrix $\mathbf{\Sigma}$, and is independent of $X$?
In words, what is the best estimator for a signal in correlated noise, what is the distribution?
I know the answer is well-known and buried somewhere in Kay's Estimation Theory text but I cannot find it at the present. (am still looking)
If $\mathbf{\Sigma}$ were diagonal (that is, uncorrelated noise), we know that this estimator is some linear combination of the $\mathbf{Y}$ and it's Gaussian distributed about the signal with variance $\sum_{n=1}^N \frac{1}{\mathbf{\Sigma}_{n,n}}$ (cf. Maximal-ratio combining)