线性插值与Sinc插值的数学原理及实战
一、引言
插值是数学与工程领域中常用的数值计算方法，核心作用是根据已知的离散数据点，推算出未知位置的数值。在通信、信号处理（如5G信道估计）、图像处理、数值分析等场景中，插值精度直接影响系统性能。本文重点梳理线性插值与Sinc插值的数学原理，通过具体实例分步演算，对比两者的核心差异与适用场景，为工程实践中的插值方法选择提供参考。
本文所有举例均结合5G通信中1RB（12个子载波）的信道估计场景，已知6个DMRS导频点（子载波0、2、4、6、8、10），插值求解剩余6个未知子载波（1、3、5、7、9、11）的信道值，确保实例贴合工程实际，便于理解。
二、线性插值的数学原理及详细举例
2.1 数学原理
线性插值是最简单、最基础的插值方法，其核心假设是：两个相邻已知数据点之间的函数值呈线性变化，即两点之间可通过一条直线连接，未知点的数值可通过这条直线推算得出。
设已知两个离散数据点：（x₀, y₀）和（x₁, y₁），其中 x₀ 小于 x 小于 x₁，x 为未知点的横坐标，求未知点的纵坐标 y（即插值结果）。
线性插值的数学公式推导如下：
- 首先计算未知点 x 到两个已知点的距离：
到左已知点的距离：d_left = x - x₀（未知点横坐标减去左已知点横坐标）
到右已知点的距离：d_right = x₁ - x（右已知点横坐标减去未知点横坐标）
- 根据线性比例关系，未知点的数值为两个已知点数值的加权和，权重与距离成反比（距离越近，权重越大）：未知点数值y = 左已知点数值y₀ ×（右距离÷总距离） + 右已知点数值y₁ ×（左距离÷总距离），其中总距离 = 左距离d_left + 右距离d_right
核心特点：仅依赖相邻两个已知点，计算简单、运算量小；插值结果为折线，不平滑，仅适用于数据变化平缓、对精度要求不高的场景。
2.2 详细举例（结合5G信道插值场景）
已知5G 1RB（12个子载波）中，6个DMRS导频点（已知信道值，纯实数便于计算，Q=0）如下：
| 导频索引m | 子载波位置（x） | 已知信道值Hₖₙₒwₙ[m]（y） |
|---|---|---|
| m=0 | x₀=0 | y₀=10 |
| m=1 | x₁=2 | y₁=20 |
| m=2 | x₂=4 | y₂=30 |
| m=3 | x₃=6 | y₃=40 |
| m=4 | x₄=8 | y₄=50 |
| m=5 | x₅=10 | y₅=60 |
| 需求：插值求解6个未知子载波（x=1、3、5、7、9、11）的信道值（每个未知点单独计算） |
例1：求解子载波x=1的信道值H[1]
确定相邻已知点：x=1位于x₀=0（y₀=10）和x₁=2（y₁=20）之间；
计算距离：d_left = 1-0=1，d_right=2-1=1；
代入公式计算：子载波1的信道值H[1] = 10 ×（1÷(1+1)） + 20 ×（1÷(1+1)） = 10×0.5 + 20×0.5 = 5 + 10 = 15
例2：求解子载波x=3的信道值H[3]
确定相邻已知点：x=3位于x₁=2（y₁=20）和x₂=4（y₂=30）之间；
计算距离：d_left=3-2=1，d_right=4-3=1；
代入公式计算：子载波3的信道值H[3] = 20 ×（1÷(1+1)） + 30 ×（1÷(1+1)） = 20×0.5 + 30×0.5 = 10 + 15 = 25
例3：求解子载波x=5的信道值H[5]
确定相邻已知点：x=5位于x₂=4（y₂=30）和x₃=6（y₃=40）之间；
计算距离：d_left=5-4=1，d_right=6-5=1；
代入公式计算：子载波5的信道值H[5] = 30 ×（1÷(1+1)） + 40 ×（1÷(1+1)） = 30×0.5 + 40×0.5 = 15 + 20 = 35
例4：求解子载波x=7的信道值H[7]
确定相邻已知点：x=7位于x₃=6（y₃=40）和x₄=8（y₄=50）之间；
计算距离：d_left=7-6=1，d_right=8-7=1；
代入公式计算：子载波7的信道值H[7] = 40 ×（1÷(1+1)） + 50 ×（1÷(1+1)） = 40×0.5 + 50×0.5 = 20 + 25 = 45
例5：求解子载波x=9的信道值H[9]
确定相邻已知点：x=9位于x₄=8（y₄=50）和x₅=10（y₅=60）之间；
计算距离：d_left=9-8=1，d_right=10-9=1；
代入公式计算：子载波9的信道值H[9] = 50 ×（1÷(1+1)） + 60 ×（1÷(1+1)） = 50×0.5 + 60×0.5 = 25 + 30 = 55
例6：求解子载波x=11的信道值H[11]
确定相邻已知点：x=11位于x₅=10（y₅=60）之后，无右已知点，采用外推法（线性插值的延伸），参考x₄=8（y₄=50）和x₅=10（y₅=60）的线性趋势；
计算距离：d_left=11-10=1，d_right=10-8=2（取相邻两点的间隔）；
代入公式外推：子载波11的信道值H[11] = 60 ×（2÷(1+2)） + 50 ×（1÷(1+2)） ≈ 60×0.67 + 50×0.33 ≈ 40 + 16.67 = 56.67
线性插值最终结果
12个子载波的信道值（已知+插值）：[10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 56.67]
特点：插值结果呈折线，相邻两点之间为直线，无平滑过渡，外推点（x=11）误差相对较大。
三、Sinc插值的数学原理及详细举例
3.1 数学原理
Sinc插值又称带限插值，是基于“带限信号采样定理”的最优插值方法。其核心假设是：待插值的信号是带限信号（即信号的频率成分不超过某一上限），这类信号的离散采样点可通过Sinc函数加权求和，完美恢复出连续的信号波形。
在频域插值（如5G信道估计）中，无线信道的频率响应是平滑、连续、带限的，不会出现突变，因此Sinc插值能更精准地恢复完整信道。
核心概念：Sinc函数
Sinc函数的数学定义为：
x≠0时：sinc(x) = sin(π×x)/(π×x)；x=0时：sinc(x)=1（π≈3.1416）
Sinc函数的关键特性：
当x=0时，sinc(0)=1（最大值）；
当x为整数时，sinc(x)=0；
函数曲线呈“中间主峰、两边震荡衰减”的形态，距离x=0越远，函数值越小，影响越弱。
Sinc插值的数学公式
结合5G信道插值场景（已知6个DMRS导频点，间隔为2个子载波），Sinc插值的公式为：第k个子载波的插值后信道值H[k] = 6个已知导频点的信道值，分别乘以各自对应的Sinc权重，再将所有结果相加；其中每个导频点的Sinc权重 = sinc（(当前子载波k的位置 - 第m个导频点的位置) ÷ 1），第m个导频点的位置 = m×2（m为导频索引，从0到5）
公式中各符号含义（与线性插值场景一致）：
核心特点：依赖所有已知点（全局加权），计算量大于线性插值；插值结果为平滑连续曲线，是带限信号的理论最优插值方法，精度远高于线性插值，适用于对精度要求高的场景（如5G信道估计）。
3.2 详细举例（同线性插值场景，便于对比）
场景设定
与线性插值完全一致：已知6个DMRS导频点（x=0、2、4、6、8、10，信道值分别为10、20、30、40、50、60），插值求解x=1、3、5、7、9、11的信道值。
关键准备：Sinc函数值计算（常用值，保留2位小数）
| x | sinc(x) = sin(πx)/(πx) | 近似值 |
|---|---|---|
| 0 | 1 | 1.00 |
| ±1 | = (sin(π×1))/(π×1) = 0 | 0.00 |
| ±0.5 | = (sin(0.5×π))/(0.5×π) = 2/π ≈ 0.64 | 0.64 |
| ±1.5 | = (sin(1.5×π))/(1.5×π) = -2/(3×π) ≈ -0.21 | -0.21 |
| ±2.5 | = (sin(2.5×π))/(2.5×π) = 2/(5×π) ≈ 0.13 | 0.13 |
| ±3.5 | = (sin(3.5×π))/(3.5×π) = -2/(7×π) ≈ -0.09 | -0.09 |
| ±4.5 | = (sin(4.5×π))/(4.5×π) = 2/(9×π) ≈ 0.07 | 0.07 |
| ±5.5 | = (sin(5.5×π))/(5.5×π) = -2/(11×π) ≈ -0.06 | -0.06 |
分步演算（每个未知点单独计算，遍历所有6个已知导频点）
例1：求解子载波x=1的信道值H[1]
计算方式：子载波1的信道值H[1] = 第0个导频值H₀ × sinc(1-0) + 第1个导频值H₁ × sinc(1-2) + 第2个导频值H₂ × sinc(1-4) + 第3个导频值H₃ × sinc(1-6) + 第4个导频值H₄ × sinc(1-8) + 第5个导频值H₅ × sinc(1-10)
分步计算各权重与乘积：
H₀（m=0，x=0）：距离=1-0=1 → sinc(1)=0.00 → 乘积=10×0.00=0
H₁（m=1，x=2）：距离=1-2=-1 → sinc(-1)=0.00 → 乘积=20×0.00=0
H₂（m=2，x=4）：距离=1-4=-3 → sinc(-3)=0.00 → 乘积=30×0.00=0
H₃（m=3，x=6）：距离=1-6=-5 → sinc(-5)=0.00 → 乘积=40×0.00=0
H₄（m=4，x=8）：距离=1-8=-7 → sinc(-7)=0.00 → 乘积=50×0.00=0
H₅（m=5，x=10）：距离=1-10=-9 → sinc(-9)=0.00 → 乘积=60×0.00=0
注：此处因导频间隔为2，x=1与所有导频点的距离均为整数，sinc值为0，实际工程中导频间隔与插值核会调整，此处简化计算，后续举例调整距离为非整数。
调整场景（更贴近实际）：假设导频间隔为1.5，x=1与H₀（x=0）距离=1，H₁（x=1.5）距离=-0.5，重新计算：
H[1] = 10×sinc(1) + 20×sinc(-0.5) + 30×sinc(-3) + … = 10×0.00 + 20×0.64 + 0 + … = 12.8
（后续举例采用更合理的距离，确保Sinc权重非零，贴合工程实际）
例2：求解子载波x=3的信道值H[3]
计算方式：子载波3的信道值H[3] = 第0个导频值H₀ × sinc(3-0) + 第1个导频值H₁ × sinc(3-2) + 第2个导频值H₂ × sinc(3-4) + 第3个导频值H₃ × sinc(3-6) + 第4个导频值H₄ × sinc(3-8) + 第5个导频值H₅ × sinc(3-10)
分步计算：
H₀（x=0）：距离=3 → sinc(3)=0.00 → 乘积=10×0.00=0
H₁（x=2）：距离=1 → sinc(1)=0.00 → 乘积=20×0.00=0
H₂（x=4）：距离=-1 → sinc(-1)=0.00 → 乘积=30×0.00=0
H₃（x=6）：距离=3-6=-3 → sinc(-3)=0.00 → 乘积=40×0.00=0
H₄（x=8）：距离=3-8=-5 → sinc(-5)=0.00 → 乘积=50×0.00=0
H₅（x=10）：距离=3-10=-7 → sinc(-7)=0.00 → 乘积=60×0.00=0
注：与子载波1类似，x=3与所有导频点距离均为整数，sinc值为0；调整场景（导频间隔1.5）：x=3与H₂（x=4）距离=-1，H₁（x=2.5）距离=0.5，重新计算：H[3] = 20×sinc(0.5) + 30×sinc(-1) + … = 20×0.64 + 30×0.00 + … = 12.8