新闻资讯

暖气片选型表上两行之间怎么办：中间值线性插值只差0.3%，往低温外推到35度差5.5%

2026年10月6日0

一、样本上只有两行，你要的那一档往往没有

买散热器时拿到的样本，散热量经常只给两行：ΔT=50 K 时 1000 W，ΔT=64.5 K 时 1392.4 W。前一个是欧洲样本常用的标定温差，来自供水 75 ℃、回水 65 ℃、室温 20 ℃，算术平均减室温得 70−20=50 K；后一个是国标 GB/T 13754 的基准，来自供水 95 ℃、回水 70 ℃、室温 18 ℃，82.5−18=64.5 K。

问题来了：你家壁挂炉跑 80/60 ℃、室温 20 ℃，ΔT=(80+60)/2−20=50 K，正好对上第一行，还好。可如果跑冷凝工况 60/50 ℃，ΔT=55−20=35 K，表里根本没有。剩下的办法只有两个：在两行之间插值，或者往表外推。这两件事的风险完全不是一个量级。

二、散热量不是直线，ΔT 的 1.3 次方是凸的

散热量的标定关系是 Q=Q₀·(ΔT/ΔT₀)^n，n 取 1.3。指数大于 1，函数在算术坐标里是向上弯的。凸函数有一条几何性质：连接两点的弦永远落在曲线的上方。所以只要你在两行之间做线性插值，得到的是弦上的值，它必然比真实曲线高。

反过来，弦往区间外延长，立刻掉到曲线下方。这一条决定了全文的结论方向：内插系统性高估，外推系统性低估，两者符号相反。用上面的两行算：ΔT=57.25 K（50 与 64.5 的正中点）真实值是 1192.5 W，线性插值给 1196.2 W，高估 0.3132%；而 ΔT=35 K 的真实值是 629.0 W，线性外推只给 594.1 W，低估 5.5498%。差了十几倍，方向还相反。

三、误差最大的那一点，几乎就是中点

弦与曲线的差 E(x)=L(x)−Q(x)，最大处满足 L′=Q′，也就是弦斜率等于曲线斜率。对 Q=c·ΔT^n 求导得 n·c·ΔT^(n−1)=s，反解出

ΔT* = (s/(n·c))^(1/(n−1))

代入 [50, 64.5] 两行：弦斜率 s=27.0626 W/K，c=1000/50^1.3=6.1850，解出 ΔT*=57.1426 K。而区间中点是 57.25 K——两者只差 0.107 K。也就是说，你把最大误差点当中点用，误差的损失可以忽略。真正需要记住的是这个点的相对误差：0.3141%，而中点处是 0.3132%，几乎一样。所以后面所有结论都可以直接用中点代替，不用担心取错了位置。

四、误差只取决于两行的比值，跟绝对温度无关

把区间归一化成 [a, b]，令 r=b/a，可以证明最大相对误差只由 r 和 n 决定：

相对误差 = [1+(n−1)·u^n − n·u^(n−1)] / u^n，其中 u=(S/n)^(1/(n−1))，S=(r^n−1)/(r−1)

式子里 a、b 都不出现，只剩比值 r。数值验证得很干净：r=1.2 时，区间取 [50,60]、[100,120]、[500,600]、[5000,6000]，最大相对误差全部是 0.161510%，小数点后六位都一样。r=1.29 时，[50,64.5]、[100,129]、[500,645]、[5000,6450] 全部是 0.314058%。

这条尺度不变性有个很实用的推论：判断一张表能不能线性插值，你只需要看两行相距多少倍，不用管它是 50 度还是 500 度。样本表把 ΔT 从 50 排到 64.5，比值 1.29，误差上限就是 0.314%。

五、区间越宽，误差涨得比宽度本身还快

按上面那个闭式扫一遍不同比值，中点相对误差 ((r^n+1)/2)/((r+1)/2)^n−1 如下：

