信号与系统 连续时间信号卷积运算
连续时间系统卷积的计算
实验报告实验名称：连续时间系统卷积的数值计算班级：120241姓名：彭壮学号：12021327一、实验目的：1、加深对卷积概念及原理的理解；2、掌握借助计算机计算任意信号卷积的方法。
二、实验原理：卷积积分不仅可以通过直接积分或查表的方法来求解，还可以用积分的数值计算方法来求解。
在线性系统的分析过程中，有时会遇到复杂的激励信号，或者有时只是一组测试数据或曲线，冲激响应也可能出现同样的情况。
显然，此时直接计算积分或查表都有困难，而采用近似的数值计算方法可以解决这个问题，求得卷积积分。
1、卷积的定义卷积积分可以表示为2卷积计算的几何算法卷积积分的计算从几何上可以分为四个步骤：翻转→平移→相乘→叠加。
3卷积积分的应用卷积积分是信号与系统时域分析的基本手段，主要用于求系统零状态响应，它避开了经典分析方法中求解微分方程时需要求系统初始值的问题。
设一个线性零状态系统，已知系统的单位冲激响应为h(t)，当系统的激励信号为e(t)时，系统的零状态响应为由于计算机技术的发展，通过编程的方法来计算卷积积分已经不再是冗繁的工作，并可以获得足够的精度。
因此，信号的时域卷积分析法在系统分析中得到了广泛的应用。
卷积积分的数值运算实际上可以用信号的分段求和来实现，即：如果我们只求当t = nΔt （n为正整数,nΔt 记为t ）时r(t)的值，则由上式可以得到：当Δt 足够小时，r(t )就是e(t)和h(t)卷积积分的数值近似，由上面的公式可以得到卷积数值计算的方法如下：1 将信号取值离散化，即以 Ts 为周期，对信号取值，得到一系列宽度间隔 为 Ts 的矩形脉冲原信号的离散取值点，用所得离散取值点矩形脉冲来表示原来的连续时间信号；2 将进行卷积的两个信号序列之一反转，与另一信号相乘，并求积分，所得为 t=0 时的卷积积分的值。
以 Ts 为单位左右移动反转的信号，与另一信号相乘求积分，求的t<0和t>0时卷积积分的值；3 将所得卷积积分值与对应的t 标在图上，连成一条光滑的曲线，即为所求卷积积分的曲线。
信号与系统第二章第一讲
则相应于1的k阶重根，有k项：
( A1t k 1 A2t k 2 Ak 1t Ak )e1t ( Ai t k i )e1t
i 1
k
例2-3
信 号 与 系 统
求如下所示的微分方程的齐次解。
Hale Waihona Puke d3 d2 d r (t ) 7 2 r (t ) 16 r (t ) 12r (t ) e(t ) 3 dt dt dt
等式两端各对应幂次的系数应相等，于是有：
信 号 与 系 统
特解为： 联立解得：
3B1 1 4 B1 3B2 2 2 B 2 B 3 B 0 2 3 1
统
线性时不变系统
线性的常系数微分方程
按照元件的约束特性及 系统结构的约束特性
也即：
具体系统物理模型
常系数微分方程建立
（1）元件端口的电压与电流约束关系
iR (t ) R
信 号 与 系 统
vR (t )
C
vR (t ) iR (t ) R
dvC (t ) iC (t ) C dt
vR (t ) Ri R (t )
与
时域经典法就是直接求解系统微分方程的方法。这种方 系 法的优点是直观，物理概念清楚，缺点是求解过程冗繁，应 用上也有局限性。所以在20世纪50年代以前，人们普遍喜欢 统 采用变换域分析方法(例如拉普拉斯变换法)，而较少采用时 域经典法。20世纪50年代以后，由于δ(t)函数及计算机的普 遍应用，时域卷积法得到了迅速发展，且不断成熟和完善， 已成为系统分析的重要方法之一。时域分析法是各种变换域 分析法的基础。
信 号 与 系 统
is (t )
信号与系统王明泉版本~第二章习题解答
第2章 线性时不变连续系统的时域分析2.1 学习要求（1）会建立描述系统激励与响应关系的微分方程；（2）深刻理解系统的完全响应可分解为：零输入响应与零状态响应，自由响应与强迫响应，瞬态响应与稳态响应；（3）深刻理解系统的零输入线性与零状态线性，并根据关系求解相关的响应； （4）会根据系统微分方程和初始条件求解上述几种响应； （5）深刻理解单位冲激响应的意义，并会求解；（6）深刻理解系统起始状态与初始状态的区别，会根据系统微分方程和输入判断0时刻的跳变情况； （7）理解卷积运算在信号与系统中的物理意义和运算规律，会计算信号的卷积。
； 2.2 本章重点（1）系统（电子、机械）数学模型（微分方程）的建立； （2）用时域经典法求系统的响应； （3）系统的单位冲激响应及其求解；（4）卷积的定义、性质及运算，特别是()t δ函数形式与其它信号的卷积； （5）利用零输入线性与零状态线性，求解系统的响应。
2.3 本章的知识结构2.4 本章的内容摘要2.4.1系统微分方程的建立电阻：)(1)(t v Rt i R R =电感：dtt di L t v L L )()(= )(d )(1)(0t i v Lt i L tL L +=⎰∞-ττ 电容：dtt dv C t i C C )()(= ⎰+=tt L C C t i i Ct v 0)(d )(1)(0ττ 2.4.2 系统微分方程的求解 齐次解和特解。
齐次解为满足齐次方程t n t t h e c e c e c t y 32121)(λλλ+⋅⋅⋅++=当特征根有重根时，如1λ有k 重根，则响应于1λ的重根部分将有k 项，形如t k t k t k t k h e c te c e t c e t c t y 111112211)(λλλλ++⋅⋅⋅++=--- 当特征根有一对单复根，即bi a +=2,1λ，则微分方程的齐次解bt e c bt e c t y at at h sin cos )(21+= 当特征根有一对m 重复根，即共有m 重ib a ±=2,1λ的复根，则微分方程的齐次解bt e t c bt te c bt c t y at m m at h cos cos cos )(121-+⋅⋅⋅++= bt e t d bt te d bt e d at m m at at sin sin sin 121-+⋅⋅⋅+++ 特解的函数形式与激励函数的形式有关。
《信号与系统》课程实验报告
《信号与系统》课程实验报告《信号与系统》课程实验报告一图1-1 向量表示法仿真图形2．符号运算表示法若一个连续时间信号可用一个符号表达式来表示，则可用ezplot命令来画出该信号的时域波形。
上例可用下面的命令来实现（在命令窗口中输入，每行结束按回车键）。
t=-10:0.5:10;f=sym('sin((pi/4)*t)');ezplot(f,[-16,16]);仿真图形如下：图1-2 符号运算表示法仿真图形三、实验内容利用MATLAB实现信号的时域表示。
三、实验步骤该仿真提供了7种典型连续时间信号。
用鼠标点击图0-3目录界面中的“仿真一”按钮，进入图1-3。
图1-3 “信号的时域表示”仿真界面图1-3所示的是“信号的时域表示”仿真界面。
界面的主体分为两部分：1) 两个轴组成的坐标平面（横轴是时间，纵轴是信号值）；2) 界面右侧的控制框。
控制框里主要有波形选择按钮和“返回目录”按钮，点击各波形选择按钮可选择波形，点击“返回目录”按钮可直接回到目录界面。
图1-4 峰值为8V，频率为0.5Hz，相位为180°的正弦信号图1-4所示的是正弦波的参数设置及显示界面。
在这个界面内提供了三个滑动条，改变滑块的位置，滑块上方实时显示滑块位置代表的数值，对应正弦波的三个参数：幅度、频率、相位；坐标平面内实时地显示随参数变化后的波形。
在七种信号中，除抽样函数信号外，对其它六种波形均提供了参数设置。
矩形波信号、指数函数信号、斜坡信号、阶跃信号、锯齿波信号和抽样函数信号的波形分别如图1-5～图1-10所示。
图1-5 峰值为8V，频率为1Hz，占空比为50%的矩形波信号图1-6 衰减指数为2的指数函数信号图1-7 斜率=1的斜坡信号图1-8 幅度为5V，滞后时间为5秒的阶跃信号图1-9 峰值为8V，频率为0.5Hz的锯齿波信号图1-10 抽样函数信号仿真途中，通过对滑动块的控制修改信号的幅度、频率、相位，观察波形的变化。
《信号与系统》课程讲义1-2
ii)抽样特性： (t ) f (t )dt f (0)
证明： (t ) f (t )dt ( ) f ( )d ( ) ( ) f 0 d f 0
iv)延时抽样： v)关系：
t t f t dt f (t )
1 t
-1 0 f(-t-2) 1 -3 -2 0 t 2 t
0 1
1 -1
2 3
f(-3t-2)
0
t
§1.3信号的运算
②已知f(t)定义域为[-1,4]，求f(-2t+5)的定义域 解：
i)方法一：f(t)→f(-t) [-4,1]；f(-t)→f(-t+5) [1,6]；
ii)方法二： 1 2t 5 4 6 2t 1
f (t ) f 1 ( t ) f 2 ( t )
§1.3信号的运算
7．信号相乘 ① f (t ) f1 (t ) f 2 (t )
②常用在调制解调中 8．卷积
f (t ) f1 (t ) f 2 (t )
f1 ( ) f 2 (t )d
9．相关
a
Ke at (a 0)
③特性：微积分后仍为指数信号
§1.2 信号描述分类和典型示例
2．正弦信号 ①表达式：
f (t ) K sin(t )
②参数：K振幅， 角频率， 初相位 f(t) ③特性 i)周期信号， 0 2 1 T f ii)微积分后仍为正弦信号
3 8
t
t
f(t)
t
0 ln 2 2 ln 2 3 ln 2
3
练习
信号与系统信号的时域分解与卷积积分
28
三、卷积的性质及卷积计算
(2) (t-t0 ) 是卷积的延迟器
y(t) f (t) (t t0 )=f (t t0 )
物理意义
f (t)
有用推论
(t t0 )
f (t t0 )
f (t t1) (t t2 ) f (t t1 t2 )
若：f1(t) f2 (t) y(t) 则： f1(t t1) f2(t t2) y(t t1 t2)
s 平面和z平面的对应关系
×
衰减振荡信号
j
×虚指数信号 ×
增长振荡信号
指数×衰减信号
×
直流信号
×
指数增长信号
jIm[z]
z esT rej r eT , T
× 虚指数信号
衰减振荡信号
×
×
× 指×数增长
指数衰减信号 直流 Re[z]
增长振荡信号
× 2
温故知新，上讲回顾
信号波形的翻转、展缩与平移
)
f3 (t
)]d
f1( )
f2 (t
)d
f1 (
)
f3 (t
)d
f1(t) f2 (t) f1(t) f3 (t)
物理意义：两个LTI系统并联，其总的单位冲激响应等
于各个子系统的单位冲激响应之和。也可通过交换律/
线性系统性质证明
f1 (t )
f2 (t) f3 (t)
f1(t) [ f2 (t) f3 (t)]
f1(t) f2 (t ) f3 (t) yzs (t) f1 (t) [ f2 (t) f3 (t)]
表明:两个LTI系统级联时，系统总的单位冲激响 应等于各个子系统单位冲激响应的卷积。
信号与系统实验教程只有答案
信号与系统实验教程（只有答案））（实验报告目录实验一信号与系统的时域分析 (2)三、实验内容及步骤 (2)实验二连续时间信号的频域分析 (14)三、实验内容及步骤 (14)实验三连续时间LTI系统的频域分析 (35)三、实验内容及步骤 (35)实验四通信系统仿真 (42)三、实验内容及步骤 (42)实验五连续时间LTI系统的复频域分析 (51)三、实验内容及步骤 (51)实验一信号与系统的时域分析三、实验内容及步骤实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。
实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。
并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。
实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。
Q1-1：修改程序Program1_1，将dt改为0.2，再执行该程序，保存图形，看看所得图形的效果如何？dt = 0.01时的信号波形dt = 0.2时的信号波形这两幅图形有什么区别，哪一幅图形看起来与实际信号波形更像？答：Q1-2：修改程序Program1_1，并以Q1_2为文件名存盘，产生实指数信号x(t)=e-0.5t。
要求在图形中加上网格线，并使用函数axis()控制图形的时间范围在0~2秒之间。
然后执行该程序，保存所的图形。
修改Program1_1后得到的程序Q1_2如下：信号x(t)=e-0.5t的波形图clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.2; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = exp(-0.5*t); % Generate the signalplot(t,x)grid on;axis ([0 2 0 1 ])title('Sinusoidal signal x(t)')xlabel('Time t (sec)')Q1-3：修改程序Program1_1，并以Q1_3为文件名存盘，使之能够仿真从键盘上任意输入的一个连续时间信号，并利用该程序仿真信号x(t)=e-2t。
信号的卷积
计算机与信息工程学院实验报告专业：通信工程年级/班级：2012级通信工程2013—2014学年第二学期课程名称计算机网络实验指导教师本组成员学号姓名实验地点实验时间项目名称信号的卷积实验类型一、实验目的1. 理解卷积的物理意义；2. 掌握运用计算机进行卷积运算的原理和方法；3. 熟悉卷积运算函数conv 的应用；二、实验仪器或设备一台安装MATLAB的计算机一台三、实验原理1.卷积的定义连续时间和离散时间卷积的定义分别如下所示：=[n-k]2.卷积的计算由于计算机技术的发展，通过编程的方法来计算卷积积分和卷积和已经不再是冗繁的工作，并可以获得足够的精度，因此信号的时域卷积分析法在系统分析中得到了广泛的应用。
卷积积分的数值运算可以应用信号的分段求和来实现，即：数值运算只求当t = nΔ时的信号值 f (nΔ)，则由上式可以得到：上式中实际上就是连续信号f1(t ) f 2(t )等间隔均匀抽样的离散序列f1(nΔ) f 2(nΔ)的卷积和当Δ足够小的时候 f (nΔ)就是信号卷积积分的数值近似。
因此，在利用计算机计算两信号卷积积分时，实质上是先将其转化为离散序列，再利用离散卷积和计算原理来计算。
3.卷积的应用3. 1 求解系统响应卷积是信号与系统时域分析的基本手段，主要应用于求解系统响应，已知一 LTI系统的单位冲激响应和系统激励信号则系统响应为激励与单位冲激响应的卷积。
四、实验步骤给定如下因果线性时不变系统：y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3=0.9x[n]-0.45x[n-1]+0.35x[n-2]+0.002x[n-3] (1)不用impz 函数，使用filter 命令，求出以上系统的单位冲激响应h[n]的前20个样本;clear all;N=20;num=[2.24 2.49];den=[1 -0.4];y=impz(num,den,N);stem(y);xlabel(‘时间序号’);ylabel(‘振幅’);title(‘冲激响应’);grid;(2)得到h[n]后，给定x[n],计算卷积输出y[n]；并用滤波器h[n]对输入x[n]滤波，求得y1[n];x=[1 -2 3 -4 3 2 1];%输入序列y=conv(h,x);%h 由(1)中filter 命令求出n=0:25;subplot(2,1,1);stem(n,y);xlabel(‘时间序号n’);ylabel(‘振幅’);title(‘用卷积得到的输出’);grid;x1=[x zeros(1,19)];y1=filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel(‘时间序号n’);ylabel(‘振幅’);title(‘用滤波得到的输出’);grid;年月日。
信号与系统知识点汇总总结
信号与系统知识点汇总总结一、信号与系统概念1. 信号的定义和分类2. 系统的定义和分类3. 时域和频域分析二、连续时间信号与系统1. 连续时间信号与系统的性质2. 连续时间信号的基本操作3. 连续时间系统的性质4. 连续时间系统的特性方程和驻点三、离散时间信号与系统1. 离散时间信号与系统的性质2. 离散时间信号的基本操作3. 离散时间系统的性质4. 离散时间系统的特性方程和驻点四、傅里叶分析1. 傅里叶级数2. 傅里叶变换3. 傅里叶变换的性质4. 傅里叶变换的逆变换五、拉普拉斯变换1. 拉普拉斯变换的定义2. 拉普拉斯变换定理3. 拉普拉斯变换的性质4. 拉普拉斯变换的逆变换六、Z变换1. Z变换的定义2. Z变换的性质3. Z变换与拉普拉斯变换的关系4. Z变换在离散时间系统分析中的应用七、系统的时域分析1. 系统的冲击响应2. 系统的单位脉冲响应3. 系统的阶跃响应4. 系统的时域性能指标八、系统的频域分析1. 系统的频率响应2. 系统的幅频特性3. 系统的相频特性4. 系统的频域性能指标九、信号与系统的稳定性1. 连续时间系统的稳定性2. 离散时间系统的稳定性3. 系统的相对稳定性十、线性时不变系统1. 线性系统的性质2. 时不变系统的性质3. 线性时不变系统的连续时间性能分析4. 线性时不变系统的离散时间性能分析十一、激励响应系统1. 激励响应系统的特性2. 激励响应系统的连续时间分析3. 激励响应系统的离散时间分析十二、卷积运算1. 连续时间信号的卷积运算2. 离散时间信号的卷积运算3. 卷积的性质和应用结语信号与系统是电子信息专业的重要基础课程，掌握好这门课程的知识对学生日后的学习和工作都有重要的帮助。
通过本文的知识点汇总总结，相信读者对信号与系统这门课程会有更深入的理解和掌握，希望对大家的学习有所帮助。
信号与系统中卷积计算方法探讨
ｊ
图２ 竖式乘法计算离散信号卷积
即有 （） （） ６５ ２ ，２ ２ ，｝ ＝｛，，３１ ，１５ 。用竖式乘法求离散信号卷
积运算 的优点在于 ， 对两 个较长的有限长序列求 卷积和时 ， 运算 简便 、 易行且不易 出错 ， 缺点是只对有 限长序列有效 。 ２ 用 图表法计算离散信号卷积和 ． ６ 用 图表法计算离 散信号卷积 的具体步骤 如下 ： １ ： 表 格 ， 第 步 列 如 图 ３ 所示 。第 ２ ： ａ 步 将表格 中每一行与每一列的两组数字相乘 ， 并将 乘积 填写在相 应的空格 中 ， 图３ 所示 。第 ３ ： 如 ｂ 步 在填写乘积的表格 中 作一组 自 下方指向右上方的平行 斜线 ， 左 如图 ３ 所示 。第 ４ ： ｃ 步 将每组 平行斜线 内的数字相加 ， 得到 的 白上而下的一组数 字就是离散卷积在
延 时性 ； ，是卷积 的单 位元 ； ，是 ｔ秒的延时器 ； （ 是 微分器 ； （ ） （ ） ｏ ’） ｕｔ是积分器。 （ ）
２ 算子计算 （ ． ４ 微分冲激法）
定 ： ＝ （，。（ 去 （，称Ｐ 分 为 分 义 ｆ ）ｐｔ Ｊ厂）＝ 厂） ｆ）一 ｒ 。 ｆ则 和 别 微
３
一
１
０
１
２
ｔ
（ ａ ）
（ ｂ ）
ｋ ， ，， ，， 的结果。即 （） （ ：｛， ，３１ ，１５ ＝Ｏ １２ ３４ ５ ６ ５２ ，２ ２ ，｝ ） ２ Nhomakorabea＇ ４ １
ｆ＾ ＾ 一 ． ＾ ＾ 一 １∽ ０ ｆ ∽ ０
。 囡 ．
信号与系统【 卷积计算方法舔讨 ｌ 】
衡 阳师范学院物理与电子信息科学系 唐建锋 杨 辉 罗湘南
信号与系统-连续信号和离散信号的表示与卷积实验报告
实验一：连续信号和离散信号的表示与卷积一.实验目的1. 学习MATLAB 软件产生信号和实现信号的可视化2. 学习和掌握连续和离散信号的时域表示方法3. 学习和掌握连续信号和离散信号卷积方法二．实验原理1. 信号的表示方法● 常用信号：➢ 连续函数()θω+=t t f sin )(, atAe t f =)(，ttt Sa sin )(=➢ 离散信号()n n f 0sin ][ω=，njw e n f 0][=，][][n u a n f n=● 奇异信号：➢ 连续函数：冲激函数)(t δ，阶跃函数)(t u ，斜坡函数)(t R ➢ 离散信号：冲激函数][n δ，阶跃函数][n u ，斜坡函数][n R2．卷积连续函数的卷积：⎰∞∞--=τττd t f f t g )()()(21离散函数的卷积：∑∞-∞=-=m m n f m f n g ][][][21三．实验内容1. 熟悉matlab 工作环境（1） 运行matlab.exe ，进入matlab 工作环境，如图（1）所示。
图1 matlab工作环境（2）matlab工作环境由Command Window（命令窗口）、Current Direcroty（当前目录）、workspace （工作空间）、command History（历史命令）和Editor（文件编辑器）5部分组成。
其中所有文件的编辑和调试、运行在Editor编辑窗口下进行。
程序的运行也可以在命令窗口进行。
程序调试的信息显示在命令窗口。
（3）程序文件的产生：点击菜单file下的New下的M_files，进入编辑器界面，如图2。
图2 M 文件编辑器（4） 在m 文件编辑器下键入程序代码，保存程序文件（文件命名规则同C 语言）。
如果所定义的是函数文件，则要求函数名为M 文件名。
（5） 程序运行需要给定义的函数参数赋值。
切换到命令窗口下运行例如指数函数定义格式 [t,y]=exp1_exp(t1,t2,dt,A,a)指数函数文件调用方式：[t,y]=exp1_exp(-10,10,0.1,3,-1,1)2 连续和离散信号的时域表示方法(1)单边指数信号 )()(t u Ae t y tα=；function y=exp1_exp(t1,t2,dt,A,a,options)%指数函数，其中t1,t2,dt 分别为起始时间、终止时间和时间间隔 %A,a 为常数 y(t)=Aexp(a*t)%options 参数等于1时为单边指数函数，其他时为双边指数函数 %函数调用的格式 y=exp1_exp(-10,10,0.1,3,-1,1) if options==1t=0:dt:t2;%单边指数函数时间范围 elset=t1:dt:t2;%双边指数函数时间范围endy=A*exp(a*t);%指数函数plot(t,y)%画图grid onxlabel('t')%X轴坐标ylabel('y(t)')%Y轴坐标if options==1title(' 单边指数信号')%标题elsetitle(' 双边指数信号')%标题end实验要求：1）在同一张图上画出a>0,a=0,a<0时指数函数波形，如图3所示. 注意：a的取值范围要适中，不要导致纵坐标相差太大。
信号与系统教案第2章
不难求得其齐次解为Czs1e-t + Czs2e-2t + 3
代入初始值求得
yzs(t)= – 4e-t + e-2t + 3 ，t≥0
第2-11页
■
©江西科技师范大学通信与电子学院
信号与系统 电子教案
2.2 冲激响应和阶跃响应
2.2 冲激响应和阶跃响应
一、冲激响应
由单位冲激函数δ(t)所引起的零状态响应称为单位冲 激响应，简称冲激响应，记为h(t)。h(t)=T[{0},δ(t)]
例1 描述某系统的微分方程为 y”(t)+5y’(t)+6y(t)=f(t) 求其冲激响应h(t)。
解 根据h(t)的定义 有 h”(t) + 5h’(t) + 6h(t) = δ(t) h’(0-) = h(0-) = 0
先求h’(0+)和h(0+)。
第2-12页
■
©江西科技师范大学通信与电子学院
代入初始条件求得C1=1,C2=-1, 所以 h(t)=( e-2t - e-3t)ε(t)
第2-13页
■
©江西科技师范大学通信与电子学院
信号与系统 电子教案
2.2 冲激响应和阶跃响应
例2 描述某系统的微分方程为
y”(t)+5y’(t)+6y(t)= f”(t) + 2f’(t) + 3f(t) 求其冲激响应h(t)。
齐次解的函数形式仅与系统本身的特性有关，而与激励 f(t)的函数形式无关，称为系统的固有响应或自由响应； 特解的函数形式由激励确定，称为强迫响应。
例 描述某系统的微分方程为
y”(t) + 5y’(t) + 6y(t) = f(t) 求（1）当f(t) = 2e-t，t≥0；y(0)=2，y’(0)= -1时的全解；
《信号与系统》第二版第二章：LTI连续时间系统的时域分析
零状态（zero state）响应 yzs (t ) ：不考虑起始时刻系统储能的作用，即Y(0-) ≡0，由系统的外加激励信号 v (t ) = v (t )u (t ) ≠ 0 所产生的响应。
零输入响应 yzi (t ) ：
5
《信号与系统》
第二章：LTI 连续时间系统的时域分析
∏(p −αi )
i =1
（αi 为互异特征根）
= N (p) ⎡⎣eαnt ∗ ∗ eα1t ∗ v (t )⎤⎦
（2-19）
n
∑ yzs (t ) = 齐次解 Aieαit ＋特解 B (t ) i =1
（2-20）
特解 B (t ) 反映系统输入对输出的强迫。
非零状态线性系统： 定义（非零状态线性系统）：系统 T 的初始状态为X(0-)≠0
令： D (p) pn + an−1pn−1 + ... + a1p + a0
N (p) bmpm + ... + b1p + b0
4
《信号与系统》
有：
第二章：LTI 连续时间系统的时域分析
y
(t)
=
N (p) D(p)
v(t
)
H (p)v(t)
（2-13）
其中，
H
(p)
=
N (p) D(p)
称为系统算子。
≤ ∫ ∫ f (τ ) g (t −τ ) dτ dt ΩΩ
= ∫ f (τ ) ∫ g (t −τ ) dtdτ
Ω
Ω
=∫
f (τ )
g (t ) dτ = 1
f (t) 1
g (t ) 1
信号与系统实验2：信号与卷积
Matlab 的函数 conv 函数计算卷积积分 x(t) h(t) ，并画出图形。
(2) 画出函数 f1(t) eatu(t) 和 f2 (t) sin tu(t) 的图形，并利用 Matlab 的函数 conv 函数计算卷积 积分 f1(t) f2 (t) ，并画出图形。
(3) 画出教材 P131 例 3-45 中 x[k] 1, 2,3, 4;k 0,1, 2,3, y[k] 1,1,1,1,1;k 0,1, 2,3, 4 的图形，
4.3 step 函数： 计算并画出系统阶跃响应曲线 调用格式：该函数与函数 impulse() 一样，也有相似的调用格式。
三、实验内容及步骤 对书中的例题进行仿真：
1、连续时间系统的冲激响应、阶跃响应
(1)
利用 impulse 函数画出教材 P127 例 3-42:
LTI
系统
d
2
yzs dt
(t
)
stem(f1,'fill'); title('f1=[1 2 3 4];'); xlabel('t'); ylabel('f1(t)'); subplot(1,3,2); stem(f2,'fill') title('f2=[1 1 1 1]'); xlabel('t'); ylabel('f2(t)'); subplot(1,3,3); stem(f,'fill') title('卷积结果'); xlabel('t(sec)'); ylabel('f1(t)');
调用格式： impulse(b,a) 该调用格式以默认方式绘出向量 a 和 b 定义的连续系统的冲激响应的时域波形。
信号与系统实验报告
实验三常见信号的MATLAB 表示及运算一、实验目的1．熟悉常见信号的意义、特性及波形2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法二、实验原理根据MATLAB 的数值计算功能和符号运算功能,在MATLAB 中,信号有两种表示方法,一种是用向量来表示,另一种则是用符号运算的方法;在采用适当的MATLAB 语句表示出信号后,就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了;1.连续时间信号从严格意义上讲,MATLAB 并不能处理连续信号;在MATLAB 中,是用连续信号在等时间间隔点上的样值来近似表示的,当取样时间间隔足够小时,这些离散的样值就能较好地近似出连续信号;在MATLAB 中连续信号可用向量或符号运算功能来表示; ⑴ 向量表示法对于连续时间信号()f t ,可以用两个行向量f 和t 来表示,其中向量t 是用形如12::t t p t =的命令定义的时间范围向量,其中,1t 为信号起始时间,2t 为终止时间,p 为时间间隔;向量f 为连续信号()f t 在向量t 所定义的时间点上的样值; ⑵ 符号运算表示法如果一个信号或函数可以用符号表达式来表示,那么我们就可以用前面介绍的符号函数专用绘图命令ezplot 等函数来绘出信号的波形; ⑶ 常见信号的MATLAB 表示 单位阶跃信号单位阶跃信号的定义为：10()0t u t t >⎧=⎨<⎩方法一： 调用Heavisidet 函数首先定义函数Heavisidet 的m 函数文件,该文件名应与函数名同名即;%定义函数文件,函数名为Heaviside,输入变量为x,输出变量为y function y= Heavisidety=t>0; %定义函数体,即函数所执行指令%此处定义t>0时y=1,t<=0时y=0,注意与实际的阶跃信号定义的区别;方法二：数值计算法在MATLAB 中,有一个专门用于表示单位阶跃信号的函数,即stepfun 函数,它是用数值计算法表示的单位阶跃函数()u t ;其调用格式为：stepfunt,t0其中,t 是以向量形式表示的变量,t0表示信号发生突变的时刻,在t0以前,函数值小于零,t0以后函数值大于零;有趣的是它同时还可以表示单位阶跃序列()u k ,这只要将自变量以及取样间隔设定为整数即可; 符号函数符号函数的定义为：10sgn()1t t t >⎧=⎨-<⎩在MATLAB 中有专门用于表示符号函数的函数sign ,由于单位阶跃信号 t 和符号函数两者之间存在以下关系：1122()sgn()t t ε=+,因此,利用这个函数就可以很容易地生成单位阶跃信号;2.离散时间信号离散时间信号又叫离散时间序列,一般用()f k 表示,其中变量k 为整数,代表离散的采样时间点采样次数;在MATLAB 中,离散信号的表示方法与连续信号不同,它无法用符号运算法来表示,而只能采用数值计算法表示,由于MATLAB 中元素的个数是有限的,因此,MATLAB 无法表示无限序列；另外,在绘制离散信号时必须使用专门绘制离散数据的命令,即stem 函数,而不能用plot 函数; 单位序列()k δ单位序列()k δ的定义为10()0k k k δ=⎧=⎨≠⎩单位阶跃序列()u k单位阶跃序列()u k 的定义为10()0k u k k ≥⎧=⎨<⎩3.卷积积分两个信号的卷积定义为：MATLAB 中是利用conv 函数来实现卷积的;功能：实现两个函数1()f t 和2()f t 的卷积;格式：g=convf1,f2说明：f1=f 1t,f2=f 2t 表示两个函数,g=gt 表示两个函数的卷积结果;三、实验内容1.分别用MATLAB 的向量表示法和符号运算功能,表示并绘出下列连续时间信号的波形： ⑴ 2()(2)()tf t e u t -=- ⑵[]()cos()()(4)2tf t u t u t π=--1 t=-1::10;t1=-1::; t2=0::10;f1=zeros1,lengtht1,ones1,lengtht2;f=2-exp-2t.f1; plott,faxis-1,10,0, syms t;f=sym'2-exp-2theavisidet'; ezplotf,-1,10;2t=-2::8;f=0.t<0+cospit/2.t>0&t<4+0.t>4; plott,f syms t;f=sym'cospit/2heavisidet-heavisidet-4 '; ezplotf,-2,8;2.分别用MATLAB 表示并绘出下列离散时间信号的波形： ⑵ []()()(8)f t k u k u k =-- ⑶()sin()()4k f k u k π= 2 t=0:8; t1=-10:15;f=zeros1,10,t,zeros1,7; stemt1,faxis-10,15,0,10; 3 t=0:50; t1=-10:50;f=zeros1,10,sintpi/4; stemt1,faxis-10,50,-2,23.已知两信号1()(1)()f t u t u t =+-,2()()(1)f t u t u t =--,求卷积积分12()()()g t f t f t =*,并与例题比较;t1=-1::0; t2=0::1; t3=-1::1;f1=onessizet1; f2=onessizet2; g=convf1,f2;subplot3,1,1,plott1,f1; subplot3,1,2,plott2,f2; subplot3,1,3,plott3,g;与例题相比较,gt 的定义域不同,最大值对应的横坐标也不同;4.已知{}{}12()1,1,1,2,()1,2,3,4,5f k f k ==,求两序列的卷积和 ;N=4; M=5; L=N+M-1; f1=1,1,1,2;f2=1,2,3,4,5; g=convf1,f2; kf1=0:N-1; kf2=0:M-1; kg=0:L-1;subplot1,3,1,stemkf1,f1,'k';xlabel'k'; ylabel'f1k';grid onsubplot1,3,2,stemkf2,f2,'k';xlabel'k'; ylabel'f2k';grid onsubplot1,3,3;stemkg,g,'k';xlabel'k'; ylabel'gk';grid on 实验心得：第一次接触Mutlab 这个绘图软件,觉得挺新奇的,同时 ,由于之前不太学信号与系统遇到一些不懂的问题,结合这些图对信号与系统有更好的了解;实验四 连续时间信号的频域分析一、实验目的1．熟悉傅里叶变换的性质 2．熟悉常见信号的傅里叶变换3．了解傅里叶变换的MATLAB 实现方法二、实验原理从已知信号()f t 求出相应的频谱函数()F j ω的数学表示为：()F j ω()j t f t e dt ω∞--∞=⎰傅里叶反变换的定义为：1()()2j t f t F j e d ωωωπ∞-∞=⎰在MATLAB 中实现傅里叶变换的方法有两种,一种是利用MATLAB 中的Symbolic Math Toolbox 提供的专用函数直接求解函数的傅里叶变换和傅里叶反变换,另一种是傅里叶变换的数值计算实现法;1.直接调用专用函数法①在MATLAB 中实现傅里叶变换的函数为：F=fourier f 对ft 进行傅里叶变换,其结果为Fw F =fourierf,v 对ft 进行傅里叶变换,其结果为Fv F=fourier f,u,v 对fu 进行傅里叶变换,其结果为Fv ②傅里叶反变换f=ifourier F 对Fw 进行傅里叶反变换,其结果为fx f=ifourierF,U 对Fw 进行傅里叶反变换,其结果为fu f=ifourier F,v,u 对Fv 进行傅里叶反变换,其结果为fu 注意：1在调用函数fourier 及ifourier 之前,要用syms 命令对所有需要用到的变量如t,u,v,w 等进行说明,即要将这些变量说明成符号变量;对fourier 中的f 及ifourier 中的F 也要用符号定义符sym 将其说明为符号表达式;2采用fourier 及fourier 得到的返回函数,仍然为符号表达式;在对其作图时要用ezplot 函数,而不能用plot 函数;3fourier 及fourier 函数的应用有很多局限性,如果在返回函数中含有δω等函数,则ezplot 函数也无法作出图来;另外,在用fourier 函数对某些信号进行变换时,其返回函数如果包含一些不能直接表达的式子,则此时当然也就无法作图了;这是fourier 函数的一个局限;另一个局限是在很多场合,尽管原时间信号ft 是连续的,但却不能表示成符号表达式,此时只能应用下面介绍的数值计算法来进行傅氏变换了,当然,大多数情况下,用数值计算法所求的频谱函数只是一种近似值;2、傅里叶变换的数值计算实现法严格说来,如果不使用symbolic 工具箱,是不能分析连续时间信号的;采用数值计算方法实现连续时间信号的傅里叶变换,实质上只是借助于MATLAB 的强大数值计算功能,特别是其强大的矩阵运算能力而进行的一种近似计算;傅里叶变换的数值计算实现法的原理如下： 对于连续时间信号ft,其傅里叶变换为：其中τ为取样间隔,如果ft 是时限信号,或者当|t|大于某个给定值时,ft 的值已经衰减得很厉害,可以近似地看成是时限信号,则上式中的n 取值就是有限的,假定为N,有： 若对频率变量ω进行取样,得： 通常取：02k k k MM ωπωτ==,其中0ω是要取的频率范围,或信号的频带宽度;采用MATLAB 实现上式时,其要点是要生成ft 的N 个样本值()f n τ的向量,以及向量k j n eωτ-,两向量的内积即两矩阵的乘积,结果即完成上式的傅里叶变换的数值计算;注意：时间取样间隔τ的确定,其依据是τ必须小于奈奎斯特Nyquist 取样间隔;如果ft 不是严格的带限信号,则可以根据实际计算的精度要求来确定一个适当的频率0ω为信号的带宽;三、 实验内容1.编程实现求下列信号的幅度频谱1 求出1()(21)(21)f t u t u t =+--的频谱函数F 1jω,请将它与上面门宽为2的门函数()(1)(1)f t u t u t =+--的频谱进行比较,观察两者的特点,说明两者的关系;2 三角脉冲21||||1()0||1t t f t t -≤⎧=⎨>⎩3 单边指数信号3()()tf t e t ε-=4 高斯信号23()t f t e -=1 syms t w Gt=sym'Heaviside2t+1-Heaviside2t-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi 0与()(1)(1)f t u t u t =+--的频谱比较,1()(21)(21)f t u t u t =+--的频谱函数F 1jω最大值是其的1/2; 2syms t w;Gt=sym'1+tHeavisidet+1-Heavisidet+1-tHeavisidet-Heavisidet-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi 0 3syms t w Gt=sym'exp-tHeavisidet';Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi -1 2 4syms t w Gt=sym'exp-t^2';Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; ezplotFFw,-30 30;grid; axis-30 30 -1 22.利用ifourier 函数求下列频谱函数的傅氏反变换122()16F j j ωωω=-+ 222()58()()65j j F j j j ωωωωω+-=++1syms t w Fw=sym'-i2w/16+w^2'; ft=ifourierFw,w,t; ft运行结果： ft =-exp4theaviside-t+exp-4theavisidet 2syms t wFw=sym'iw^2+5iw-8/iw^2+6iw+5'; ft=ifourierFw,w,t; ft运行结果： ft =diract+-3exp-t+2exp-5theavisidet实验心得matlab 不但具有数值计算能力,还能建模仿真,能帮助我们理解不同时间信号的频域分析;实验五 连续时间系统的频域分析一、实验目的1. 学习由系统函数确定系统频率特性的方法;2. 学习和掌握连续时间系统的频率特性及其幅度特性、相位特性的物理意义;3.通过本实验了解低通、高通、带通、全通滤波器的性能及特点;二、实验原理及方法频域分析法与时域分析法的不同之处主要在于信号分解的单元函数不同;在频域分析法中,信号分解成一系列不同幅度、不同频率的等幅正弦函数,通过求取对每一单元激励产生的响应,并将响应叠加,再转换到时域以得到系统的总响应;所以说,频域分析法是一种变域分析法;它把时域中求解响应的问题通过 Fourier 级数或 Fourier 变换转换成频域中的问题；在频域中求解后再转换回时域从而得到最终结果;在实际应用中,多使用另一种变域分析法：复频域分析法,即 Laplace 变换分析法;所谓频率特性,也称频率响应特性,是指系统在正弦信号激励下稳态响应随频率变化的情况,包括幅度随频率的响应和相位随频率的响应两个方面;利用系统函数也可以确定系统频率特性,公式如下：幅度响应用()ωj H 表示,相位响应用)(ωϕH 表示;本实验所研究的系统函数Hs 是有理函数形式,也就是说,分子、分母分别是m 、n 阶多项式; 要计算频率特性,可以写出为了计算出()ωj H 、)(ωϕH 的值,可以利用复数三角形式的一个重要特性： 而⎥⎦⎤⎢⎣⎡+=2sin 2cosππωωj j ,则()⎥⎦⎤⎢⎣⎡+=2sin 2cos ππωωn j n j n n利用这些公式可以化简高次幂,因此分子和分母的复数多项式就可以转化为分别对实部与虚部的实数运算,算出分子、分母的实部、虚部值后,最后就可以计算出幅度()ωj H 、相位)(ωϕH 的值了;三、实验内容a)sm m ms H )(1)(2-+=,m 取值区间 0,1,绘制一组曲线 m=,,,,； b) 绘制下列系统的幅频响应对数曲线和相频响应曲线,分析其频率特性; a %figurealpha=,,,,;colorn='r' 'g' 'b' 'y' 'k'; % r g b y m c k 红,绿,蓝,黄,品红,青,黑 for n=1:5b=0 alphan; % 分子系数向量a=alphan-alphan^2 1; % 分母系数向量 printsysb,a,'s' Hz,w=freqsb,a; w=w./pi; magh=absHz;zerosIndx=findmagh==0; maghzerosIndx=1; magh=20log10magh; maghzerosIndx=-inf; angh=angleHz;angh=unwrapangh180/pi; subplot1,2,1plotw,magh,colornn;hold onsubplot1,2,2plotw,angh,colornn;hold onendsubplot1,2,1hold offxlabel'特征角频率\times\pi rad/sample' title'幅频特性曲线 |Hw| dB';subplot1,2,2hold offxlabel'特征角频率 \times\pi rad/sample' title'相频特性曲线 \thetaw degrees';b1 %b=1,0; % 分子系数向量a=1,1; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';2 %b=0,1,0; % 分子系数向量a=1,3,2; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';3 %b=1,-1; % 分子系数向量a=1,1; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';实验心得：虽然之前用公式转换到频域上分析,但是有时会觉得挺抽象的,不太好理解;根据这些图像结合起来更进一步对信号的了解;同时,这个在编程序时,虽然遇到一些问题,但是总算解决了;实验六离散时间系统的Z域分析一、 实验目的1. 学习和掌握离散系统的频率特性及其幅度特性、相位特性的物理意义;2. 深入理解离散系统频率特性和对称性和周期性;3. 认识离散系统频率特性与系统参数之间的系统4.通过阅读、修改并调试本实验所给源程序,加强计算机编程能力; 二、 实验原理及方法对于离散时间系统,系统单位冲激响应序列)(n h 的 Fourier 变换)(ωj e H 完全反映了系统自身的频率特性,称)(ωj eH 为离散系统的频率特性,可由系统函数)(z H 求出,关系式如下：ωωj j e z z H e H ==)()( 6 – 1由于ωj e是频率的周期函数,所以系统的频率特性也是频率的周期函数,且周期为π2,因此研究系统频率特性只要在πωπ≤≤-范围内就可以了;∑∑∑∞-∞=∞-∞=∞-∞=--==n n n j j n n h j n n h en h e H )sin()()cos()()()(ωωωω6 – 2容易证明,其实部是ω的偶函数,虚部是ω的奇函数,其模ωj e H (的ω的偶函数,相位[])(arg ωj e H 是ω的奇函数;因此研究系统幅度特性)(ωj e H 、相位特性[])(arg ωj e H ,只要在πω≤≤0范围内讨论即可;综上所述,系统频率特性)(ωj eH 具有周期性和对称性,深入理解这一点是十分重要的;当离散系统的系统结构一定,它的频率特性)(ωj e H 将随参数选择的不同而不同,这表明了系统结构、参数、特性三者之间的关系,即同一结构,参数不同其特性也不同; 例如,下图所示离散系统,其数学模型由线性常系数差分方程描述：)()1()(n x n ay n y +-=系统函数：a z az z z H >-=,)(系统函数频率特性：ωωωωωsin )cos 1(1)(ja a a e e e H j j j +-=-=幅频特性：ωωcos 211)(2a a eH j -+=相频特性：[]ωωωcos 1sin arctan)(arg a a eH j --= 容易分析出,当10<<a 时系统呈低通特性,当01<<-a 时系统呈高通特性；当0=a 时系统呈全通特性;同时说明,在系统结构如图所示一定时,其频率特性随参数a 的变化而变化;三、 实验内容a 2281.011)(----=z z z H ;b 1.04.06.01.03.03.01.0)(2323+++-+-=z z z z z z z Hc 2181.011)(--+-=zz z H a %b=1,0,-1; % 分子系数向量a=1,0,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';带通b %b=,,,; % 分子系数向量a=1,,,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';高通c %b=1,-1,0; % 分子系数向量a=1,0,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';带通实验心得：本来理论知识不是很强的,虽然已经编出程序得到相关图形,但是不会辨别相关通带,这让我深刻地反省;。
信号与系统(教案) 第二章
二、图解机理
用图形方式理解卷积运算过程，包括以下6个步骤： Step1：换元。画出f1(t)与f2(t)波形，将波形图中的t轴 改换成τ轴，分别得到f1(τ)和f2(τ)。 Step2：翻转。将f2(τ)波形以纵轴为中心轴翻 180°，得 到f2(-τ)波形。 4
信号与系统
2.2
卷积积分
Step3:平移。给定t值,将f2(-τ)波形沿τ轴平移|t|。
卷积积分是一种数学运算，它有许多重要的性质 （或运算规则），灵活地运用它们能简化卷积运算。 下面讨论均设卷积积分是收敛的（或存在的）。
性质1.卷积代数 满足乘法的三律： 1. 交换律： f1(t)* f2(t) =f2(t)* f1(t) 2. 分配律： f1(t)*[ f2(t)+ f3(t)] =f1(t)* f2(t)+ f1(t)* f3(t) 3. 结合律： [f1(t)* f2(t)]* f3(t)] =f1(t)*[ f2(t) * f3(t)]
1.奇异信号
单位冲激信号 (t), 单位阶跃信号 (t).
2.正弦信号
也称为虚指数信号。 f (t ) A cos( t ) A [e j (t ) e j (t ) ] 2
式 中A、和分 别 为 正 弦 信 号 的 振 幅 角 频 率 和 初 相 。 、 f ( t )是 周 期 信 号 ， 其 周 期 2 T＝
1 0
f 1(t)
2
t
14
信号与系统 例：f1(t), f2(t)如图，求f1(t)* f2(t) 解： f1(t) = 2ε (t) –2ε (t –1) f2(t) = ε (t+1) –ε (t –1)
2.2 卷积积分 2.2 卷积积分
《信号与系统》实验报告
信号与系统实验报告班级:姓名:信息与通信工程学院实验一 系统的卷积响应实验性质：提高性 实验级别：必做 开课单位：信息与通信工程学院 学 时：2一、实验目的：深刻理解卷积运算，利用离散卷积实现连续卷积运算；深刻理解信号与系统的关系，学习MATLAB 语言实现信号通过系统的仿真方法。
二、实验设备： 计算机，MATLAB 软件 三、实验原理： 1、 离散卷积和： 调用函数：conv （）∑∞-∞=-==i i k f i f f f conv S )()(1)2,1(为离散卷积和，其中，f1(k), f2 (k) 为离散序列，K=…-2, -1, 0 , 1, 2, …。
但是，conv 函数只给出纵轴的序列值的大小，而不能给出卷积的X 轴序号。
为得到该值，进行以下分析：对任意输入：设)(1k f 非零区间n1~n2，长度L1=n2-n1+1；)(2k f 非零区间m1~m2，长度L2=m2-m1+1。
则：)(*)()(21k f k f k s =非零区间从n1+m1开始，长度为L=L1+L2-1，所以S （K ）的非零区间为：n1+m1~ n1+m1+L-1。
2、 连续卷积和离散卷积的关系：计算机本身不能直接处理连续信号，只能由离散信号进行近似： 设一系统（LTI ）输入为)(t P ∆，输出为)(t h ∆，如图所示。
)t)()(t h t P ∆∆→)()(lim )(lim )(0t h t h t P t =→=∆→∆∆→∆δ若输入为f(t):∆∆-∆=≈∑∞-∞=∆∆)()()()(k t P k f t f t f k得输出：∆∆-∆=∑∞-∞=∆∆)()()(k t hk f t y k当0→∆时：⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==ττδτd t f k t P k f t f t f k )()()()(lim)(lim )(0⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==τττd t h f k t hk f t y t y k )()()()(lim)(lim )(0所以：∆∆-∆=-==∑⎰→∆)()(lim)()()(*)()(212121k t f k fd t f f t f t f t s τττ如果只求离散点上的f 值)(n f ∆])[()()()()(2121∑∑∞-∞=∞-∞=∆-∆∆=∆∆-∆∆=∆k k k n f k f k n f k fn f所以，可以用离散卷积和CONV （）求连续卷积，只需∆足够小以及在卷积和的基础上乘以∆。
信号与系统实验报告实验一 信号与系统的时域分析
实验一信号与系统的时域分析一、实验目的1、熟悉与掌握常用的用于信号与系统时域仿真分析的MA TLAB函数;2、掌握连续时间与离散时间信号的MA TLAB产生,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MA TLAB编程;3、牢固掌握系统的单位冲激响应的概念,掌握LTI系统的卷积表达式及其物理意义,掌握卷积的计算方法、卷积的基本性质;4、掌握利用MA TLAB计算卷积的编程方法,并利用所编写的MA TLAB程序验证卷积的常用基本性质;掌握MA TLAB描述LTI系统的常用方法及有关函数,并学会利用MATLAB求解LTI系统响应,绘制相应曲线。
基本要求:掌握用MA TLAB描述连续时间信号与离散时间信号的方法,能够编写MATLAB程序,实现各种信号的时域变换与运算,并且以图形的方式再现各种信号的波形。
掌握线性时不变连续系统的时域数学模型用MA TLAB描述的方法,掌握卷积运算、线性常系数微分方程的求解编程。
二、实验原理信号(Signal)一般都就是随某一个或某几个独立变量的变化而变化的,例如,温度、压力、声音,还有股票市场的日收盘指数等,这些信号都就是随时间的变化而变化的,还有一些信号,例如在研究地球结构时,地下某处的密度就就是随着海拔高度的变化而变化的。
一幅图片中的每一个象素点的位置取决于两个坐标轴,即横轴与纵轴,因此,图像信号具有两个或两个以上的独立变量。
在《信号与系统》课程中,我们只关注这种只有一个独立变量(Independent variable)的信号,并且把这个独立变量统称为时间变量(Time variable),不管这个独立变量就是否就是时间变量。
在自然界中,大多数信号的时间变量都就是连续变化的,因此这种信号被称为连续时间信号(Continuous-Time Signals)或模拟信号(Analog Signals),例如前面提到的温度、压力与声音信号就就是连续时间信号的例子。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
连续时间信号的卷积运算的MATILAB实现
薛皓20091453
例1：已知两连续时间信号如图9-3所示，试用matlab求f(t)=f1(t)*f2(t)，并绘出f（t）的时域波形图。
图1-1 连续时间信号波形图示例
实现上述过程的matlab命令如下：
p=0.5;
k1=0:p:2;
f1=0.5*k1;
k2=k1;
f2=f1;
[f,k]=sconv(f1,f2,k1,k2,p)
上述命令绘制的波形图也在图9-3中示出。
图9-3中给出了抽样时间间隔p=0.5时的处理效果。
而图9-4给出了抽样时间间隔p=0.01时的处理效果。
图1-2 例1的连续时间信号波形图
习题1：已知f1(t)=1(2t 1≤≤)，f2(t)=1(3t 2≤≤)，用matlab 实现其卷积并绘制出卷积曲线。
解：程序代码如下：
>> p=0.01;
k1=1:p:2;
f1=ones(size(k1)).*(k1>1);
k2=2:p:3;
f2=ones(size(k2)).*(k2>2);
f=conv(f1,f2);
f=f*p;
k0=k1(1)+k2(1);
k3=k1(length(k1))+k2(length(k2));
subplot(2,2,1)
plot(k1,f1)
title('f1(t)')
xlabel('t')
ylabel('f1(t)')
subplot(2,2,2)
plot(k2,f2)
title('f2(t)')
xlabel('t')
ylabel('f2(t)')
subplot(2,2,3)
plot(k,f);
h=get(gca,'position');
h(3)=2.5*h(3); 0
set(gca,'position',h)
title('f(t)=f1(t)*f2(t)')
xlabel('t')
ylabel('f(t)')
绘制图形如图2-1所示。
图2-1
习题2：)1()2/1t ()t (2f ),1t ()t ()t (1f δ-+δ=-ε-ε=，求其卷积。
程序代码：
p=0.01;
t1=0:p:1;
f1=ones(size(t1)).*(t1>0);
t2=-0.5:p:1;
f2=(t2==-0.5)-(t2==1);
f=conv(f1,f2);
f=f*p;
t=-0.5:p:2;
k1=t1;k2=t2;k=t;
subplot(2,2,1)
plot(k1,f1)
title('f1(t)')
xlabel('t')
ylabel('f1(t)')
subplot(2,2,2)
plot(k2,f2)
title('f2(t)')
xlabel('t')
ylabel('f2(t)')
subplot(2,2,3)
plot(k,f);
h=get(gca,'position');
h(3)=2.5*h(3);
set(gca,'position',h)
title('f(t)=f1(t)*f2(t)')
xlabel('t')
ylabel('f(t)')
绘制波形如图2-2示：
图2-3
实验心得：感觉知识太贫乏了，知识尚且不够，何谈智慧呢，知耻后勇吧！。