试图逼近 - Woflram Alpha 给出但我得到?乙[ f( X) ]E[f(X)]乙[ f( X) ] ≈13√E[f(X)]≈13乙[ f( X) ] ≈ 0E[f(X)]≈0

机器算法验证 随机变量 期望值 近似 增量法 泰勒系列
2022-03-10 03:09:55

我想近似XN(μX,σX2)=N(0,1)f(x)=ex2E[f(X)]

Wolfram Alpha给出

E[f(X)]13.

使用泰勒展开方法,并注意到,我得到 \begin{align} E[f(X)] &\approx f(\mu_X) + \frac{f''(\ mu_X)}{2} \sigma_X^2 \\ & = f(0) + \frac{f''(0)}{2} \\ & = 1 + \frac{-2}{2} \\ & = 0. \end{对齐}f(0)=2

E[f(X)]f(μX)+f(μX)2σX2=f(0)+f(0)2=1+22=0.

为什么我的近似值与 Wolfram Alpha 结果不匹配?可以做些什么来修复它?

2个回答

这里不需要近似值。使用矩生成函数的属性,是标准正态,因此一个 df进行卡方,矩生成函数(对于 .) 然后注意 是定义,所以 我们可以在 R 中通过快速模拟(总是进行模拟检查的好主意):XX2MX2(t)=112tt<1/2

MX(t)=EetX
EeX2=MX²(1)=112(1)=13

 mean( exp(-rnorm(1E6)^2) )
[1] 0.5774847
 1/sqrt(3)
[1] 0.5773503

在评论中回答:

如果 𝑋 不是标准正态而是正态且均值 和方差您的方法是否仍然可以使用,或者它是否特定于标准正态分布的情况?𝜇𝑋𝜎𝑋2

它仍然可以使用。我不会提供完整的细节。首先,简单的案例然后标准正常,所以在上面的参数中,你得到参数代替用于 mgf(矩生成函数)。对于完全一般的情况,请参阅例如非中心卡方分布的矩生成函数 (MGF) 并从那里开始工作。XN(0,σ2)X=(σZ)2Zσ21

的确切值时,无需“近似” 让我们应用无意识统计学家定律 (LoTUS) 获得: E[f(X)]

E[f(X)]=+ex212πexp(x22) dx=20+12πexp(3x22) dx=20+12π16zez dz=223π0+ezz121 dz=13πΓ(12)=13ππ=13

希望这可以帮助。:)