多采样率信号处理 信号的抽取与插值
第六部分：多速率信号处理
12
n
X D (e ) =
jw
k =−∞
∑
+∞
xD [k ]e− jwk
1
X (e jw )
因为
xD [n] = x p [nD ]
X D (e jw ) =
k =−∞
∑ x [kD]e
p
+∞
− jwk
−2π
−π
−ωM ωM
1 D
π
2π
ω
X p (e jw )
如果令n=kD，上式等效为 ， 如果令
D为周期的 为周期的 脉冲串采样
D倍抽取 倍抽取 表示、传输和存储这 个已采样序列是很不 经济的，因为在采样 点之间明知都是零
脉冲串采样过程
p[n] =
k =−∞
∑ δ [n − kD]
xp[n]
+∞
x[n]
x p [n] = x[n] p[n]
=
k =−∞
∑ x[kD]δ [n − kD]
+∞
− DωM
DωM π
2π
ω
由上图可知，已采样序列xp[n] 和抽取序列xD[n] 的频谱差别只是频率尺度上的或归一化上 抽取的效果是将原来序列的频谱扩展到一个较宽 的频带部分，这也反映了频域和时域之间的关系。 抽取相当于时域压缩，故频域会扩展 同时可以看出，如果要避免混叠，则：
DωM < π
即
ωM < π / D
取样率变换的多级实现
前面所讨论的取样率变换（抽取和内插），都是按 单级实现考虑的，即内插和抽取都一次完成。但 是实际中，当抽取倍数D和内插倍数I很大时，所需 的低通滤波器h[n]的阶数将非常高，乃至无法实现。 所以一个简单的想法就是通过多次小倍数的抽取和 内插完成
数字信号处理知识点汇总
数字信号处理知识点汇总数字信号处理是一门涉及多个领域的重要学科，在通信、音频处理、图像处理、控制系统等众多领域都有着广泛的应用。
接下来，让我们一同深入了解数字信号处理的主要知识点。
一、数字信号的基本概念数字信号是在时间和幅度上都离散的信号。
与模拟信号相比，数字信号具有更强的抗干扰能力和便于处理、存储等优点。
在数字信号中，我们需要了解采样定理。
采样定理指出，为了能够从采样后的信号中完全恢复原始的连续信号，采样频率必须至少是原始信号最高频率的两倍。
这是保证数字信号处理准确性的关键原则。
二、离散时间信号与系统离散时间信号可以通过序列来表示，常见的有单位脉冲序列、单位阶跃序列等。
离散时间系统则是对输入的离散时间信号进行运算和处理，产生输出信号。
系统的特性可以通过线性、时不变性、因果性和稳定性等方面来描述。
线性系统满足叠加原理，即多个输入的线性组合产生的输出等于各个输入单独作用产生的输出的线性组合。
时不变系统的特性不随时间变化，输入的时移会导致输出的相同时移。
因果系统的输出只取决于当前和过去的输入，而稳定系统对于有界的输入会产生有界的输出。
三、Z 变换Z 变换是分析离散时间系统的重要工具。
它将离散时间信号从时域转换到复频域。
通过 Z 变换，可以方便地求解系统的差分方程，分析系统的频率特性和稳定性。
Z 变换的收敛域决定了其特性和应用范围。
逆 Z 变换则可以将复频域的函数转换回时域信号。
四、离散傅里叶变换（DFT）DFT 是数字信号处理中的核心算法之一。
它将有限长的离散时间信号转换到频域。
DFT 的快速算法——快速傅里叶变换（FFT）大大提高了计算效率，使得在实际应用中能够快速处理大量的数据。
通过 DFT，可以对信号进行频谱分析，了解信号的频率成分和能量分布。
五、数字滤波器数字滤波器用于对数字信号进行滤波处理，分为有限冲激响应（FIR）滤波器和无限冲激响应（IIR）滤波器。
FIR 滤波器具有线性相位特性，稳定性好，但设计相对复杂。
数字信号处理实验(民航无线电监测关键技术研究)
《数字信号处理》实验报告实验名称数字信号处理实验（民航无线电监测关键技术研究）实验时间一、实验目的：通过实验，理解和掌握民航无线电监测关键技术中调制解调、FIR 数字滤波器、多采样率数字信号处理、FFT、语音数字信号处理、静噪等技术，培养学生对数字信号处理技术的兴趣，并提高学生基于数字信号处理技术的工程应用能力。
二、实验环境：Matlab三、实验原理、内容与分析（包括实验内容、MATLAB程序、实验结果与分析）实验总体框图如上图所示，主要实现民航无线电监测关键技术中调制解调、FIR 数字滤波器、多采样率数字信号处理、FFT、语音数字信号处理、静噪等技术。
1．有限长单位脉冲（FIR）滤波器的设计FIR 数字滤波器是一种非递归系统，其冲激响应h(n)是有限长序列,其差分方程表达式为：系统传递函数可表达为：N-1 为FIR 滤波器的阶数。
在数字信号处理应用中往往需要设计线性相位的滤波器，FIR 滤波器在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性。
为了使滤波器满足线性相位条件，要求其单位脉冲响应h(n)为实序列，且满足偶对称或奇对称条件，即h(n)=h(N-1-n)或h(n)=-h(N-1-n)。
这样，当N 为偶数时，偶对称线性相位FIR 滤波器的差分方程表达式为：由上可见FIR 滤波器不断地对输入样本x(n)延时后，再做乘法累加算法，将滤波器结果y(n)输出，因此，FIR 实际上是一种乘法累加运算。
而对于线性相位FIR 而言，利用线性相位FIR 滤波器系数的对称特性，可以采用结构精简的FIR 结构将乘法器数目减少一半。
2．AM 调制解调AM 调制解调过程如下：3．多采样率数字信号处理一般认为，在满足采样定理的前提下，首先将以采样率F1 采集的数字信号进行D/A 转换, 变成模拟信号，再按采样率F2 进行A/D 变换，从而实现从F1 到F2 的采样率转换。
但这样较麻烦，且易使信号受到损伤，所以实际上改变采样率是在数字域实现的。
数字信号处理综合实验
通信与信息工程学院数字信号处理课程设计班级：电子信息工程13级03班姓名：学号：指导教师：张释如、李国民、张龙妹、王瑜设计时间：2015.12.28 --- 2016.1.8成绩：评语：通信与信息工程学院二〇一五年数字信号处理课程设计报告一、课程设计时间2015年12月28日至2016年1月8日二、课程设计目的数字信号处理主要研究如何对信号进行分析、变换、综合、估计与识别等加工处理的基本理论和方法。
通过课程设计，使学生巩固所学基本理论，掌握最基本的数字信号处理的理论和方法，提高综合运用所学知识，提高计算机编程的能力。
进一步加强学生独立分析问题、解决问题的能力、综合设计及创新能力的培养，同时注意培养学生实事求是、严肃认真的科学作风和良好的实验习惯，为今后的工作打下良好的基础。
三、课程设计任务及要求1、掌握数字信号处理IIR滤波器设计及FIR滤波器设计原理和实现，能根据不同的应用设计合理的滤波器；2、掌握多频率采样的原理，并能分析其频谱特性；3、了解语音信号处理的原理，并能根据实际情况设计合理的滤波器进行除燥处理；3、编程实现以下实验内容：（1）数字信号的基本运算（2）多采样率数字信号处理（3）数字滤波器的设计及仿真（4）语音信号滤波处理。
一、数字信号的基本运算一、实验目的：(1) 掌握数字信号的时间翻转、上采样、下采样等基本运算；(2) 学会用MATLAB对数字信号进行时间翻转、上采样、下采样等运算；二、设计内容：(1) 利用Windows下的录音机以采样频率8000Hz录制语音“新年好”和“好”，在Matlab 软件平台下，利用wavread函数得到两个语音数据（信号长度不够时信号补零使其长度为8000）；(2) 对采样得到的语音数据x(k)分别进行处理模仿回音效果，演示回声的效果，数据处理如下式：x(k)=x(k)+a*x(k-d)其中d为时延，a为时延信号的衰减幅度。
(参数：时延d=0.4秒，衰减幅度a=0.5 对上述语音信号进行时间反转x(-k)、上采样x(k/2)、下采样x(2k)操作，并演示运算效果。
插值fir滤波器工作原理
插值fir滤波器工作原理
插值FIR滤波器是一种数字滤波器，它通过插入零值来增加输入信号的采样率，并且以此增加输出信号的频率分辨率。
它的工作原理如下：
1. 零值插入：首先，插值FIR滤波器会将输入信号的每个采样值之间插入多个零值，以增加采样率。
这样可以在频域中扩展输入信号，使其能够更好地与滤波器的频域响应匹配。
2. 频域响应设计：根据所需的滤波特性，设计一个合适的频域响应，例如低通、带通、高通或带阻滤波器等。
这通常通过选择合适的滤波器系数来实现。
3. 卷积运算：将插值后的输入信号序列与滤波器的系数序列进行卷积运算，得到输出信号序列。
4. 采样返回：最后，将输出信号进行抽取，以还原到原始的采样率，得到插值FIR滤波器的输出信号。
总的来说，插值FIR滤波器通过零值插入、频域响应设计、卷积运算和采样返回等步骤，实现对输入信号的插值和滤波处理。
它可以用于信号重构、频率转换、滤波等应用。
[离散时间信号处理学习笔记]14.多采样率信号处理
![[离散时间信号处理学习笔记]14.多采样率信号处理](https://img.taocdn.com/s3/m/201a0352a9956bec0975f46527d3240c8447a1a2.png)
[离散时间信号处理学习笔记]14.多采样率信号处理多采样率信号处理⼀般是指利⽤增采样、减采样、压缩器和扩展器等⽅式来提⾼信号处理系统效率的技术（These multirate techniques refer in general to utilizing upsampling, downsampling, compressors, and expanders in a variety of ways to increase the efficiency of signal-processing systems. ）本⽂章主要讨论多采样率技术中的两个研究成果：滤波与压缩器/扩展器的互换；多相分解。
尽管上⼀篇⽂章中已经讨论过这部分内容，不过由于这部分是理解本⽂所必须的关键知识点，这⾥将在时域与频域展开更详细的分析。
压缩器假设压缩器的压缩率为M，那么压缩器在时域上的表⽰为x_d[n] = x[nM]x[n]的采样频率为T，那么x_d[n]的采样频率为T_d = MT，按照，有\begin{align*} X(e^{j\omega}) &= \frac{1}{T}\sum_{k=-\infty}^{\infty}X_c\left[ j\left(\frac{\omega}{T}-\frac{2\pi k}{T}\right)\right ]\\ X_d(e^{j\omega}) &= \frac{1}{MT}\sum_{r=-\infty}^{\infty}X_c\left[ j\left(\frac{\omega}{MT}-\frac{2\pi r}{MT}\right)\right ] \end{align*}压缩前的序列频谱X(e^{j\omega})与压缩后的序列频谱X_d(e^{j\omega})之间有如下关系\begin{align*} X_d(e^{j\omega}) &= \frac{1}{MT}\sum_{r=-\infty}^{\infty}X_c\left[ j\left(\frac{\omega}{MT}-\frac{2\pi r}{MT}\right)\right ] \\ & = \frac{1}{MT}\left\{\cdot\cdot\cdot+X_c\left[j\left(\frac{\omega}{MT}-\frac{-2\pi}{MT} \right ) \right ] +X_c\left[ j\left(\frac{\omega}{MT}-\frac{0} {MT}\right)\right ] + X_c\left[j\left(\frac{\omega}{MT}-\frac{2\pi}{MT} \right ) \right ]+\cdot\cdot\cdot \right \}\\ & = \frac{1}{MT}\left\{\cdot\cdot\cdot+X_c\left[j\left(\frac{\omega}{MT}-\frac{0}{MT} \right ) \right ]+\cdot\cdot\cdot +X_c\left[ j\left(\frac{\omega}{MT}-\frac{2(M-1)\pi} {MT}\right)\right ]\right.\\ &\quad\qquad\qquad\left.+ X_c\left[j\left(\frac{\omega}{MT}-\frac{2M\pi}{MT} \right ) \right ]+\cdot\cdot\cdot+X_c\left[j\left(\frac{\omega}{MT}-\frac{2M\pi}{MT} -\frac{2(M-1)\pi}{MT}\right ) \right ]+\cdot\cdot\cdot \right \}\\ \end{align*} \begin{align*} \qquad\quad\ &= \frac{1}{MT}\left\{\cdot\cdot\cdot+\sum_{i=0}^{M-1}X_c\left[j\left(\frac{\omega}{MT}-\frac{2i\pi}{MT} \right ) \right ]+\sum_{i=0}^{M-1}X_c\left[j\left(\frac{\omega}{MT}-\frac{2i\pi}{MT}-\frac{2\pi}{T} \right ) \right ]+\cdot\cdot\cdot\right\}\\ &= \frac{1} {MT}\sum_{k=-\infty}^{\infty} \sum_{i=0}^{M-1}X_c\left[j\left(\frac{\omega}{MT}-\frac{2\pi i}{MT}-\frac{2\pi k}{T} \right ) \right ] \\ &=\frac{1} {M}\sum_{i=0}^{M-1}\left\{\frac{1}{T}\sum_{k=-\infty}^{\infty}X_c\left[j\left(\frac{\omega-2\pi i}{MT}-\frac{2\pi k}{T} \right ) \right ]\right\}\\&=\frac{1}{M}\sum_{i=0}^{M-1}X(e^{j(\omega-2\pi i)/M}) \end{align*}如下图所⽰扩展器假设扩展器的扩展率为L，那么扩展器在时域上的表⽰为x_e[n] = \left\{\begin{matrix} x[n/L], &n=0,\pm L,\pm 2L,\cdot\cdot\cdot \\ 0, &else \end{matrix}\right.扩展前的序列频谱X(e^{j\omega})与扩展后的序列频谱X_e(e^{j\omega})之间有如下关系\begin{align*} X_e(e^{j\omega}) &= \sum_{n=-\infty}^{\infty}x_e[n]e^{-j\omega n}\\ &=\sum_{n=-\infty}^{\infty}x[n/L]e^{-j\omega n}\quad n=0,\pm L,\pm 2L,\cdot\cdot\cdot\\ &=\sum_{k=-\infty}^{\infty}x[k]e^{-j\omega kL}\quad letting\ n=kL\\ &=X(e^{j\omega L}) \end{align*}如下图所⽰滤波器与压缩器互换如上⼀篇⽂章所描述的减采样就是⼀个滤波器与压缩器的级联系统。
数字信号处理的常见问题及解决方法总结
数字信号处理的常见问题及解决方法总结数字信号处理在科学、工程领域中的应用越来越广泛。
在实际应用过程中，我们常常会遇到一些问题。
本文总结了一些常见的问题及其解决方法，以帮助读者更好地理解和应用数字信号处理技术。
问题一：信号滤波数字信号往往包含噪声和干扰，需要进行滤波处理以提取有效信息。
常见的信号滤波问题包括滞后滤波器、移动平均滤波器、低通滤波器等。
解决这些问题的方法通常包括设计合适的滤波器参数、选择适当的滤波器类型，并进行滤波器性能评估。
问题二：信号采样率选择在数字信号处理中，采样率的选择对信号重构和频谱分析等方面具有重要影响。
选择过低的采样率会导致信号失真，选择过高的采样率会浪费存储和计算资源。
解决这个问题的方法包括根据信号的带宽和特性选择合适的采样率，并根据需要进行抽取或插值处理。
问题三：频谱分析频谱分析是数字信号处理中的重要步骤，用于研究信号的频域特性。
常见的频谱分析问题包括功率谱密度估计、傅里叶变换等。
解决这些问题的方法包括选择合适的频谱分析方法（如快速傅里叶变换）、处理频谱分辨率问题，并进行频谱分析结果的解释和应用。
问题四：数字滤波器设计数字滤波器的设计是数字信号处理领域的关键问题之一。
常见的数字滤波器设计问题包括低通滤波器设计、高通滤波器设计、带通滤波器设计等。
解决这些问题需要根据滤波器的要求和性能指标，选择适当的设计方法（如窗函数法、频率抽样法），并进行滤波器参数调整和性能评估。
问题五：数字信号压缩数字信号压缩是在保证信号质量的前提下，减少信号数据量的一种技术。
常见的数字信号压缩问题包括有损压缩和无损压缩。
解决这些问题的方法通常包括选择适当的压缩算法（如哈夫曼编码、小波变换压缩），根据压缩效率和信号质量要求进行参数调整。
以上是数字信号处理中常见问题的一些总结及解决方法。
希望能够帮助读者更好地应用数字信号处理技术，解决实际应用中的问题。
采样率转换程序
采样率转换程序采样率是指在数字信号处理中，对模拟信号进行采样的频率。
它是一个非常重要的参数，决定了数字信号的质量和精度。
在实际应用中，我们经常需要对采样率进行转换，以适应不同的需求和系统要求。
本文将介绍采样率转换程序的原理和实现方法。
一、采样率转换的原理采样率转换主要是通过插值和抽取的方式实现的。
插值是指在已有的采样点之间插入新的采样点，从而增加采样率；抽取是指在已有的采样点之间舍弃部分采样点，从而降低采样率。
这样就可以通过插值和抽取的组合，实现不同采样率之间的转换。
二、采样率转换的实现方法1. 线性插值法：线性插值法是最简单的插值方法之一。
它通过连接相邻的采样点，按照一定的比例插入新的采样点。
这种方法简单易懂，计算量小，但是会引入一定的误差。
2. 多项式插值法：多项式插值法是一种更精确的插值方法。
它通过拟合已有的采样点，构造一个多项式函数，并在函数中插入新的采样点。
这种方法计算量较大，但精度更高。
3. 快速傅里叶变换（FFT）：FFT是一种高效的频域分析方法，也可以用于采样率转换。
它将信号从时域转换到频域，然后通过修改频域上的采样点，再通过反变换将信号恢复到时域。
这种方法计算量较大，但可以实现高质量的采样率转换。
三、采样率转换程序的实现采样率转换程序的实现需要根据具体的需求选择合适的算法和方法。
一般情况下，我们可以使用现成的信号处理库或软件包来实现采样率转换。
这些库和软件包中已经实现了各种采样率转换算法，并提供了相应的接口和函数，方便我们进行调用和使用。
在使用采样率转换程序时，我们需要提供两个参数：输入信号的采样率和目标采样率。
程序会根据这两个参数，选择合适的算法和方法进行转换。
转换后的信号可以保存到文件中，或者直接输出到音频设备中。
四、采样率转换的应用采样率转换在很多领域都有广泛的应用。
例如，在音频处理中，我们经常需要将不同采样率的音频进行转换，以适应不同的播放设备或系统要求。
在图像处理中，采样率转换可以用于图像的放大和缩小，以及图像的压缩和解压缩等。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
5
H ( z) =
记 则 若再记 则
l =0
信号的多相表示
M −1
给定序列h(n)，令 n = 0 ~ ∞ ，有 给定序列 ，
z −l ∑ h ( Mn + l ) z − Mn ∑
∞
n =0
∞
E l ( z ) = ∑ h ( Mn + l ) z − n
n =0 =0
H ( z) =
NCEPUBD
时域
2.2 先滤波再抽取
1 H (e ) = 0
jω
h(n)为一理想低通滤波器 为一理想低通滤波器: 为一理想低通滤波器
| ω |≤ π M else
∞
滤波后的输出为
υ (n) =
∞
k = −∞
∑ h(k ) x (n − k )
∞
抽取后的序列为y(n) : 对v(n)抽取后的序列为 抽取后的序列为
NCEPUBD
1
1.3
引 言
研究内容
核心内容:信号抽样率的转换及滤波器组。 核心内容:信号抽样率的转换及滤波器组。
信号的“抽取 减少抽样率以 信号的“抽取(decimatiom) ” ：减少抽样率以 去掉过多数据 信号的“插值 增加抽样率以 信号的“插值(interpolation) ” :增加抽样率以 增加 增加数据 分析滤波器组和综合滤波器组 滤波器组:分析滤波器组和综合 滤波器组 分析滤波器组和综合滤波器组
jω jω jω jω jω L
所以应取c=L以保证 以保证y(n)=x(n/L) 所以应取 以保证
NCEPUBD
时域
3.3 先插值再滤波
y ( n) = υ ( n) * h( n) = ∑ υ ( k ) h( n − k )
k
= ∑ x ( k L ) h( n − k )
k
即
y (n) =
Ql ( z ) = ∑ h ( Mn − l ) z −n
n =0
∞
H ( z) =
显然
M −1 l =0
∑ z Q (z
l l
M
)
Ql ( z ) = z E M −l ( z )
NCEPUBD
−1
类 型 III 多 相 表 示
6
几个重要的恒等关系
多采样率信号 处理
——信号的抽取与插值
孙正 2007.9
NCEPUBD
主要内容
• • • • • • • • 引言 信号的抽取 信号的插值 抽取与插值相结合的抽样率转换 信号的多相表示 几个重要的恒等关系 抽取和插值的滤波器实现 抽取与插值的编程实现
学习要求：掌握数字域升、降采样的基本原理， 学习要求：掌握数字域升、降采样的基本原理， 降采样过程中滤波器的设计、特性和作用， 升、降采样过程中滤波器的设计、特性和作用， 以及插值和抽取前后信号频谱变化。 以及插值和抽取前后信号频谱变化。
j ( ω − 2 πk )
) H (e
M
)
在(-π/M～π/M)内, 抽取后信号的频谱与原信号频 ～ 内 谱只是幅度相差M倍 谱只是幅度相差 倍。
NCEPUBD
3
• 插值的概念
信号的插值
• 插零后的信号及其频谱 • 先插值再滤波
NCEPUBD
3.1 插值的概念
设 x(n) = x(t ) |t = nTs ，将x(n)中每两个点之中补 中每两个点之中补 L-1个0，组成一个新的序列 个 ，组成一个新的序列v(n)，即 ，
π π
该滤波器既去除了 插值后的映像又防 止了抽取后的混迭
使用2 使用2个低通滤波器
使用1 使用1个低通滤波器
NCEPUBD
时域上x(n)和y(n)的关系 4.1 时域上 和 的关系
y(n ) =
k = −∞
∑ x(k )h( Mn − Lk )
∞
y(n)正是单独抽取和单独插值时时域关系的结合 正是单独抽取和单独插值时时域关系的结合 因为h(n)是因果滤波器，所以 是因果滤波器， 因为 是因果滤波器 Mn − Lk ≥ 0 k ≤ ( M L) n 记
g (n, m) = h(nL + < Mm > L )
因为
−∞ < n, m < ∞
g (n, m + kL) = h(nL + < Mm + kML > L ) = h(nL + < Mm > L )
所以g(n,m)是变量 的周期函数，周期为 。 是变量m的周期函数 周期为L。 所以 是变量 的周期函数，
jω
)
插零后信号的频谱V(ejω)在(-π/L~π/L)内等于 jω),相 内等于X(e 相 插零后信号的频谱 在 内等于 当于将X(e 作了周期压缩 换句话说，就是V(e 在 作了周期压缩。 当于将 jω)作了周期压缩。换句话说，就是 jω)在 (-π~π)内包含了 个X(ejω)的压缩样本。 内包含了L个 的压缩样本。 内包含了 的压缩样本
M −1 l =0
∑z
− ( M −1− l )
Rl ( z )
M
式中 Rl ( z ) = E M −1−l ( z ) = ∑ h ( Mn + M − 1 − l ) z − n
n =0
∞
类 型 II 多 相 表 示
NCEPUBD
5
信号的多相表示
代替类型I中的 用-l代替类型 中的 ，则有 代替类型 中的l，
NCEPUBD
频域上x(n)和y(n)的关系 4.2 频域上 和 的关系
L M −1 1 Y (e jω ) = U (e j ( ω − 2π k ) M ) = M ∑ M k =0
M −1
X (e j (ω L − 2π k ) M ) | ω |≤ min( π , π ) ∑ k =0 M L 其它 0
NCEPUBD
1
1.2
应用举例： 应用举例：
引 言
研究目的
既可传输一般的语音信号， 既可传输一般的语音信号，也可传输播视频信 号的数字传输系统； 号的数字传输系统； 在音频世界，存在着多种抽样频率； 在音频世界，存在着多种抽样频率； 当需要将数字信号在两个具有独立时钟的数字 系统之间传递时； 系统之间传递时； 对信号（如语音，图象）作谱分析或编码时； 对信号（如语音，图象）作谱分析或编码时； 对一个信号抽样时，若抽样率过高， 对一个信号抽样时，若抽样率过高，必然会造 成数据的冗余； 成数据的冗余；
k = [ Mn L] − m
取整
Mn − [ Mn L]L = Mn mod L = Mn
L
则
y (n) = ∑ x([ Mn L] − m)h(mL + Mn L )
m =0
∞
Mn对模 求余 对模L求余 对模
NCEPUBD
时域上x(n)和y(n)的关系 4.1 时域上 和 的关系
y(n)可以看作是将 可以看作是将x(n)通过一个时变滤波器所 通过一个时变滤波器 可以看作是将 通过一个时变滤波器所 得到的输出。 得到的输出。记该时变系统的单位抽样响应为
NCEPUBD
1
• 研究背景 • 研究目的 • 研究内容
引
言
NCEPUBD
1
1.1
引 言
研究背景
至今， 至今，我们讨论的数字系统中只有一个 抽样率。 抽样率。 但是，在实际应用中， 但是，在实际应用中，各系统之间的采 样率往往是不同的
NCEPUBD
1
1.2
引 言
研究目的
要求一个数字系统能工作在“ 要求一个数字系统能工作在“多抽样率 （multirate）”状态，以适应不同抽样 ） 状态， 信号的需要。 信号的需要。 对一个数字信号， 对一个数字信号，能在一个系统中以不 同的抽样频率出现。 同的抽样频率出现。
NCEPUBD
2
信号的抽取
• 抽取对信号频谱的影响 • 先滤波再抽取
NCEPUBD
2.1 抽取对信号频谱的影响
中每M个点中抽取一个 设 x (n) = x (t ) |t = nTs ，将x(n)中每 个点中抽取一个， 中每 个点中抽取一个， 依次组成一个新的序列y(n)，即 ， 依次组成一个新的序列
NCEPUBD
5
意义： 意义：
信号的多相表示
使用多相表示可在抽样率转换的过 程中去掉许多不必要的计算， 程中去掉许多不必要的计算，因而大大 提高运算速度。 提高运算速度。
NCEPUBD
5
∞ n =0
信号的多相表示
给定序列h(n)，令 n = 0 ~ ∞ ，假定M=4 给定序列 ， 假定
H ( z ) = ∑ h( n) z − n = h0 + h4 z −4 + h8 z −8 + h12 z −12 + L
NCEPUBD
2.1 抽取对信号频谱的影响
x (t )
x (n )
抽样
x (n ) y ( n)
保证 f s ≥ 2 f c 不会发生频谱的混迭
M倍抽取 倍抽取
保证 f s ≥ 2 Mf c 不会发生频谱的混迭
是可变的， 出现混迭， 若M是可变的，为防止抽取后在Y (e jω )出现混迭，应对 是可变的 x(n)抽取前先作低通滤波，压缩其频带。 抽取前先作低通滤波 抽取前先作低通滤波，压缩其频带。
y (n) = υ ( Mn) =
k =−∞
∑ h(k ) x( Mn − k ) = ∑ x(k )h( Mn − k )
k =−∞
NCEPUBD
2.2 先滤波再抽取
NCEPUBD
频域
2.2 先滤波再抽取
M −1 k =0