北交大数字信号处理4多速率要点.doc

《数字信号处理》课程研究性学习报告 试点班专用 多速率专题 姓名 学号 同组成员 指导教师 陈后金 时间 多速率信号处理专题研讨 【目的】 (1) 掌握序列抽取运算与内插运算的频谱变化规律。 (2) 掌握确定抽取滤波器与内插滤波器的频率指标。 (3) 掌握有理数倍抽样率转换的原理及方法。 (4) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 基本题 1．抽取、内插信号特征的时域/频域分析 对于给定的单频模拟信号y(t)=sin(1000?t)，确定一个合适的采样频率fsam，获得离散信号y[k]，试进行以下问题的分析： (1) 对离散信号y[k]进行M=2倍抽取，对比分析y[k]和y[Mk]在时域/频域的关系； (2) 对离散信号y[k]进行L=2倍内插，对比分析y[k]和y[k/L]在时域/频域的关系。 【温磬提示】 在多速率信号分析中，离散序列的抽取和内插是多速率系统的基本运算，抽取运算将降低信号的抽样频率，内插运算将提高信号的抽样频率。两种运算的变换域描述中，抽取运算可能出现频谱混叠线性，而内插运算将出现镜像频谱。 【设计步骤】 1.对所给定信号进行抽样，由于所给定信号为y(t)=sin(1000?t)，其频率f=500Hz,且实验过程中我们发现，当T取较大值时，Matlab 绘出的图像会有较大误差，经查询发现这是由于Matlab的数值计算特点所决定的，于是我们减小了T的值，同时加大了采样频率fsam, ，使其能产生理想的图像。 2. 对离散信号进行M=2倍的抽取时，我们通过xD=x(1:M:end); 进行MATLAB的仿真运算。 3.对离散信号进行L=2倍内插时，我们通过xL=zeros(1,L*length(x));xL(1:L:end)=x; 进行MATLAB的仿真运算。 4.对于各离散信号的频谱计算，我们都采用了FFT进行计算。 【仿真结果】 对离散信号y[k]抽取和内插的时域/频域对比分析结果。 【结果分析】 抽取、内插运算在频域的描述和性质。 通过图像我们可以发现，在时域内，离散序列的抽取表示每隔M-1点抽取原序列中的一点；离散信号的内插是在原序列每两个样本点之间插入L-1个零值样本点。 在频域内，M倍抽取后序列的频谱为将原序列频谱扩展M倍，得到，再将其周期化后的 函数相加并乘以因子1/M，即可得到抽取后序列的频谱。 由于原序列频谱在区间[-π，π]范围内满足，所以抽取后信号的频谱不会发生混叠。 信号L倍内插的频谱为将原序列频谱压缩L倍得到，除了与原序列频谱相差一个尺度因子外，两个频谱的形状保持不变。由于的周期为2π，所以的周期为2π/L。内插序列在区间[-π/L，π/L]内的频谱将在区间[-π，π]内重复L次，此重复的部分为镜像频谱 【自主学习内容】 【阅读文献】 [1] 陈后金.数字信号处理[M].北京：高等教育出版社.2008.11 【发现问题】 (专题研讨或相关知识点学习中发现的问题)： 当连续信号频率过大时，T值较大时MATLAB绘图会产生较大误差 【问题探究】 【仿真程序】 抽取： T=0.004; fs=8000; Ts=1/fs; N=T/Ts; x=zeros(1,N); t=0:N-1; x=sin(1000*pi*(t*Ts)); subplot(2,2,1); stem(t*Ts,x); xD=x(1:2:end); subplot(2,2,3); stem(xD); X=fft(x,32); X=fftshift(X); subplot(2,2,2); plot(abs(X)); Y=fft(xD,16); Y=fftshift(Y); subplot(2,2,4); plot(abs(Y)); 内插： T=0.004; fs=8000; Ts=1/fs; N=T/Ts; x=zeros(1,N); t=0:N-1; x=sin(1000*pi*(t*Ts)); subplot(2,2,1); stem(t*Ts,x); xL=zeros(1,2*length(x)); xL(1:2:end)=x; subplot(2,2,3); stem(xL); X=fft(x,32); X=fftshift(X); subplot(2,2,2); plot(abs(X)); Y=fft(xD,64); Y=fftshift(Y); subplot(2