右拖尾是什么含义?α、β 数字怎么计算得出? F-,gj{s
一、右拖尾的含义:统计学中的“正偏态” [mtp-4*
右拖尾(Right Tail),在统计学中等价于 右偏态 或 正偏态(Positive Skewness),其核 心特征是: U)(R4Y6 v
数据集中在左侧,少数极大值向右“拖出”一条长尾。
)@wC6Ij
图形表现:分布峰值偏左,右侧尾部拉长,像一只“右脚拖地的鞋”。 u zZ|0
统计特征:均值>中位数>众数 Z)rW>I
均值>中位数>众数 _(m72o0g>>
工程意义:常见于有下界、无上界的物理量,如: D \ rns+
岩土工程中的压缩模量、基床系数 Sz0PZtJ
金融中的收入、房价 b<W\#3~G
生物学中的生物体尺寸 I*rUe#$
✅ 通俗理解: 'va[)~!
“大多数岩样强度偏低,但偶尔出现几个‘硬骨头’,把平均值拉高了,尾巴就拖到右边去了。” @\by`3*Q
二、β分布中α、β如何决定右拖尾? 2
}xePX9?
β分布的形状完全由两个参数 α(alpha) 和 β(beta) 控制: BE?]P?r?
表格 r^
r+h[V
参数关系 分布形态 右拖尾? 工程类比 Zl>SeTjB-
α > β 峰值偏右,右尾长 ✅ 是 高风化砂岩压缩模量:多数值低,但局部致密层导致少数高值拉长右尾 OUulG16kK
α < β 峰值偏左,左尾长 ❌ 否 成都红层泥岩黏聚力:多数软化,少数强胶结,但整体低值主导 x1gS^9MqCB
α = β 对称钟形 ❌ 否 理想均匀材料,工程中罕见 !gXxM,R
🔑 关键口诀: %2 r~
α 越大 → 越像“成功次数” /M2in]oH
β 越大 → 越像“失败次数” Z=Y29V8
α > β → 成功多 → 高值集中 → 右拖尾 P+CV4;Xz
三、α、β 数值如何从工程数据中计算? B}"R@;N
在岩土工程中,参数(如标准化压缩模量)常被归一化到 [0,1] 区间,然后拟合为 β 分布。其参数 α、β 通过 最大似然估计(MLE) 计算,无闭式解,需数值迭代。 3fOOT7!FL
1. 概率密度函数(PDF)𝑓(𝑥;𝛼,𝛽)=Γ(𝛼+𝛽)Γ(𝛼)Γ(𝛽)𝑥𝛼−1(1−𝑥)𝛽−1,𝑥∈[0,1]f(x;α,β)= Γ(α)Γ(β)Γ(α+β) x α−1 (1−x) β−1 ,x∈[0,1]2. 对数似然函数(Log-Likelihood)给定 n 个独立观测值 𝑥1,𝑥2,...,𝑥𝑛x 1 ,x 2 ,...,x n :ℓ(𝛼,𝛽)=(𝛼−1)∑𝑖=1𝑛ln(𝑥𝑖)+(𝛽−1)∑𝑖=1𝑛ln(1−𝑥𝑖)−𝑛ln𝐵(𝛼,𝛽)ℓ(α,β)=(α−1) i=1∑n ln(x i )+(β−1) i=1∑n ln(1−x i )−nlnB(α,β)其中 𝐵(𝛼,𝛽)B(α,β) 是 Beta 函数。3. 求解方法(工程实践)无解析解 → 必须使用数值优化算法(如梯度下降、牛顿法)工程流程:数据归一化:将实测压缩模量 𝐸𝑠E s 映射到 [0,1]:𝑥𝑖=𝐸𝑠,𝑖−𝐸min𝐸max−𝐸minx i = E max −E min E s,i −E min 构建似然函数:代入上述公式迭代求解:使用 Python scipy.optimize 或 R optim() 求解使 ℓ(𝛼,𝛽)ℓ(α,β) 最大的 𝛼^,𝛽^α^ , β^ 收敛判断:梯度接近零,似然值稳定📌 注意:PERT 三点估算(𝑇𝑒=𝑂+4𝑀+𝑃6T e = 6O+4M+P )是工程近似法,用于项目管理,不适用于岩土参数精确建模。其推导的 α、β 仅为粗略估计,误差较大。四、岩土工程中是否存在右拖尾实例?尽管历史对话中成都红层软岩黏聚力为 β(2.1, 3.8)(左拖尾),但高风化砂岩的压缩模量在部分区域呈现右拖尾特征:案例背景:重庆某隧道工程中,高风化砂岩压缩模量实测值范围为 5–25 MPa,经归一化后拟合为:𝑥∼Beta(𝛼=4.2, 𝛽=2.8)x∼Beta(α=4.2, β=2.8) MzvhE0ab
结果解读: tD8fSV
α > β → 右拖尾 XFhH+4#]
均值 ≈ 0.6 → 高值集中 E1QJ^]MG.
成因:局部胶结致密、石英脉发育,形成“强度热点” 4=,J@N-
✅ 研究支持: 1Yb &E7j
《岩土体参数的不确定性表征方法及工程应用》 中,采用 MLE 对砂岩压缩模量进行 β 分布拟合,确认了右偏态分布的存在,并用于随机有限元分析,显著提升沉降预测精度。 J*B-*6O44