1、 1.1 离散时间信号序列 1.2 线性移不变系统 1.3 常系数线性差分方程 1.4 连续时间信号的抽样 第1章 离散时间信号和系统天津师范大学计算机与信息工程学院 1.1 离散时间信号序列 1.1.1 序列的定义 1.1.2 序列的基本运算 1.1.3 常用的基本序列 1.1.4 序列的周期性 1.1.5 用单位脉冲表示任意序列 天津师范大学计算机与信息工程学院1.1.1 序列的定义 信号在数学上定义为一个函数，这个函 数表示一种信息，通常是关于一个物理系统的状态或特性的。信号的函数表示是关于一个或几个独立变量的，关于一个独立变量的信号称为一维信号，关于多个独立变量的信号称为多维信号。在本
2、中，主要 的信号是一信号x(t), 一般情况下x(t)随化的信号，称 信号或域信号。天津师范大学计算机与信息工程学院 若t是定义在时间上的连续变量，称x(t)为连续时间信号，也就是模拟信号；若t仅在时间的离散点上取值，称x(t)为离散时间信号或时域离散信号。离散时间信号可以通过对连续时间信号的采样得到,这种情况下把信号记为x(nT) ，T 表示的是采样点之间的时间间隔，n是一个整数。 离散时间信号可以表示成下列形式： x(nT) n=0,1,2,3,. 天津师范大学计算机与信息工程学院 在大多数DSP系统中，x(nT)的存放是按n下标来放置的，不同的x(nT)只要靠n就可区别。因此，将x(nT
3、)表示为x(n)，这是一种数学的抽象。所以一个离散时间信号定义为： x(n) n=0,1,2,3,.x(n)定在n等于整数点上，在n不等于整数点上，x(n)没有定，但并不表示信号 零。从数学的角度看，上面的定式表示一个序列，所以也把离散时间信号称作离散时间序列，常常化x(n) 。 天津师范大学计算机与信息工程学院 序列除了数学表达式外，还常常采用图形方式来表示，如图1.1所示。虽然横坐标画成一条连续的直线，但x(n)仅仅对于整数的n值才有意义。图1.1 离散时间信号的图形表示 天津师范大学计算机与信息工程学院 离散时间信号在幅度上定义成连续 的，如果将幅度进行量化，一般为等间隔量化。在时间和幅
4、度上都取离散值的信号称为“数字信号”。因此，离散时间信号并不等于数字信号，但由于数字信号是幅度量化得到的，在数学表示和推导中不如序列形式方便和容易，所以一般都采用离散时间信号来讨论数字信号处理的理论和算法。天津师范大学计算机与信息工程学院1.1.2 序列的基本运算 和和 积积 移位移位 标乘标乘 翻转翻转 累加累加 差分差分 时间尺度变换时间尺度变换 序列的能量序列的能量 卷积和卷积和天津师范大学计算机与信息工程学院基本运算序列的和 设序列为设序列为x x( (n n) )和和y y( (n n) )，则序列，则序列z z( (n n) )= x= x( (n n) )+ y+ y( (n n
5、) ) 表示两个序列的和，定义为表示两个序列的和，定义为同序号同序号的序的序列值列值逐项对应逐项对应相加。相加。天津师范大学计算机与信息工程学院例：序列的和例例1.1.11.1.1 设序列设序列计算序列的和计算序列的和x x( (n n)+ )+ y y( (n n) )。解：解：天津师范大学计算机与信息工程学院例：序列求和图示天津师范大学计算机与信息工程学院基本运算序列的积 设序列为设序列为x x( (n n) )和和y y( (n n) )，则序列，则序列z z( (n n) )= x= x( (n n) ) y y( (n n) ) 表示两个序列的积，定义为表示两个序列的积，定义为同序号
6、同序号的序的序列值列值逐项对应逐项对应相乘。相乘。天津师范大学计算机与信息工程学院例：序列的积例例1.1.21.1.2 设序列设序列计算序列的积计算序列的积x x( (n n) ) y y( (n n) )。解：解：天津师范大学计算机与信息工程学院例：序列求积图示x(n)天津师范大学计算机与信息工程学院基本运算序列的移位 设序列为设序列为x x( (n n) )，则序列，则序列y y( (n n)= )= x x( (n n- -mm) ) 表示将序列表示将序列x x( (n n) )进行移位。进行移位。 mm为正时为正时x x( (n n - -mm) )：x x( (n n) )逐项依次逐
7、项依次延时延时( (右移右移) )mm位位x x( (n n+ +mm) )：x x( (n n) )逐项依次逐项依次超前超前( (左移左移) )mm位位 m m为负时，则为负时，则相反相反。天津师范大学计算机与信息工程学院例：序列的移位例例1.1.31.1.3 设序列设序列计算序列的移位序列计算序列的移位序列x x( (n n+1)+1)。解：解：天津师范大学计算机与信息工程学院例：序列移位图示x(n)天津师范大学计算机与信息工程学院基本运算序列的标乘 设序列为设序列为x x( (n n) )，a a为常数为常数(a 0)(a 0)，则序列，则序列y y( (n n)= a)= ax x(
8、(n n) ) 表示将序列表示将序列x x( (n n) )的标乘，定义为的标乘，定义为各序列值各序列值均乘以均乘以a a，使新序列的，使新序列的幅度幅度为原序列的为原序列的a a倍。倍。天津师范大学计算机与信息工程学院例：序列的标乘例例1.1.41.1.4 设序列设序列计算序列计算序列4 4x x( (n n) )。解：解：天津师范大学计算机与信息工程学院基本运算序列的翻转 设序列为设序列为x x( (n n) )，则序列，则序列y y( (n n)= )= x x(- (-n n) ) 表示以表示以n n= 0= 0的的纵轴为对称轴纵轴为对称轴将序列将序列x x( (n n) )加加以翻转
9、。以翻转。天津师范大学计算机与信息工程学院例：序列的翻转例例1.1.51.1.5 设序列设序列计算序列计算序列x x(- (-n n) )。解：解：天津师范大学计算机与信息工程学院基本运算序列的累加 设序列为x(n)，则序列 定义为对x(n)的累加，表示将n 以前的所有x(n)值求和。天津师范大学计算机与信息工程学院例: 序列的累加 设序列为则其累加序列即 y(0)=x(0)=1, y(1)=x(0)+x(1)=y(0)+x(1)=3, y(2)=y(1)+x(2)=7天津师范大学计算机与信息工程学院基本运算序列的差分 前向差分：将序列前向差分：将序列先先进行进行左移左移，再相减，再相减x(n
10、) = x(n+1)- x(n) x(n) = x(n+1)- x(n) 后向差分：将序列后向差分：将序列先先进行进行右移右移，再相减，再相减x(n) = x(n)- x(n-1) x(n) = x(n)- x(n-1) 由此容易得出由此容易得出x(n) = x(n-1)x(n) = x(n-1)天津师范大学计算机与信息工程学院多阶差分运算 二阶前向差分二阶前向差分 二阶后向差分二阶后向差分 单位延迟算子单位延迟算子D D，有，有 DyDy( (n n)= )= y y( (n-n-1) 1) y y( (n n)= )= y y( (n n)- )- y y( (n-n-1)= 1)= y
11、y( (n n)- )- Dy Dy( (n n)= (1- )= (1- D D) )y y( (n n) )= 1-= 1-D D k k 阶后向差分阶后向差分(按二项式定理展开)二阶后向差分二阶后向差分天津师范大学计算机与信息工程学院例: 差分运算例1.1.6 设序列 求x(n)x(n)和和x(n)x(n)。解：解：前向差分前向差分 天津师范大学计算机与信息工程学院例: 差分运算后向差分天津师范大学计算机与信息工程学院基本运算时间尺度(比例)变换 设序列为设序列为x x( (n n) )，mm为正整数，则序列为正整数，则序列 抽取序列：抽取序列：y y( (n n)= )= x x( (
12、mnmn) ) x x( (mnmn) ) 和和x x( (n/mn/m) )定义为对定义为对x(n)x(n)的时间尺度变的时间尺度变换。换。 插值序列：插值序列： 天津师范大学计算机与信息工程学院抽取序列 x x( (mnmn) )：对：对x x( (n n) )进行抽取运算进行抽取运算 不是简单在时间轴上按比例增加到不是简单在时间轴上按比例增加到mm倍倍 以以1/1/mm倍的取样频率倍的取样频率每隔每隔mm-1-1个点个点抽取抽取1 1点。点。 保留保留 x x(0)(0)天津师范大学计算机与信息工程学院插值序列 x x( (n/mn/m) ) ：对：对x x( (n n) )进行插值运算
13、进行插值运算 表示在原序列表示在原序列x x( (n n) )相邻两点相邻两点之间插入之间插入mm-1-1个零个零值点值点 保留保留 x x(0)(0)天津师范大学计算机与信息工程学院基本运算序列的能量 设序列为设序列为x x( (n n) )，则序列，则序列 定义为序列的能量，表示序列各取样值的定义为序列的能量，表示序列各取样值的平方平方之和；之和； 若为复序列，取若为复序列，取模值模值后再求平方和。后再求平方和。天津师范大学计算机与信息工程学院基本运算序列的卷积和 设序列为设序列为x x( (n n) )和和z z( (n n) )，则序列，则序列 定义为序列定义为序列x x( (n n)
14、 )和和z z( (n n) )的的卷积和卷积和。卷积和卷积和又称为又称为离散卷积离散卷积或或线性卷积线性卷积，是，是很重要很重要的公的公式。式。天津师范大学计算机与信息工程学院卷积和计算的四个步骤 翻转翻转：x x( (mm) ) ，z z( (mm) ) z z(- (-mm) ) 移位移位：z z(- (-mm) ) z z( (n n- -mm) ) n n为正数时，右移为正数时，右移n n位位 n n为负数时，左移为负数时，左移n n位位 相乘相乘：z z( (n n- -mm) ) x x( (mm) ) （mm值相同）值相同） 相加相加：y y( (n n) =) =z z( (
15、n n- -mm) ) x x( (mm) )天津师范大学计算机与信息工程学院对应点相乘！对应点相乘！例：卷积和计算例例1.1.71.1.7 设序列设序列求求y y( (n n)= )= x x( (n n)* )*z z( (n n) ) 。解：解：n n n n0 0时，时，x(m)x(m)与与z(n-m)z(n-m)没有重叠没有重叠，得，得y(n)=0y(n)=0。n n 00n n44时，时，对应点相乘！对应点相乘！天津师范大学计算机与信息工程学院例：卷积和计算n n 4 4n6n6时，时，n n 6 6n10n10时，时，n nn n1010时，时，x x( (mm) )与与z z(
16、 (n-mn-m) )没有重叠，得没有重叠，得y y( (n n)= 0)= 0。 天津师范大学计算机与信息工程学院1.1.3 几种常用序列 单位脉冲单位脉冲( (抽样抽样) )序列序列 单位阶跃序列单位阶跃序列 矩形序列矩形序列 实指数序列实指数序列 正弦序列正弦序列 复指数序列复指数序列 天津师范大学计算机与信息工程学院单位脉冲序列 ( (n n) )只在只在n n= 0= 0时取确定时取确定值值1 1，其它均为零，其它均为零 ( (n n) )类似于类似于 ( (t t) ) ( (n n- -mm) )只有在只有在n n= = mm时时取确定值取确定值1 1，而其余点，而其余点取值均为
17、零取值均为零 天津师范大学计算机与信息工程学院单位阶跃序列 u u( (n n) )类似于类似于u u( (t t) ) u u( (t t) )在在t t= 0= 0时常不定义，时常不定义，u u( (n n) )在在n n= 0= 0时为时为u u(0)= 1(0)= 1 ( (n n) )和和u u( (n n) )的关系：的关系：(n) = u(n)-u(n-1) 天津师范大学计算机与信息工程学院单位矩形序列 N N 为矩形序列的为矩形序列的长度长度 和和u u( (n n) )、 ( (n n) )的关系的关系 ：天津师范大学计算机与信息工程学院实指数序列 a a为实数为实数 当当|
18、a|a|1 1时序列时序列收敛收敛 当当|a|a|1 1时序列时序列发散发散 天津师范大学计算机与信息工程学院正弦序列 A A为幅度为幅度 为数字域角频率为数字域角频率 为起始相位为起始相位 x x( (n n) )由由x x( (t t)= sin)= sint t 取样取样得到x(n)= Asin(n+) 归一化归一化: : =T=T =/fs (与与线性关系 )天津师范大学计算机与信息工程学院复指数序列 为为数字域数字域角频率角频率 用用实部实部与与虚部虚部表示表示 用极坐标极坐标表示 =0=0时，序列具有以时，序列具有以22为周期的为周期的周期性周期性 复指数序列在实际中不存在，它是为
19、了数学上的表示和分析方便而引入的，它的特性和正弦或余弦序列的特性基本一致。天津师范大学计算机与信息工程学院1.1.4 序列的周期性 对于序列对于序列x x( (n n) )，如果对所有，如果对所有n n 存在一个最小存在一个最小的正整数的正整数N N，满足，满足x x( (n n)= )= x x( (n+Nn+N) )则序列则序列x x( (n n) )是周期序列是周期序列 ，最小周期最小周期为为N N 。 以以正弦序列正弦序列 为例讨论周期性 设设 x x( (n n)=)= A Asin(sin(nn+ + ) ) 则有则有 x x( (n+Nn+N) =) =A Asinsin( (n
20、+Nn+N)+)+ = =A Asin(sin(NN+ +n+n+) ) 若若满足条件满足条件N= N= 2 2k k，则，则x x( (n+Nn+N)=)= A Asinsin( (n+Nn+N)+)+ = = A Asin(sin(n+n+) ) = x= x( (n n) )天津师范大学计算机与信息工程学院周期性讨论 N N、k k 为整数，为整数，k k 的取值满足条件，且的取值满足条件，且保证保证N N 是最小正整数。其周期为是最小正整数。其周期为 2/2/为为整数整数时，取时，取k k = 1= 1，保证为最小正整，保证为最小正整数。此时为周期序列，周期为数。此时为周期序列，周期为
21、2/2/。 例例1.1.81.1.8 序列序列 ，因为，因为2/2/= 8= 8，所以，所以是一个周期序列，其周期是一个周期序列，其周期N N= 8= 8。 例例1.1.81.1.8 序列序列 ，因为，因为2/2/= 8= 8，所以，所以是一个周期序列，其周期是一个周期序列，其周期N N= 8= 8。 例例1.1.81.1.8 序列序列 ，因为，因为2/2/= 8= 8，所以，所以是一个周期序列，其周期是一个周期序列，其周期N N= 8= 8。 天津师范大学计算机与信息工程学院周期性讨论2/2/为为有理数而非整数有理数而非整数时，仍然是周期序时，仍然是周期序列，周期大于列，周期大于2/2/。例
22、例1.1.91.1.9 序列序列 ，2/2/= 8/3= 8/3是有理数是有理数，所以是周期序列，取，所以是周期序列，取k k= 3= 3，得到周期，得到周期N N= 8= 8。 2/2/为为无理数无理数时，时，任何k 都不能使N 为正整数，这时正弦序列不是周期序列。 指数为纯虚数指数为纯虚数的复指数序列的周期性的复指数序列的周期性与正弦序列的情况与正弦序列的情况相同相同。 例例1.1.101.1.10 序列序列 , 2/ , 2/= 8/3= 8/3是无是无理数，所以不是周期序列。理数，所以不是周期序列。例例1.1.91.1.9 序列序列 ，2/2/= 8/3= 8/3是有理数是有理数，所以
23、是周期序列，取，所以是周期序列，取k k= 3= 3，得到周期，得到周期N N= 8= 8。 天津师范大学计算机与信息工程学院周期性讨论 判断一个正弦序列是否是周期序列的方法是： 用2除以它的数字频率，若得出的是整数或有理数，则序列为周期序列；若得出的是无理数，序列就不是周期序列。 但无论序列是否为周期序列，仍把称作序列的数字频率。天津师范大学计算机与信息工程学院下面来说明模拟频率和数字频率之间的关系。设模拟正弦信号为 对该 以T为采样间隔进行采样离散,得 将离散后的信号表示成离散正弦序列，即 天津师范大学计算机与信息工程学院可知 其中， 称为采样频率。该式即为数字频率和模拟角频率0、模拟频率
24、f 0之间关系式，它们是依靠采样间隔 T 或采样频率 f s 进行关联的。整理后可得可以看出： 是一个相对频率，它是连续正弦信号的频率f0 对抽样频率 fs的相对频率乘以2，或说是连续正弦信号的角频率0对抽样频率 fs 的相对频率。天津师范大学计算机与信息工程学院 数字频率的特点：（1）是一个连续取值的量；（2）的量纲为一种角度的量纲单位：弧度（rad）。它表示序列在采样间隔T内正弦信号变化的角度，表示了信号相对变化的快慢程度；（3） 序列对于是以2为周期的，或者说，的独立取值范围为0,2)或-,)。天津师范大学计算机与信息工程学院正弦型序列是周期序列的条件为：（有理数）则当 N 个抽样间隔等
25、于 k 个连续时间信号的周期时，由正弦信号抽样得到的正弦序列是周期序列。天津师范大学计算机与信息工程学院1.1.5 用单位脉冲序列表示任意序列 任何序列都可以用单位脉冲序列的移位加权和来表示，任何序列都可以用单位脉冲序列的移位加权和来表示，即即式中式中天津师范大学计算机与信息工程学院例如如图序列，可以表示成天津师范大学计算机与信息工程学院 x x( (n n) ) 可看成是可看成是x x( (n n) )和和 ( (n n) )的卷积和，即的卷积和，即 例例1.1.111.1.11 天津师范大学计算机与信息工程学院1.2 线性移不变系统1.2.1 系统定义 数字信号处理的任何处理都是依靠系统来
26、完成的,所以系统是数字信号处理的核心，系统一般包括系统硬件和系统所完成的处理算法。 系统在数学上定义为将输入序列x(n)映射成输出序列y(n)的唯一性变换或运算。这种映射是广义的，实际上表示的是一种具体的处理，或是变换，或是滤波。天津师范大学计算机与信息工程学院系统可以表示为 其中，符号T 表示系统的映射或处理，可以把T 简称为系统。 系统的图形表示如下图所示，输入x(n)称为系统的激励，输出y(n)称为系统的响应。由于它们均为离散时间信号，将系统T 称为离散时间系统或时域离散系统。天津师范大学计算机与信息工程学院1.2.2 线性离散时间系统 满足叠加原理的系统，或满足齐次性和可加性的系统称为
27、线性系统。 设y1(n)=T x1(n)， y2(n)=T x2(n)对任意常数a,b，若 T ax1(n)+bx2(n)=aT x1(n)+bT x2(n) =a y1(n)+b y2(n)则称T 为线性离散时间系统。天津师范大学计算机与信息工程学院 推广到一般情况，设 yk(n) = T xk(n) ， k=1，2，.N线性系统满足 1kN 线性系统的特点是多个输入的线性组合的系统输出等于各输入单独作用的输出的线性组合。天津师范大学计算机与信息工程学院例1.2.1 证明由线性方程表示的系统是非线性系统。证明 设 所以，该系统是非线性系统。天津师范大学计算机与信息工程学院1.2.3 非时变离
28、散时间系统 若满足下列条件，系统称为非时变（非移变）系统，或时不变（移不变）系统。 设 y(n) = T x(n) 对任意整数k,有 y(n-k)=T x(n-k) 即系的映射T 不随 化，只要入x(n)是相同的，无何 行激励，出y(n)是相同的，正是系非 性的特征。 天津师范大学计算机与信息工程学院下图形象说明了系统非时变性的概念。天津师范大学计算机与信息工程学院例1.2.2 设系统的映射 y = T x(n) = nx(n) ，判断系统的线性和时不变性。解 设 y1(n) = nx1(n), y2(n) =nx2(n) x(n) = a1x1(n)+a2x2(n) 则 T x(n) = n
29、x(n) = na1x1(n)+na2x2(n) = a1y1(n)+ a2y2(n) 所以，系统为线性系统。设 y(n) = nx(n), x1(n) = x(n-k) y1(n) = nx1(n) = nx(n-k) 而 y(n-k) = (n-k)x(n-k)y1(n)所以，系统为时变系统。天津师范大学计算机与信息工程学院 1.2.4 线性时不变离散系统 定义 同时具备线性和时不变性的系统称作线 性非时变系统或线性时不变系统。它的重要意义在于，系统的处理过程可以统一采用这种系统的特征描述之一单位取样响应，以一种相同的运算方式卷积运算，进行统一的表示。任何一个信号可以表示成单位取样序列的线
30、性组合，即天津师范大学计算机与信息工程学院 系统对 的响应为 设系统对单位取样序列 的响应为 ， 即 称为系统的“单位取样响应”，它是描述系统的一个非常重要的信号。天津师范大学计算机与信息工程学院 根据时不变性，有 则系统输出y(n)可表示为 上式表明：当线性非时变系统的单位取样响应h(n)确定时，系统对任何一个输入x(n)的响应y(n)就确定了， y(n)可以表示成x(n)和h(n)之间的一种简单的运算形式。将上式的运算方式称作“离散卷积”，简称“卷积”，采用符号“*”表示，即天津师范大学计算机与信息工程学院 1.2.5 离散卷积的运算规律 (1) 交换律 h(n)*x(n) = x(n)*
31、h(n) 它的意义可以解释为，如果互换系统的单位取样响应h(n)和输入x(n)，系统的输出保持不变。 x(n)h(n)y(n)=h(n)x(n)y(n)天津师范大学计算机与信息工程学院交换律证明：令n-m=k（m=n-k），则天津师范大学计算机与信息工程学院 (2) 结合律 x(n)*h1(n)*h2(n)=x(n)*h2(n)*h1(n) =x(n)*h2(n)*h1(n) 它的意义可以解释为一种级联系统结构，级联顺序可以交换，或系统级联可以等效为一个系统，输出保持不变。x(n)y(n)h1(n)h2(n)x(n)y(n)h2(n)h1(n)x(n)y(n)h1(n)*h2(n)天津师范大学
32、计算机与信息工程学院结合律证明：令k-m=r（k=r+m），则天津师范大学计算机与信息工程学院 (3) 分配律 x(n)*h1(n)+x(n)*h2(n) = x(n)*h1(n)+ h2(n) 它的意义可以解释为一个并联系统结构，或并联系统可以等效为一个系统，输出保持不变。 x(n)y(n)h1(n)+h2(n)h1(n)h2(n)x(n)y(n)天津师范大学计算机与信息工程学院分配律证明：天津师范大学计算机与信息工程学院 (4) 与(n)卷积的不变性 x(n)*(n) = x(n) 它的意义可以解释为输入通过一个零相位的全通系统。 (5) 与(n-k)卷积的移位性 x(n)*(n-k) =
33、 x(n-k) 它的意义可以解释为输入通过一个线性相位的全通系统。天津师范大学计算机与信息工程学院 1.2.6 离散卷积的计算 卷积的计算一般采用两种方法：解析法和图解法，或是两种方法的结合。 例1.2.3 设线性时不变系统的单位脉冲响应和输入序列如下图所示，画出输出的波形。天津师范大学计算机与信息工程学院解：（1）采用图解法。 解法的程如1.2所示。天津师范大学计算机与信息工程学院 图1.2 例1.2.3图解法天津师范大学计算机与信息工程学院（2）采用解析法。 因为 所以 将x(n)的表达式代入上式，得到 两种方法结果一致。天津师范大学计算机与信息工程学院1.2.7 系统的稳定性和因果性 一
34、、 稳定性 稳定系统是有界输入产生有界输出的系统。 若 则 线性时不变离散系统是稳定系统的充要条件（稳定性定理）：即，系统的单位抽样响应绝对可和。天津师范大学计算机与信息工程学院证明：充分条件 若系统满足条件 ，且输入x(n)有界，有，对所有n，M是一个任意大的有限数，此时系统的输出为两边取绝对值，得即输出y(n)有界，故系统是稳定的。 必要条件 利用反证法，已知系统稳定，假设 ，可以找到一个有界的输入则即输出无界，这不符合稳定的假设，因而假设不成立，所以是稳定的必要条件。天津师范大学计算机与信息工程学院 二、 因果性 若系统 n时刻的输出，只取决于n时刻以及n时刻以前的输入序列，而与n时刻以
35、后的输入无关，则称该系统为因果系统。 线性时不变离散系统是因果系统的充要条件（因果性定理）：天津师范大学计算机与信息工程学院证明：充分条件 若n0时，h(n)=0，根据卷积和公式 因为只有当n-m 0时，h( n-m)才有值，所以mn，这就证明了y(n)的值只取决于x(n)在nn时的值，因此系统是因果的。 必要条件 利用反证法，已知因果系统，假设当nn时的值有关，也就是y(n)值与n以后的x(n)有关，所以该系统不是因果系统，与已知条件矛盾，因而假设不成立。可见要使y(n)与nn时的x(n)无关，则必须使 天津师范大学计算机与信息工程学院 结论:因果稳定的线性时不变系统的单位取样响应是因果的,
36、且是绝对可和的,即天津师范大学计算机与信息工程学院 例1.2.4 某线性时不变离散系统，其单位取样响应为 试讨论其是否是因果的、稳定的。 解： 因果性: 该系统是非因果系统。 稳定性: 当 时系统稳定,当 时系统不稳定。天津师范大学计算机与信息工程学院例1.2.5 设系统输入输出关系为，判断其线性，移不变性，因果性和稳定性。解：因而所以此系统为线性系统而因而所以此系统不是移不变系统，也就是系统是移变的。天津师范大学计算机与信息工程学院若x(n)有界，即 ，则而 ，所以 。即有界的输入产生有界的输出，因此系统是稳定的。只与x(n)的当前值有关，而与未来值无关，所以系统是因果的。天津师范大学计算机
37、与信息工程学院.常系数线性差分方程 连续线性时不变系统的输入输出关系常用常系数线性微分方程表示，而离散线性移不变系统的输入输出关系常用常系数线性差分方程表示，即 或者常系数是指决定系统特征的系数是常数，若系数中含有n，则称为“变系数”。差分方程的阶数等于y(n)的变量序号的最高值与最低值之差，例如上式就是N阶差分方程。线性是指各y(n-i)项和各x(n-i)项都只有一次幂而且不存在它们的相乘项，否则就是非线性。天津师范大学计算机与信息工程学院求解差分方程有如下几种方法：递推法、时域经典法、卷积法、变换域法等等递推解法比较简单，适合计算机求解，但是只能得到数值解，不易直接得到闭合形式（公式）解答
38、。时域经典法和微分方程的解法比较类似，比较麻烦，实际应用中很少采用。卷积法则必须知道系统的单位抽样响应h(n) ，这样利用卷积和就能得到任意输入时的输出响应。变换域法是利用Z变换的方法求解差分方程。 当系统的初始状态为零，单位抽样响应h(n)就能完全代表系统，那么对于线性移不变系统，任意输入下的系统输出就可以利用卷积和求得。差分方程在给定输入和边界条件下，可用迭代的方法求系统的响应，当输入为(n)时，输出（响应）就是单位抽样响应h(n)。天津师范大学计算机与信息工程学院例1.3.1 常系数差分方程（）初始条件为n0时，y(n)=0，求其单位抽样响应；（）初始条件为n0时，y(n)=0，求其单位
39、抽样响应。解：（）设 ，且，必有依次迭代所以单位抽样响应为天津师范大学计算机与信息工程学院（）设，由初始条件知，必有将原式该写为另一种递推关系则所以单位抽样响应为由本例看出，差分方程相同，但是初始条件不同，得到的单位抽样响应不同，也就是对应着不同的系统天津师范大学计算机与信息工程学院1.4 连续时间信号的采样1.4.1 采样的基本概念 从原理上说，采样器就是一个开关， 通过控制开关的接通和断开来实现信号的 采样，它的概念如图1.3所示。 天津师范大学计算机与信息工程学院图1.3 采样过程天津师范大学计算机与信息工程学院 采样在数学上等效为下列运算： 式中s(t)是一个开关函数， 是原信号， 是
40、采样后的信号理想采样情况下，s(t)是无 限多项单位冲击信号 等间隔构成的一个单 位冲击串，即 式中T是采样间隔 。天津师范大学计算机与信息工程学院则式中， 只在 时不为零，因而 只在这些点上才有定义的值，为 ，可见采样的结果是使原来的模拟信号变成为在 这些点上的离散信号，这就是采样的简单原理。天津师范大学计算机与信息工程学院1.4.2 采样过程中频谱的变化 周期信号T(t)可以进行傅里叶级数展开，如下式： 可以求解出 式中， ，是 的基波频率，同时也是采样频率。 天津师范大学计算机与信息工程学院令 ，求得Ak为T(t)等效为 因此天津师范大学计算机与信息工程学院 上式表示 是无限多个载波 被
41、 调制之和，从频域变化来看， 的频谱被搬移到无限多个频率点，这些频率点是 ，所以的频谱就变成了周期 函数，周期等于 。所以 天津师范大学计算机与信息工程学院 分析上式，xS(t)与xa(t)的频谱比较，主 要的变化是：它的频谱变成了周期的，即 是周期函数，周期为 ，也就是说， 离散时间信号的频谱是连续时间信号频谱以 采样频率为周期进行无限项周期延拓的结 果，这是信号采样带来的最重要的变化。另 一点变化是频谱幅度变为原来幅度的1/T。 图1.4表示了这种频谱的变化。天津师范大学计算机与信息工程学院图1.4 理想采样信号的频谱天津师范大学计算机与信息工程学院1.4.3 低通信号采样定理 设xa(t
42、)表示一个带限的低通模拟信号， 最高频率分量为fmax，它的频谱为Xa(j)， 如图1.5所示。图1.5 带限的低通模拟信号天津师范大学计算机与信息工程学院 对该信号以采样频率fs进行采样，采样后的离散时间信号的频谱Xs(j)变成了以fs为周期的周期频谱。显然，在这种情况下，Xs(j)和Xa(j)包含的信息是相同的，或者说，采样后的离散信号能完全表示原来的模拟信号。 若fs2fmax，这时周期频谱的各周期出现了混叠，造成实际的周期频谱的一个周期不等于原信号的频谱，也就是说，采样以后，信号出现了失真。天津师范大学计算机与信息工程学院 |Xa(j)| |Xs(j)| |T(j)| |Xs(j)|S
43、S-S-SS/2-S/2S-S(a)(b)(c)(d)-hhh天津师范大学计算机与信息工程学院 采样定理 对一个低通带限信号进行均匀理想采样，如果采样频率大于等于信号最高频率的两倍，采样后的信号可以精确地重建原信号。可以表示为 或 式中 ， 是信号的最高频率。 当 时的采样频率为临界采样频率或 称为“奈奎斯特采样频率”。天津师范大学计算机与信息工程学院 例1.4.1 模拟信号 其中 。 (1) 求 的周期,采样频率应为多少?采样间隔应为多少? (2) 若选采样频率 ,采样间隔为多少?写出采样信号 的表达式。天津师范大学计算机与信息工程学院解 (1) 由 ,得 的周期为: 采样频率应为: 采样间
44、隔应为: (2) 选 则采样间隔为天津师范大学计算机与信息工程学院所以天津师范大学计算机与信息工程学院1.4.4 信号恢复 当满足采样定理的条件时，可以推导出从离散时间信号恢复原来的模拟信号的内插公式。 先从频域分析入手，已知采样后的信号的频谱在一个周期里可以表示为 因此，只要设计一个截止频率为 的理想低通滤波器，就可以恢复原信号的频谱 。 天津师范大学计算机与信息工程学院设该理想低通滤波器的频率响应为根据模拟系统的频域描述理论，有所以， 将等于原信号的频谱 。 T 0 天津师范大学计算机与信息工程学院 设滤波器的单位冲击响应为天津师范大学计算机与信息工程学院 滤波器的输出为这就是信号恢复的内
45、插公式。其中 称为内插函数，它是一个关于t的连续函数，关于k的离散函数。天津师范大学计算机与信息工程学院 有一个重要的特点： 在取样点n = k时， =1, 在其他取样点nk， =0， 在当前的取样点上， 和 完全相等，在取样点之间， 是由无限多个 被相应的采样值 作加权系数的线性组合构成，图1.6是这种过程的示意图。 在每一个抽样点上，由于只有该抽样值所对应的内插函数不为零，抽样内插公式保证了各采样点上信号值不变，而采样点之间的信号值则是由各抽样值对应的内插函数的波形延伸叠加而成 天津师范大学计算机与信息工程学院图1.6 采样信号的恢复天津师范大学计算机与信息工程学院1.4.5 窄带信号采样
46、定理 所谓窄带信号就是信号带宽远远小于它的中心频率的信号。 设窄带信号的数学模型为 其中， 是低频信号，其最高频率远远小于 ，它们通常携带有信息，分别被调制在频率为 的载波的幅度和相位上。 天津师范大学计算机与信息工程学院 当信号 的最高频率是带宽的整数倍， 即 结论 :对窄带信号 ，当上限频率正好是带宽的整数倍时，采样频率 只要等于2倍的带宽频率 ，可以由 重 建 。( k正整数) 天津师范大学计算机与信息工程学院对一般情况，即 不是带宽 的整数倍时，令此时， 不是整数，在这种情况下，可以考虑在保持 不变的情况下，增加带宽 ，使之为 ，使得结论:如果 不是 的整数倍，若想重建 ，应增加采样频
47、率为 ， 是原 的m倍, 。 天津师范大学计算机与信息工程学院9、静夜四无，荒居旧 。三月-21三月-21Monday, March 22, 202110、雨中黄叶，灯下白人。17:34:4617:34:4617:343/22/2021 5:34:46 PM11、以我独沈久，愧君相 。三月-2117:34:4617:34Mar-2122-Mar-2112、故人江海，几度隔山川。17:34:4617:34:4617:34Monday, March 22, 202113、乍翻疑梦，相悲各年。三月-21三月-2117:34:4617:34:46March 22, 202114、他生白，旧国青山。22
48、 三月 20215:34:46 下午17:34:46三月-2115、比不了得就不比，得不到的就不要。三月 215:34 下午三月-2117:34March 22, 202116、行出成果，工作出富。2021/3/22 17:34:4617:34:4622 March 202117、做前，能 四周；做，你只能或者最好沿着以脚起点的射向前。5:34:46 下午5:34 下午17:34:46三月-219、没有失，只有 停止成功！。三月-21三月-21Monday, March 22, 202110、很多事情努力了未必有果，但是不努力却什么改也没有。17:34:4617:34:4617:343/22/
49、2021 5:34:46 PM11、成功就是日复一日那一点点小小努力的累。三月-2117:34:4617:34Mar-2122-Mar-2112、世成事，不求其 ，留一份不足，可得无限完美。17:34:4617:34:4617:34Monday, March 22, 202113、不知香寺，数里入云峰。三月-21三月-2117:34:4617:34:46March 22, 202114、意志强的人能把世界放在手中像泥一任意揉捏。22 三月 20215:34:46 下午17:34:46三月-2115、楚塞三湘接， 九派通。三月 215:34 下午三月-2117:34March 22, 20211
50、6、少年十五二十，步行得胡 。2021/3/22 17:34:4617:34:4622 March 202117、空山新雨后，天气晚来秋。5:34:46 下午5:34 下午17:34:46三月-219、柳散和，青山澹吾。三月-21三月-21Monday, March 22, 202110、 一切好如同和去最杰出的人 。17:34:4617:34:4617:343/22/2021 5:34:46 PM11、越是没有本的就越加自命不凡。三月-2117:34:4617:34Mar-2122-Mar-2112、越是无能的人，越喜挑剔人的儿。17:34:4617:34:4617:34Monday, Ma
