信号的插值计算

信号的插值计算

---

前言

最近在学习信号的插值，学到了一些常用的方法，遂在此博客记录之

---

一、二次插值

二次函数插值是一种常见的插值方法，如下图：

假设在$x(n),x(n-1),x(n-2)$三个点已知，为了拟合估算$x(n-u)$的值，可以使用二次插值得到。假设二次曲线为$x(t)=a{t}^2+bt+c$，将$x(n),x(n-1),x(n-2)$三个点代入计算即可得到$a、b、c$的值。但这样计算略显繁琐，可以使用牛顿法直接插值得到这个二次曲线。
$$\begin{array}{rl}x(t)& =x(n)+\frac{x(n)-x(n-1)}{n-(n-1)}(t-n)+\left[\frac{x(n)-x(n-1)}{n-(n-1)}-\frac{x(n-1)-x(n-2)}{(n-1)-(n-2)}\right]\frac{(t-n)[t-(n-1)]}{n-(n-2)}\\ & =x(n)+\left[x(n)-x(n-1)\right](t-n)+\frac{1}{2}\left[x(n)+x(n-2)-2x(n-1)\right](t-n)(t-n+1)\end{array}$$

所以$x(n-u)$的插值结果为
$$\begin{array}{rl}x(n-u)& =\left(\frac{1}{2}{u}^2-\frac{3}{2}u+1\right)x(n)+(-u^2+2u)x(n-1)+\left(\frac{1}{2}{u}^2-\frac{1}{2}u\right)x(n-2)\\ & =c_{0}x(n)+c_{1}x(n-1)+c_{2}x(n-2)\end{array}$$
其中，
$$\{\begin{array}{ll}{c}_0& =\frac{1}{2}{u}^2-\frac{3}{2}u+1\\ c_1& =-u^2+2u\\ c_2& =\frac{1}{2}{u}^2-\frac{1}{2}u\end{array}$$
很明显这是一个滤波器的结构，所以对任意的三个已知点$x(n),x(n-1),x(n-2)$可以用滤波器插值得到$x(n-u)$的值，当$0<u<1$时，该滤波器实现了信号的小数倍延时。

二、分段式二次插值

以前在学习插值公式时，看到一个用4个点实现分段式二次插值的方法。假设$x(n+1),x(n),x(n-1),x(n-2)$四个点已知，利用这四个点得到类似于二次插值的系数。这个方法我查了一些资料^{1 , 2}，得到了一些点拨³，得到了计算方法。已知的4个点如下图所示。

这里的插值分为3段。
第1段利用$x(n),x(n-1),x(n-2)$三个点插值得到第1条二次曲线$x_1(t)$，同上一节类似，可以利用牛顿插值得到。
$$\begin{array}{r}{x}_1(t)=x(n)+\left[x(n)-x(n-1)\right](t-n)+\frac{1}{2}\left[x(n)+x(n-2)-2x(n-1)\right](t-n)(t-n+1)\end{array}$$
并且计算出$x_1(n-u)$，
$$\begin{array}{r}{x}_1(n-u)=\left(\frac{1}{2}{u}^2-\frac{3}{2}u+1\right)x(n)+(-u^2+2u)x(n-1)+\left(\frac{1}{2}{u}^2-\frac{1}{2}u\right)x(n-2)\end{array}$$
第2段利用$x(n+1),x(n),x(n-1)$三个点插值得到第2条二次曲线$x_2(t)$，
$$\begin{array}{r}{x}_2(t)=x(n+1)+\left[x(n+1)-x(n)\right](t-n-1)+\frac{1}{2}\left[x(n+1)+x(n-1)-2x(n)\right](t-n-1)(t-n)\end{array}$$
并且计算出$x_2(n-u)$，
$$\begin{array}{r}{x}_2(n-u)=\left(\frac{1}{2}{u}^2-\frac{1}{2}u\right)x(n+1)+(-u^2+1)x(n)+\left(\frac{1}{2}{u}^2+\frac{1}{2}u\right)x(n-1)\end{array}$$
第3段利用$x(n),x(n-1)$两个点线性插值得到第3条曲线$x_3(t)$，
$$x_3(t)=(t-n+1)x(n)+(-t+n)x(n-1)$$
可以计算出$x_3(n-u)$，
$$x_3(n-u)=(1-u)x(n)+ux(n-1)$$
这里，定义一个变量$\alpha ，0<\alpha <1$，将三段插值结果一起来看，即
$$x(t)=2\alpha x_1(t)+2\alpha x_2(t)+(1-4\alpha )x_3(t)$$
所以可以计算出$x(n-u)$的值
$$\begin{array}{rl}x(n-u)& =(\alpha u^2-\alpha u)x(n+1)+[-\alpha u^2-(1-\alpha )u+1]x(n)+[-\alpha u^2+(1+\alpha )u]x(n-1)+(\alpha u^2-\alpha u)x(n-2)\\ & =c_{-1}x(n+1)+c_{0}x(n)+c_{1}x(n-1)+c_{2}x(n-2)\end{array}$$
其中，
$$\{\begin{array}{ll}{c}_{-1}& =\alpha u^2-\alpha u\\ c_0& =-\alpha u^2-(1-\alpha )u+1\\ c_1& =-\alpha u^2+(1+\alpha )u\\ c_2& =\alpha u^2-\alpha u\end{array}$$
这种结构的插值利用了三段插值的加权相加的结果，具有更好的稳定性，能更好的避免边缘点插值时的插值结果和真实值相差太大的现象(被称为龙格现象)。
接下来可以看看在u取不同值时的滤波器响应如何，如下的结果是在$\alpha =0.5$的条件下仿真的。
幅度响应和群延迟如下所示。

我们希望滤波器的幅度响应为0dB，群延迟$\frac{d\varphi }{d\omega }=-u$，但是很显然，这种插值滤波器的响应和我们所期望的结果存在一定的误差。尤其是大约当$\omega >1.7$后，随着$\omega$逐渐增大，误差也越来越大。如此看来这种滤波器适合工作在采样率和信号带宽比值较大的条件下。

---

参考

---

1. 邱健伟. 小数倍多速率信号处理应用研究[D]. 北京邮电大学, 2017. ↩︎

2. https://www.doc88.com/p-1334885958480.html?r=1 ↩︎

3. https://www.douban.com/group/topic/6231956/ ↩︎