降采样

✍ dations ◷ 2025-04-02 13:43:24 #信号处理

在数位信号处理领域中,降采样,又作减采集,是一种多速率数字信号处理的技术或是降低信号采样率的过程,通常用于降低数据传输速率或者数据大小。跟插值互补,插值是用来增加取样频率。降采样的过程中会运用滤波器降低混叠造成的失真,因为降采样会有混叠的情形发生,系统中具有降采样功能的部分称为降频器。

降采样因子(常用表示符号为"M")一般是大于1的整数或有理数。这个因子表达采样周期变成原来的M倍,或者等价表示采样率变成原来的1/M倍。采样率的降低会造成频谱的压缩,因此需要利用滤波器确保在较低的采样频率下不发生混叠,确保奈奎式采样定理依旧成立。

降采样降低数据传输速率,降采样因子(速度降低的倍率M)通常为整数或是大于一的分数,这个数字为取样周期的倍数或是取样频率的因数。举例来说一个取样率为44,100 赫兹的16位元数字音乐讯号若被降采样到22,050 赫兹,此时降采样因子等于2,位元率同时降低了一半,假设降采样后每个取样长度依旧为16位元,则位元率从1,411,200 位元/秒(44100 (样本) × 16 (bit/样本) × 2 (声道))降低至705,600位元秒。

降采样整数M倍的过程可以被分解成两个部分:

若仅少了第一个步骤,讯号的高频成分会在降低采样速率的过程中被混入低频讯号中,这样的失真效应我们称为混叠,因此第一个步骤是为了降低混叠影响至可以接受的程度,降采样中的滤波器称为抗混叠滤波器,下面讨论如何设计抗混叠滤波器。

若用无限脉冲响应的方法设计抗混叠滤波器,则滤波的过程需要在降低采样速率前从输出端回传资讯到输入端,若用有限脉冲响应的方法设计抗混叠滤波器,滤波的过程较为简单因为我们只需要考虑间隔为M的采样点,滤波的过程表示成:

序列h是数位滤波器的脉冲响应、K是它的长度,序列x为要被降采样的原始讯号的取样点,一般来说在算出y后计算y最简单的方法是让x的索引增加M然后重新计算卷积,在大部分的情况下M=2,h可以被设计成半频带滤波器,序列h中将近一半的系数为0,因此卷积的计算过程会被大大简化。

每隔M个取样点取样脉冲响应的系数成为一个次序列,总共有M个降采样的脉冲响应次序列,内积的结果是每个次序列.跟对应的取样点x 内积的总和。此外,由于采样周期变成M倍,每个降采样后的次序列会独立在每个内积中,数列x 的采样点只会出现在一组内积之中而不会出现在其他组的内积里,因此可以用M组有限长度频率响应数位滤波器同时多工地算出M组内积,M组输出值最后平行加总。这个观点提供不同的硬件实作方法,很可能适合多处理器的架构。换句话说,输入的串流经过分工并送去M组滤波器并将输出加总,上述的架构称为多项位滤波器。

我们可以从图一的三组波型中了解为何降采样的过程需要抗混叠滤波器,注意三组波型图除了横轴的单位不同以外其余皆相同。每一组波形图的第一张图是取样讯号x(t),经过傅立叶转换得到的周期性频谱,X(f);第二章图为取样讯号x(t)以1/3倍采样速率降取样后经过傅立叶转换得到的周期性频谱。对于这三组波型,混叠不发生的条件同为: B < 1 M 1 2 T , {\displaystyle B<{\tfrac {1}{M}}\cdot {\tfrac {1}{2T}},}   T为采样的时间间格,1/T为采样速率,1/2T为奈奎氏取样频率,抗混叠滤波器借由设计截止频率小于1/M倍的奈奎氏取样频率确保混叠不会发生。

最上面那组图片的上面那张图的横轴代表离散时间傅立叶转换,也就是X(f)的傅立叶级数表示方式:

n = x ( n T ) x   e i 2 π f n T DTFT = 1 T k = X ( f k / T ) . {\displaystyle \underbrace {\sum _{n=-\infty }^{\infty }\overbrace {x(nT)} ^{x}\ e^{-i2\pi fnT}} _{\text{DTFT}}={\frac {1}{T}}\sum _{k=-\infty }^{\infty }X(f-k/T).}

 

 

 

 

(1)

T的单位为秒,f的单位为赫兹。将MT取代公式中的T我们可以得到以1/M倍速率采样后的数列的离散时间傅立叶转换,x:

