Matlab环境下采样定理的验证Word下载.docx

Matlab环境下采样定理的验证Word下载.docx

《Matlab环境下采样定理的验证Word下载.docx》由会员分享，可在线阅读，更多相关《Matlab环境下采样定理的验证Word下载.docx（15页珍藏版）》请在冰豆网上搜索。

计算机与信息工程学院

2014年1月3日

课程设计任务书

2012—2013学年第1学期

计算机与信息工程学院电子信息工程专业11电信1班班级

课程设计名称：

数字信号处理

设计题目：

Matlab环境下采样定理的验证

完成期限：

自2014年12月23日至2014年1月3日共2周

设计依据、要求及主要内容：

一．课程设计依据

时域采样定理和频域采样定理是数字信号处理中的重要理论，在掌握采样定理内容及原理的基础上，编写Matlab程序验证采样定理。

二．课程设计内容

1.连续信号画出连续信号的时域波形及频谱特性曲线

2.对信号进行采样得到采样序列，画出采样频率分别是200Hz，100Hz，60Hz时的采样序列波形；

3.对不同采样频率下的采样序列进行频谱分析，绘制出幅频曲线，对比各频率下采样序列的幅频曲线有无区别；

4.由采样序列恢复出连续信号，画时域波形，对比原连续时间信号波形；

5.信号，编写程序分别对=FT[x（n）]在上等间隔采样32点和16点，得到，再分别对进行32点和16点IFFT得到，分别画出，的频谱图，并画出x（n）,的波形，进行对比。

三．课程设计要求

1.要求独立完成设计任务。

2.课程设计说明书封面格式要求见《天津城市建设学院课程设计教学工作规范》附表1

3.课程设计的说明书要求简洁、通顺，计算正确，图纸表达内容完整、清楚、规范。

4.运行程序，观察并保存程序运行结果，能够对运行结果进行结果分析。

5.课设说明书要求：

1）说明题目的设计原理和思路、采用方法及程序。

2）详细说明调试方法和调试过程，并给程序加注释。

3）给出程序运行结果，并对其进行说明和分析。

指导教师（签字）：

系/教研室主任（签字）：

批准日期：

2013年12月19日

目录

第1章数字信号处理介绍

背景知识

数字信号处理是20世纪60年代，随着信息学科和计算机学科的高速发展而迅速发展起来的一门新兴学科。

它的重要性日益在各个领域的应用中表现出来。

数字信号处理是把信号用数字或符号表示成序列，通过计算机或通用（专用）信号处理设备，用数字的数值计算方法处理（例如：

滤波、变换、增强、估计、识别等），达到提取有用信息便于应用的目的。

数字信号处理（DigitalSignalProcessing,简称DSP）是一门设计许多学科而又广泛应用于许多领域的新兴学科。

DSP有两种含义：

DigitalSignalProcessing（数字信号处理）、DigitalSignalProcessor（数字信号处理器）。

我们常说的DSP指的是数字信号处理器。

数字信号处理器是一种适合完成数字信号处理运算的处理器。

广义来说，数字信号处理是研究用数字方法对信号进行分析、变换、滤波、检测、、以及快速算法的一门技术学科。

但很多人认为：

数字信号处理主要是研究有关数字滤波技术、离散变换快速算法和谱分析方法。

随着技术以及计算机技术的发展，数字信号处理技术也相应地得到发展，其应用领域十分广泛。

MATLAB软件介绍

MATLAB是矩阵实验室（MatrixLaboratory）之意。

除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。

经过不断完善MATLAB已经发展成为适合多学科，多种工作平台的功能强大大大型软件。

成为线性代数,自动控制理论，数理统计，数字信号处理，时间序列分析，动态系统仿真等高级课程的基本教学工具。

MATLAB的语言特点：

（1）语言简洁紧凑，使用方便灵活，库函数极其丰富。

（2）运算符丰富。

（3）MATLAB既具有结构化的控制语句（如for循环，while循环，break语句和if语句），又有面向对象编程的特性。

（4）程序限制不严格，程序设计自由度大。

（5）MATLAB的图形功能强大。

（6）MATLAB的缺点是，它和其他高级程序相比，程序的执行速度较慢。

由于MATLAB的程序不用编译等预处理，也不生成可执行文件，程序为解释执行，所以速度较慢。

第2章设计目的及要求

设计目的

1理解并掌握序列的采样定理的条件，并进行验证；

2学会运用MATLAB验证函数的采样定理；

3学会运用MATLAB画出函数的幅频曲线；

设计要求

1掌握利用MATLAB在数字信号处理中的基本应用，并会对结果用所学知识进行分析；

2对连续信号进行采样，在满足采样定理和不满足采用定理两种情况下对连续信号和采样信号进行FFT频谱分析；

3从采样信号中恢复原信号，对不同采样频率下的恢复信号进行比较分析。

第3章设计原理

整体设计原理

利用Matlab对连续信号进行采样，对比采样频率，，三种情况下的频谱，并与理想频谱进行分析，进而从实验上验证时域采样定理。

由时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：

“时域采样频谱周期延拓，频域采样时域信号周期延拓”。

因此放在一起进行实验。

时域采样定理

模拟信号经过（A/D）变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率fs，重复出现一次。

