2
$\begingroup$

I want to find the normalization constant for the following pdf

$ f(x_1,x_2,\ldots,x_n) = C\cdot\Pi_{m=1}^n\exp\left\{-\frac{1}{2a^2}(x_m-\sum_{l=1}^L\alpha_lx_{m-l})^2+y_m\sum_{l=0}^L\beta_lx_{m-l}\right\} $ where $n>L$ and $x_l=0$ for $l\leq0$. Namely, I want to find $C$ for which

$ \int_{\mathbb{R}^n}f(x_1,x_2,\ldots,x_n)dx_1\ldots dx_n = 1. $

Any idea?

  • 0
    @joriki Actually, the outer product is what bothers me.2012-12-15

1 Answers 1

1

From your comment, it seems that you may not be aware that the outer product is equivalent to a sum in the exponent. Using $\exp(a)\exp(b)=\exp(a+b)$ and transforming the indices as suggested in a comment, you can rewrite your function as

$ f(x_1,x_2,\ldots,x_n) = C\exp\left\{\sum_{m=1}^n\left(-\frac{1}{2a^2}\left(x_m-\sum_{k=m-L}^{m-1}\alpha_{m-k}x_k\right)^2+y_m\sum_{l=0}^L\beta_{m-k}x_k\right)\right\}\;. $

You can see from this representation that the matrix of the quadratic form (whose determinant you need) is a band matrix with bandwidth $2L+1$. The diagonal entries are constant ($1+\sum_i\alpha_i^2$) up to $m=n-L$, but decrease for greater $m$ (because successive $\alpha_i^2$ go missing) down to $1$ for $m=n$. Because of this asymmetry, I doubt that you'll find a nice closed form for the determinant.

  • 0
    Yes, my previous comment is not clear (I was aware that the outer product is equivalent to a sum in the exponent). Anyway, it really looks that I will not find a nice solution. Thx.2012-12-15