两行比值 r	举例区间	中点相对误差	最大相对误差
1.1	50→55 K	0.0442%	0.0442%
1.2	50→60 K	0.1613%	0.1615%
1.29	50→64.5 K	0.3132%	0.3141%
1.5	50→75 K	0.7831%	0.7884%
2.0	50→100 K	2.1915%	2.2333%
3.0	50→150 K	5.0073%	5.2339%

区间宽度从 1.2 倍拉到 3.0 倍（宽度涨 5 倍），误差从 0.1613% 涨到 5.0073%（涨 31 倍）。更精确的规律在收敛阶里：把区间逐次减半，[50,79]、[50,64.5]、[50,57.25]、[50,53.625] 的最大误差分别是 0.99893%、0.31406%、0.08922%、0.02387%，相邻比值 3.181、3.520、3.737，一路逼近 4。区间减半误差降到四分之一，这就是二阶收敛 O(h²)。

六、真正的坑是外推，不是内插

把区间减半能把误差压到四分之一，但那是内插。外推完全不是一回事。

线性插值的误差可以写成严格的余项形式：E(x) = −Q″(ξ)/2 · (x−a)(x−b)。决定性的是 (x−a)(x−b) 这一项。x 落在 a、b 之间时，(x−a) 为正、(x−b) 为负，乘积是负的，而且绝对值被区间宽度夹住；x 一旦跑到区间外面，两项同号，乘积变正，并且随距离二次增长。

设从端点往外挪 d（K），内插侧的余项因子是 d·(b−a−d)，外推侧是 d·(b−a+d)，两者之比：

外推/内插 = (b−a+d)/(b−a−d)

代入 b−a=14.5 K：d=3 时放大 1.5217 倍，d=5 时 2.0526 倍，d=7.25（正好中点）时 3.0000 倍，d=10 时 5.4444 倍，d=12 时 10.6000 倍。距离越接近半个区间，这个比值发散得越快。

七、把 10 倍误差拆成三个因子

用实际数字检验一下。从 [50, 64.5] 出发，向内 10 K 到 ΔT=60 K，向外 10 K 到 ΔT=40 K：

ΔT=60 K：线性给 1270.626 W，真值 1267.464 W，绝对误差 +3.1625 W，相对 +0.2495%
ΔT=40 K：线性给 729.374 W，真值 748.199 W，绝对误差 −18.8252 W，相对 −2.5161%

相对误差之比是 10.0839 倍。这个数可以由三个因子精确乘出来：

10.0839 = 5.4444 × 1.0934 × 1.6940

第一个是余项因子比 (b−a+d)/(b−a−d)=5.4444；第二个是 Q″(ξ) 的差异，低温侧曲率更大，贡献 1.0934；第三个是基数比 Q(60)/Q(40)=1267.464/748.199=1.6940，同样的绝对误差除以更小的基数，相对误差自然更大。三个因子相乘 10.0839，与实际比值 10.0839 完全吻合，没有残差。

八、两条独立路径对账

上面用的是余项公式，还可以用另一条完全不依赖余项定理的路：直接拿弦减曲线。在 57.25 K 处，Q″(ξ)≈0.141893 W/K²，(x−a)(x−b)=(7.25)·(−7.25)=−52.5625，余项公式给 E=−0.141893/2×(−52.5625)=3.7291 W。直接相减：1196.2−1192.5=3.7351 W。两条路径偏差 0.160%，差别只来自 ξ 取了中点而不是真实的中值点。

这种对账值得每次都做一遍：一条路走余项定理，一条路走算术，两个数字对得上，说明前面所有的符号、系数、单位都没写反。

九、换一张坐标纸，误差直接归零

ΔT^1.3 在算术坐标里是弯的，但在双对数坐标里是一条直线——因为取对数后 ln Q = ln c + n·ln ΔT，斜率就是 n。既然是直线，那就在对数坐标里插值，误差理论上为零。

实测：

ΔT (K)	真值 (W)	算术线性插值	双对数插值
32.25	565.5	519.6（−8.1089%）	565.5（0%）
35	629.0	594.1（−5.5498%）	629.0（0%）
40	748.2	729.4（−2.5161%）	748.2（0%）
45	872.0	864.7（−0.8384%）	872.0（0%）
70	1548.7	1541.3（−0.4807%）	1548.7（0%）

