MATLAB辅助现代工程数字信号处理 第二版 教学课件 作者 李益华 第1 5章 第5章.pptVIP

第5章 多采样率数字信号处理 　　　　5.1 信号的抽取与内插5.1.1 信号的整数倍抽取　　在时域中，若原始信号为x(n)，抽取因子为D，抽取后得到的信号为y(n)，则整个信号的抽取过程可表示为 　　　　　　　　y(n)＝x(Dn) 　　　　　　　 (5.1) 　 信号的抽取过程可分为两个步骤: 第一步，将x(n)与一个周期为D的采样脉冲序列相乘，即每D个点中保留一个点，其他(D－1)个点为零，以便得到采样信号ω(n); 第二步，去掉ω(n)中的零点后，便得到一个低速率的信号y(n)。中间采样信号ω(n)为 　 　　抽取后的信号为　　　　　　　　　　　　y(n)=ω(Dn)　　　　　　　(5.3)　　【例 5.1】 设x(n)=2 sin(0.089πn)，抽取因子D为3，求抽取输出信号y(n)。　　MATLAB程序如下: 　　　%MATLAB PROGRAM 5－1　　　clc; 　　　n=[0: 49]; 　　　m=[0: 50*3－1]; 　　　x=2*sin(0.089*pi*m); 　　　y=x([1: 3: length(x)]); 　　　subplot(211); 　　stem(n，x(1: 50)); xlabel(′n′); ylabel(′x(n)′); 　　subplot(212); 　　stem(n，y); xlabel(′n′); ylabel(′y(n)′); 程序运行结果如图5.1所示。 图 5.1 信号的整数倍抽取时域分析 　　为分析下行抽样信号的频谱，先计算ω(n)的频谱为 　　【例 5.2】 利用MATLAB的函数fir2构造一个带限输入序列: f=[0， 0.45， 0.49，1]; 幅值mag=[0， 1， 0， 0]。 试分析信号抽取因子D=2的频域特性。　　MATLAB程序如下: 　　　　%MATLAB PROGRAM 5-2　　　　freq=[0， 0.45， 0.49， 1]; 　　　　mag=[0， 1， 0， 0]; 　　　　x=fir2(101，freq，mag); 　　　　%求取并绘制输入谱 　　　　[Xz，w]=freqz(x，1，512); 　　　　subplot(211); 　　　　plot(w/pi，abs(Xz)); 　　　　xlabel(′＼omega/＼pi′); ylabel(′|Xz(w)|′); grid; 　　　　y=x([1: 2: length(x)]); %产生抽取输出谱　　　　%求取并绘制抽取输出谱　　　　[Yz，w]=freqz(y，1，512); 　　　　subplot(212); 　　　　plot(w/pi，abs(Yz)); 　　　　xlabel(′＼omega/＼pi′); ylabel(′|Yz(w)|′); grid; 　　程序运行结果如图5.2所示。 图 5.2 整数倍抽取频域分析 5.1.2 信号的整数倍内插　　将采样率增加到I倍可以通过在信号的两个连续值之间内插(I－1)个样本来实现。内插过程可以由许多方式实现，下面介绍的是一个保持信号序列x(n)的谱形状不变的方式。在时域中，若已知输入信号x(n)，内插因子为I，插值后得到的信号为y(n)，则整个信号的内插过程可表示为 　　【例 5.3】 设x(n)=2 sin(0.26πn)，内插因子I为3，求内插输出信号y(n)。　　MATLAB程序如下:  　　%MATLAB PROGRAM 5-3　　　n=[0: 50]; 　　　x=2*sin(0.26*pi*n); 　　　y=zeros(1，3*length(x)); 　　　y([1: 3: length(y)])=x; 　　　subplot(211); 　　　stem(n，x); 　　　xlabel(′n′); ylabel(′ x(n) ′); 　　　subplot(212); 　　　stem(n，y(1: length(x))); 　　　xlabel(′n′); ylabel(′ y(n) ′); 　　程序运行结果如图5.3所示。 图 5.3 信号的整数倍内插时域分析 　　对输出信号进行傅里叶变换，得到信号y(n)的频谱为 　　　　[Xz，w]=freqz(x，1，512， ′whole′); 　　　　subplot(211); 　　　　plot(w/pi，abs(Xz)); 　　　　xlabel(′＼omega/＼pi′); ylabel(′|Xz(w)|′); grid; 　　　　%产生抽取输出谱　　　　y=zeros(1，