降采样

✍ dations ◷ 2025-10-14 07:35:14 #信号处理

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

降采样因子(常用表示符号为"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}}} (周期数/取样)会是较低的截止频率。。

相关

  • 杰夫·贝佐斯杰弗里·普雷斯顿·“杰夫”·贝索斯(英语:Jeffrey Preston "Jeff" Bezos,/ˈbeɪzoʊs/;1964年1月12日-),生于美国新墨西哥州,美国互联网巨头亚马逊公司创始人及现任董事长兼CEO,华
  • 南海路南海路是中国青岛市市南区的一条东南-西北方向临海道路,也是全市最早建成的一批道路之一。道路东南接八大关风景区的正阳关路,与汇泉路、武昌路、文登路交汇,西北接莱阳路。路
  • 蛸亚纲见内文蛸亚纲(学名:Coleoidea)是头足纲下的一类软体动物。它们不像鹦鹉螺亚纲般有硬壳,最多只有用来控制浮沉的内骨。一些物种甚至完全没有骨头,一些则以软骨来支撑。具有墨囊和
  • 玛雅历玛雅历是一套以不同历法与年鉴所组成的系统,为前哥伦布时期中部美洲的玛雅文明所使用。现仍使用于墨西哥的瓦哈卡州、恰帕斯州、韦拉克鲁斯州以及危地马拉高地的玛雅社群中。
  • 欧巴宾海蝎欧巴宾海蝎(属名:Opabinia,意为“(来自)欧巴宾山口的”。欧巴宾山口(Opabin Pass)是本种的化石发现地),是一种存活在寒武纪的生物。此种生物只有一个种。另一种学名为Myoscolex的同一
  • 苏格兰教会(1689年彻底分开)苏格兰教会(英语:Church of Scotland;低地苏格兰语:The Scots Kirk,常称 the Kirk ;苏格兰盖尔语:Eaglais na h-Alba)或称苏格兰长老会,是苏格兰的一个长老宗教会,是苏
  • 相对性状相对性状(英语:Relative character)是指一种生物的同一种性状的不同表现类型,这些性状是能够稳定地遗传给后代的。
  • 边大绶边大绶,字素一,号长白,北直隶任丘县(今河北省任丘市)人。明末政治人物。崇祯十二年(1639年)中举人,曾任陕西米脂县知县。陕西总督汪乔年下令边大绶在崇祯十五年(1642年)正月初八,掘毁李
  • 图像志图像志(德语:Ikonographie;英语:Iconography),又叫图像演变学、图像学,是艺术史的一个方法学。图像志早在希腊罗马时代的学术研究中意指“肖像学”,十六世纪以来则以“基督图像志”
  • 春晓油气田春晓油气田(日本称之为“白桦油气田”(白樺ガス田))是中国在东海大陆架盆地西湖凹陷中开发的一个大型油气田,位于浙江宁波东南约350公里,探明的天然气储量达700亿立方米以上,由中国