梳理一下各大平台使用的sample rate convert算法

梳理一下各大平台使用的sample rate convert算法梳理一下各大平台使用的resample算法

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

前言

转采样属于数字信号重建的范畴,整数倍的升降采样可以通过插值抽取+带限滤波的方法【1】,并且通过级联来实现任意比值的采样率转换。这里需要对奈奎斯特采样定理,以及数字信号时频关系有比较清晰的理解。还有一种方法是插值,插值的方法比较多,在参考里有一阶保持FOH零阶保持ZOH三次样条函数splinesinc函数,引用的博文中有比较直观的分析和事例。那么在各大平台,最后都是用什么实现的resample重建呢?

Smarc对Interpolation和Decimation的执著

来自法国的Smarc的内核是经典的数字信号差值抽取方法,能够支持任意采样率,众所周知,48k到44.1k的互相转采样是个挑战,在SMARC method详细的讲述的算法的特点,通过因式分解和多相、多阶实现降低了FIR的阶数,提高了算法效率。简单的翻译一些文中的思想如下:
经典的香农插值抽取法虽然便于理解,但碍于为了实现非整数(有理数)比的转采样,滤波器级数可能长的无法用于实现,文中举了一个例子,可以这样理解 44100 48000 = 147 160 \frac{44100}{48000}=\frac{147}{160} 4800044100=160147那么两个频率的最小公倍数为: 44100 ∗ 160 = 48000 ∗ 147 = 7056000 = 7.056 M 44100*160=48000*147=7056000=7.056M 44100160=48000147=7056000=7.056M
为了满足48000->44100的变换,满足截止频率22.05k的设计要求,需要FIR的阶数38560,这个抽头和延时是无法接受,这个计算是利用smarc内置的remez方法。smarc则采样了多级滤波器的设计思想,规避了长延时的问题,每一级采用很小的系数比,同时结合多相滤波的实现,进一步优化算法。对于降采样抽取:

1

1

1/M

x(n)

h(n)

w(n)

M decimation

y(n)

经过低通滤波消除混叠频率后就可以抽取除想要的新采样。

CCRMA & speex

搜了一下外文网站,能够完整讲述插值重建的就是斯坦福大学的这篇Digital Audio Resampling Home Page,文中讲述了带限插值的原理。简单的翻译一下,根据奈奎斯特和傅里叶变换的知识,时域信号可以被重建:
x ( t ) ^ ≜ ∑ n = − ∞ ∞ x ( n T s ) h s ( t − n T s ) ≡ = x ( t ) \hat{x(t)} \triangleq\sum_{n=-\infty}^{\infty}x(nT_s)h_s(t-nT_s)\equiv=x(t) x(t)^n=x(nTs)hs(tnTs)=x(t)
这里的 h s ( t ) ≜ s i n c ( F s t ) ≜ s i n ( π F s t ) π F s t h_s(t)\triangleq sinc(F_st)\triangleq \frac{sin(\pi F_st)}{\pi F_st} hs(t)sinc(Fst)πFstsin(πFst)
s i n c sinc sinc函数的神奇就不赘述,但这个公式这么干是不行滴,因为正负无穷谁也受不了,所以还得加个窗,减小纹波和过冲。文中利用凯撒窗来评估对截取5个过零点的sinc函数卷积的频响比较,截取拼接了一下,效果变化很明显。
在这里插入图片描述如果要利用这个公式进行新采样率 F s ′ F_s^\prime Fs的抽取,可以改写: x ( n T s ′ ) = ∑ n = − ∞ ∞ x ( n T s ) h s ( n T s ′ − n T s ) = ∑ n = − ∞ ∞ x ( n T s ) h s ( n ( T s ′ − T s ) ) \begin{aligned} x(nT_s^\prime)&=\sum_{n=-\infty}^{\infty}x(nT_s)h_s(nT_s^\prime-nT_s)\\ &=\sum_{n=-\infty}^{\infty}x(nT_s)h_s(n(T_s^\prime-T_s)) \end{aligned} x(nTs)=n=x(nTs)hs(nTsnTs)=n=x(nTs)hs(n(TsTs))据此,首先准备sinc函数的样点,找了一张图:
在这里插入图片描述
如果能对准备好的 h s ( t ) h_s(t) hs(t)样点进行插值,那么就可以根据公式求出累加后的新采样点 x ( n T s ′ ) x(nT_s^\prime) x(nTs),这里假设 ρ = F s ′ F s \rho=\frac{F_s^\prime}{F_s} ρ=FsFs ρ \rho ρ小于1意味着降采样, ρ \rho ρ大于1则是升采样。这相当于设计一个低通滤波器,那么这个滤波器实际上只需要存储半边的系数。那么剩下来的问题就是这个插值如何计算了。原文写的很节略,刚开始看的云山雾绕,不如先回归公式,假设就使用一个矩形窗,截断无限长的冲激响应。 x ′ ( n ) ≜ ∑ n = − N N x ( n ) h s ( n T s ′ − n T s ) ≜ ∑ n = − N N x ( n ) h s ( n ( T s ′ − T s ) ) ≜ ∑ n = − N N x ( n ) h s ( n + ( T s ′ − T s ) T s ) \begin{aligned} x^\prime(n)&\triangleq\sum_{n=-N}^{N}x(n)h_s(nT_s^\prime-nT_s)\\ &\triangleq\sum_{n=-N}^{N}x(n)h_s(n(T_s^\prime-T_s))\\ &\triangleq\sum_{n=-N}^{N}x(n)h_s(n+\frac{(T_s^\prime-T_s)}{T_s}) \end{aligned} x(n)n=NNx(n)hs(nTsnTs)n=NNx(n)hs(n(TsTs))n=NNx(n)hs(n+Ts(TsTs))
上面这个卷积不是因果系统,这点跟插值抽取一样的特性。只要输入 2 N + 1 2N+1 2N+1个输入序列,同时能求得偏离 h s ( n ) h_s(n) hs(n)的插值,即可以算出来转采样之后的点。具体实现细节不追了。