这个周期性的讯号以较低速率采样过后,频率跟振幅都被降低了1/M倍,如同图一的第二张图所示,当X(f)与附近的波型重叠,混叠将会发生,抗混叠滤波器的目的就是在于确保降采样后周期降低但不会发生波型重叠,也就是混叠。

中间那组图频率f的单位改为标准化频率,这时周期为1而0.5的位置对应到奈奎氏取样频率,一般抗混叠滤波器的设计会以图形的频率单位为标准,因此截止频率也要改由以标准化频率作为单位,此时截止频率 B m a x = 1 M 1 2 T , {\displaystyle B_{max}={\tfrac {1}{M}}\cdot {\tfrac {1}{2T}},} 被标准化为 T B m a x = 1 M 1 2 = 0.5 M . {\displaystyle TB_{max}={\tfrac {1}{M}}\cdot {\tfrac {1}{2}}={\tfrac {0.5}{M}}.}   这个值的单位为(秒/取样)*(周期数/秒)=周期数/取样。

最下面那组图形为原本数列以及以低速率取样后数列的Z转换,受到复数变数 z = e i ω . {\displaystyle z=e^{i\omega }.}   的限制,x的数列的形式为傅立叶级数,与1比较,我们推论:

这个方程式描述了图一中第五张图,同理,第六章图为:

令 M/L代表降采样因子,M,L都是整数,M>L,分数倍降采样可以被分解成两个步骤

升采样需要低通滤波器过滤资料数率增加的讯号,降采样需要低通滤波器过滤输入讯号,因此这两个滤波过程可以被和而为一,借由用单一个低通滤波器取代,此单一低通滤波器的截止频率为两者低通滤波器的低者。当M > L, 抗混叠滤波器的截止频率   0.5 M {\displaystyle {\tfrac {0.5}{M}}} (周期数/取样)会是较低的截止频率。。

相关

  • 高等动物高等动物即“低等动物”的对称。在动物学中,与低等动物无明确的界限。一般指体制复杂,组织及器官分化显著,并具有脊椎的动物,即脊椎动物。
  • 最小二乘法最小二乘法(英语:least squares method),又称最小平方法,是一种数学优化方法。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便的求得未知的数据,并使得
  • 第509混合飞行大队第509混合飞行大队(英语:509th Composite Group)是美国陆军航空军在二战期间的一个空军编队,1944年12月17日组建,负责执行运送核武器的任务。该编队在1945年8月将两枚原子弹投放
  • 粪蝇科粪蝇科(学名:Scathophagidae),又名拟花蝇科,是家蝇总科下的一个小科。这个科下的黄粪蝇是北半球最为普遍及大量的苍蝇。粪蝇科下的物种在幼虫阶段较为多样化,包括有吃植物的、水中
  • 高溴酸高溴酸,或称过溴酸,化学式为HBrO4,是溴的含氧酸之一,其中溴的氧化态为+7。与高氯酸和高碘酸不同,高溴酸不稳定,不能通过溴酸盐热分解反应或卤素间的置换反应制备。首次制得利用的
  • 革命广场站 (莫斯科地铁)革命广场站(俄语:Пло́щадь Револю́ции,罗马化:Ploshchad Revolyutsii)是莫斯科地铁的一个车站,位于莫斯科市中心。车站名称来自于革命广场,是阿尔巴特-波克罗夫卡
  • 玩命鸳鸯《玩命鸳鸯》(英语:,发音: /zəˈzɪzɪks/ )是一出于2006年上映的惊悚美国电影,由约翰·彭尼编剧和导演,由利奥·格里洛、凯瑟琳·海格和汤姆·西斯摩尔主演。此电影因为其极低的
  • 美丽男孩《美丽男孩》(英语:)是一部2018年美国传记片,由菲利克斯·范·葛罗尼根(英语:Felix Van Groeningen)执导,本片为导演首部的英语长片。葛罗尼根和卢克·戴维斯(英语:Luke Davies)合力撰
  • 约翰·贝鲁西约翰·亚当·贝鲁西(英语:John Adam Belushi,1949年1月24日-1982年3月5日)是一名美国男演员、喜剧演员和歌手,同时也是NBC喜剧小品节目《周六夜现场》的七位元老之一。在芝加哥的
  • 栃木山守也栃木山守也,(2/2/1892-3/10/1959),本名中田守也,日本栃木县下都贺郡赤麻村(现在栃木县栃木市藤冈町赤麻)出身的前大相扑力士。第27代横纲。他身高172cm,重104公斤,是历代横纲中体重