时域抽样与频率抽样精讲.doc

本科学生实验报告 学号 姓名 学院 专业、班级 实验课程名称 数字信号处理 教师及职称 开课学期 2014 至 2015 学年 下 学期 填报时间 2015 年 3 月 28 日 云南师范大学教务处编印 实验设计方案 实验序号 二 实验名称 时域抽样与频域抽样 实验时间 2015年3月25日 实验室 同析312 实验目的 加深理解连续时间信号离散化的过程中的数学概念和物理概念，掌握时域抽样定理的基本内容。掌握由抽样序列重建原连续信号的基本原理与实现方法，理解其工程概念。加深理解频谱离散的过程中的数学概念和物理概念，掌握频谱抽样定理的基本内容。 2． 实验原理、实验流程或装置示意图 （1）离散系统在处理信号时，信号必须是离散的序列。因此，在利用计算机等离散系统分析处理连续时间信号时必须对信号进行离散化处理。时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件：对于基带信号，信号抽样频率大于等于2倍的信号最高频率。信号的重建是信号抽样的逆过程。 （2）非周期离散化信号的频率是连续谱。计算机在分析离散子信号的频谱时，必须将其连续频谱离散化。频域抽样定理给出了连续频谱抽样过程中信号不失真的约束条件。 【例题1.2.1】 利用MATLAB实现对信号的抽样。 [解] 绘出该连续信号 t0=0:0.001:0.1; x0=cos(2*pi*20*t0); plot(t0,x0,r) hold on%信号最高频率为20Hz，绘出按100Hz抽样频率进行得到的离散序列 Fs=100; t=0:1/Fs:0.1; x=cos(2*pi*20*t); stem(t,x); hold off title(连续信号函数及其抽样信号) 图1.2.1 连续信号及其抽样信号 【例1.2.2】 已知序列 ，对其频谱进行抽样，分别取N=2,3，10，观察频域抽样造成的混叠现象。 [解] x=[1,1,1];L=3; N=256; omega=[0:N-1]*2*pi/N; X0=1+exp(-j*omega)+exp(-2*j*omega); figure(1); plot(omega./pi,abs(X0)); x label(Omega/PI);hold on N=2;omegam=[0:N-1]*2*pi/N; Xk=1+exp(-j*omegam)+exp(-2*j*omegam); stem(omegam./pi,abs(Xk),r,o);hold off 运行结果如图1.2.2所示，可以看到的频谱及其频率抽样点 图1.2.2 序列频谱及其频率抽样点N=2 当N=3时，运行结果如图1.2.3所示 图1.2.3 序列频谱及其频率抽样点 N=3 当N=10时，运行结果如图1.2.4所示 图1.2.4 序列频谱及其频率抽样点 N=10 由频域抽样点恢复时域信号： x=[1,1,1];L=3; N=256; omega=[0:N-1]*2*pi/N; X0=1+exp(-j*omega)+exp(-2*j*omega); figure(1); plot(omega./pi,abs(X0)); xlabel(Omega/PI);hold on N=2;omegam=[0:N-1]*2*pi/N; Xk=1+exp(-j*omegam)+exp(-2*j*omegam); stem(omegam./pi,abs(Xk),r,o);hold off x1=real(ifft(Xk)); figure(2);stem(x1,r) 这里用到ifft函数，实现由频域抽样点计算其对应时域序列。可以得到时域序列为 ，产生混叠导致失真。运行结果如图1.2.5所示，可以得到由频域抽样点恢复时域信号 图1.2.5 频域抽样点 N=2恢复时域信号 3．实验设备及材料 装有Matlab的计算机 一台 4．实验方法步骤及注意事项 利用Matlab中的函数分析并绘出常用基本信号的波形。 注意事项： 在使用MATLAB时应注意中英输入法的切换，在中文输入法输入程序时得到的程序是错误的； MATLAB中两个信号相乘有好多批量乘,即.*； 使用MATLAB编写程序时，应新建一个.m文件。