第1章 时域离散时间信号和时域离散系统

第1章时域离散信号和时域离散系统本章主要内容时域离散信号的基本概念及典型序列时域离散系统的定义及其性质线性时不变系统的输入/输出δ(n
)求解法时域离散系统的输入输出法：线性常系数差分方程模拟信号数字处理方法Matlab实现1.1引言信号的分类系统的分类信号的分类时
域连续信号（模拟信号）：信号的自变量和函数值都取连续值，例如语言信号、温度信号等；时域离散信号：如果自变量取离散值，而函数值取连续
值，这种信号通常来源于对模拟信号的采样；数字信号：信号的自变量和函数值均取离散值。模拟信号采样间隔T=0.005s进行等间隔采样
，得时域离散信号x(n)，={…，0.0，0.6364，0.9，0.6364，0.0，-0.6364，-0.9，-0.63
64，…}显然,时域离散信号是时间离散化的模拟信号。如果用四位二进制数表示该时域离散信号，得到相应的数字信号x［n］={…，
0.000，0.101，0.111，0.101，0.000，1.101，1.111，1.101，…}数字信号是幅度、时间均离散化
的模拟信号，或者说是幅度离散化的时域离散信号。系统的分类模拟系统时域离散系统数字系统模拟网络和数字网络构成的混合系统1.2
时域离散信号—序列序列的定义及表示常用的典型序列序列的周期性用单位脉冲序列表示任意序列序列的基本运算1.2.1序列的定义及
表示序列的定义数字序列：离散时间信号{-2,5,-6,8,3,-7}一般只在均匀间隔的离散时间nT上给出数值{…,x(
-2T),X(-1T),X(0),X(T),X(2T),…}序列的表示用集合符号表示用公式表示用图形表示序列表示用集合符号
表示x(n)={x(n)}，-∞＜n＜+∞x(n)={……,x(-2),x(-1),x(0),x(1),x(2),….
}n代表nTT采样时间间隔nT指均匀间隔的离散时间点n为非整数时没有定义，不能认为此时x(n)的值是零用公式表示序列
表示用图形表示1.2.1常用的典型序列单位脉冲序列单位阶跃序列矩形序列实指数序列正弦序列复指数序列周期序列任意序列表示单
位脉冲序列δ(n-m)只有在n=m时取确定值1，而其余点取值均为零δ(n)只在n=0时取确定值1，其它均为零δ(n)类似
于δ(t)，注意二者的定义与区别单位阶跃序列u(n)类似于u(t)u(t)在t=0时常不定义u(n)在n=0时为u(0)=1
δ(n)和u(n)的关系：δ(n)=u(n)-u(n-1)矩形序列N为矩形序列的长度和u(n)、δ(n)的关系：
实指数序列a为实数当|a|＞1时序列发散当|a|＜1时序列收敛正弦序列A为幅度ω为数字域频率φ为起始相位x(n)=Asin
(ωn+φ)设x(n)由x(t)=sinΩt取样得到（A、φ与频率无关不考虑）ω=ΩT=Ω/fs,ω与Ω线性关
系,ω的单位为rad复指数序列ω为数字域频率用实部与虚部表示用极坐标表示只考虑频率令σ=0，序列频率ω呈现以2π为周期的
周期性后续研究中频率域只考虑或就够了周期序列对于序列x(n)，如果对所有n存
在一个最小的正整数N，对任意整数m满足x(n)=x(n+mN)则序列x(n)是周期序列，最小周期为N。以正弦
序列为例讨论周期性设x(n)=Asin(ωn+φ)则有x(n+N)=Asin[ω(n+N)+φ]=Asin(ωN+
ωn+φ)若满足条件ωN=2kπ，则x(n+N)=Asin[ω(n+N)+φ]=Asin(ωn+φ)=x(n)周期
序列N、k为整数，k的取值满足条件，且保证N最小正整数。其周期为2π/ω为整数时，取k=1，保证为最小正整数。此时为周
期序列，周期为2π/ω。例1.4序列，因为2π/ω=8，所以是
一个周期序列，其周期N=8。周期序列2π/ω为有理数而非整数时，仍然是周期序列，周期大于2π/ω。例1.5序列
，2π/ω=8/3是有理数，所以是周期序列，取k=3，得到周期N=8。2π
/ω为无理数时，任何k都不能使N为正整数，这时正弦序列不是周期序列。例序列指数为纯虚数的复指数序列的周期性与正弦序列的情况
相同。用单位脉冲序列表示任意序列任何序列都可以用单位脉冲序列的移位加权和来表示，即x(n)可看成是x(n)和δ(n)的卷积
和，式中例1.61.2.2序列的基本运算和积移位标乘翻转累加差分时间尺度变换序列能量卷积和基本运算—序列的和设序列为x(n
)和y(n)，则序列z(n)=x(n)+y(n)表示两个序列的和，定义为同序号的序列值逐项对应相加。例：序列的和例：设序
列计算序列的和x(n)+y(n)。解：例：序列求和图示基本运算—序列的积设序列为x(n)和y(n)，则序列z(n)=x(n
)?y(n)表示两个序列的积，定义为同序号的序列值逐项对应相乘。例：序列的积例：设序列计算序列的和x(n)?y(n)。
解：例：序列求积图示x(n)基本运算—序列的移位设序列为x(n)，则序列y(n)=x(n-m)表示将序列x(n)进行移位。
m为正时x(n-m)：x(n)逐项依次延时(右移)m位x(n+m)：x(n)逐项依次超前(左移)m位m为负时，则相反。例：序列
的移位例：设序列计算序列的和x(n+1)。解：例：序列移位图示x(n)基本运算—序列的标乘设序列为x(n)，a为常数(a≠0
)，则序列y(n)=ax(n)表示将序列x(n)的标乘，定义为各序列值均乘以a，使新序列的幅度为原序列的a倍。例：序列的
标乘例：设序列计算序列4x(n)。解：基本运算—序列的翻转设序列为x(n)，则序列y(n)=x(-n)表示以n=0的纵
轴为对称轴将序列x(n)加以翻转。例：序列的翻转例：设序列计算序列x(-n)。解：基本运算—序列的累加设序列为x(n)，则序列
定义为对x(n)的累加，表示将n以前的所有x(n)值求和。基本运算—序列的差分前向差分：将序列先进行左移，再相减Δx(n)
=x(n+1)-x(n)后向差分：将序列先进行右移，再相减▽x(n)=x(n)-x(n-1)由此，容易得出▽x(
n)=Δx(n-1)基本运算—时间尺度(比例)变换设序列为x(n)，m为正整数，则序列抽取序列y(n)=x(mn)插
值序列x(mn)和x(n/m)定义为对x(n)的时间尺度变换。插值序列x(n/m)：对x(n)进行零值内插运算表示在原
序列x(n)相邻两点之间插入m-1个零值点保留x(0)基本运算—序列的能量设序列为x(n)，则序列定义为序列的能量，表示
序列各取样值的平方之和；若为复序列，取模值后再求平方和。1.3时域离散系统时域离散系统的定义及表示线性时不变系统线性时
不变系统h(n)与I/O关系线性时不变系统的性质系统的因果性和稳定性时域离散系统的定义及表示时域离散系统定义为将输入序列x(
n)映射成输出序列y(n)的惟一变换或运算。以T[·]表示这种运算y(n)=T[x(n)]对变换T[·]加以不同的约束条件，
所定义的系统就具有不同的特性和功能。线性时不变系统:最重要、最常用，可表征许多物理过程。1.3.1、1.3.2线性时不变系统线
性系统满足叠加原理叠加原理包含可加性和齐次性两方面性质时不变系统系统的响应与输入信号施加于系统的时刻无关运算关系在整个运算过程中
不随时间而变化线性时不变系统既满足叠加原理，又满足时不变性的系统线性系统设系统的输入序列与输出分别为可加性:如果系统
的输入之和与输出之和满足齐次性(或比例性):设a为常数，系统的输入增大a倍，输出也增大a倍线性系统、非线性系统（不满足可加性与齐
次性）例：证明一个线性系统注意：必须证明系统同时满足可加性和齐次性，且信号及比例常数都可以是复数。例：试分析下列系统是否是线性系
统(1)y(n)=2x(n)-3，(2)y(n)=x(Mn)，其中M为正整数。满足叠加原理，线性系统不满足叠加原理，非
线性系统时不变系统输入序列x(n)移动任意m位后，输出序列y(n)也移动m位，数值却保持不变。m为任意常整数时不
变系统也称为移不变系统例：证明一个时不变系统例：试分析下列系统的时不变性(1)y(n)=2x(n)-3，(2)y(n)
=x(Mn)，其中M为正整数。时变系统二者相等，具有时不变性1.3.3线性时不变系统h(n)与I/O关系时域离散系统：
完全响应=零输入响应+零状态响应单位脉冲响应(单位取样响应)h(n)=T[δ(n)]线性时不变系统输入为δ(n)时的零状态
响应线性时不变系统特性都可以用它的单位脉冲响应h(n)来表征已知h(n)可得到线性时不变系统对任意输入的输出I/O关系推导
用δ(n)表示x(n)系统输出叠加原理时不变性I/O关系:线性时不变系统的输出等于输入序列和单位脉冲响应h(n)的卷积。
线性时不变系统的性质交换律结合律分配律可以推广到多个系统的情况，由卷积和的定义可以很容易加以证明。序列的卷积和设序列为x
(n)和z(n)，则序列定义为x(n)和z(n)的卷积和。卷积和又称为离散卷积或线性卷积，是很重要的公式。线性时不变系统的I
/O关系：就是序列卷积和的运算！卷积和计算的四个步骤(1)翻转：x(m)，x(m)→x(-m)(2)移位：z(m)
→z(n-m)n为正数时，右移n位n为负数时，左移n位(3)相乘：x(m)z(n-m),（m值相同）(4)相加
：y(n)=∑x(m)z(n-m)例：卷积和计算例设序列求y(n)=x(n)z(n)。对应点相乘！解：n＜0时，x(m
)与z(n-m)没有重叠，得y(n)=0。对应点相乘！0≤n≤4时，例：卷积和计算4＜n≤6时，6＜n≤10时，n＞10
时，x(m)与z(n-m)没有重叠，得y(n)=0。线性时不变系统的I/O求解例：已知x(n)=R4(n),h(n)
=R4(n),求y(n)=x(n)h(n)。解计算卷积的基本运算是翻转、移位、相乘和相加。计算方法：图解法、解析法图解
法首先将h(n)用h(m)表示，并将波形翻转，得到h(－m)，然后将h(－m)移位n，得到h(n－m)，n>0,序列右移；n
<0，序列左移。如n=1，得到h(1-m)，接着将h(m)和h(n－m)相乘后，再相加,得到y(n)的一个值。对所有的n重复这种
计算,最后得到卷积结果，如图1.3.2(f)所示,y(n)表达式为y(n)={1,2,3,4,3,2,1}例
1.3线性卷积表1.3图解法（列表法）解析法例：设x(n)=anu(n)，h(n)=R4(n)，求y(n)=x(
n)h(n)。解关键：根据求和号内的两个信号乘积的非零值区间确定求和的上、下限。因为n≥m时，u(n-m)才能取非零值；0
≤m≤3时，R4(m)取非零值；所以，求和区间中m要同时满足下面两式：m≤n0≤m≤3这样求和限与n有关系，必须将n进行
分段然后计算。n<0时，y(n)=00≤n≤3时，乘积的非零值范围为0≤m≤n，因此n≥4时，乘积的非零区间为0≤m≤3，因此写
成统一表达式1.3.4系统的因果性和稳定性一般因果系统定义系统某时刻的输出y(n)只取决于此时刻x(n)和以前的输入x(
n-1)，x(n-2)，…，而和此时刻以后的输入x(n+1)，x(n+2)，…无关。先因后果因果系统的响应不会出现于外加输入之前
。非因果系统当前的输出还取决于未来的输入，不符合因果关系。因果性的充分必要条件线性时不变系统具有因果性的充要条件h(n)=
0，n＜0证明充分条件若n＜0时，h(n)=0，则m>n时，h(n)=0因而n0时刻的输出可见，y(n0)只与
m≤n0时的x(m)有关，因而是因果系统。因果条件证明证明利用反证法证明必要条件假设因果系统，n＜0时h(n)≠0，则在所
设条件下，第二个求和式中至少有一项不为零，y(n)将至少和m＞n时的某一个x(n)值有关，这不符合因果性，假设不成立。例：判断因果
系统例：判断差分系统的因果性。(1)前向差分系统:y(n)=x(n+1)-x(n);(2)后向差分系统:y(n)=
x(n)-x(n-1)。解因为前向差分系统的y(n)决定于x(n+1)，故系统为非因果的。而后向差分系统定义为y(n)=x
(n)-x(n-1)，显然是因果的。稳定系统一般稳定系统定义系统的每个有界输入，对应产生的输出都有界。如果输入满足|x(n)|
≤M＜+∞(M为正常数)，有输出|y(n)|≤P＜+∞(P为正常数)。判断系统不稳定只要找出一个特别的有界输入，对应的输出是
无界的，则该系统就是不稳定的。判断系统稳定必须证明所有有界输入，其输出都是有界的。稳定性的充分必要条件线性时不变系统具有稳定性
的充要条件是其单位脉冲响应绝对可和，即证明充分条件若式成立，对于所有n都有|x(n)|≤M，得即输出y(n)有界，系统稳定
。稳定条件证明假设系统稳定，但单位脉冲响应不绝对可和证明利用反证法证明必要条件定义一个有界输入计算输出，令n=0则有即y(
0)无界，系统不稳定，因此假设不成立。例：判断稳定系统例：判断累加器系统的稳定性解考虑有界输入x(n)=u(n)，累加器的
输出为虽然n为有限值时，系统输出也为有限值，但对于所有n值(包括+∞)不存在有限值P，使得(n+1)≤P＜+∞，故系统输出无界。
系统不稳定例：判断因果稳定系统例：已知线性时不变系统的单位脉冲响应式中a为实常数，讨论其因果性和稳定性。解因为n＜0
时，u(-n-1)=1，所以h(n)≠0，故系统是非因果系统。因为所以|a|＞1时系统稳定，|a|≤1时不稳定。收敛序列：如
|a|＞1时，h(n)模值随n加大而减小发散序列：如|a|≤1时，h(n)模值随n加大而加大1.4时域离散系统的输入输出描述法
——线性常系数差分方程描述一个系统可以不管系统内部的结构如何，将系统看成一个黑盒子，只描述系统的输出与输入之间的关系，这种描述法被
称为输入输出描述法。微分方程模拟系统差分方程时域离散系统状态变量描述法线性时不变系统线
性常系数差分方程T[·]y(n)x(n)时域离散系统用方程来描述两种不同的描述方法返回1.4.1线性常系数差分方程一个N阶线
性常系数差分方程用下式描述：或，a0=1式中
，x(n)和y(n)分别表示系统的输入和输出，系数ai和bi均为常系数，且x(n-i)和y(n-i)只有次幂，也没有相互交叉的线性
相乘项，故称为线性常系数差分方程。1.4.2线性常系数差分方程的求解已知系统的输入信号和描述系统的线性常系数差分方程，求解系统
的输出一般有三种方法：经典解法:和求解微分方程解法类似，齐次解＋特解递推解法：由初始值和输入值递推解出系统以后输出值Z变换解法：
适合计算机求解线性常系数差分方程的递推解法递推解法：观察上式，如果已知输入信号x(n)，求n时刻的输出，需要知道输入信号x(n
)，以及n时刻以前的N个输出信号值：y(n-1)，y(n-2)，y(n-3)，…，y(n-N)。这N个输出信号值就构成初始条件。可
以看到，上式是一个递推方程。如果已知输入信号x(n)和N个初始条件，就可以求出n个时刻的输出；如果将这公式中的n用n+1代替，就可
求出n+1时刻的输出，依此类推，可求出各个时刻的输出。递推法求差分方程【例2.14】设系统用差分方程y(n)=ay(n－1)+x
(n)描述，输入序列x(n)=δ(n)，求输出序列y(n)。解：该系统差分方程是一阶差分方程，需要一个初始条件。（2）设初始
条件:（1）设初始条件:1.5模拟信号数字处理方法1.5.1采样定理及A/D变换式中δ(t)是单位冲激信号，在上式
中只有当t=nT时，才可能有非零值，因此写成下式：对进行傅里叶
变换，得到式中，Ωs=2π/T，称为采样角频率，单位是rad/s理想采样信号的频谱是原模拟信号的频谱沿频率轴，每间隔采样角频率Ω
s重复出现一次，或者说理想采样信号的频谱是原模拟信号的频谱以Ωs为周期，进行周期性延拓而成的。图1.5.3采样信号的频谱
采样恢复图1.5.4采样恢复采样恢复设xa(t)是带限信号，最高频率为Ωc，其频谱Xa(jΩ)如图1.5.3(a)所示。p
δ(t)的频谱Pδ(jΩ)如图1.5.3(b)所示，那么按照（1.5.5）式，的频谱如图1.5.3(c)所示，图中原模
拟信号的频谱称为基带频谱。如果满足Ωs≥2Ωc，或者用频率表示该式，即满足Fs≥2fc，基带谱与其它周期延拓形成的谱不重叠，如图1
.5.3(c)所示情况，可以用理想低通滤波器G(jΩ)从采样信号中不失真地提取原模拟信号，如图1.5.4所示。但如果选择采样频率太
低，或者说信号最高截止频率过高，使Fs<2fc,Xa(jΩ)按照采样频率Fs周期延拓时，形成频谱混叠现象，用图1.5.3(d)表
示。这种情况下，再用图1.5.4所示的理想低通滤波器对Xa(t)进行滤波，得到的是失真了的模拟信号。折叠频率Fs/2这里需要
说明的是，一般频谱函数是复函数，相加应是复数相加，图1.5.3和图1.5.4仅是示意图。一般称Fs/2为折叠频率，只有当信号最高频
率不超过Fs/2时，才不会产生频率混叠现象，否则超过Fs/2的频谱会折叠回来而形成混叠现象，因此频率混叠在Fs/2附近最严重。
采样定理（1）对连续信号进行等间隔采样形成采样信号，采样信号的频谱是原连续信号的频谱以采样频率Ωs为周期进行周期性的延拓形成的
，用公式（1.5.5）表示。（2）设连续信号xa(t)属带限信号，最高截止频率为Ωc，如果采样角频率Ωs≥2Ωc，那么让采样
信号通过一个增益为T、截止频率为Ωs/2的理想低通滤波器，可以唯一地恢复出原连续信号xa(t)。否则,Ωs<2Ωc会造
成采样信号中的频谱混叠现象，不可能无失真地恢复原连续信号。实际中对模拟信号进行采样，需根据模拟信号的截止频率，按照采样定理的要求
选择采样频率，即Ωs≥2Ωc，但考虑到理想滤波器G（jΩ)不可实现，要有一定的过渡带，为此可选Ωs=(3~4)Ωc。另外，可以在采
样之前加一抗混叠的低通滤波器，滤去高于Ωs/2的一些无用的高频分量，以及滤除其它的一些杂散信号。这就是在图1.5.1中采样之前加预
滤的原因。A/DC模/数转换器模/数转换器原理框图例如：模拟信号xa(t)=sin(2πft+π/8)，式中f=50
Hz,选采样频率Fs=200Hz，将t=nT代入xa(t)中，得到采样数据：当时，得到序列x(n)如下：x(n)={
,0.382683,0.923879,－0.382683,－0.923879,}如果A/D
C按照M=6进行量化编码，即上面的采样数据均用6位二进制码表示，其中一位为符号位，则数字信号用表示：={,0
.01100,0.11101,1.01100,1.11101,}用十进制数表示的为={
,0.37500,0.90625,-0.37500,-0.90625,}显然量化编码以后的
和原x(n)不同。这样产生的误差称为量化误差，这种量化误差的影响称为量化效应，这部分内容将在第9章介绍。1.5.2将数
字信号转换成模拟信号我们已经知道模拟信号xa(t)经过理想采样，得到采样信号,xa(t)和之间的关系用（1.5.2）式描述
。如果选择采样频率Fs满足采样定理，的频谱没有频谱混叠现象，可用一个传输函数为G（jω)的理想低通滤波器不失真地将原模拟信号x
a(t)恢复出来，这是一种理想恢复。下面先分析推导该理想低通滤波器的输入和输出之间的关系，以便了解理想低通滤波器是如何由采样信号恢
复原模拟信号的，然后再介绍在实际中数字信号如何转换成模拟信号。由低通滤波器的传输函数G(jΩ）推导其单位冲激响应g(t)：因为
Ωs=2πFs=2π/T，因此g(t)也可以用下式表示：理想低想滤波器的输入、输出分别为和ya(t)，将（1.5.7）式
表示的g(t)和（1.5.2）式表示的代入上式，得到：（1.5.8）理想恢复由于满足采样定理，ya(t)=xa(t)，因
此得到：（1.5.9）式中，当n=,－1,0,1,2,时，xa(nT)是一串离散的采样值，而xa(t)是模拟信号，
t取连续值，g(t)的波形如图1.5.6所示。其特点是:t=0时，g(0)=1;t=nT(n≠0)时，g(t)=0。在（1
.5.9）式中，g(t)保证了在各个采样点上，即t=nT时，恢复的xa(t)等于原采样值，而在采样点之间，则是各采样值乘以g(t－
nT)的波形伸展叠加而成的。图1.5.6内插函数g(t)波形理想恢复理想恢复这种伸展波形叠加的情况如图1.5.7所示。g
(t)函数所起的作用是在各采样点之间内插，因此称为内插函数，而（1.5.9）式则称为内插公式。这种用理想低通滤波器恢复的模拟信号完
全等于原模拟信号xa(t)，是一种无失真的恢复。但由于g(t)是非因果的，因此理想低通滤波器是非因果不可实现的。下面介绍实际的数字
信号到模拟信号的转换。图1.5.7理想恢复图1.5.8D/AC方框图D/AC数模转换器实际中采用D/AC（Digi
tal/AnalogConverter）完成数字信号到模拟信号的转换。D/AC包括三部分，即解码器、零阶保持器和平滑滤波器，D/
AC方框图如图1.5.8所示。解码器的作用是将数字信号转换成时域离散信号xa(nT)，零阶保持器和平滑滤波器则将xa(nT)变成模
拟信号。由时域离散信号xa(nT)恢复模拟信号的过程是在采样点内插的过程。理想低通滤波的方法是用g(t)函数作内插函数，还可以
用一阶线性函数作内插。零阶保持器是将前一个采样值进行保持，一直到下一个采样值来到，再跳到新的采样值并保持，因此相当于进行常数内插。
零阶保持器的单位冲激函数h(t)以及输出波形如图1.5.9所示。对h(t)进行傅里叶变换，得到其传输函数：（1.5.10）图1
.5.9零阶保持器的输出波形其幅度特性和相位特性如图1.5.10所示。由该图看到,零阶保持器是一个低通滤波器，能够起到将时域
离散信号恢复成模拟信号的作用。图中虚线表示理想低通滤波器的幅度特性。零阶保持器的幅度特性与其有明显的差别，主要是在|Ω|>π/T区
域有较多的高频分量，表现在时域上，就是恢复出的模拟信号是台阶形的。因此需要在D/AC之后加平滑低通滤波器，滤除多余的高频分量，对时
间波形起平滑作用，这也就是在图1.5.1模拟信号数字处理框中，最后加平滑滤波器的原因。虽然这种零阶保持器恢复的模拟信号有些失真，但
简单、易实现，是经常使用的方法。实际中，将解码器与零阶保持器集成在一起，就是工程上的D/AC器件。图1.5.10零阶保持器的
频率特性1.6Matlab实现单位脉冲序列单位阶跃序列矩形序列实指数序列正弦序列复指数序列常用序列的Matlab实现序
列运算的Matlab实现Matlab求解离散系统的差分方程翻转序列的能量卷积和单位脉冲序列δ(n-1)n=[-3:3];
%生成位置向量x=[(n-1)==0];
%生成单个脉冲序列stem(n,x);axis([-3,3,0,1.5]);
%标示坐标单位阶跃序列u(n+1)n=[-3:3];%生成位置向量x=[(n+1)>=0]
;%生成阶跃序列stem(n,x);axis([-3,3,0,1.5]);矩形序列生成函数function[x,n]=
rectseq(n0,n1,n2,N)%单位矩形序列生成函数%调用方式[x,n]=rectseq(n0,n1,n2,N)
n=[n0:n2];%生成位置向量x=[(n-n1)>=0&((n1+N-1)-n)>=0];%
生成矩形脉冲序列矩形序列[x,n]=rectseq(-3,-1,4,5);stem(n,x);axis([-3,5,0,1
.5]);实指数序列n=[0:10];%生成位置向量x=(0.6).^n;
%生成实指数序列stem(n,x);axis([0,10,0,1.5]);正弦序列3sin(0.1
πn+π/3)n=[0:1:20];%生成位置向量x=3si
n(0.1pin+pi/3);%生成正弦序列stem(n,x);axis([0,20,-4,4]);复指
数序列n=[-2:10];x=exp((0.2-0.5j)n);%复
指数序列subplot(1,2,1),stem(n,real(x));%用空心圆画点line([-5,10],[0,
0]);%画横坐标subplot(1,2,2),stem(n,imag(
x),''filled'');%用实心圆画点%line([-5,10],[0,0])翻转:调用fliplrn=[-3:3]
;%生成一个序列x=[0,0,1,0.5,0.25,0.125,0];stem(n,x);x
=fliplr(x);%x排列次序左右翻转n=-fliplr(n);%向量n对
n=0翻转stem(n,x);序列的能量conj求共轭复数sum求总和E=sum(x.conj(x));abs求
幅值sum求总和E=sum(abs(x).^2);卷积和：调用convx=[3,-3,7,0,-1,5,2];%
序列x的非零区间-4≤n≤2h=[2,3,0,-5,2,1];%序列x的非零区间-1≤n≤4%调用conv计算卷积和y
=conv(x,h);运行结果：无位置信息y=635619-313018-27-19
2卷积和函数：convextd.mfunction[y,ny]=convextd(x,nx,h,nh)%序列y为序列
x和序列h的卷积%ny，nx，nh分别为序列y，x和h的位置向量%调用方式[y,ny]=convextd(x,nx,h
,nh)ny1=nx(1)+nh(1);%计算卷积后的起点位置ny_end=nx(end)
+nh(end);%计算卷积后的终点位置y=conv(x,h);%计算卷积和
序列的数值ny=[ny1:ny_end];%计算卷积和序列的位置向量卷积和：包含位置向量x=[3,-3,7,0,-1,5,2];nx=[-4:2];%给定输入序列h=[2,3,0,-5,2,1];nh=[-1:4];%给定脉冲响应序列[y,ny]=convextd(x,nx,h,nh);%带位置序列的卷积结果运行结果：有位置信息y=635619-313018-27-192ny=-5-4-3-2-10123456解差分方程：调用filter函数的调用方式为y=filter(b,a,x);输入参数b、a为差分方程的系数，b=[b0，b1，…，bM]a=[a0，a1，…，aN]输入参数x是输入序列求得的输出序列y和输入x的长度一样系数a0必须不为零。例：解差分方程例1.15线性常系数差分方程y(n)-y(n-1)+0.75y(n-2)=x(n)，求输入x(n)=δ(n)时系统的输出序列。(1)求单位脉冲响应h(n)b=1;a=[1,-1,0.75];x=impseq(-10,0,50);%生成单位脉冲序列h=filter(b,a,x);%计算单位脉冲响应n=[-10:50];stem(n,h);%脉冲响应曲线axis([-10,50,-1,1.5])%标出坐标title(''ImpulseResponse'');xlabel(''n'');ylabel(''h(n)'');例：判断系统稳定(2)求得单位脉冲响应的和sum(abs(h));%计算单位脉冲响应的和程序的运行结果为ans=6.1718绝对可和，说明系统是稳定的。