Skip to content
OdorAssist
OdorAssist

10. 附录:核心算法参照说明

为彰显实验室测试数据的合法合规性与权威性,本系统摒弃了所有人工修约或估算的可能。本系统后台统计分析引擎中内嵌的全部计算模型,均完全采纳并对齐《环境空气和废气 臭气的测定 三点比较式臭袋法(HJ 1262-2022)》标准原文提供的官方公式。

以下列举了系统核心自动核算环节所依据的主要公式参数,供实验室质量审计与报告查验参考。

10.1 环境空气和无组织排放监控点空气

Section titled “10.1 环境空气和无组织排放监控点空气”

此类样品的实验依据 MM 值(正解率)判断实验进程,并在终止后利用内插法计算最终浓度。

系统实时收集 6 名嗅辨员在同一稀释倍数下重复 3 次实验的共 18 个结果,代入公式计算:

M=1.00×a+0.33×b+0×c18M = \frac {1.00 \times a + 0.33 \times b + 0 \times c}{18}
  • aa:答案正确且肯定的总人次数(记为 1.00 权重)
  • bb:答案正确但猜测/不明确的总人次数(记为 0.33 权重)
  • cc:答案错误的总人次数(记为 0 权重)
  • 系统终止拦截条件:当计算得 M0.58M \le 0.58 时,该样品的嗅辨实验宣告结束。

利用临界值前后相邻两级的正解率,通过内插推导:

  • 幂参数 α\alphaα=M10.58M1M2\alpha = \frac {M_1 - 0.58}{M_1 - M_2} (其中 M2M_20.58\le 0.58 的正解率,M1M_1 为上一级 >0.58> 0.58 的正解率)
  • 幂参数 β\betaβ=lgt2t1\beta = \lg \frac {t_2}{t_1} (其中 t1t_1t2t_2 分别对应 M1M_1M2M_2 时的稀释倍数)
  • 最终臭气浓度数值生成: Y=t1×10αβY = t_1 \times 10^{\alpha \cdot \beta}

此类样品依据两次或三次的平行检验以及显著性分析模型计算。

系统依据单人判定的“最高正确稀释倍数”和“初次失误稀释倍数”的对数进行平均计算:

Xi=lga1+lga22X_i = \frac {\lg a_1 + \lg a_2}{2}
  • a1a_1:个人正解最大稀释倍数
  • a2a_2:个人误解稀释倍数

2. 小组平均嗅觉阈值 (Xˉ\bar{X})

Section titled “2. 小组平均嗅觉阈值 (Xˉ\bar{X}Xˉ)”

在单次平行实验中,系统汇总 nn 名(不少于 4 名)嗅辨员的个人嗅觉阈值并求其算术平均值:

Xˉ=i=1nXin\bar{X} = \frac {\sum_{i=1}^{n} X_i}{n}
  • nn:参与该次平行实验的有效嗅辨员人数

为了客观评估嗅辨小组在两次平行试验中状态的离散程度,系统调用 tt 检验公式:

t=Xˉ1Xˉ2SX12+SX222γSX1SX2n1t = \frac {\left| \bar{X}_1 - \bar{X}_2 \right|}{\sqrt{\frac{S^2_{X_1} + S^2_{X_2} - 2\gamma S_{X_1} S_{X_2}}{n-1}}}
  • Xˉ1\bar{X}_1Xˉ2\bar{X}_2:两次测试的平均阈值
  • SX12S^2_{X_1}SX22S^2_{X_2}:两次测试的小组方差
  • γ\gamma:两次测试结果的相关系数
  • 自动判定动作:若系统算得的统计量 tt 大于对应的 tt 分布临界值(如自由度 3、α=0.05\alpha=0.05 时为 3.182),表明存在显著性差异,系统强制开启第三次平行测定。

选取通过 t 检验的个人阈值求平均(Xˉ\bar{X}),生成该固定源样品浓度:

Y=10XˉY = 10^{\bar{X}}

10.3 嗅觉灵敏度准入核算(正丁醇测试)

Section titled “10.3 嗅觉灵敏度准入核算(正丁醇测试)”

系统在日常建档与每日实验前执行的“正丁醇在线标定任务”,其核心准入拦截逻辑建立在重复性标准偏差的分析上:

S=i=1n(xixˉ)2n1S = \sqrt{\frac{\sum_{i=1}^{n}(x_i - \bar{x})^2}{n-1}}
  • nn:取时间最近的 10 次有效测试结果
  • xˉ\bar{x}:10 次测试阈值的算术平均值
  • 准入通关条件:系统强制要求平均阈值浓度处于 20×10380×103 μmol/mol20 \times 10^{-3} \sim 80 \times 10^{-3} \ \mu \mathrm{mol/mol} 区间,且: 10S2.310^S \le 2.3 否则触发合规熔断并对账号状态采取长达 2 小时的业务禁入锁定。