2. 土壤湿度传感器多点分段线性标定
为什么两点标定不够用
电容式土壤湿度传感器在干土端与饱和端之间往往呈现明显的S形响应：低含水量区灵敏度偏低，中间区近似线性，高含水量区又趋于饱和。用两点（烘干土、饱和土）拟合一条直线，虽然端点误差为零，但中段可能出现数个百分点体积含水量的系统性偏差。多点分段线性标定的思路是：在五个已知含水量梯度上采集响应，先看整体线性拟合的残差分布，再在残差出现结构性拐点处插入分段节点，用分段折线逼近真实响应曲线。
五个梯度点的制备与采集
本例采用体积含水量 θv（m³/m³）作为基准量，传感器原始输出取归一化计数 N（0–1000，随含水量单调上升）。五个目标梯度设为 0.05、0.12、0.20、0.28、0.38。制备方法：取同一批次风干土，按目标质量含水量加水并闷置 24 h 使水分均匀，装入固定容积环刀（100 cm³）压实至同一容重，每个梯度制备 3 个平行样。传感器插入深度保持一致，静置 2 min 待读数稳定后记录，取 3 个平行样的中位数作为该梯度的响应值。
下表为教学示例数据（示例值，非实测）：
| 梯度点 | θv 基准 | N 中位数 | 平行样极差 |
|---|---|---|---|
| P1 | 0.050 | 182 | 9 |
| P2 | 0.120 | 341 | 11 |
| P3 | 0.200 | 512 | 8 |
| P4 | 0.280 | 648 | 13 |
| P5 | 0.380 | 742 | 15 |
平行样极差均在 15 计数以内，说明制样均匀性可接受；若某梯度极差超过 30，应先排查装填密实度或传感器贴合问题，再进入拟合。
整体线性拟合与残差诊断
对五点做最小二乘直线拟合 θ = a·N + b。用示例数据计算，得到 a ≈ 5.72×10⁻⁴，b ≈ −0.0546。将各点 N 代入反算 θ̂，残差 r = θ − θ̂ 如下：
| 梯度点 | θ 基准 | θ̂ 线性预测 | 残差 r |
|---|---|---|---|
| P1 | 0.050 | 0.0495 | +0.0005 |
| P2 | 0.120 | 0.1405 | −0.0205 |
| P3 | 0.200 | 0.2383 | −0.0383 |
| P4 | 0.280 | 0.3161 | −0.0361 |
| P5 | 0.380 | 0.3699 | +0.0101 |
残差呈“两端小、中间负偏大”的弧形分布，最大绝对残差 0.0383 m³/m³（约 3.8% 体积含水量），远超田间常用的 ±0.02 精度要求。这种中间段系统性低估正是单直线无法描述S形响应的典型表现，说明必须引入分段。
分段方案与节点确定
残差在 P2 之后迅速跌出精度带、在 P4 之后回升，说明响应曲线在 N≈341 与 N≈648 附近存在斜率变化。据此把量程分为三段：低湿段 [P1,P2]、中湿段 [P2,P4]、高湿段 [P4,P5]，节点取 P2（N=341）与 P4（N=648）。分段线性插值在每个节点上残差为零，段内误差由折线与真实曲线的偏离决定。用示例数据计算各段中点（N=261、N=580、N=695）的插值残差，得到约 +0.006、−0.009、+0.004，全部落在 ±0.02 以内。
对比结论：整体线性最大绝对残差 0.0383，分段线性最大绝对残差约 0.009，精度提升约 4 倍，满足田间 ±0.02 要求。分段节点即 N=341 与 N=648，对应 θ 基准 0.120 与 0.280。
段1 [182,341]：θ = 0.120 + (N−341)×(0.120−0.050)/(341−182) = 0.120 + (N−341)×4.403×10⁻⁴
段2 [341,648]：θ = 0.280 + (N−648)×(0.280−0.120)/(648−341) = 0.280 + (N−648)×5.212×10⁻⁴
段3 [648,742]：θ = 0.380 + (N−742)×(0.380−0.280)/(742−648) = 0.380 + (N−742)×1.064×10⁻³
验证与落地要点
用未参与标定的 3 个检验点（θ 约 0.09、0.24、0.33）复测，分别代入所在段公式反算 θ，与烘干法基准比较，若三点绝对误差均 ≤0.02 则方案通过。落地时注意：分段节点应写入采集端固件或上位机标定表，节点处采用左闭右开区间避免重复计算；土壤类型或容重变化后需重新制样标定，不可直接沿用系数。