信号与系统实验1
更新时间:2023-12-01 00:50:01 阅读量: 教育文库 文档下载
- 信号与系统推荐度:
- 相关推荐
实验篇 信号与系统实验指导
实验一、MATLAB编程基础及典型实例
一、实验目的
(1) 熟悉MATLAB软件平台的使用; (2) 熟悉MATLAB编程方法及常用语句; (3) 掌握MATLAB的可视化绘图技术;
(4) 结合《信号与系统》的特点,编程实现常用信号及其运算。
二、实验原理
连续信号是指自变量的取值范围是连续的,且对于一切自变量的取值,除了有若干个不连续点以外,信号都有确定的值与之对应。严格来说,MATLAB并不能处理连续信号,而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时,这些离散的样值就能较好地近似连续信号。
矩阵是MATLAB进行数据处理的基本单元,矩阵运算是MATLAB最重要的运算。通常意义上的数量(也称为标量)在MATLAB系统中是作为1×1的矩阵来处理的,而向量实际上是仅有一行或者一列的矩阵。
通常用向量表示信号的时间取值范围,如n = -5:5,但信号x(n)、向量n本身的下标都是从1开始的,因此必须用一个与向量x等长的定位时间变量n,以及向量x,才能完整地表示序列x(n)。这一点详情可参考预备篇示例7的程序说明。
三、实验内容与步骤
(1) 新建一个文件夹,以自己的汉语名字命名,以后就用该文件夹专门存放自己所编制的M文件和产生的图形;将该文件夹设置成当前工作目录。
2(2) 绘制信号x(t)?e?2tsin(t)的曲线,t的范围在0 ~ 30s,取样时间间隔为0.1s.
3?2n,?3?n?3(3) 在n = [-10:10] 范围产生离散序列:x(n)?? ,并绘图。
0,其余n? 1
四、实验报告要求
整理并给出“实验内容与步骤”(2)、(3)的程序代码与产生的图形;并回答下面的问题。
(1) 在调用某一函数文件时,该文件中除了输入、输出变量外的其它变量在调用
函数结束后是否还存在?这些变量是全局还是局部变量?
(2) 设n = -10:0.2:20,你可以通过哪些方法查看向量n的维数?经过关系运
算y = (n >= 3)以后,y的维数是多少?y又等于什么?
(3) 通过MATLAB的帮助系统,学习fliplr函数的功能和使用方法。再此基础上,
写出能够产生如下图形的程序,其中x(n)?2n,?4?n?4。
1614121086420-416141210x(-n)-2024x(n)86420-4-2024
2
实验二、 连续时间系统的时域分析
一、实验目的
(1) 深刻理解卷积运算,掌握连续线性卷积的计算方法;
(2) 加深对线性时不变系统中零状态响应概念的理解,掌握其求解方法; (3) 掌握给定连续系统的冲激响应和阶跃响应。
二、实验原理
(1)
线性时不变 (LTI) 连续时间系统用常系数线性微分方程进行描述,系统的零状态响应就是在系统初始状态为零条件下微分方程的解。MATLAB控制系统工具箱提供了一个lsim函数来求解连续时间系统的零状态响应。设系统方程为:
a3y(3)(t)?a2y(2)(t)?a1y(1)(t)?a0y(t)?b3f(3)(t)?b2f(2)(t)?b1f(1)(t)?b0f(t),
该方程左边、右边的系数向量分别为a?[a3,a2,a1,a0],b?[b3,b2,b1,b0],所对应的系统模型sys可借助MATLAB中的tf函数得到:sys = tf(b, a) .
这样,系统的零状态响应为:y = lsim(sys, f, t) ,其中f是输入信号向量,t是与f对应的时间变量。 (2)
连续系统的冲激响应、阶跃响应分别是输入信号为?(t)和u(t)所对应的零状态响应。MATLAB控制系统工具箱专门提供了两个函数求解连续系统的冲激响应和阶跃响应。
冲激响应:y = impulse(sys, t) ; 阶跃响应:y = step(sys, t) . 其中sys, t的含义同上。 (3)
卷积是信号与系统中一个最基本、也是最重要的概念之一。在时域中,对于LTI 连续时间系统,其零状态响应等于输入信号与系统冲激响应的卷积;而利用卷积定理,这种关系又对应频域中的乘积。如实验一所述,我们用离散卷积来代替连续卷积,只要取样时间间隔足够小时,就可得到满意的效果。
MATLAB信号处理工具箱提供了一个计算两个离散序列卷积和的函数
3
conv。设向量a、b代表待卷积的两个序列,则c = conv(a, b)就是a与b卷积后得到的新序列。
我们知道两个序列卷积以后,一般而言所得新序列的时间范围、序列长度都会发生变化。例如设f1(n)长度为5,-3≤n≤1;f2(n)长度为7,2≤n≤8;则卷积后得到的新序列长度为11,-1≤n≤9。但是用conv函数求出卷积后没有给出新序列所对应的时间变量。为此,我们在下面的程序示例中给出了一个函数文件dconv,它在完成conv函数功能的同时,还产生了一个对应新序列的时间变量。 (4)
对于连续卷积
?f(t)?f1(t)?f2(t)????f(?)?f12(t??)d??lim?f1(k?)?f2(t?k?)??,
??0k????令t?n?(n为整数),则
f(n?)?k????f(k?)?f1?2(n??k?)?????f1(k?)?f2[(n?k)?] (*)
k????由(*)式,连续卷积积分可由离散卷积和近似代替,只要取样时间间隔?足够小,就可以得到高精度卷积积分的数值计算。在示例3中给出了一个函数文件cconv来完成该功能。
三、程序示例
示例1:已知系统的微分方程为
y(2)(t)?4y(1)(t)?4y(t)?f(1)(t)?3f(t),f(t)?e?t?(t)。求零状态响应y(t)。
a = [1 4 4]; b = [1 3]; sys = tf(b, a); td = 0.01; t = 0 : td : 10; f = exp(-t);
y = lsim(sys, f, t); plot(t, y);
xlabel('t(sec)'); ylabel('y(t)'); grid on
程序运行结果见下图。
4
0.350.30.250.2y(t)0.150.10.050012345t(sec)678910
示例2:利用conv函数,编制一个函数文件dconv,其输出为两个序列卷积后的新序列以及与该新序列对应的时间变量。
function [f, k] = dconv(f1, f2, k1, k2)
% 计算 f1与f2的卷积,并返回与得到的新序列相对应的时间变量 f = conv(f1, f2); k_start = k1(1) + k2(1);
k_end = length(f1) + length(f2) - 2; k = k_start : (k_start + k_end);
示例3:在dconv函数和(*)式的基础上,编制一个函数文件cconv,利用离散卷积和来近似计算连续卷积积分。
function [f, k] = cconv(f1,f2,k1,k2,td)
% 计算 f1与f2的连续卷积,并返回与得到的新序列相对应的时间变量 f = td*conv(f1,f2); % 实验讲义中的(*)式计算 k_start = k1(1) + k2(1);
k_end = length(f1) + length(f2) - 2; k = k_start :td: (k_start + k_end*td);
5
四、实验内容与步骤
(1) 已知系统的微分方程为
y(2)(t)?2y(1)(t)?2y(t)?f(1)(t),f(t)??(t)。计算系统的零状态响应y(t)、
冲激响应?(t)和阶跃响应g(t),并画出相应的图形。
(2) 编程实现如下图所示的两个波形;并利用cconv函数计算这两个信号的卷
积、画出卷积后的波形。
f1(t)21f2(t)-1 0 1t-2 0 2t
五、实验报告要求
整理并给出“实验内容与步骤”(1)、(2)、(3)中的程序代码与产生的图形;并回答下面的问题。
(1) 在“实验内容与步骤”(1),零状态响应y(t)和阶跃响应g(t)是否相同?为
什么?
(2) 两序列进行卷积后得到新的序列,说明新序列在时域长度、时域区间上与原
来两序列的关系。
6
四、实验内容与步骤
(1) 已知系统的微分方程为
y(2)(t)?2y(1)(t)?2y(t)?f(1)(t),f(t)??(t)。计算系统的零状态响应y(t)、
冲激响应?(t)和阶跃响应g(t),并画出相应的图形。
(2) 编程实现如下图所示的两个波形;并利用cconv函数计算这两个信号的卷
积、画出卷积后的波形。
f1(t)21f2(t)-1 0 1t-2 0 2t
五、实验报告要求
整理并给出“实验内容与步骤”(1)、(2)、(3)中的程序代码与产生的图形;并回答下面的问题。
(1) 在“实验内容与步骤”(1),零状态响应y(t)和阶跃响应g(t)是否相同?为
什么?
(2) 两序列进行卷积后得到新的序列,说明新序列在时域长度、时域区间上与原
来两序列的关系。
6
正在阅读:
信号与系统实验112-01
蚁族不应该坚守大城市11-25
2019七年级数学上册第一章1.1正数和负数1.1.1正数和负数备课资料03-26
SMT设备操作五项注意事项06-02
国际金融期中考总结资料09-20
研究生英语读写 下册翻译06-08
高中地理必修一第一章单元测试题08-05
第2章 流体静力学10-06
我的写真集作文600字07-01
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 信号
- 实验
- 系统
- 安徽省建设项目涉及国家安全事项管理规定
- 小学生名言警句大全
- 泡沫灭火系统-计算实例 - 图文
- 中小学教师信息技术应用能力培训心得体会篇
- 五上易错题1 - 图文
- db2命令实例大全
- 21克拉情侣主题餐厅创业计划书 - 图文
- 焊接作业安全控制措施
- 五年级下期数学专项复习 - 分数加减法应用题
- 小学阳光体育运动实施方案
- 网格员管理制度
- 城市设计经典案例总结
- 专升本生物化学试卷(b)
- ××考勤管理专项检查自查报告(上传)
- 小米公司战略分析报告
- 重庆市2008建筑工程定额说明和计算规则
- 三年级数学上册第六单元蚂蚁做操说课稿北师大版
- 八年级英语下册-Unit-7-10知识点详解 - 图文
- 水质工程学二100道简答题
- 山东省对外贸易现状、存在问题及对策