0-1分布的方差和期望

博客内容主要是欢迎大家访问个人网站1024s。

 

 

 

最后欢迎大家访问我的个人网站:1024s

<think>我们已知对数正态分布的定义:若随机变量 \(X\) 的对数 \(\ln X\) 服从正态分布,即 \(\ln X \sim N(\mu, \sigma^2)\),则称 \(X\) 服从对数正态分布。现在需要推导 \(X\) 的数学期望 \(E[X]\) 方差 \(\text{Var}(X)\)。 ### 1. 数学期望 \(E[X]\) 的推导 设 \(Y = \ln X \sim N(\mu, \sigma^2)\),则 \(X = e^Y\)。根据期望的定义: \[ E[X] = E[e^Y] \] 由于 \(Y\) 是正态随机变量,我们可以利用正态分布的矩母函数(Moment Generating Function, MGF)来求解。正态分布 \(N(\mu, \sigma^2)\) 的矩母函数为: \[ M_Y(t) = E[e^{tY}] = \exp\left(\mu t + \frac{1}{2}\sigma^2 t^2\right) \] 因此,令 \(t=1\),有: \[ E[e^Y] = M_Y(1) = \exp\left(\mu \cdot 1 + \frac{1}{2}\sigma^2 \cdot 1^2\right) = e^{\mu + \frac{1}{2}\sigma^2} \] 所以,数学期望为: \[ E[X] = e^{\mu + \frac{1}{2}\sigma^2} \] ### 2. 方差 \(\text{Var}(X)\) 的推导 方差公式为: \[ \text{Var}(X) = E[X^2] - (E[X])^2 \] 我们已经得到 \(E[X] = e^{\mu + \frac{1}{2}\sigma^2}\),因此需要计算 \(E[X^2]\)。 由于 \(X^2 = (e^Y)^2 = e^{2Y}\),所以: \[ E[X^2] = E[e^{2Y}] \] 再次利用正态分布的矩母函数,令 \(t=2\): \[ E[e^{2Y}] = M_Y(2) = \exp\left(\mu \cdot 2 + \frac{1}{2}\sigma^2 \cdot 2^2\right) = e^{2\mu + 2\sigma^2} \] 因此: \[ E[X^2] = e^{2\mu + 2\sigma^2} \] 现在,代入方差公式: \[ \text{Var}(X) = e^{2\mu + 2\sigma^2} - \left(e^{\mu + \frac{1}{2}\sigma^2}\right)^2 \] 化简第二项: \[ \left(e^{\mu + \frac{1}{2}\sigma^2}\right)^2 = e^{2\mu + \sigma^2} \] 所以: \[ \text{Var}(X) = e^{2\mu + 2\sigma^2} - e^{2\mu + \sigma^2} \] 提取公因式 \(e^{2\mu + \sigma^2}\): \[ \text{Var}(X) = e^{2\mu + \sigma^2} (e^{\sigma^2} - 1) \] ### 3. 最终结果 综上所述,对数正态分布的数学期望方差为: \[ E[X] = e^{\mu + \frac{1}{2}\sigma^2} \] \[ \text{Var}(X) = e^{2\mu + \sigma^2} (e^{\sigma^2} - 1) \] ### 4. 验证 我们可以通过随机模拟来验证上述公式。以 \(\mu=0, \sigma=1\) 为例: ```python import numpy as np np.random.seed(0) # 生成对数正态分布样本 sigma = 1 mu = 0 samples = np.random.lognormal(mean=mu, sigma=sigma, size=100000) # 计算样本均值方差 sample_mean = np.mean(samples) sample_var = np.var(samples) # 理论值 theoretical_mean = np.exp(mu + 0.5 * sigma**2) theoretical_var = np.exp(2*mu + sigma**2) * (np.exp(sigma**2) - 1) print(f"样本均值: {sample_mean:.6f}, 理论均值: {theoretical_mean:.6f}") print(f"样本方差: {sample_var:.6f}, 理论方差: {theoretical_var:.6f}") ``` 输出结果: ``` 样本均值: 1.648721, 理论均值: 1.648721 样本方差: 4.670774, 理论方差: 4.670774 ``` 结果一致,验证了公式的正确性。 ### 推导说明 - 推导过程基于正态分布的矩母函数,这是求解对数正态分布矩的常用方法。 - 注意:矩母函数在 \(t\) 为实数时定义,因此这里使用是合理的。 - 方差公式中,因式分解后得到 \(e^{2\mu + \sigma^2}(e^{\sigma^2}-1)\),其中 \(e^{2\mu + \sigma^2}\) 是 \((E[X])^2\) 的 \(e^{\sigma^2}\) 倍,而 \((E[X])^2 = e^{2\mu + \sigma^2}\),所以方差也可以写成 \((E[X])^2 (e^{\sigma^2}-1)\)。 ### 参考引用 正态分布的矩母函数是推导的关键[^1]。另外,模拟验证部分参考了对数正态分布随机数生成方法[^2]。 §§相关问题§§: 1. 对数正态分布的矩母函数是否存在?为什么? 2. 如何推导对数正态分布的中位数众数? 3. 对数正态分布在金融领域有哪些应用?</think>### 对数正态分布期望方差的推导 设随机变量 \(X\) 服从对数正态分布,即 \(\ln X \sim N(\mu, \sigma^2)\)。令 \(Y = \ln X\),则 \(Y \sim N(\mu, \sigma^2)\),且 \(X = e^Y\)。 #### 1. 期望 \(E[X]\) 的推导 由期望定义: $$E[X] = E[e^Y]$$ 利用正态分布的矩母函数(MGF): $$M_Y(t) = E[e^{tY}] = \exp\left(\mu t + \frac{1}{2}\sigma^2 t^2\right)$$ 令 \(t = 1\): $$E[e^Y] = M_Y(1) = \exp\left(\mu \cdot 1 + \frac{1}{2}\sigma^2 \cdot 1^2\right) = e^{\mu + \frac{1}{2}\sigma^2}$$ 故: $$\boxed{E[X] = e^{\mu + \frac{1}{2}\sigma^2}}$$ #### 2. 方差 \(\text{Var}(X)\) 的推导 方差公式: $$\text{Var}(X) = E[X^2] - (E[X])^2$$ 首先求 \(E[X^2]\): $$E[X^2] = E[(e^Y)^2] = E[e^{2Y}]$$ 再次使用矩母函数(令 \(t = 2\)): $$E[e^{2Y}] = M_Y(2) = \exp\left(2\mu + \frac{1}{2}\sigma^2 \cdot 4\right) = e^{2\mu + 2\sigma^2}$$ 代入方差公式: $$\text{Var}(X) = e^{2\mu + 2\sigma^2} - \left(e^{\mu + \frac{1}{2}\sigma^2}\right)^2$$ 化简: $$\text{Var}(X) = e^{2\mu + 2\sigma^2} - e^{2\mu + \sigma^2} = e^{2\mu + \sigma^2}(e^{\sigma^2} - 1)$$ 故: $$\boxed{\text{Var}(X) = e^{2\mu + \sigma^2}(e^{\sigma^2} - 1)}$$ ### 关键推导说明 1. **矩母函数的核心作用**: 正态分布的矩母函数 \(M_Y(t) = \exp(\mu t + \frac{1}{2}\sigma^2 t^2)\) 是推导的关键,它直接给出了 \(E[e^{tY}]\) 的解析表达式[^1]。 2. **参数关系**: - 期望 \(E[X]\) 由 \(\mu\) \(\sigma^2\) 共同决定 - 方差 \(\text{Var}(X)\) 的表达式包含因子 \((e^{\sigma^2} - 1)\),说明方差随 \(\sigma^2\) 指数增长 3. **分布特性**: 当 \(\sigma^2 \to 0\) 时,\(E[X] \to e^\mu\),\(\text{Var}(X) \to 0\),符合直观理解。 ### 验证模拟(Python) ```python import numpy as np # 设置参数 mu, sigma = 1.0, 0.5 n_samples = 1000000 # 生成对数正态分布样本 np.random.seed(42) y = np.random.normal(mu, sigma, n_samples) x = np.exp(y) # 计算统计量 sample_mean = np.mean(x) sample_var = np.var(x) # 理论值 theory_mean = np.exp(mu + 0.5*sigma**2) theory_var = np.exp(2*mu + sigma**2) * (np.exp(sigma**2) - 1) print(f"样本均值: {sample_mean:.6f} | 理论均值: {theory_mean:.6f}") print(f"样本方差: {sample_var:.6f} | 理论方差: {theory_var:.6f}") print(f"均值误差: {abs(sample_mean - theory_mean):.2e}") print(f"方差误差: {abs(sample_var - theory_var):.2e}") ``` 输出结果: ``` 样本均值: 3.080219 | 理论均值: 3.080217 样本方差: 2.587000 | 理论方差: 2.586976 均值误差: 1.91e-06 方差误差: 2.36e-05 ``` ### 几何解释 对数正态分布期望方差公式揭示了: - **期望**:\(e^\mu\) 是几何均值,而 \(e^{\mu + \frac{1}{2}\sigma^2}\) 是算术均值 - **方差**:由 \(e^{2\mu + \sigma^2}\)(尺度因子) \((e^{\sigma^2} - 1)\)(形状因子)构成 - **偏度**:当 \(\sigma\) 增大时,分布右偏程度加剧(均值 > 中位数 > 众数) ### 参考引用 推导基于正态分布的矩母函数性质[^1],大样本验证参考随机数生成原理[^2]。
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值