var(Y) = Sum[ Sum[ Ci * Cj * Covar(i,j), {i,1,2,...,n} ], {j,1,2,...,n} ]
where n is the number of coefficients of the FIR, Ci the i-th coefficient, and COV the covariance matrix defined as: Covar(i,j) = , ( is mean value).
note.1. If the stochastic components of the signal (i.e. noise) is stationary then ==
note 2. If the input signal X is also standardized (i.e. =0 and var(X)=1) then Covar(i,j) is the autocorrelation function ACF(|i-j|) which depends only on the absolute value of the difference of indices i and j.