为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。

（1）对模拟信号以间隔T进行时域等间隔理想采样，形成的采样信号的频谱是原模拟信号频谱以采样角频率（）为周期进行周期延拓。

公式为：

（2）采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。

理想采样信号和模拟信号之间的关系为

对上式进行傅立叶变换，得到：

在上式的积分号内只有当时，才有非零值，因此：

采样信号的恢复

可用传输函数的理想低通滤波器不失真地将原模拟信号恢复出来，只是一种理想恢复。

因为

理想低通滤波器的输入输出和，

=*=

采样信号的混叠现象

对连续信号进行等间隔采样时，如果不能满足采样定理，采样后信号的频率就会重叠，即高于采样频率一半的频率成分将被重建成低于采样频率一半的信号。

这种频谱的重叠导致的失真称为，而重建出来的信号称为原信号的混叠替身，因为这两个信号有同样的样本值。

频域采样定理

频域采样点数N必须大于等于时域离散信号的长度M（即N≥M），才能使时域不产生混叠，则N点IDFT[]得到的序列就是原序列x（n）,即=x（n）。

如果N>

M，比原序列尾部多N-M个零点；

如果N<

M，z则=IDFT[]发生了时域混叠失真，而且的长度N也比x（n）的长度M短，因此。

与x（n）不相同。

频域采样定理公式的推导

对信号x（n）的频谱函数X（ejω）在[0，2π]上等间隔采样N点，得到

则N点IDFT[]得到的序列就是原序列x（n）以N为周期进行周期延拓后的主值区序列，公式为：

第4章设计过程的实现

设计步骤

1.连续信号画出连续信号的时2.域波形及频谱特性曲线；

2对信号进行采样得到采样序列，画出采样频率分别是200Hz，100Hz，60Hz时的采样序列波形；

设计步骤框图

绘制时域波形及幅频特性

采样

200Hz

100Hz

60Hz

频谱分析

恢复信号计算误差

图1采样定理示意

第5章设计结果及分析

用MATLAB设计程序对应图形

图2时域波形与幅频特性曲线

图3信号恢复后的时域波形图

图4频域采样的频谱图

设计分析

时域采样定理的分析

由公式知原信号的最高截止频率为100Hz,由上图可知，在三种不同的采样频率下，只有采样频率为200Hz时信号恢复较好。

具体分析如下：

（1）频率为60Hz即<

时，为原信号的欠采样信号和恢复，采样频率不满足时域采样定理，那么频移后的各相临频谱会发生相互重叠，这样就无法将他们分开，因而也不能再恢复原信号。

频谱重叠的现象被称为混叠现象。

（2）频率为100Hz即=时，为原信号的临界采样信号和恢复，上图有其采样的离散波形和频谱，从恢复后信号和原信号先对比可知，恢复情况很差，前后时域波形差别很大。

（3）频率为200Hz即>

时，为原信号的过采样信号和恢复，由上图采样信号离散波形和频谱，可以看出采样信号的频谱是原信号频谱进行周期延拓形成的。

此时的采样是成功的，它能够恢复原信号，从时域波形可看出，从上图采样恢复后的波形可看出与原信号误差很小了，并且比上面采样所得的冲激脉冲串包含的细节要多，说明恢复信号的精度已经很高，这样我们可以利用低通滤波器得到无失真的重建。

频域采样定理的分析

由频域采样定理可知，采样后的信号频谱是原信号频谱以采样频率为周期进行周期延拓形成的，周期性在上面两个图中都有很好的体现。

但是从16点和32点采样后的结果以及与员连续信号频谱对比可以看出，16点对应的频谱出现了频谱混叠而并非原信号频谱的周期延拓。

这是因为N取值过小导致采样角频率，因此经周期延拓出现了频谱混叠。

而N取32时，其采样角频率，从而可以实现原信号频谱以抽样频率为周期进行周期延拓，并不产生混叠.

本例中，x（n）的长度M=26。

从图中可以看出，当采样点数N=32大于等于序列长度M=26时，则x32（n）=IDFT[X（k）]=x（n），时域无混叠；

当采样点数N小于序列长度M时，对信号x（n）的频谱函数X（ejω）在[0，2π]上等间隔采样N=16时，N点IDFT[]得到的序列正是原序列x（n）以16为周期进行周期延拓后的主值区序列。

由于存在时域混叠，x316（n）=IDFT[X（k）]≠x（n）。

由此可知频域采样定理的条件是：

频域采样点数N≧M。

总结

这次的数字信号处理的课设，我们做的是基于MATLAB的采样定理的验证。

采样定理在数字信号处理中是一个重要的内容。

以前的学习当中，有很多的不了解。

通过这次的课设，我清楚地掌握了采样定理的内容，学到了不少。

首先，通过本次课设，我掌握了MATLAB软件的使用，了解了其优缺点，学会了利用MATLAB实现函数的波形图及频谱图的产生。

在设计过程中，有好几次图都无法产生，通过查阅资料及网上查询，最终将所要求的图形全都画了出来。

通过这个过程，我明白了所有的软件程序调试过程中，我们只有自己亲自动手，亲自编写程序，才能学会解决问题的方法，提高自己发现问题、分析