简介本资源是一份系统、精炼的概率论与数理统计核心公式汇编文档面向高校数学、统计、计算机及相关专业本科生与考研学生解决课程复习、考前速记与公式快速检索等实际需求。文档以清晰逻辑分章组织覆盖随机事件与概率基础含加法/乘法原理、古典与几何概型、事件关系与运算、条件概率与贝叶斯公式、独立性判定、全概率公式、离散与连续型随机变量分布律及密度函数、分布函数性质等关键内容第二章起已延伸至随机变量分布层面具备完整知识链路。资源为单文件Word文档.doc格式体积轻量仅1.41MB便于移动查阅与打印背诵。目前已有1052人学习下载内容排版规范、公式书写严谨、符号标注统一特别适合作为教材补充笔记、期末冲刺手册或机器学习前置数学工具包使用。1. 这份《概率论与数理统计公式整理.doc》不是复习提纲而是你调参、推导、debug时伸手就抄的“公式扳手”你有没有过这种时刻在写贝叶斯网络的后验采样代码时卡在联合分布分解上调试一个异常检测模块发现 p 值校准总对不上理论阈值甚至只是想快速确认 t 分布自由度 n−1 是怎么来的却要翻三本教材、查五次维基、再比对两篇论文附录——最后发现它们写的还不一致这份名为《概率论与数理统计公式整理.doc》的文档本质不是给期末划重点用的它是一线数据工程师、算法研究员、量化策略开发人员在真实项目推进中反复打开、高频粘贴、直接套用的“可执行公式集”。它不讲“什么是期望”但会明确写出当 X∼Gamma(α,β)尺度参数化时E[X]αβVar[X]αβ²并标注「注意此处 β 是尺度参数非率参数若你用 scipy.stats.gamma(aα, scaleβ)则完全匹配若用 rateβ 的写法请先转为 scale1/β」。它不罗列全部定理但每条公式都带使用边界中心极限定理适用条件独立同分布 有限二阶矩、大数定律弱/强形式对收敛性的不同要求、χ² 拟合优度检验中期望频数 ≥5 的硬约束。它面向的是那个正在改模型评估脚本、正在重写 AB 实验置信区间计算逻辑、正在把旧版 R 代码迁移到 Python 的你——不是学生是持工具上岗的实践者。2. 公式整理不是抄书而是按“使用动线”重构知识图谱从随机变量出发到假设检验收口2.1 为什么必须放弃“章节顺序整理”而采用“操作流驱动”的结构设计传统教材按“随机事件→随机变量→数字特征→极限定理→统计推断”线性展开这符合教学逻辑但严重违背工程实践路径。真实场景中你几乎不会从“样本空间”开始建模。更常见的是你拿到一批用户停留时长数据直方图右偏 → 思考是否用对数正态或 Gamma 建模 → 需立刻查 Gamma 的 PDF 形式、MLE 估计量、与指数分布的关系你在做多臂老虎机策略需计算 Beta-Binomial 共轭更新 → 要同时看到 Beta 先验参数更新规则、后验预测分布、以及如何从 posterior sampling 得到 Thompson Sampling 动作你部署了一个新特征上线监控需实时计算滑动窗口内均值的置信区间 → 必须能 3 秒内定位到“单样本 t 区间”公式、自由度修正逻辑、以及当 n30 时为何不能用 z 替代。因此这份 .doc 的骨架不是按教材章回而是按“你面对什么数据 → 你想做什么 → 你需要哪条公式 → 它依赖哪些前提 → 它输出什么可落地结果”四层动线组织。例如“连续型随机变量”一节下不列所有分布只放四类高频实战分布正态系Normal / Log-Normal / Truncated Normal含 PDF/CMF、矩生成函数、截断后的归一化常数计算式计数系Poisson / Negative Binomial / Zero-Inflated Poisson含过离散诊断方法、NB 的两种参数化r/p vs μ/α、ZIP 的 EM 初始化技巧等待系Exponential / Gamma / Weibull含失效率函数、生存函数、与泊松过程的双向映射共轭系Beta / Dirichlet / Inverse-Gamma含先验-似然-后验完整更新链、超参数物理意义如 Beta(α,β) 中 α 可解释为“成功先验计数”β 为“失败先验计数”。提示所有分布表格首列必含「Scipy/PyTorch/TensorFlow 对应实现名及关键参数名」。例如 Gamma 行写明scipy.stats.gamma(ashape, scalescale)torch.distributions.Gamma(concentrationshape, rate1/scale)tfp.distributions.Gamma(concentrationshape, scalescale)。避免因参数命名差异导致的 silent bug。2.2 公式不是孤立存在必须绑定“可验证数值示例”和“典型误用反例”光写 E[X] ∫x f(x)dx 没用。真正救命的是数值锚点当 X∼N(μ2,σ²4)则 P(X0) Φ((0−2)/2) Φ(−1) ≈ 0.1587。你复制该行进 Pythonfrom scipy.stats import norm; print(norm.cdf(0, loc2, scale2))输出0.15865525393145707误差 1e−12 —— 这就是可信公式的刻度。反例警示写完“若 X,Y 独立则 Var(XY)Var(X)Var(Y)”后紧跟一句“⚠️ 但若仅知 Cov(X,Y)0不相关不能推出独立反例X∼U(−1,1), YX²则 Cov(X,Y)0但 Y 完全由 X 决定显然不独立。”我们为每个核心公式配一个最小可运行验证块Python NumPy/SciPy并强制要求输入值取自真实业务范围如用户留存率不用 0.5而用 0.03~0.12订单金额不用 100而用 log-normal(μ8.2, σ0.9) 模拟输出包含解析解、数值解、相对误差三列误差 1e−10 则标红并触发检查常因浮点精度、积分限截断、特殊函数实现差异引起所有反例用assert not显式断言如assert not (np.corrcoef(X, Y)[0,1] 0 and np.allclose(Y, X**2))让误用当场暴露。# 示例验证 Beta-Binomial 共轭更新用于 A/B 实验贝叶斯分析 import numpy as np from scipy.stats import beta, binom # 先验Beta(α2, β8) → 期望转化率 2/(28)0.2反映较保守的初始信念 alpha_prior, beta_prior 2, 8 # 观测1000 次曝光230 次转化即成功数 s230, 失败数 f770 s, f 230, 770 # 后验参数α_post α_prior s, β_post β_prior f alpha_post alpha_prior s # 232 beta_post beta_prior f # 778 # 验证后验均值应等于 (α_prior s) / (α_prior β_prior s f) posterior_mean_analytic alpha_post / (alpha_post beta_post) # 232/1010 ≈ 0.2297 posterior_mean_numeric beta.mean(alpha_post, beta_post) # scipy 计算值 print(f解析均值: {posterior_mean_analytic:.6f}) print(f数值均值: {posterior_mean_numeric:.6f}) print(f绝对误差: {abs(posterior_mean_analytic - posterior_mean_numeric):.2e}) # 输出解析均值: 0.229703数值均值: 0.229703绝对误差: 1.11e-16 → 可信这段代码的价值不在“能跑”而在它把抽象的共轭更新具象成两个整数加法、一次除法、一次 scipy 调用并用误差量化告诉你这个公式在此参数域下是数值稳定的。当你下次在 AB 平台后台写贝叶斯置信区间时可以直接复用alpha_post alpha_prior convs这一行而不是再推导一遍。3. 把公式变成“可调试单元”每个公式块自带参数说明、边界检查与降级方案3.1 公式不是黑匣子必须拆解为“输入参数 → 计算步骤 → 输出含义 → 失效条件”以最常用的威尔科克森秩和检验Mann-Whitney U Test为例教科书只写 U n₁n₂ n₁(n₁1)/2 − R₁。但这对工程师毫无指导意义。我们的整理方式是字段内容说明用途检验两独立样本是否来自相同分布非参数不假设正态适用于转化率、停留时长等常不服从正态的指标输入样本 An₁ 个观测、样本 Bn₂ 个观测要求 n₁,n₂ ≥ 5小样本需查精确表若 n₁3,n₂4拒绝使用正态近似应调用scipy.stats.mannwhitneyu(..., methodexact)计算步骤① 合并两样本并排序得秩次② 对 A 中每个值统计 B 中小于它的个数之和 → U₁③ U min(U₁, U₂)步骤②可用np.searchsorted向量化实现避免双循环正态近似公式Z (U − μᵤ) / σᵤ其中 μᵤ n₁n₂/2σᵤ √[n₁n₂(n₁n₂1)/12]关键边界当 n₁,n₂ 20 且无大量结ties时可用若有 10% 数据为结必须用correctionTrue连续性校正输出解读p 值 0.05 仅表明分布位置不同P(XY)≠0.5不意味均值差异显著常见翻车用 U 检验结论替代 t 检验结论误判“提升均值”这种结构让公式从静态文本变为动态决策树。你不再问“U 公式是什么”而是问“我当前 n₁15,n₂18数据有 3 个重复值该用哪个分支”——答案就在表格里。3.2 每个公式必须配套“失效降级方案”否则就是埋雷没有万能公式。当条件不满足时硬套只会产出错误结论。我们为每个主流公式预设三条退路参数预警在公式旁用⚠️标出硬性阈值。如卡方检验中写“期望频数 5 的格子数占比 20% → 禁用卡方改用 Fisher 精确检验”。自动降级代码提供带 guard 的封装函数内部自动判断并切换算法def safe_chi2_test(observed, alpha0.05): 安全卡方检验自动检测期望频数不足时降级为 Fisher from scipy.stats import chi2_contingency, fisher_exact import numpy as np # 计算期望频数 chi2_stat, p_chi2, dof, expected chi2_contingency(observed) # 检查期望频数小于 5 的格子数占比 low_exp_count np.sum(expected 5) total_cells observed.size if low_exp_count / total_cells 0.2: print(f⚠️ 期望频数过低{low_exp_count}/{total_cells}降级为 Fisher 精确检验) # Fisher 只支持 2x2故对更大表先做合并此处省略合并逻辑实际需实现 if observed.shape (2, 2): _, p_fisher fisher_exact(observed) return {method: Fisher, p_value: p_fisher} else: raise ValueError(Fisher 仅支持 2x2 表请手动合并稀疏列/行) return {method: Chi2, p_value: p_chi2, statistic: chi2_stat} # 使用传入你的混淆矩阵即可无需人工判断 result safe_chi2_test(np.array([[12, 8], [3, 17]])) print(result) # 自动选 Fisher 或 Chi2物理意义兜底当所有统计方法失效如 n1回归业务直觉。例如“单次实验无法做假设检验应结合历史基线波动范围如过去 30 天转化率标准差为 ±0.008判断本次 0.023 是否超出 3σ 异常”。4. 避坑那些让模型上线前夜崩溃的“公式幻觉”与参数陷阱4.1 现象t 检验 p 值在 Python 和 Excel 里差 10 倍原因Excel 的T.TEST()默认计算双侧 p 值而部分 Python 教程错误地用scipy.stats.ttest_1samp(...).pvalue也是双侧但业务需求常是单侧检验如“新策略是否显著提升转化率”只关心上侧。更致命的是Excel 的T.TEST(array1,array2,1,2)中第三个参数1表示单侧第四个2表示双样本异方差——而很多人复制 Excel 公式时漏掉参数或误读文档。解决统一用scipy.stats.ttest_ind(a, b, alternativegreater)显式指定greater上侧或less下侧并注释“此处 alternativegreater 对应 H₁: μ₁ μ₂拒绝域为 t t_{α,df}”。4.2 现象用np.cov(X,Y)计算协方差结果与理论值偏差 5%原因np.cov默认biasFalse无偏估计除以 n−1但某些教材公式写的是有偏估计除以 n。当 n 较小时如 n10n/(n−1)1.11直接导致 11% 系统性偏差。解决在公式旁强制标注分母。例如协方差定义写为Cov(X,Y) \frac{1}{n-1}\sum_{i1}^n (x_i-\bar{x})(y_i-\bar{y}) \quad \text{样本无偏估计numpy 默认}Cov(X,Y) \frac{1}{n}\sum_{i1}^n (x_i-\bar{x})(y_i-\bar{y}) \quad \text{总体矩估计需 np.cov(..., biasTrue)}并在代码块中显式写出np.cov(X, Y, biasFalse)[0,1]。4.3 现象贝叶斯后验预测分布采样结果始终偏高AB 实验结论矛盾原因混淆了 Gamma 分布的两种参数化。某团队用scipy.stats.gamma(aα, scaleβ)建模但先验设定时参考了某论文的 “Gamma(α,β) with rate β”误将 rate 当作 scale 传入导致先验均值从 α/β 错为 αβ放大 100 倍。解决所有分布公式块顶部加粗声明参数化约定。Gamma 行写【参数约定】本文所有 Gamma(α,β) 均指 shapeα, scaleβ即 PDF∝x^{α−1}e^{−x/β}对应 scipy.stats.gamma(aα, scaleβ)NOT rateβ。并附验证代码assert abs(scipy.stats.gamma.mean(a2,scale3) - 6) 1e-10。4.4 现象KL 散度计算报错invalid value encountered in true_divide原因KL(P∥Q) Σ pᵢ log(pᵢ/qᵢ)当 qᵢ0 且 pᵢ0 时log(∞) 导致 NaN。但很多教程忽略此细节直接写求和公式。解决提供鲁棒实现用np.where屏蔽非法除法def kl_divergence_safe(p, q, eps1e-15): 安全 KL 散度p,q 为概率向量自动处理 q_i0 且 p_i0 的情况 p np.clip(p, eps, 1-eps) # 防 p_i0 q np.clip(q, eps, 1-eps) # 防 q_i0 return np.sum(p * np.log(p / q)) # 关键业务中若 Q 来自模型输出必须确保其 softmax 后无零值若 Q 是经验分布需平滑如加 Dirichlet prior4.5 现象用scipy.stats.norm.ppf(0.975)计算 95% 置信区间但线上服务返回 NaN原因ppf分位函数在输入超出 [0,1] 时返回nan。而业务代码中若上游传入的置信水平alpha因浮点误差变成 1.0000000000000002ppf(1.0000000000000002)即崩。解决所有分位函数调用前加裁剪def safe_ppf(dist, q, eps1e-12): q_clipped np.clip(q, eps, 1-eps) return dist.ppf(q_clipped) # 使用safe_ppf(norm(loc0, scale1), 0.975) → 永不崩这些不是“理论细节”而是你凌晨两点排查线上报警时真正决定能否在晨会前修复的命门。它们被塞进 .doc 的每一个公式角落不是为了炫技是为了让你少一次重启服务。5. 进阶把公式文档变成“可执行知识库”——用 Python 自动生成校验报告与参数看板5.1 用脚本驱动公式验证每次修改 .doc 前先跑通所有数值测试公式文档的生命力在于持续可信。我们用 Python 脚本将 .doc 中所有公式块解析为可执行测试用例。流程如下结构化解析用python-docx读取 .doc按标题层级识别“公式块”含公式 LaTeX、参数说明、数值示例代码块代码提取正则匹配python ... 间的代码提取输入变量如alpha2, beta8和断言如assert abs(...)1e-10批量执行用exec()在隔离 namespace 中运行每段代码捕获AssertionError、RuntimeWarning、nan输出生成报告HTML 报告列出✅ 全部通过 / ⚠️ 3 个警告如RuntimeWarning: divide by zero/ ❌ 1 个失败AssertionError: 0.2297 ! 0.2298这样当你修改“负二项分布的方差公式”时脚本会立刻告诉你修改前Varμμ²/r测试通过修改后若误写为Varμ²/r报告立即标红失败并显示具体哪一行assert不成立。这不再是“人肉校对”而是把公式文档变成了有 CI 的代码仓库。5.2 构建“参数影响看板”交互式观察公式输出对输入的敏感度公式不是静态等式而是参数空间的曲面。我们用plotly为关键公式生成交互看板。以t 分布置信区间宽度为例X 轴样本量 n5→200Y 轴置信水平 1−α0.9→0.999Z 轴区间半宽t_{α/2, n−1} × s/√n滑块可拖动调整样本标准差 s看板直观揭示当 n30 时α 从 0.05 降到 0.01区间宽度暴增 2.3 倍而 n100 后再降 α 影响微弱。这直接指导你小样本实验必须严控 α大样本可放宽。此类看板导出为 HTML嵌入团队 Wiki新人点开即懂“为什么我们 AB 实验要求 n≥50”。5.3 公式即 API封装为 PyPI 包让同事pip install probstats-formula直接调用最终形态不是一份文档而是一个生产级 Python 包probstats.formula.ttest带自动单/双侧选择、异方差检验、效应量Cohens d计算probstats.formula.confint统一接口输入数据methodt,bootstrap,wilcoxon自动选最优probstats.formula.distributionGamma、Beta 等分布的 PDF/CMF/PPF全部预编译 jit 加速包内每个函数 docstring 严格对应 .doc 中的公式块含数学公式 LaTeX、参数表、使用示例、避坑提示。当同事写代码时from probstats.formula import confint他调用的不是黑盒而是你亲手验证过的、带血泪经验的公式实体。我坚持这样做十年每次在模型评审会上当有人质疑“这个置信区间怎么算的”我不翻 PPT而是打开终端python -c from probstats.formula import confint; print(confint.__doc__)—— 文档即代码代码即文档。它不华丽但每次都能让争论戛然而止因为所有人看到的是同一份经过数值锤炼、边界清晰、可执行可验证的公式。希望帮到你。本文还有配套的精品资源点击获取
