基于sinc的音频重采样（二）：实现

上篇（基于sinc的音频重采样（一）：原理）讲了基于sinc方法的重采样原理，并给出了数学表达式，如下：

                     (1)

本文讲如何基于这个数学表达式来做软件实现。软件实现的细节很多，这里主要讲核心部分。函数srcUD()和filterUD()就是实现的主要函数（这两个函数是在源码基础上作了一定的改动，核心思想没变）。srcUD()是实现一帧中点的重采样，一个点一个点的做。filterUD()被srcUD()调用。数学表达式就体现在函数filterUD()里。粗看肯定会懵，怎么也跟上面的数学表达式联系不起来。下面就讲讲实现细节，让代码和数学表达式联系起来。

 

先看函数srcUD( )，下图是其实现，是将一帧中原采样点的值转换成新的采样率下的值，图中对输入参数的意思已做了解释。

 

主要看3个细节。标1的红框处是算新的采样率下两个点之间的采样间隔。设定Tx为原采样率下的采样间隔，Ty为新的采样率下的采样间隔，factor = Tx / Ty，所以 Ty = Tx / factor。软件实现时对采样间隔做归一化处理，设定Tx = 1，即原采样点x(n)的采样间隔归一化到1，则新采样点y(m)的采样间隔为Ty = 1/factor。把