《应用 MATLAB实现连续信号的采样与重构仿真.》由会员分享，可在线阅读，更多相关《应用 MATLAB实现连续信号的采样与重构仿真.（18页珍藏版）》请在金锄头文库上搜索。
1、课 程 设 计 报 告课程名称 信号与系统 系 别： 机电工程系 专业班级： 自动化1002班 学 号： 1009101022 姓 名： 乔 垒 垒 课程题目： LTI连续系统分析仿真 完成日期： 2013年6月10日 指导老师： 权 宏 伟 目录第一章 绪论21.1 信号与系统的背景21.2 MATLAB软件简介2第二章 连续信号的采样与重构仿真32.1、课程设计的目的32.2、课程设计的内容及要求32.3、课程设计的原理42.3.1连续信号的采样定理42.3.2信号采样52.3.3信号重构7第三章 应用MATLAB仿真93.1 MATLAB设计的思路93.2 详细设计过程93.2.1Sa(
2、t)的临界采样及重构93.2.2 Sa(t)的过采样及重构113.2.3Sa(t)的欠采样及重构132.5设计方案优缺点15第四章 收获和体会16参考文献17第一章 绪论1.1 信号与系统的背景人们之间的交流是通过消息的传播来实现的，信号则是消息的表现形式，消息是信号的具体内容。信号与系统课程是一门实用性较强、涉及面较广的专业基础课，该课程是将学生从电路分析的知识领域引入信号处理与传输领域的关键性课程,对后续专业课起着承上启下的作用. 该课的基本方法和理论大量应用于计算机信息处理的各个领域,特别是通信、数字语音处理、数字图像处理、数字信号分析等领域,应用更为广泛。近年来，计算机多媒体教序手段的
3、运用逐步普及，大量优秀的科学计算和系统仿真软件不断涌现，为我们实现计算机辅助教学和学生上机实验提供了很好的平台。通过对这些软件的分析和对比，我们选择MATLAB语言作为辅助教学工具，借助MATLAB强大的计算能力和图形表现能力，将信号与系统中的概念、方法和相应的结果，以图形的形式直观地展现给我们，大大的方便我们迅速掌握和理解老师上课教的有关信号与系统的知识。1.2 MATLAB软件简介MATLAB 是MathWork 公司于1984 年推出的一套面向工程和科学运算的高性能软件。它具有强大的矩阵计算能力和良好的图形可视化功能，为用户提供了非常直观和简洁的程序开发环境，因此被称为第四代计算机语言。
4、MATLAB 强大的图形处理功能及符号运算功能，为我们实现信号的可视化及系统分析提供了强有力的工具。MATLAB 强大的工具箱函数可以分析连续信号、连续系统，同样也可以分析离散信号、离散系统，并可以对信号进行各种分析域计算，如相加、相乘、移位、反折、傅里叶变换、拉氏变换、Z 变换等等多种计算。此次课程设计是在MATLAB软件下进行LTI连续系统的分析仿真，有助于我对该连续信号的分析和理解。MATLAB 强大的功能为此次求连续信号冲激阶跃响应、系统零输入、零状态响应，及幅频相频等各种信号求解提供很好的视觉效果，对我们有很大的学习帮助。第二章 连续信号的采样与重构仿真2.1、课程设计的目的1. 掌
5、握利用MATLAB分析系统频率响应的方法，增加对仿真软件MATLAB的感性认识，学会该软件的操作和使用方法。2. 掌握利用MATLAB实现连续信号采用与重构的方法，加深理解采样与重构的概念。3 . 初步掌握线性系统的设计方法，培养独立工作能力。4. 学习MATLAB中信号表示的基本方法及绘图函数的调用，实现对常用连续时间信号的可视化表示，加深对各种电信号的理解。5. 加深理解采样对信号的时域和频域特性的影响；验证信号与系统的基本概念、基本理论，掌握信号与系统的分析方法。6. 加深对采样定理的理解和掌握，以及对信号恢复的必要性；掌握对连续信号在时域的采样与重构的方法。2.2、课程设计的内容及要求
6、离散正弦序列的MATLAB表示与连续信号类似，只不过是用stem函数而不是用plot函数来画出序列波形。命令窗口没打开时，从“Desktop”菜单中选择“Command Window”选项可以打开它。“”符号是输入函数的提示符，在提示符后面输入数据和运行函数。退出MATLAB时，工作空间中的内容随之清除。可以将当前工作中的部分或全部变量保存在一个MAT文件中，它是一种二进制文件，扩展名为.mat。然后可在以后使用它时载入它。用MATLAB的当前目录浏览器搜索、查看、打开、查找和改变MATLAB路径和文件。在MATLAB桌面上，从“Desktop”菜单中选择“Current Directory”
7、选项，或者在命令窗口键入“file browser”，打开当前目录浏览器。使用当前目录浏览器可以完成下面的主要任务：查看和改变路径；创建、重命名、复制和删除路径和文件；打开、运行和查看文件的内容；由于函数不是严格的带限信号，其带宽可根据一定的精度要求做一近似。根据以下三种情况用MATLAB实现采样信号及重构并求出两者误差，分析三种情况下的结果。（1）的临界采样及重构：，；（2）的过采样及重构：，；（3）的欠采样及重构：，。2.3、课程设计的原理2.3.1连续信号的采样定理模拟信号经过 (A/D) 变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率 fs，重复出现
8、一次。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。时域采样定理从采样信号恢复原信号必需满足两个条件：(1) 必须是带限信号，其频谱函数在 各处为零；（对信号的要求，即只有带限信号才能适用采样定理。）(2) 取样频率不能过低，必须 2 （或 2）。（对取样频率的要求，即取样频率要足够大，采得的样值要足够多，才能恢复原信号。）如果采样频率大于或等于，即（为连续信号的有限频谱），则采样离散信号能无失真地恢复到原来的连续信号 。一个频谱在区间（- ，）以外为零的频带有限信号，可唯一地由其在均匀间隔 （ ）上的样点值所确定。根据时域与频域的对称性，可以由
9、时域采样定理直接推出频域采样定理。一个时间受限信号，它集中在（）的时间范围内，则该信号的频谱在频域中以间隔为的冲激序列进行采样，采样后的频谱可以惟一表示原信号的条件为重复周期，或频域间隔（其中）。采样信号 的频谱是原信号频谱 的周期性重复，它每隔 重复出现一次。当2 时，不会出现混叠现象，原信号的频谱的形状不会发生变化，从而能从采样信号 中恢复原信号 。（注：2 的含义是：采样频率大于等于信号最高频率的2倍；这里的“不混叠”意味着信号频谱没有被破坏，也就为后面恢复原信号提供了可能！）(a)(b) (c)图2.1 抽样定理a) 等抽样频率时的抽样信号及频谱（不混叠）b) 高抽样频率时的抽样信号及
10、频谱（不混叠）c) 低抽样频率时的抽样信号及频谱（混叠）2.3.2信号采样 如图2.2所示，给出了信号采样原理图图2.2 信号采样原理由图2.2可见，其中，冲激采样信号的表达式为：其傅立叶变换为，其中。设，分别为，的傅立叶变换，由傅立叶变换的频域卷积定理，可得若设是带限信号，带宽为， 经过采样后的频谱就是将在频率轴上搬移至处（幅度为原频谱的倍）。因此，当时，频谱不发生混叠；而当时，频谱发生混叠。一个理想采样器可以看成是一个载波为理想单位脉冲序列的幅值调制器，即理想采样器的输出信号，是连续输入信号调制在载波上的结果，如图2.3所示。图2.3信号的采样用数学表达式描述上述调制过程，则有理想单位脉冲
11、序列可以表示为其中是出现在时刻，强度为1的单位脉冲。由于的数值仅在采样瞬时才有意义，同时，假设所以又可表示为2.3.3信号重构设信号被采样后形成的采样信号为，信号的重构是指由经过内插处理后，恢复出原来信号的过程。又称为信号恢复。若设是带限信号，带宽为，经采样后的频谱为。设采样频率，则由式（9）知是以为周期的谱线。现选取一个频率特性（其中截止频率满足）的理想低通滤波器与相乘，得到的频谱即为原信号的频谱。显然，与之对应的时域表达式为而将及代入式（10）得式（11）即为用求解的表达式，是利用MATLAB实现信号重构的基本关系式，抽样函数在此起着内插函数的作用。例：设，其为：即的带宽为，为了由的采样信
12、号不失真地重构，由时域采样定理知采样间隔，取（过采样）。利用MATLAB的抽样函数来表示，有。据此可知：通过以上分析，得到如下的时域采样定理：一个带宽为wm的带限信号f(t)，可唯一地由它的均匀取样信号fs(nTs)确定，其中，取样间隔Ts/wm, 该取样间隔又称为奈奎斯特间隔。 根据时域卷积定理，求出信号重构的数学表达式为： 式中的抽样函数Sa(wct)起着内插函数的作用，信号的恢复可以视为将抽样函数进行不同时刻移位后加权求和的结果，其加权的权值为采样信号在相应时刻的定义值。利用MATLAB中的抽样函数来表示Sa(t)，有，于是，信号重构的内插公式也可表示为： 第三章 应用MATLAB仿真3
13、.1 MATLAB设计的思路连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值与之对应。严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。时域对连续时间信号进行采样，是给它乘以一个采样脉冲序列，就可以得到采样点上的样本值，信号被采样前后在频域的变化，可以通过时域频域的对应关系分别求得了采样信号的频谱。 在一定条件下，一个连续时间信号完全可以用该信号在等时间间隔上的瞬时值来表示，并且可以用这些样本值把信号完全恢复过来。这样，抽样定理为连续时间信号
14、与离散时间信号的相互转换提供了理论依据。通过观察采样信号的频谱，发现它只是原信号频谱的线性重复搬移，只要给它乘以一个门函数，就可以在频域恢复原信号的频谱，在时域是否也能恢复原信号时，利用频域时域的对称关系，得到了信号。3.2 详细设计过程3.2.1Sa(t)的临界采样及重构1实现程序代码当采样频率小于一个连续的同信号最大频率的2倍，即时，称为临界采样. 修改门信号宽度、采样周期等参数，重新运行程序，观察得到的采样信号时域和频域特性，以及重构信号与误差信号的变化。Sa(t)的临界采样及重构程序代码；wm=1;%升余弦脉冲信号带宽wc=wm; %频率Ts=pi/wm; %周期ws=2.4*pi/Ts; %理想低通截止频率n=-100:100; %定义序列的长度是201nTs=n*Ts %采样点f=sinc(nTs/pi); %抽样信号Dt=0.005;t=-20:Dt:20;fa=f*Ts*wc/pi*sinc(wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t); %信号重建t1=-20:0.5:20;f1=sinc(t1/pi);subplot(211);stem(t1,f1);xlabel(kTs);ylabel(f(kTs);title(sa(t)=sinc(t/pi)的临界采样信号
