第九章多采样率数字信号处理
- 格式:ppt
- 大小:584.02 KB
- 文档页数:46
实验七 多采样率数字信号处理一.实验目的1.掌握信号抽取和插值的基本原理和实现2.掌握信号的有理数倍速率转换 二.实验内容例7.1 对信号进行抽取,使采样率为原来的1/4倍。
MATLAB 程序: t=0:.00025:1;x=sin(2*pi*30*t)+sin(2*pi*60*t); y=decimate(x,4);figure,subplot(2,2,1),stem(x(1:120)); title('原始信号时域图'),xlabel('(a)');subplot(2,2,2),plot(abs(fft(x))),title('原始信号频域图'),xlabel('(b)'); subplot(2,2,3),stem(y(1:30))title('抽样后的信号时域图'),xlabel('(c)'); subplot(2,2,4),plot(abs(fft(y)));title('抽样后的信号频域图'),xlabel('(d)');050100150-2-1012原始信号时域图(a)2000400060000100020003000原始信号频域图(b)0102030-2-1012抽样后的信号时域图(c)050010001500200400600抽样后的信号频域图(d)例 7.2 信号()⎪⎪⎭⎫ ⎝⎛=s f f nn x π3cos ,采样频率16/1=sf f,现将采样率提高为原来的4倍。
MATLAB 程序: n=0:30;x=cos(3*pi*n/16); y=interp(x,4);figure,subplot(2,2,1),stem(x);title('原始信号时域图'),xlabel('(a)');subplot(2,2,2),plot(abs(fft(x))),title('原始信号频域图'),xlabel('(b)'); subplot(2,2,3),stem(y(1:30))title('抽样后的信号时域图'),xlabel('(c)');subplot(2,2,4),plot(abs(fft(y)));title('抽样后的信号频域图'),xlabel('(d)');10203040-1-0.500.51原始信号时域图(a)102030405101520原始信号频域图(b)0102030-1-0.500.51抽样后的信号时域图(c)05010015020406080抽样后的信号频域图(d)例7.3 序列()()()2.12.1cos 5.0sin ++=n n n x ππ,调用resample 函数对想按因子3/7进行采样率转换,并绘出图形。
第9章 多采样率信号处理多采样率信号处理广泛应用于要求转换采样率,或要求系统工作在多采样率状态的信号处理系统中。
如多种媒体——语音、视频、数据的传输,它们的频率很不相同,采样率自然不同,必须实行采样率的转换;又如信号要在两个时钟频率的数字系统中传输时,为了便于信号的处理、编码、传输和存储,要求根据时钟频率对信号的采样率加以转换;再如一种信号处理算法在系统的不同部分采用不同的采样率(如子带编码等),使处理更加有效,等等。
本章首先介绍直接在数字域对离散时间信号进行采样率转换的抽取(Decimation )和内插(Interpolation )方法,然后讨论抽取滤波器与内插滤波器的设计与实现方法。
由不同采样率构成的系统称为多采样率系统,大部分多采样率系统使用了滤波器组,以正交镜像滤波器组(QMF )为基础的树状结构滤波器组是一典型的多采样率系统,它与离散小波(Wavelet )变换的关系密切。
本章最后由正交镜像滤波器组的概念引入小波变换的基本原理和多分辨率分析的概念,以利开拓思路,为进一步深入学习打下基础。
9.1 离散信号的抽取与内插9.1.1 抽取与内插的时域描述离散序列的抽取与内插是多采样率系统中的基本运算,抽取运算将降低信号的采样频率,内插运算将提高信号的采样频率。
离散序列()x k 的M 倍抽取定义为()(),D x k x Mk k =∈¢ (9-1-1)其中M 为一正整数。
抽取运算的框图如图9-1所示。
()x k ()D k图9-1 M倍抽取运算的框图图9-2画出了M =3时序列抽取的示意图。
由图可知,离散序列的抽取表示保留第M 个样本点,而去除两个样本之间的M-1个样本点,设原离散信号()x k 的采样周期为T ,经M 倍抽取后的信号()D x k 的采样周期为'T ,满足'T MT =。
为了强调此概念,在图9-2中,有意将抽取后的序列的间隔画为原序列的3倍。
这时新的采样频率's f 为Mf MT T f ss ===11'' (9-1-2) 式中,s f 为原有的采样频率。
目录1.背景 12.具体过程 22.1 整数因子抽取 22.2 整数因子内插 22.3 I/D的采样率转换 22.4多采样率数字信号处理的应用 23.实验过程 23.1整数倍抽取实验 23.2整数倍插值实验 23.3用有理因子I/D的采样率转换进行的实验 2 4.实验结果 24.1信号的整数倍抽取 24.2信号的整数倍插值 24.3用有理因子I/D的采样速率转换 25.结论 25.1整数因子抽取 25.2整数因子插值 25.3有理因子I/D的采样速率转换 26.心得体会与总结 21.背景现在实际系统中,经常要求一个数字系统能工作在多采样率状态,例如:在数字电视系统中,图像采集系统一般按4:4:4标准或4:2:2标准采集数字电视信号,再根据不同的电视质量要求将其转换成其它标准的数字电视信号(如4:2:2,4:1:1,2:1:1)进行处理。
在数字电话系统中,传输的信号既有语音信号又有传真信号,甚至有视频信号。
这些信号的频域成分相差甚远。
因此该系统应具有多种采样率,并能根据所传输的信号自动完成采样率转换。
对一个非平稳随机信号(如语音信号)做频谱分析或编码时,对不同的信号段可根据其频域成分的不同而采用不同的采样率,已到达既满足采样定理,又最大限度的减少数据量的目的。
如果以高采样率采集的数据存在冗余,这时就希望在该数字信号的基础上降低采样率。
多采样率数字信号处理是建立在单抽样率信号处理基础上的一类信号处理。
在传输信号时,由于语音﹑图像、视频信号的中心频率相差很大,所以需要以多种抽样频率来对信号采样来满足各种传输类型的需要。
2.具体过程2.1 整数因子抽取信号的抽取是实现频率降低的方法。
在第二章曾经讨论过,当采样频率大于信号最高频率的2倍时,不会产生混叠失真。
显然,当采样频率远高于信号最高频率时,采样后的信号就会有冗余数据。
此时,通过信号的抽取来降低采样频率,同样不会产生混叠失真。
Xd(n)整数因子抽取原理图:设x(n)=x(t)|t=nTs,欲使fs减少D倍,最简单的方法就是从x(n)中每D个点中抽取一个,依次组成一个新的序列xd(n),即xd(n)=x(Dn)因为是舍去部分点,故可引入冲激函数来进行抽样,得到xd(n)与x(n)之间的表达式:xd(n)=x(n) D(n)其中为周期单位脉冲序列,当且仅当n为D的整数倍时, D(n)的值为1,n为其他值时为零。
多采样率信号处理1.绪论随着数字信号处理的发展,信号的处理、编码、传输和存储等工作量越来越大。
为了节省计算工作量及存储空间,在一个信号处理系统中常常需要不同的采样率及其相互转换,在这种需求下,多速率数字信号处理产生并发展起来。
它的应用带来许多好处,例如:可降低计算复杂度、降低传输速率、减少存储量等。
在信号处理领域,多速率信号处理最早于20世纪70年代提出,由其引出的多速率滤波在数学领域里基于多格算法解决了大量的微分等式。
在多速率数字信号处理发展中,一个突破点是70年代两通道正交镜像滤波器组应用于语音信号的压缩。
在该方法中,信号通过分析滤波器组被分成低通和高通两个子带,每个子带经过2倍抽取和量化后再进行压缩,之后可以通过综合滤波器组近似地重建出原始信号,重建的近似误差一部分源于子带信号的压缩编码,一部分是由分析和综合滤波器组产生的误差,其中最主要的误差是混叠误差,它是由分析滤波器组不是理想带限而引起的。
在很多应用系统中,混叠误差存在一定程度的影响,因此就需要对其进行改进。
多速率系统应用于通信、语音信号处理、谱分析、雷达系统和天线系统,以及在数字音频系统、子带编码技术( 用于声音和图像的压缩) 和模拟语音个人系统(如标准电话通信) 等方面的应用。
另外还应用于多相理论和多速率系统在一些非传统领域,包括:高效率信号压缩的多速率理论;高效窄带滤波器的脉冲响应序列的编码新技术的推导;可调整的多级响应FIR滤波器的设计等。
基于上述研究的发展,从20世纪80年代初开始,多速率数字信号处理技术在工程实践中得到广泛的应用,主要用于通信系统、语音、图像压缩、数字音频系统、统计和自适应信号处理、差分方程的数值解等。
多速率信号处理在基础理论和应用领域的蓬勃发展,也促进了整个数字信号处理界的发展。
2.采样率转换基础理论实现采样率转换的方法有三个:一是若原模拟信号x (t)可以再生,或是己记录下来了的话,那么可重新抽样;二是将x(n)通过D/A变成模拟信号x(t)后,对x (t)经A/D再抽样;三是发展一套算法,对抽样后的数字信号x(n)在“数字域”作采样率转换,以得到新的抽样。