sox & deadbeef

从它github主页上的描述

 for real-time resampling, libsoxr may have a higher latency
than non-FFT based resamplers.  For example, when using the `High Quality'
configuration to resample between 44100Hz and 48000Hz, the latency is
around 1000 output samples, i.e. roughly 20ms (though passband and FFT-
size configuration parameters may be used to reduce this figure).

推测这是个频域的转采样算法,据说是被deadbeef采用的,但实时性一般,

WebRTC

webRtc也有sinc带限转采样方法,对于整数的也有抽取和插值的方法。稍微多说的是早期的webrtc是不支持多媒体44.1族的采样率。在AOSP的webrtc代码里能找到sinc框架。

ffmpeg

早期用的是libavcodec,后来改用libswresampleGitHub有源码可供分析

Secret Rabbit Code

这个平台借用一个客户的评价“完胜大部分主流”,那这个平台里的用到的就是插值方法:

  1. ZOH插值
  2. 线性插值
  3. sinc插值

其中前两种是简单的,但声音质量一般,后一种是“完美”插值重建,而且还分为fast/medium/quality三种滤波器插值系数可选,满足用户不同的需求,是一个比较成熟的实现平台。

参考

1.采样率变换和多速率filter
2.【 MATLAB 】MATLAB 实现模拟信号采样后的重建(三)一阶保持(FOH)内插
3.【 MATLAB 】MATLAB 实现模拟信号采样后的重建(三)应用三次样条函数spline实现内插
4.【 MATLAB 】MATLAB 实现模拟信号采样后的重建(二)零阶保持(ZOH)
5.【 MATLAB 】MATLAB 实现模拟信号采样后的重建(一)
【信号与系统学习笔记】—— 采样与恢复(内插重建方法解析)
数字信号处理之内插
零阶保持器(ZOH)
Resampling Audio Algorithms
Secret Rabbit Code
Resampling from hydrogenaudio
Resampling from dspguru
Digital Audio Resampling Home Page
FFmpeg/libswresample/
Planet CCRMA at home
CCRMA:libresample
shibatch SSRC
chirlu /soxr
重采样Resample 的一些研究记录
What is DeaDBeeF?
speex
Speex编码器中回声消除算法的分析与评估
[WebRTC架构分析]采样率转换

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/181211.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • 线性代数的消元法_高斯消元法例题

    线性代数的消元法_高斯消元法例题1.消元的思想针对下面的方程,我们无法直接得到方程的解。$$\begin{alignedat}{2}&x\space\space&2&y\space=\space

  • 关于Vue使用es6模板字符串没反应的问题「建议收藏」

    关于Vue使用es6模板字符串没反应的问题「建议收藏」错误示范VScode发get请求的地址及参数使用单引号”包裹时,发现${this.keyWord}没有变颜色,跟字符串一个颜色,也就是没有将this.keyWord识别成变量,当成了一般字符串,发请求时带的参数问题请求不到结果search(){ this.$axios.get(‘https://api.github.com/search/users?q=${this.keyWord}’).then( res=>{ console.log(res); }, err=&gt

  • 手把手教你制作网页导航栏

    手把手教你制作网页导航栏手把手教你制作网页导航栏众所周知,导航栏是网页的重要组成部分,本篇文章将会带你由浅入深的制作导航栏,子菜单,样式风格变化等。导航栏的重要部分——ul标签在导航栏中的文字,都是以无序列表ul和li标签实现的。下面通过几个小例子来为大家介绍如何实现网页导航栏的制作。1.用ul标签做一个简单般的菜单样式,首先在一个div盒子里创建一个无序列表,如图所示运行结果如下图所示:2.之后我们可以让列表横过来,需要用到css的浮动样式也就是float标签。我们需要让它向左浮动如图所示:运行结果如

  • 编译器警告RegisterStartupScript已经过时

    编译器警告RegisterStartupScript已经过时protectedvoidAlertMsg(stringmsg){this.Page.RegisterStartupScript(“alert”,”<scriptlanguage=\”javascript\”>alert(‘”+msg+”‘);</script>”);}使用上面的代码建立的客户端脚本块编译器会警告Reg…

  • 安捷伦(Agilent)示波器使用简介

    安捷伦(Agilent)示波器使用简介安捷伦示波器使用简介

    2022年10月12日
  • UML时序图(Squence diagram)的设计介绍[通俗易懂]

    UML时序图(Squence diagram)的设计介绍[通俗易懂]本文以starUML工具为例。一、概念UML时序图,也叫顺序图,是一种详细表示对象之间以及对象与参与者实例之间交互的图,它由一组协作的对象(或参与者实例)以及它们之间可发送的消息组成,它强调消息之间的顺序。它详细而直观地表现了一组相互协作的对象在执行一个(或少量几个)用例时的行为依赖关系,以及操作和消息的时序关系。类图是反映类、对象、接口等事物之间的依赖、继承、实现关系等,主要是静…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号