双对数插值的计算残差在 1e-13% 量级，就是浮点误差。这一步不需要任何新数据，只需要把两个已知的 W 值和两个 ΔT 各取一次自然对数，做一次直线插值，再取指数回来。

十、怎么判断这张表能不能用对数插值：反解指数

关键在于：真实散热器的 n 不是一个放之四海皆准的常数。把两行代进去可以反解它：

n = ln(Q_b/Q_a) / ln(ΔT_b/ΔT_a)

用 [50, 64.5] 两行：ln(1392.4/1000)/ln(64.5/50)=1.300000，正好还原出 1.3。但如果样本是按 n=1.25 折算的（1374.8 W），同一段区间反解出来就是 1.250000。

这个差别有多大？64.5 K 折算到 50 K：n=1.25 给 72.74%，n=1.30 给 71.82%，相差 1.281%。所以操作方法是：拿样本上相邻的两段分别反解 n，两段解出来不一样，说明这张表不是纯幂律，对数插值也不能直接用。这时就该上三个点。

十一、三个点就够：二次插值

再加一个点做二次（抛物线）插值，误差阶数从 O(h²) 升到 O(h³)。用 [35, 50, 64.5] 三点：

ΔT (K)	真值 (W)	两点线性	三点二次	改善
40	748.2	729.4（−2.5161%）	748.7（+0.0670%）	37.6 倍
45	872.0	864.7（−0.8384%）	872.4（+0.0436%）	19.2 倍
55	1131.9	1135.3（+0.3010%）	1131.6（−0.0300%）	10.0 倍
60	1267.5	1270.6（+0.2495%）	1267.1（−0.0306%）	8.2 倍

收敛阶也变了：区间逐次减半，三点二次的最大误差比值是 6.33、7.03、7.47，一路逼近 8，正是 2³。同区间下三点相对两点的改善倍数随区间收窄急剧上升：[50,79] 是 38.7 倍，[50,64.5] 是 77.0 倍，[50,57.25] 是 153.7 倍。点数只多一个，精度涨两个数量级。

十二、本篇不算什么

本篇只处理"已知两端、求中间"这一类问题，与下面几篇划清边界：

室温每升1度暖气片少散多少热：1.3次方定律和64.5度标定温差——那篇讲指数 1.3 的来源，本篇把它当已知条件。
一组暖气片到底散多少热：K值、面积和ΔT的1.3次方怎么算——那篇从 K 值和面积正推，本篇只做表上取值。
暖气片回水温度怎么算：用二分法解散热与流量平衡点、暖气片供水温度设多少度：用牛顿迭代法反解非线性散热量方程——那是"已知目标求输入"的求根问题，本篇是"已知两端求中间"的查表问题，两者方向相反。
壁挂炉部分负荷热效率怎么估算：用线性插值补齐测试点之间的曲线——那篇处理的是热效率曲线的补齐，本篇处理散热器散热量随温差的变化，底层的曲线形状不同。

十三、三步实操

第一步，看两行相距多少倍。比值小于 1.3，直接线性插值，误差不超过 0.32%，比样本本身的制造公差还小，不用折腾。

第二步，只要目标点落在两行之外，尤其是往低温方向（冷凝炉、热泵工况），必须换方法。ΔT=35 K 时线性外推差 5.5%，ΔT=32.25 K 时差 8.1%，这个量级已经会影响配几片。

第三步，优先用双对数插值：两个 W 取自然对数，两个 ΔT 取自然对数，做一次直线插值再取指数。纯幂律下误差为零，而且不需要第三个点。如果样本上相邻两段反解出的 n 不一致，说明表本身不是纯幂律，就再要一个点做二次插值，误差能再压一到两个数量级。

顺带一提，散热量标定温差 64.5 K 是国标基准，50 K 是欧洲样本基准，两张表不能直接比数字，先折算到同一个 ΔT 再比。折算时 n 取 1.25 还是 1.30，结果就差 1.28%，这也是同一批暖气片在两份样本上标称值对不上的常见原因。