利用MATLAB实现SA信号的抽样与重构

利用MATLAB实现Sa信号的抽样与重构仿真1目录中文摘要.................................................I1概论...................................................21.1MATLAB的介绍....................................21.2课程设计的目的和要求..............................32MATLAB实现连续信号采样与重构的理论基础................42.1连续时间信号......................................42.2信号的采样........................................42.3信号的重构........................................63MATLAB实现Sa信号的抽样与重构仿真程序分析.............83.1Sa信号的临界采样及重构...........................93.1.1程序实现及运行结果图..................................93.1.2程序分析.............................................113.2Sa信号的过抽样及重构............................113.2.1程序实现及运行结果图.................................113.2.2程序分析.............................................143.3Sa信号的欠抽样及重构............................143.3.1程序实现及运行结果图.................................143.3.2程序分析.............................................174总结..................................................18参考文献................................................19利用MATLAB实现Sa信号的抽样与重构仿真21概论1.1MATLAB的介绍MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple、MathCAD并称为四大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。MATLAB由一系列工具组成。这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面。包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器。随着MATLAB的商业化以及软件本身的不断升级，MATLAB的用户界面也越来越精致，更加接近Windows的标准界面，人机交互性更强，操作更简单。而且新版本的MATLAB提供了完整的联机查询、帮助系统，极大的方便了用户的使用。简单的编程环境提供了比较完备的调试系统，程序不必经过编译就可以直接运行，而且能够及时地报告出现的错误及进行出错原因分析。本文基于Matlab强大功能来分析Sa信号的抽样与重构,由于Sa信号易于生成,分析方便,故在许多实际应用与仿真中(如数字通信系统等)运用较为广泛。对于连续时间信号的处理,往往将它转换为相应的离散信号(或数字信号),并进行加工和处理,然后再将处利用MATLAB实现Sa信号的抽样与重构仿真3理后的离散信号转换为连续时间信号。而抽样定理为连续信号与离散时间信号的相互转换提供了理论依据。1.2课程设计的目的和要求本次课程设计应用MATLAB实现连续信号的采样与重构仿真。了解MATLAB软件，学习应用MATLAB软件的仿真技术。它主要侧重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。初步掌握线性系统的设计方法，培养独立工作能力。加深理解采样与重构的概念，掌握利用MATLAB分析系统频率响应的方法和掌握利用MATLAB实现连续信号采用与重构的方法。计算在临界采样、过采样、欠采样三种不同条件下重构信号的误差，并由此总结采样频率对信号重构误差的影响。离散正弦序列的MATLAB表示与连续信号类似，只不过是用stem函数而不是用plot函数来画出序列波形。命令窗口没打开时，从“Desktop”菜单中选择“CommandWindow”选项可以打开它。“”符号是输入函数的提示符，在提示符后面输入数据和运行函数。退出MATLAB时，工作空间中的内容随之清除。可以将当前工作中的部分或全部变量保存在一个MAT文件中，它是一种二进制文件，扩展名为.mat。然后可在以后使用它时载入它。用MATLAB的当前目录浏览器搜索、查看、打开、查找和改变MATLAB路径和文件。在MATLAB桌面上，从“Desktop”菜单中选择“CurrentDirectory”选项，或者在命令窗口键入“filebrowser”，打开当前目录浏览器。使用当前目录浏览器可以完成下面的主要任务：查看和改变路径；创建、重命名、复制和删除路径和文件；打开、运行和查看文件的内容；由于函数)(tSa不是严格的带限信号，其带宽m可根据一定的精度要求做一近似。根据以下三种情况用MATLAB实现采样信号及重构并求出两者误差，分析三种情况下的结果。（1）)(tSa的临界采样及重构:wm=1，wc=wm，Ts=Pi/wm；（2）)(tSa的过采样及重构：wm=1，wc=1.1*wm，，Ts=0.7*Pi/wm。（3）)(tSa的欠采样及重构：，wm=1，wc=wm，，Ts=1.5*Pi/wm。利用MATLAB实现Sa信号的抽样与重构仿真42MATLAB实现连续信号采样与重构的理论基础2.1连续时间信号连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值与之对应。严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。在一定条件下，一个连续时间信号完全可以用该信号在等时间间隔上的瞬时值来表示，并且可以用这些样本值把信号完全恢复过来。这样，抽样定理为连续时间信号与离散时间信号的相互转换提供了理论依据。通过观察采样信号的频谱，发现它只是原信号频谱的线性重复搬移，只要给它乘以一个门函数，就可以在频域恢复原信号的频谱，在时域是否也能恢复原信号时，利用频域时域的对称关系，得到了信号。2.2信号的采样模拟信号经过(A/D)变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率fs，重复出现一次。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。时域采样定理从采样信号)(tfs恢复原信号)(tf必需满足两个条件：(1))(tf必须是带限信号，其频谱函数在m各处为零；（对信号的要求，即只有带限信号才能适用采样定理。）(2)取样频率不能过低，必须ms2（或fmfs2）。（对取样频率的要求，即取样频率要足够大，采得的样值要足够多，才能恢复原信号。）如图1所示，给出了信号采样原理图利用MATLAB实现Sa信号的抽样与重构仿真5图2-1信号采样原理图由图2-1可见，)()()(ttftfsTs，其中，冲激采样信号)(tsT的表达式为：nsTnTtts)()(其傅立叶变换为nssn)(，其中ssT2。设)(jF，)(jFs分别为)(tf，)(tfs的傅立叶变换，由傅立叶变换的频域卷积定理，可得nssnsssnjFTnjFjF)]([1)(*)(21)(若设)(tf是带限信号，带宽为m，)(tf经过采样后的频谱)(jFs就是将)(jF在频率轴上搬移至,,,,,02nsss处（幅度为原频谱的sT1倍）。因此，当ms2时，频谱不发生混叠；而当ms2时，频谱发生混叠。一个理想采样器可以看成是一个载波为理想单位脉冲序列)(tT的幅值调制器，即理想采样器的输出信号)(*te，是连续输入信号)(te调制在载波)(tT上的结果，如图2所示。图2-2信号的采样用数学表达式描述上述调制过程，则有)()()(*tteteT理想单位脉冲序列)(tT可以表示为0)()(nTnTtt其中)(nTt是出现在时刻nTt，强度为1的单位脉冲。由于的数值仅在采样瞬时才有意义，同时，假设00)(tte)(te利用MATLAB实现Sa信号的抽样与重构仿真6所以)(*te又可表示为*0()()()netenTtnT2.3信号的重构设信号)(tf被采样后形成的采样信号为)(tfs，信号的重构是指由)(tfs经过内插处理后，恢复出原来信号)(tf的过程,又称为信号恢复。若设)(tf是带限信号，带宽为m，经采样后的频谱为)(jFs。设采样频率ms2，则由式（9）知)(jFs是以s为周期的谱线。现选取一个频率特性ccsTjH0)(（其中截止频率c满足2scm）的理想低通滤波器与)(jFs相乘，得到的频谱即为原信号的频谱)(jF。)()()(jHjFjFs与之对应的时域表达式为)(*)()(tfthtfs而nssnssnTtnTfnTttftf)()()()()()()]([)(1tSaTjHFthccs将)(th及)(tfs代入得nscscsccssnTtSanTfTtSaTtftf)]([)()(*)()(此式即为用)(snTf求解)(tf的表达式，是利用MATLAB实现信号重构的基本关系式，抽样函数)(tSac在此起着内插函数的作用。例：设tttSatfsin)()(，其)(jF为：101)(jF利用MATLAB实现Sa信号的抽样与重构仿真7即)(tf的带宽为1m，为了由)(tf的采样信号)(tfs不失真地重构)(tf，由时域采样定理知采样间隔msT，取7.0sT（过采样）。利用MATLAB的抽样函数tttSinc)sin()(来表示)(tSa，有)/()(tSinctSa。据此可知：nscscsccs