数字信号处理实验报告三.docVIP

数字信号处理实验报告 实验三 抽样 班级 学号：2011081418 姓名：孙明凤 日期：2013年11月30日 指导老师：田园 成绩： 实验三 抽样 实验要求：研究抽样过程，分析产生混叠效应的原因，实现不同重建方案。 实验目的： 1、加深理解连续时间信号的离散化过程中的数学概念和物理概念，掌握时域抽样定理的基本内容； 2、掌握由抽样序列重建原连续信号的基本原理与实现方法； 3、加深理解频谱离散化过程中的数学概念和物理概念，掌握频域抽样定理的基本内容。 实验原理： 本部分说明抽样过程的两个基本原理：混叠和重建，涉及正弦波和线性调频信号的混叠。然后使用DTFT在频域中展开混叠过程。最后实践几种可从其抽样样本把信号恢复的不同手段。 3.3.1抽样引起的混叠 由于在MATLAB中不能产生模拟信号，所以需要做实时t轴的仿真。因此，把仿真时的△t与所研究的抽样周期Ts明确地区分开是很重要的。 程序： a. nn=0:79; sinus=sin(2*pi*300/8000*nn); stem(nn,sinus) b. plot(nn,sinus) c.把正弦的频率从100hz变至475hz nn=0:7 9; sinus1=sin(2*pi*100/8000*nn); subplot(221) stem(nn,sinus1) sinus2=sin(2*pi*225/8000*nn); subplot(222) stem(nn,sinus2) sinus3=sin(2*pi*350/8000*nn); subplot(223) stem(nn,sinus3) sinus4=sin(2*pi*475/8000*nn); subplot(224) stem(nn,sinus4) 正如预期的那样，正弦信号频率在逐渐增加。 d.把正弦的频率从7525Hz变至7900Hz nn=0:79; sinus1=sin(2*pi*7525/8000*nn); subplot(221) stem(nn,sinus1) sinus2=sin(2*pi*7650/8000*nn); subplot(222) stem(nn,sinus2) sinus3=sin(2*pi*7775/8000*nn); subplot(223) stem(nn,sinus3) sinus4=sin(2*pi*7900/8000*nn); subplot(224) stem(nn,sinus4) 正弦信号频率在逐渐减少。 e.把正弦的频率从32100Hz变至32475Hz，每次增加125Hz nn=0:79; sinus1=sin(2*pi*32100/8000*nn); subplot(221) stem(nn,sinus1) sinus2=sin(2*pi*32225/8000*nn); subplot(222) stem(nn,sinus2) sinus3=sin(2*pi*32350/8000*nn); subplot(223) stem(nn,sinus3) sinus4=sin(2*pi*32475/8000*nn); subplot(224) stem(nn,sinus4) 如预期的那样，正弦信号频率在逐渐增加。 实验分析：当正弦波的频率逐渐接近采样频率的整数倍时，其绘制出的波形频率减少，反之则增大。 3.3.2抽样的频域视图 对连续时间信号抽样时，因为频域以抽样频率延拓，故其频谱显现出混叠效应。为了实际演示这一效应，需要使用示波器。在MATLAB中，只能仿真这一效应，而这正是本节的目的。 为了仿真模拟信号，必须用非常高的抽样率——至少是任何模拟信号所允许具有的最高频率的5倍。因此在这个问题中将会有两个“抽样率”，一个用于所研究的实际抽样中，另一个用于仿真连续信号。第二个问题是如何显示连续时间信号的傅里叶变换。同样，这也只能仿真实现。下面的M文件应该用来绘出模拟频谱。 function fmagplot(xa,dt) L=length(xa); Nfft=round(2.^round(log2(5*L))); Xa=fft(xa,Nfft); range=0:(Nfft/4); ff=range/Nfft/dt; plot(ff/1000,abs(Xa(1:length(range)))) title(CONT-TIME FOURIER T