Symbolic expectations and double integrals

Green diod
123 ●9 ●13 ●21

I'd like to compute the following expectation ( $U$ and $V$ are independent and normally distributed or Gaussian) $a_{k,p}=E(|U|^p|U+\sqrt{k-1}V|^p)$

Is there a way to directly compute those expectations in Sage ?

If I write it down, I obtain a double integral which I already tried to compute with maxima like this (ok here, you can only see the starting point with the inner integral with respect to $v$ first)

load(distrib);
n(x):=pdf_normal(x, 0, 1);
inner_integral: integrate(abs(u+sqrt(k-1)*v)^p*n(v), v, minf, inf);

but to no avail.

NB: The absolute central moment $a_{1,p/2}$ can be easily obtained with maxima with

ratsimp(integrate(abs(x)^p*n(x), x, minf, inf));

but no 'simple' expression as for the aforementioned double integral.

Any hint?

Comments

Since Sage uses Maxima for all of its symbolic integration, I find it unlikely we'll do better. We do have numerical routines of various types, but I don't think this is what you are looking for.

kcrisman ( 13 years ago )

@kcrisman I thought that sometimes Sage could use sympy or some other CAS for symbolic integration. But then, I'm stuck here then.

Green diod ( 13 years ago )

add a comment

absm(p)=integral(abs(x)^p*n(x), x, -infinity, infinity) absm p |--> 2^(1/2*p)*gamma(1/2*p + 1/2)/sqrt(pi) a(k,p)=integral(integral(abs(u)^p*abs(u+sqrt(k-1)*v)^p*n(v)*n(u), v, -infinity, infinity), u, -infinity, infinity) a(k,p) 1/2*integrate(abs(u)^p*integrate(abs(sqrt(k - 1)*v + u)^p*e^(-1/2*u^2 - 1/2*v^2), v, -Infinity, +Infinity), u, -Infinity, +Infinity)/pi ab(p)=a(1,p/2).simplify() ab p |--> 2^(1/2*p)*gamma(1/2*p + 1/2)/sqrt(pi) bool(absm == ab) True

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Symbolic expectations and double integrals

Comments

1 Answer

Your Answer

Question Tools

Stats

Related questions

Symbolic expectations and double integrals savecancel

Comments

1 Answer