信号采样与重建的编程实现

课程设计任务书学生姓名：陈凯鑫专业班级：电信1203班指导教师：阙大顺，王虹工作单位：信息工程学院题目:信号采集与重建的编程实现初始条件：1.Matlab6.5以上版本软件；2.课程设计辅导资料：“Matlab语言基础及使用入门”、“数字信号处理原理与实现”、“Matlab及在电子信息课程中的应用”等；3.先修课程：信号与系统、数字信号处理、Matlab应用实践及信号处理类课程等。要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1.课程设计时间：1周（课内实践）；2.课程设计内容：信号采样与重建的编程实现，具体包括：连续信号的时域采样、频谱混叠分析、由离散序列恢复模拟信号等；3.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结；4.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括：①目录；②与设计题目相关的理论分析、归纳和总结；③与设计内容相关的原理分析、建模、推导、可行性分析；④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结；⑤课程设计的心得体会（至少500字）；⑥参考文献；⑦其它必要内容等。时间安排：1）第1-2天，查阅相关资料，学习设计原理。2）第3-4天，方案选择和电路设计仿真。3）第4-5天，电路调试和设计说明书撰写。4）第6天，上交课程设计成果及报告，同时进行答辩。指导教师签名：年月日系主任（或责任教师）签名：年月日摘要数字信号处理是一门理论与实践紧密结合的课程。做大量的习题和上机实验，有助于进一步理解和巩固理论知识，还有助于提高分析和解决实际问题的能力。过去用其他算法语言，实验程序复杂，在有限的实验课时内所做的实验内容少。MATLAB强大的运算和图形显示功能，可使数字信号处理上机实验效率大大提高。特别是它的频谱分析和滤波器分析与设计功能很强，使数字信号处理工作变得十分简单、直观。本实验设计的题目是：信号的采样与恢复。通过产生一个连续时间信号并生成其频谱，然后对该连续信号抽样，并对采样后的频谱进行分析，实验中，原连续信号的频谱由于无法实现真正的连续，所以通过扩大采样点的数目来代替，理论上当采样点数无穷多的时候即可实现连续，基于此尽可能增加采样点数并以此来产生连续信号的频谱。信号采样过程中，通过采样点的不同控制采样频率实现大于或小于二倍最高连续信号的频率，从而可以很好的验证采样定理。关键词：信号采样信号重建MATLAB编程目录1.软件介绍...............................................12.课程设计的方案.........................................22.1课程设计的原理.....................................22.1.1连续信号的采样定理............................22.1.2信号采样......................................32.1.3信号重建......................................52.2设计的思路........................................52.3设计方案优缺点.....................................53.信号采样程序及仿真.....................................63.1连续信号x（t）及其抽样函数x（n）..................63.2采样程序及其波形..................................83.3200Hz幅频特性程序及波形..........................93.4400Hz幅频特性程序及波形..........................103.51000Hz幅频特性程序及波形.........................113.6波形分析.........................................114.信号的重建程序及仿真..................................124.1采样信号的重建程序与波形..........................124.2误差分析.........................................155心得体会..............................................166参考文献..............................................17武汉理工大学《信号分析与处理》课程设计说明书11.软件介绍MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩形计算、视化以线性动态线性系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多领域一面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB是矩阵实验室（MatrixLaboratory）之意。除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。经过不断完善MATLAB已经发展成为适合多学科，多种工作平台的功能强大大大型软件。成为线性代数,自动控制理论，数理统计，数字信号处理，时间序列分析，动态系统仿真等高级课程的基本教学工具。MTLAB的语言特点：(1)语言简洁紧凑，使用方便灵活，库函数极其丰富。(2)运算符丰富。(3)MATLAB既具有结构化的控制语句（如for循环，while循环，break语句和if语句），又有面向对象编程的特性。(4)程序限制不严格，程序设计自由度大。(5)MATLAB的图形功能强大。(6)MATLAB的缺点是，它和其他高级程序相比，程序的执行速度较慢。由于MATLAB的程序不用编译等预处理，也不生成可执行文件，程序为解释执行，所以速度较慢。武汉理工大学《信号分析与处理》课程设计说明书22.课程设计的方案2.1课程设计的原理2.1.1连续信号的采样定理模拟信号经过(A/D)变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率fs，重复出现一次。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。时域采样定理从采样信号)(tfs恢复原信号)(tf必需满足两个条件:(1))(tf必须是带限信号，其频谱函数在||＞s各处为零；（对信号的要求，即只有带限信号才能适用采样定理。）(2)取样频率不能过低，必须s＞2m（或sf＞2mf）。（对取样频率的要求，即取样频率要足够大，采得的样值要足够多，才能恢复原信号。）如果采样频率Ts/2大于或等于max2，即max2s（max2为连续信号t的有限频谱），则采样离散信号t能无失真地恢复到原来的连续信号t。一个频谱在区间（-m，m）以外为零的频带有限信号)(tf，可唯一地由其在均匀间隔ST上的样点值)(nTsf所确定。根据时域与频域的对称性，可以由时域采样定理直接推出频域采样定理。一个时间受限信号tf，它集中在（mm,）的时间范围内，则该信号的频谱jF在频域中以间隔为1的冲激序列进行采样，采样后的频谱)(1jF可以惟一表示原信号的条件为重复周期mtT21。采样信号)(tfs的频谱是原信号频谱)(jF的周期性重复，它每隔重复出现一次。当s＞2m时，不会出现混叠现象，原信号的频谱的形状不会发生变化，从而能从采样信号)(tfs中恢复原信号)(tf。（注：s＞2m的含义是：采样频率大于等于信号最高频率的2倍；这里的“不混叠”意味着信号频谱没有被破坏，也就为后面恢复原信号提供了可能）武汉理工大学《信号分析与处理》课程设计说明书3(a)(b)(c)图*抽样定理a)等抽样频率时的抽样信号及频谱（不混叠）b)高抽样频率时的抽样信号及频谱（不混叠）c)低抽样频率时的抽样信号及频谱（混叠）2.1.2信号采样如图1所示，给出了信号采样原理图信号采样原理图（a）由图1可见，)()()(ttftfsTs，其中，冲激采样信号)(tsT的表达式武汉理工大学《信号分析与处理》课程设计说明书4为：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信号的采样用数学表达式描述上述调制过程，则有)()()(*tteteT理想单位脉冲序列)(tT可以表示为0)()(nTnTtt其中)(nTt是出现在时刻nTt，强度为1的单位脉冲。由于)(te的数值仅在采样瞬时才有意义，同时，假设00)(tte所以)(*te又可表示为：*0()()()netenTtnT武汉理工大学《信号分析与处理》课程设计说明书52.1.3信号重建用时域内插公式()()()anxtxngtnT（其中sin()()sin()stTgtcFttT）完成信号的重建。2.2设计的思路连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值与之对应。严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。时域对连续时间信号进行采样，是给它乘以一个采样脉冲序列，就可以得到采样点上的样本值，信号被采样前后在频域的变化，可以通过时域频域的对应关系分别求得了采样信号的频谱。在一定条件下，一个连续时间信号完全可以用该信号在等时间间隔上的瞬时值来表示，并且可以用这些样本值把信号完全恢复过来。这样，抽样定理为连续时间信号与离散时间信号的相互转换提供了理论依据。通过观察采样信号的频谱，发现它只是原信号频谱的线性重复搬移，只要给它乘以一个门函数，就可以在频域恢复原信号的频谱，在时域是否也能恢复原信号时，利用频域时域的对称关系，得到了信号。2.3设计方案优缺点优点：MATLAB在绘图方面提供了相当高级的函数序及程序界面，即使用户没有丰富的程序设计经验，也能够快速地得到自己想要的结果，熟练的使用MATLAB的程序员或研究人员能缩短研究开发时间，从而提高竞争力，MATLAB和其他高级语言有良好的接口，可以方便地实现与其他语言的混合编程，从而进一步扩宽MATLAB的应用潜力。武汉理工大学《信号分析与处理》课程设计说明书6缺点：MATLAB占用内存空间很大，并且会因硬盘分区是NTFS格式还是FAT格式而有差异。3.信号采样程序及仿真3.1连续信号x（t）及其抽样函数x（n）cl