正态分布的奇妙性质:为什么奇数阶中心矩(odd central moments)为零?
正态分布的奇妙性质:为什么奇数阶矩为零?
正态分布(Normal Distribution)是统计学中最常见的分布之一,它的钟形曲线几乎无处不在,从身高体重到测量误差,都能看到它的影子。除了均值和方差这两个核心参数,正态分布还有一个有趣的特性:它的奇数阶中心矩(odd central moments)全部为零,比如 ( E [ x − μ ] = 0 E[x - \mu] = 0 E[x−μ]=0 ) 和 ( E [ ( x − μ ) 3 ] = 0 E[(x - \mu)^3] = 0 E[(x−μ)3]=0 )。这到底是怎么回事?今天我们就来聊聊这个性质的由来、证明,以及它背后的意义。
什么是中心矩?
在探讨奇数阶矩之前,我们先明白什么是中心矩。中心矩是描述随机变量偏离其均值 ( μ \mu μ ) 的统计量,定义为:
μ k = E [ ( x − μ ) k ] \mu_k = E[(x - \mu)^k] μk=E[(x−μ)k]
- ( k = 1 k = 1 k=1 ):一阶中心矩,( E [ x − μ ] E[x - \mu] E[x−μ] ),总是等于零(因为 ( E [ x ] = μ E[x] = \mu E[x]=μ ))。
- ( k = 2 k = 2 k=2 ):二阶中心矩,就是方差 ( σ 2 \sigma^2 σ2 )。
- ( k = 3 k = 3 k=3 ):三阶中心矩,衡量分布的偏度(skewness)。
- ( k = 4 k = 4 k=4 ):四阶中心矩,与峰度(kurtosis)相关。
对于正态分布,我们关心的是这些矩的特性,尤其是奇数阶(( k = 1 , 3 , 5 , … k = 1, 3, 5, \dots k=1,3,5,… ))的中心矩。
正态分布的奇数阶矩为零
正态分布的概率密度函数为:
p ( x ∣ μ , σ 2 ) = 1 2 π σ 2 exp ( − ( x − μ ) 2 2 σ 2 ) p(x | \mu, \sigma^2) = \frac{1}{\sqrt{2\pi \sigma^2}} \exp\left( -\frac{(x - \mu)^2}{2\sigma^2} \right) p(x∣μ,σ2)=2πσ21exp(−2σ2(x−μ)2)
它的一个显著特点是对称性:以均值 ( μ \mu μ ) 为中心,左右两侧完全对称。这种对称性直观地暗示了一个结论:奇数阶中心矩为零。为什么呢?
通俗解释
想象你在玩一个对称的跷跷板,中间是均值 ( μ \mu μ )。你把 ( x − μ x - \mu x−μ )(偏离均值的距离)拿来计算奇数次方,比如 ( ( x − μ ) 3 (x - \mu)^3 (x−μ)3 )。因为正态分布是对称的,对于每一个正的 ( x − μ x - \mu x−μ )(比如 +2),总有一个对应的负的 ( − ( x − μ ) - (x - \mu) −(x−μ) )(比如 -2),它们的概率密度相等。奇数次方会保留正负号:
- ( ( + 2 ) 3 = 8 (+2)^3 = 8 (+2)3=8 )
- ( ( − 2 ) 3 = − 8 (-2)^3 = -8 (−2)3=−8 )
当你把这些值按概率加权平均时,正负项正好抵消,结果为零。这种对称性是奇数阶矩为零的直观原因。
数学证明
现在,让我们用数学来证明这个性质。以 ( k = 3 k = 3 k=3 )(三阶中心矩)为例,证明 ( E [ ( x − μ ) 3 ] = 0 E[(x - \mu)^3] = 0 E[(x−μ)3]=0 ),其他奇数阶的证明类似。
步骤 1:定义期望
对于 ( x ∼ N ( μ , σ 2 ) x \sim N(\mu, \sigma^2) x∼N(μ,σ2) ),三阶中心矩是:
E [ ( x − μ ) 3 ] = ∫ − ∞ ∞ ( x − μ ) 3 ⋅ 1 2 π σ 2 exp ( − ( x − μ ) 2 2 σ 2 ) d x E[(x - \mu)^3] = \int_{-\infty}^{\infty} (x - \mu)^3 \cdot \frac{1}{\sqrt{2\pi \sigma^2}} \exp\left( -\frac{(x - \mu)^2}{2\sigma^2} \right) \, dx E[(x−μ)3]=∫−∞∞(x−μ)3⋅2πσ21exp(−2σ2(x−μ)2)dx
步骤 2:变量替换
为了简化计算,令 ( z = x − μ σ z = \frac{x - \mu}{\sigma} z=σx−μ ),则 ( x − μ = σ z x - \mu = \sigma z x−μ=σz ),( d x = σ d z dx = \sigma \, dz dx=σdz ),且 ( z ∼ N ( 0 , 1 ) z \sim N(0, 1) z∼N(0,1))(标准正态分布),概率密度为:
ϕ ( z ) = 1 2 π e − z 2 2 \phi(z) = \frac{1}{\sqrt{2\pi}} e^{-\frac{z^2}{2}} ϕ(z)=2π1e−2z2
代入后:
E [ ( x − μ ) 3 ] = ∫ − ∞ ∞ ( σ z ) 3 ⋅ 1 2 π e − z 2 2 ⋅ σ d z E[(x - \mu)^3] = \int_{-\infty}^{\infty} (\sigma z)^3 \cdot \frac{1}{\sqrt{2\pi}} e^{-\frac{z^2}{2}} \cdot \sigma \, dz E[(x−μ)3]=∫−∞∞(σz)3⋅2π1e−2z2⋅σdz
= σ 3 ∫ − ∞ ∞ z 3 ⋅ 1 2 π e − z 2 2 d z = \sigma^3 \int_{-\infty}^{\infty} z^3 \cdot \frac{1}{\sqrt{2\pi}} e^{-\frac{z^2}{2}} \, dz =σ3∫−∞∞z3⋅2π1e−2z2dz
步骤 3:分析被积函数
被积函数是 ( f ( z ) = z 3 ⋅ 1 2 π e − z 2 2 f(z) = z^3 \cdot \frac{1}{\sqrt{2\pi}} e^{-\frac{z^2}{2}} f(z)=z3⋅2π1e−2z2 )。我们需要判断这个积分是否为零。关键在于 ( f ( z ) f(z) f(z) ) 的性质:
- ( z 3 z^3 z3 ) 是奇函数(odd function),因为 ( ( − z ) 3 = − z 3 (-z)^3 = -z^3 (−z)3=−z3 )。
- ( 1 2 π e − z 2 2 \frac{1}{\sqrt{2\pi}} e^{-\frac{z^2}{2}} 2π1e−2z2 ) 是偶函数(even function),因为 ( e − ( − z ) 2 2 = e − z 2 2 e^{-\frac{(-z)^2}{2}} = e^{-\frac{z^2}{2}} e−2(−z)2=e−2z2 )。
奇函数乘以偶函数的结果还是奇函数:
f ( − z ) = ( − z ) 3 ⋅ 1 2 π e − ( − z ) 2 2 = − z 3 ⋅ 1 2 π e − z 2 2 = − f ( z ) f(-z) = (-z)^3 \cdot \frac{1}{\sqrt{2\pi}} e^{-\frac{(-z)^2}{2}} = -z^3 \cdot \frac{1}{\sqrt{2\pi}} e^{-\frac{z^2}{2}} = -f(z) f(−z)=(−z)3⋅2π1e−2(−z)2=−z3⋅2π1e−2z2=−f(z)
步骤 4:奇函数积分的性质
对于任意奇函数 ( f ( z ) f(z) f(z) ),在对称区间 ( [ − ∞ , ∞ ] [-\infty, \infty] [−∞,∞] ) 上积分(假设积分收敛):
∫ − ∞ ∞ f ( z ) d z = ∫ − ∞ 0 f ( z ) d z + ∫ 0 ∞ f ( z ) d z \int_{-\infty}^{\infty} f(z) \, dz = \int_{-\infty}^{0} f(z) \, dz + \int_{0}^{\infty} f(z) \, dz ∫−∞∞f(z)dz=∫−∞0f(z)dz+∫0∞f(z)dz
令 ( u = − z u = -z u=−z ),则:
∫ − ∞ 0 f ( z ) d z = ∫ ∞ 0 f ( − u ) ( − d u ) = ∫ 0 ∞ − f ( u ) d u = − ∫ 0 ∞ f ( u ) d u \int_{-\infty}^{0} f(z) \, dz = \int_{\infty}^{0} f(-u) (-du) = \int_{0}^{\infty} -f(u) \, du = -\int_{0}^{\infty} f(u) \, du ∫−∞0f(z)dz=∫∞0f(−u)(−du)=∫0∞−f(u)du=−∫0∞f(u)du
所以:
∫ − ∞ ∞ f ( z ) d z = − ∫ 0 ∞ f ( z ) d z + ∫ 0 ∞ f ( z ) d z = 0 \int_{-\infty}^{\infty} f(z) \, dz = -\int_{0}^{\infty} f(z) \, dz + \int_{0}^{\infty} f(z) \, dz = 0 ∫−∞∞f(z)dz=−∫0∞f(z)dz+∫0∞f(z)dz=0
对于 ( f ( z ) = z 3 ⋅ 1 2 π e − z 2 2 f(z) = z^3 \cdot \frac{1}{\sqrt{2\pi}} e^{-\frac{z^2}{2}} f(z)=z3⋅2π1e−2z2 ),由于 ( e − z 2 2 e^{-\frac{z^2}{2}} e−2z2 ) 衰减很快,积分收敛,奇函数性质保证:
∫ − ∞ ∞ z 3 ⋅ 1 2 π e − z 2 2 d z = 0 \int_{-\infty}^{\infty} z^3 \cdot \frac{1}{\sqrt{2\pi}} e^{-\frac{z^2}{2}} \, dz = 0 ∫−∞∞z3⋅2π1e−2z2dz=0
因此:
E [ ( x − μ ) 3 ] = σ 3 ⋅ 0 = 0 E[(x - \mu)^3] = \sigma^3 \cdot 0 = 0 E[(x−μ)3]=σ3⋅0=0
推广到所有奇数阶
对于任意奇数 ( k = 2 n + 1 k = 2n + 1 k=2n+1 )(( n = 0 , 1 , 2 , … n = 0, 1, 2, \dots n=0,1,2,… )),( ( x − μ ) 2 n + 1 (x - \mu)^{2n+1} (x−μ)2n+1 ) 是奇函数,乘以偶函数 ( 1 2 π σ 2 e − ( x − μ ) 2 2 σ 2 \frac{1}{\sqrt{2\pi \sigma^2}} e^{-\frac{(x - \mu)^2}{2\sigma^2}} 2πσ21e−2σ2(x−μ)2 ) 后仍为奇函数,积分从 ( − ∞ -\infty −∞ ) 到 ( ∞ \infty ∞ ) 为零。所以,所有奇数阶中心矩都为零。
补充信息
1. 为什么偶数阶矩不为零?
偶数次方(如 ( ( x − μ ) 2 (x - \mu)^2 (x−μ)2 ) 或 ( ( x − μ ) 4 (x - \mu)^4 (x−μ)4 ))是偶函数,乘以偶函数后仍是偶函数,积分不会抵消。例如:
- ( E [ ( x − μ ) 2 ] = σ 2 E[(x - \mu)^2] = \sigma^2 E[(x−μ)2]=σ2 )(方差)
- ( E [ ( x − μ ) 4 ] = 3 σ 4 E[(x - \mu)^4] = 3\sigma^4 E[(x−μ)4]=3σ4 )(四阶矩)
这些值反映了分布的宽度和形状。
2. 偏度的含义
三阶中心矩 ( E [ ( x − μ ) 3 ] E[(x - \mu)^3] E[(x−μ)3] ) 与偏度相关。偏度为零意味着分布没有左偏或右偏,正态分布的对称性恰好保证了这一点。如果分布不对称(例如指数分布),奇数阶矩就不为零。
3. 在统计中的应用
奇数阶矩为零在统计推断中很有用。例如,在计算Fisher信息矩阵时:
I 12 = E [ x − μ σ 2 ⋅ ( − 1 2 σ 2 + ( x − μ ) 2 2 ( σ 2 ) 2 ) ] I_{12} = E\left[ \frac{x - \mu}{\sigma^2} \cdot \left( -\frac{1}{2\sigma^2} + \frac{(x - \mu)^2}{2(\sigma^2)^2} \right) \right] I12=E[σ2x−μ⋅(−2σ21+2(σ2)2(x−μ)2)]
因为 ( E [ x − μ ] = 0 E[x - \mu] = 0 E[x−μ]=0 ) 和 ( E [ ( x − μ ) 3 ] = 0 E[(x - \mu)^3] = 0 E[(x−μ)3]=0 ),交叉项为零,说明 ( μ \mu μ ) 和 ( σ 2 \sigma^2 σ2 ) 信息正交。
具体参考参考笔者的另一篇博客:Fisher信息矩阵(Fisher Information Matrix,简称FIM)
4. 其他分布呢?
并非所有分布的奇数阶矩都为零。例如:
- 指数分布 ( f ( x ) = λ e − λ x f(x) = \lambda e^{-\lambda x} f(x)=λe−λx )(( x ≥ 0 x \geq 0 x≥0 )),( E [ ( x − μ ) 3 ] ≠ 0 E[(x - \mu)^3] \neq 0 E[(x−μ)3]=0 ),因为它右偏。
- 对称的均匀分布也有奇数阶矩为零,但范围有限。
正态分布的无限对称支持和指数衰减共同造就了这个特性。
总结
正态分布的奇数阶中心矩为零源于其完美的对称性。奇函数与偶函数相乘后,积分在对称区间上抵消,数学上严谨地证明了这一点。这个性质不仅让正态分布更加“优雅”,还在统计估计、信息理论中简化了计算,比如保证参数间的正交性。下次看到正态分布的钟形曲线,不妨想想它隐藏的这些奇妙特性!
后记
2025年2月24日22点07分于上海,在Grok3大模型辅助下完成。
相关文章:
正态分布的奇妙性质:为什么奇数阶中心矩(odd central moments)为零?
正态分布的奇妙性质:为什么奇数阶矩为零? 正态分布(Normal Distribution)是统计学中最常见的分布之一,它的钟形曲线几乎无处不在,从身高体重到测量误差,都能看到它的影子。除了均值和方差这两个…...
架构——Nginx功能、职责、原理、配置示例、应用场景
以下是关于 Nginx 的功能、职责、原理、配置示例、应用场景及其高性能原因的详细说明: 一、Nginx 的核心功能 1. 静态资源服务 功能:直接返回静态文件(如 HTML、CSS、JS、图片、视频等)。配置示例:server {listen 80…...
涉密载体管控系统革新:RFID技术引领,信息安全新境界
行业背景 文件载体管控系统DW-S402是用于对各种SM载体进行有效管理的智能柜(智能管理系统),实现对载体的智能化、规范化、标准化管理,广泛应用于保密、机要单位以及企事业单位等有载体保管需求的行业。 随着信息化技术发展&…...
基于 SpringBoot 的 “电影交流平台小程序” 系统的设计与实现
大家好,今天要和大家聊的是一款基于 SpringBoot 的 “电影交流平台小程序” 系统的设计与实现。项目源码以及部署相关事宜请联系我,文末附上联系方式。 项目简介 基于 SpringBoot 的 “电影交流平台小程序” 系统设计与实现的主要使用者分为 管理员 和…...
【Rust中级教程】2.9. API设计原则之显然性(obvious) :文档与类型系统、语义化类型、使用“零大小”类型
喜欢的话别忘了点赞、收藏加关注哦(加关注即可阅读全文),对接下来的教程有兴趣的可以关注专栏。谢谢喵!(・ω・) 2.9.1. 文档与类型系统 用户可能不会完全理解API的所有规则和限制。所以你写的API应该让你…...
git branch
文章目录 1.简介2.格式3.选项4.示例参考文献 1.简介 git branch 用于管理分支,包括查看、创建、删除、重命名和关联。 git branch 是 Git 版本控制系统中用于管理分支的命令。分支是 Git 的核心功能之一,允许开发者在同一个代码库中并行开发不同的功能…...
【网络编程】广播和组播
数据包发送方式只有一个接受方,称为单播。如果同时发给局域网中的所有主机,称为广播。只有用户数据报(使用UDP协议)套接字才能广播: 广播地址以192.168.1.0 (255.255.255.0) 网段为例,最大的主机地址192.168.1.255代表该网段的广…...
运维Crontab面试题及参考答案
Crontab 文件的六个域分别是什么?顺序如何? Crontab 文件用于设置定时执行任务,其六个域及顺序从左到右依次为:分钟(Minute)、小时(Hour)、日期(Day of month)…...
Lecture 1 - AI Systems (Overview)
一、Machine Learning Approach标准机器学习流程 • Train ML algorithm(训练机器学习算法):基于收集的数据训练机器学习模型。 二、Machine Learning for Adaptation(适应性机器学习) 加入了数据更新和自动化的部分…...
Ansible 学习笔记
这里写自定义目录标题 基本架构文件结构安装查看版本 Ansible 配置相关文件主机清单写法 基本架构 Ansible 是基于Python实现的,默认使用22端口, 文件结构 安装 查看用什么语言写的用一下命令 查看版本 Ansible 配置相关文件 主机清单写法...
设计模式-结构型-代理模式
1. 代理模式概述 代理模式(Proxy Pattern) 是一种结构型设计模式,它允许通过代理对象来控制对目标对象的访问。代理模式主要用于以下场景: 控制对象访问:限制某些对象的访问权限,例如权限控制。 延迟实例…...
FCC CE SRRC MIC是什么意思?
1.FCC CE SRRC MIC是什么意思? 2.4000 GHz 至 2.4835 GHz:<33 dBm(FCC),<20 dBm(CE/SRRC/MIC) 5.150 GHz 至 5.250 GHz(CE:5.170 GHz 至 5.250 GHz)&a…...
springboot005学生心理咨询评估系统(源码+数据库+文档)
源码地址:学生心理咨询评估系统 文章目录 1.项目简介2.部分数据库结构与测试用例3.系统功能结构4.包含的文件列表(含论文)后台运行截图 1.项目简介 使用旧方法对学生心理咨询评估信息进行系统化管理已经不再让人们信赖了,把现…...
Apache Doris:一款高性能的实时数据仓库
Apache Doris 是一款基于 MPP 架构的高性能、实时分析型数据库。它以高效、简单和统一的特性著称,能够在亚秒级的时间内返回海量数据的查询结果。Doris 既能支持高并发的点查询场景,也能支持高吞吐的复杂分析场景。 Apache Doris 最初是百度广告报表业务…...
使用Vue-Flow创建一个流程图可视化节点坐标查询器
在开发中遇到这样一个需求,需要后端返回数据前端网页生成流程图,由于流程图使用了Vue-Flow,所以需要坐标来辅助后端生成数据。 首先引入方法并定义添加节点数据 const { updateEdge, addEdges, addNodes} useVueFlow() const add_nodes …...
面试基础--Java 集合框架详解
Java 集合框架详解:从 ArrayList 到 HashMap 的底层原理 引言 在 Java 开发中,集合框架(Collection Framework)是处理数据存储和操作的核心工具。无论是日常开发还是大厂面试,对集合框架的理解都是考察的重点之一。本…...
轻量级日志管理平台Grafana Loki
文章目录 轻量级日志管理平台Grafana Loki背景什么是Loki为什么使用 Grafana Loki?架构Log Storage Grafana部署使用基于 Docker Compose 安装 LokiMinIO K8s集群部署Loki采集Helm 部署方式和案例 参考 轻量级日志管理平台Grafana Loki 背景 在微服务以及云原生时…...
回文串
长度为偶数的串,重排连续字串变成回文串。 Problem - D - Codeforces 代码: #include <bits/stdc.h> #define fi first #define se second using namespace std; typedef long long LL; typedef pair<int,int> PII; typedef pair<LL,L…...
《跟李沐学 AI》AlexNet论文逐段精读学习心得 | PyTorch 深度学习实战
前一篇文章,使用 AlexNet 实现图片分类 | PyTorch 深度学习实战 本系列文章 GitHub Repo: https://github.com/hailiang-wang/pytorch-get-started 本篇文章内容来自于学习 9年后重读深度学习奠基作之一:AlexNet【下】【论文精读】】的心得。 《跟李沐…...
【电机控制器】FU6832S——持续更新
【电机控制器】FU6832S——持续更新 文章目录 [TOC](文章目录) 前言一、ADC二、UART三、PWM四、参考资料总结 前言 使用工具: 提示:以下是本篇文章正文内容,下面案例可供参考 一、ADC 二、UART 三、PWM 四、参考资料 总结 本文仅仅简…...
Flutter屏幕适配终极方案:flutter_screenutil深度解析
在跨平台应用开发中,屏幕适配始终是开发者面临的核心挑战。Flutter虽然自带响应式布局体系,但面对复杂的设计稿标注时,手动计算比例效率低下。今天我们将深度解析目前Flutter社区最受欢迎的屏幕适配方案——flutter_screenutil,手…...
计算机视觉算法实战——产品分拣(主页有源码)
✨个人主页欢迎您的访问 ✨期待您的三连 ✨ ✨个人主页欢迎您的访问 ✨期待您的三连 ✨ ✨个人主页欢迎您的访问 ✨期待您的三连✨ 1. 领域简介✨✨ 产品分拣是工业自动化和物流领域的核心技术,旨在通过机器视觉系统对传送带上的物品进行快速识别、定位和分类&a…...
可视化报表
根据你的需求,以下是一些可以实现报表可视化的开源项目,这些项目提供了类似阿里巴巴 FBI 报表的功能,支持数据可视化、报表设计、仪表盘和大屏展示等功能: 1. DataEase DataEase 是一个开源的 BI 工具,帮助用户快速分…...
基于模块联邦的微前端架构:重构大型前端应用的模块化边界
引言:企业级前端的模块化困境 字节跳动广告系统采用Webpack 5模块联邦后,主应用构建时间从14分钟降至38秒,微应用独立发布频率提升至每天50次。在动态加载机制下,首屏资源加载体积减少79%,跨团队组件复用率达到92%。其…...
Android之图片保存相册及分享图片
文章目录 前言一、效果图二、实现步骤1.引入依赖库2.二维码生成3.布局转图片保存或者分享 总结 前言 其实现在很多分享都是我们自定义的,更多的是在界面加了很多东西,然后把整个界面转成图片保存相册和分享,而且现在分享都不需要第三方&…...
Linux放行端口
8080这个端口测试看telnet是不通的,您服务器内是否有对应的业务监听了这个端口呢?您到服务器内执行下: netstat -nltp |grep 8080 同时服务器内执行下: systemctl status firewalld iptables -nL 截图反馈下,我看下防火…...
Spring Boot延迟执行实现
说明:本文介绍如何在Spring Boot项目中,延迟执行某方法,及讨论延迟执行方法的是事务问题。 搭建Demo 首先,创建一个Spring Boot项目,pom.xml如下: <?xml version"1.0" encoding"UTF-…...
npm i 失败权限问题
安装完node之后, 测试全局安装一个最常用的 express 模块进行测试 失败,但是用管理员权限打开cmd 安装就成功。 报错如下: npm ERR! If you believe this might be a permissions issue, please double-check the npm ERR! permissions of the file and …...
uniapp 微信小程序打包之后vendor.js 主包体积太大,解决办法,“subPackages“:true设置不生效
现在是打包的时候,vendor.js 的内容全部打到了主包里面, 说一下我的方法: 1. 通过发行 小程序打包 这样打包的体积是最小的,打包之后打开微信开发工具,然后再上传 2.manifest.json,在“mp-weixin”里添加代码 "…...
23.2、云计算安全机制与案例分析
目录 云计算安全保护机制与技术方案云计算安全保护机制与技术方案常见云计算网络安全机制云计算安全管理与运维云计算安全综合应用案例分析 - 阿里云云计算安全综合应用案例分析 - 腾讯云云计算安全综合应用案例分析 - 华为云 云计算安全保护机制与技术方案 首先针对云计算&am…...
