1.1 离散时间信号——序列
1.1.1 序列
离散时间信号是时间上不连续的序列,用x(n)来表示。x(n)可以看成对模拟信号xa(t)的等间隔时间抽样,即
x(n)=xa(t)t=nT=xa(nT) 序列的三种表示方法:
函数表示法,例如x(n)=anu(n)
数列表示法,例如x(n)=...,−5,−3,−1,0,2,7,9,n=0时,x(0)的值要用下划线标注
1.1.2 序列的运算
3个基本运算单元:加法器、乘法器和延时单位。
序列的三类运算:
对幅度进行运算
加法:z(n)=x(n)+y(n)
乘法:w(n)=x(n)y(n),当x(n)或y(n)是常数c时,称为标度运算,即w(n)=cx(n)
累加:y(n)=∑k=−∞nx(k)
序列的绝对可和:S=∑n=−∞∞∣x(n)∣,当S=B<∞时,称序列x(n)为绝对可和序列
序列的能量:E[x(n)]=∑n=−∞∞∣x(n)∣2,若E[x(n)]=A<∞,则称x(n)为能量信号有限信号,简称能量信号
序列的平均功率:P[x(n)]=limN→∞2N+11∑n=−NN∣x(n)∣2,若极限存在,即P[x(n)]=C<∞,则称x(n)为功率有限信号,简称功率信号,对于周期信号,只需取一个周期N的平均功率即可,即P[x(n)]=N1∑n=0N−1∣x(n)∣2
基于变量的运算
移位,序列x(n)的移位序列是x(n−m),m为正数时,表示向右移位(延时),反之向左移位(超前)
翻褶,序列x(n)的翻褶序列为x(−n),x(n−m)的翻褶序列为x(−n−m)
时间尺度变换
抽取(上抽样变换),为了减少抽样频率,xd(n)=x(Dn),D为整数
插值(下抽样变换),为了增加抽样频率,例如插入零值,
xI′(n)=⎩⎨⎧x(n/I),0,n=mI,I 为整数,m=0,±1,±2,...其他 n
既对幅度运算又对变量运算
差分运算
前向差分:Δx(n)=x(n+1)−x(n)
后向差分:∇x(n)=x(n)−x(n−1)
其中:∇x(n)=Δx(n−1)
卷积和运算
y(n)=x(n)∗h(n)=m=−∞∑∞x(m)h(n−m)=m=−∞∑∞x(n−m)h(m) 相关运算
rxy(m)=n=−∞∑∞x(n)y(n−m)
1.1.3 序列的卷积和
y(n)=x(n)∗h(n)=m=−∞∑∞x(m)h(n−m) 卷积和的运算步骤:
翻褶:选哑变量为 m,作 x(m),h(m),将 h(m)以 m=0 的垂直轴为对称轴翻褶成h(−m)。
移位:将h(−m)移位n,得h(n−m),n>0 时,右移n 位,n<0 时,左移∣n∣位。
相乘:将h(n−m)与x(m)在相同m 处的对应值相乘。
相加:将以上所有m处的乘积值叠加,就得到这一个n值下的 y(n)值。依上法取 n=⋅⋅⋅,−2,−1,0,1,2,⋅⋅⋅各值,即可得到全部y(n)值。
常用的卷积和求解方法:
图解加上卷积方法,一个结论:
若x(n)在N3⩽n⩽N4范围有非零值,h(n)在N1⩽n⩽N2范围有非零值。 y(n)=x(n)∗h(n)存在的范围为:N1+N3⩽n⩽N2+N4
列表方法,适用于两个有限长序列的卷积和
已知:x(n)=1,2,4,3,h(n)=2,3,5,取 y(n)=∑m=02h(m)x(n−m),则m应为0⩽m⩽2,列表如下
故y(n)=2,7,19,28,29,15
对位相乘相加法:将序列排成两行,且将其各自n最大的序列值对其(即按右端对其),然后作乘法运算,但是不要进位,最后将同一列的乘积值相加即得到卷积和结果,如下图
卷积和的序列长度:x(n)为N点长序列,h(n)为M点长序列,则$$y(n)=x(n)*h(n)$$为L=N+M−1点长列。
1.1.4 序列的相关性
互相关函数序列
定义
rxy(m)=n=−∞∑∞x(n)y(n−m) 互相关函数没有翻褶这一步骤
性质
rxy(m)与ryx(m)互为偶对称关系,即rxy(m)=ryx(−m)
rxy(m)不是偶对称函数,即rxy(−m)=rxy(m)
当x(n),y(n)是绝对可和信号时,则有limm→∞rxy(m)=0
rxy中m的有值范围
设有限长序列x(n),y(n)有值的范围分别为
x(n):Nx1⩽n⩽Nx2;y(n):Ny1⩽n⩽Ny2 则rxy(m)的取值范围为
−(Ny2−Nx1)⩽m⩽(Nx2−Ny1) 若 Ny2⩽Nx1,则 rxy(m)只在 m⩾0 时有值。
若 Ny1⩾Nx2,则 rxy(m)只在 m⩽0 时有值。
用卷积运算表示相关运算,即
rxy(m)=n=−∞∑∞x(n)y(n−m)=n=−∞∑∞x(n)y(−(m−n))=x(m)∗y(−m)
自相关函数序列
定义
rxx(m)=∑x(n)x(n−m)=x(m)∗x(−m) 性质
rxx(m)是实偶序列,即rxx(m)=rxx(−m)
m=0时自相关序列取最大值,即rxx(0)=∑n=−∞∞x2(n)>∣rxx(m)∣
当x(n)是绝对可和信号时,则有limm→∞rxx(m)=0
对于功率信号,互相关函数及自相关函数定义为
rxy(m)=M→∞lim2M+11n=−M∑Mx(n)y(n−m)rxx(m)=M→∞lim2M+11n=−M∑Mx(n)x(n−m) 若为周期信号,可写成
rxy(m)=N1n=0∑N−1x(n)y(n−m)rxx(m)=N1n=0∑N−1x(n)x(n−m)
1.1.5. 几种常用的典型序列
单位抽样(单位冲激,单位脉冲)序列
δ(n)=⎩⎨⎧1,n=00,n=0 单位阶跃序列
u(n)=⎩⎨⎧1,0,n⩾0n<0 矩形序列
RN(n)=⎩⎨⎧1,0,0⩽n⩽N−1其他 n 实指数序列
x(n)=anu(n),a 为实数 复指数序列
x(n)=e(σ+jω0)n=eσn[cosω0n+jsin(ω0n)] 当指数为纯虚数时,复指数序列为x(n)=ejω0n
正弦型序列
x(n)=Asin(ω0n+φ) 其中 A 为幅度,ω0为数字频率,φ为起始相位
注意:
δ(n)与u(n)的关系为
δ(n)=u(n)−u(n−1)u(n)=k=−∞∑nδ(k)=⎩⎨⎧1,0,n⩾0n<0 δ(n)和RN(n)的关系为
RN(n)=m=0∑N−1δ(n−m)=δ(n)+δ(n−1)+⋯+δ(n−N+ u(n)与RN(n)的关系为
RN(n)=u(n)−u(n−N) 因果序列x(n)与单位阶跃序列u(n)卷积,得到的是该因果序列的累加序列,即
u(n)∗x(n)=m=0∑nx(m)u(n−m)=m=0∑nx(m) 对于anu(n),∣a∣<1时序列收敛,∣a∣>1时序列发散
对于e(σ+jω0)n,可以将它分解为实部序列与虚部序列或分解为模序列与相角序列
x(n)=e(σ+jω0)n=eσnejω0n 数字角频率(ω)和模拟角频率(Ω)之间的关系为
ω=ΩT=Ω/fs=2πf/fs(rad) 其中fs=T1表示抽样频率,所以数字角频率事模拟角频率被抽样频率fs归一化后的弧度,关系可用下表表示
1.1.6 序列的周期性
周期性序列:
x(n)=x(n+Nr),r=0,±1,±2,⋅⋅⋅ 正弦型模拟信号xa(t)=Asin(Ω0t+φ)中,Ω0越大,则xa(t)变化越快;但由于x(n)=sin(ω0n+φ)=sin[(ω0+2πm)n+φ],当ω0变化时,x(n)是以 2π 为周期的,并不是 ω0 越大$,x(n)变化越快。
正弦型模拟信号xa(t)=Asin(Ω0t+φ)一定是周期性信号,周期为T0=2π/Ω0
正弦序列x(n)=Asin(nω0+φ)成为周期性序列的条件:
2π/ω为整数时,周期为2π/ω
2π/ω不是整数但是有理数时,周期将大于2π/ω
2π/ω是无理数时,不存在周期
如果正弦型序列时有一个连续正弦信号抽样得到的,抽样频率为fs,抽样时间间隔为T,连续信号周期为T,则有
ω02π=Ω0T2π=f0T1=TT0=f0fs 抽样时间间隔T和连续时间信号周期T0之间的关系如下
当 2π/ω0 是整数时,例如 2π/ω0=N,则正弦序列的周期为 N,有
T0=NT或fs=Nf0 即一个正弦信号周期(T0)中有N个抽样周期(T),或说抽样频率(fs)是正弦信号频率(f0)的整数(N)倍。
当2π/ω0是有理数时,即2π/ω0=N/M,其中,N、M是互为素数的正整数,则正弦序列的周期为 N,有
MT0=NT或Mfs=Nf0 即M个正弦信号周期(T0)中应有 N 个抽样周期(T),或说N倍正弦信号频率(f0)应等于M倍抽样频率。
1.1.7 用单位抽样序列表示任意序列
δ(n)的选择性
x(m)⋅δ(n−m)=⎩⎨⎧x(n),0,当 m=n其他 m 任意序列x(n)与单位抽样序列δ(n)的卷积和等于该序列x(n)本身
x(n)=m=−∞∑∞x(m)δ(n−m)=x(n)∗δ(n) 任 意 序 列x(n)与单位抽样序列的移位序列 δ(n−n0)的卷积和就得到此序列作相同移位的序列x(n−n0)
x(n)∗δ(n−n0)=m=−∞∑∞x(m)δ(n−n0−m)=x(n−n0)
1.2 线性移不变系统
一个离散时间系统是将输入序列按照所需要的目的变换成输出序列的一种运算,若以$T[\cdot]$表示这种运算,则有
y(n)=T[x(n)] 如图
下面讨论线性移不变系统
1.2.1 离散时间线性系统
满足叠加原理
a1y1(n)+a2y2(n)=a1T[x1(n)]+a2T[x2(n)]=T[a1x1(n)+a2x2(n)] 证明一个系统是线性系统需对所有常系数(包括复数)以及所有输入(包括复数)满足叠加性(可加性和比例性)的两个条件。
在线性系统中,在全部时间上,零输入一定产生零输出。
增量线性系统(例如y(n)=4x(n)+6),任意两个输入的响应之差与两个输入之差呈线性关系
1.2.2 离散时间移不变系统
若系统的输人输出关系不随时间而变化,则称它为移(时)不变系统,即若
T[x(n)]=y(n) 则
T[x(n−n0)]=y(n−n0) 移不变系统的输出序列随输入序列的移位而作相同的移位,且保持输出序列的形状是不变的
若系统有一个移变增益,例如y(n)=nx(n),则一定是移变系统
若一个系统在时间轴上有压缩或扩展,例如y(n)=x(Mn),则一定是移变系统
1.2.3 离散时间线性移不变系统(LSI系统)
单位冲激响应
单位抽样响应(也称单位冲激响应或单位脉冲响应)是指输入为单位抽样序列δ(n)时,LSI 系统的输出序列(或称输出响应),一般用h(n)表示,即
h(n)=T[δ(n)] LSI系统输出序列与输入序列在时域中的关系——卷积和关系
y(n)=m=−∞∑∞x(m)h(n−m)=x(n)∗h(n) 即LSI系统的输出序列是输入序列与系统单位抽样响应的卷积和。
LSI系统卷积和运算的性质
交换律
x(n)∗h(n)=m=−∞∑∞x(m)h(n−m)=m=−∞∑∞x(n−m)h(m)=h(n)∗x(n) 结合律
x(n)∗h1(n)∗h2(n)=[x(n)∗h1(n)]∗h2(n)=x(n)∗[h1(n)∗h2(n)]交换律[x(n)∗h2(n)]∗h1(n) 分配律
x(n)∗[h1(n)+h2(n)]=x(n)∗h1(n)+x(n)∗h2(n)
1.2.4 因果系统
因果系统是指系统的输出不发生在输入之前的系统,即
y(n0) 只取决于 x(n)∣n⩽n0 对于因果系统,n<n0的输入相同则n<n0的输出也相同
考察系统因果性时,只看输入x(n)和输出y(n)的关系,而不讨论其他以n为变量函数的影响,例
y(n)=(n+2)x(n) 是因果系统y(n)=x(n)sin(n+4)也是因果系统 LSI系统是因果系统的充要条件是单位冲激响应h(n)是因果序列,即
h(n)=0,n<0
1.2.5 稳定系统
若系统满足稳定输入产生稳定输出,就是因果系统,即
若∣x(n)∣⩽M<∞则∣y(n)∣⩽P<∞ LSI系统是稳定系统的充要条件是单位抽样响应绝对可和,即
n=−∞∑∞∣h(n)∣=P<∞ 1.3 常系数线性差分方程
离散时间系统的输入和输出关系可用常系数线性差分方程表示,即
k=0∑Naky(n−k)=m=0∑Mbmx(n−m) 求解线性常系数差分方程的三种方法:
常系数差分方程表示的系统,只有当边界条件是使系统松弛的,则该系统是线性,移不变系统,所谓松弛状态是指
x(n)∣n<n0=0 则边界条件必须满足:
y(n)∣n<n0=0 1.4 连续时间信号的采样
1.4.1 模拟信号的连续采样
理想抽样信号
设 xu(t)为模拟信号,x^a(t)为理想抽样信号,则有
x^a(t)=xa(t)⋅p(t)=xa(t)⋅δT(t) 其中,抽样信号 p(t)就是周期性的单位冲激信号
δT(t)=m=−∞∑∞δ(t−mT) 因而有
x^a(t)=xa(t)m=−∞∑∞δ(t−mT)=m=−∞∑∞xa(t)δ(t−mT)=m=−∞∑∞xa(mT)δ(t−mT) T为采样周期,fs=1/T为抽样频率,Ωs=2πfs=2π/T是抽样角频率
理想抽样信号的频谱
X^a(jΩ)=T1k=−∞∑∞Xa[j(Ω−kT2π)] 理想抽样信号的频谱X^a(jΩ)是被抽样的模拟信号的频谱 Xa(jΩ)的周期延拓,在角频率Ω轴上其延拓周期为
Ωs=T2π=2πfs 即频率轴f上以抽样频率fs为周期而周期延拓。
1.4.2 时域采用定理
如果模拟信号是带限信号(频带有限信号),信号的最高频率分量为 fh,设抽样频率为 fs,由于抽样后信号的频谱等于模拟信号频谱按抽样频率 fs作周期延拓,只有当
fh<2fs 时,周期延拓的频谱分量才不会产生交叠,称抽样频率之半(fs/2)为折叠频率。即
2fs=2T1或2Ωs=Tπ 奈奎斯特抽样定理:若 xa(t)是频带宽度有限的信号(称带限信号),要想抽样后的信号能够不失真地还原出原信号,则必须抽样频率fs大于信号最高频率分量fh 的两倍,即
fs>2fh 数字频率ω是模拟角频率对抽样频率fs的归一化频率,因而抽样频率fs对于的数字频率为
ωs=ΩsT=2πfs/fs=2π 因此数字域的抽样频率ωs等于2π,也就是延拓周期等于2π,折叠频率fs/2所对应的数字域频率为
2ωs=π 按照采样定理,信号数字域最高频率ωh应小于π。
1.4.3 模拟信号的实际采样
在实际中,抽样脉冲串不是冲激函数串δT(t),而是一点宽度τ的矩形周期脉冲串,这时奈奎斯特抽样定理仍然有效
1.4.4 带通信号的抽样
带通信号的频谱存在某一频段范围,而不是零频范围,如图
其最高频率为fh,带宽为Δfo,算术中心频率为fa=fh−2Δfo,一般来说fa≫Δfo,即通带中心频率远大于通带带宽。
当fh=rΔfo,即带通信号的最高频率是其通带宽度的整数倍时,则选抽样频率fs为
fs=2Δfo 即所取抽样频率为带通信号通带宽度的两倍,其抽样后的频谱是带通信号频谱以此fs的整数倍而周期延拓后的频谱,如图
图中 r=5,显然没有频谱混叠现象。
当 fh=r′Δfo,r′=整数。这是最一般的情况,即带通信号最高频率不等于其带宽的整数倍时,这时,可保持fh不变,将通带下端延伸到使其带宽为Δfo′,使得满足
fh=rΔf0′,r=⌊r′⌋(⌊⌋表示取整数部分) 此时有 Δfo′>Δfo,这时选抽样频率为
fs=2Δfo′=2fh/⌊r′⌋=2fh/⌊fh/Δfo⌋ 后面的情况同上述分析。
带通信号抽样频率fs的取值范围为(其中Δfo为带宽)
2Δfo⩽fs<4Δfo fs的下限是Δfofh=整数时的情况,而其上限则对应于Δfofh=整数,且是最不利的情况。由于上面讨论的带通信号的抽样频率fs不满足fs>2fh的要求,因而可称其为亚奈奎斯特抽样频率。
1.4.5 连续时间信号
xa(t)、理想抽样信号
x^a(t)以及抽样序列
x(n)
x^a(t)在工程应用中无法实现,故工程中采样抽样序列x(n)。
x^a(t)和 x(n)的本质差别:x^a(t)本质上是连续时间信号,是在 t=nT 时的冲激串, 在t=nT 时,x^a(t)=0,其每个冲激的幅度都是无穷大,存在时间(宽度)为无穷小,x^a(t)的大小是以冲激的积分面积表示的;抽样序列 x(n)则是整数变量 n 的函数,这里时间已归一化,x(n)本身已没有抽样率的信息,当 n=整数时,x(n)=xa(nT),即抽样点上抽样序列的幅值是确定数值,在 n=整数时,x(n)无定义,不是零值。
连续时间信号 xa(t)、理想抽样信号x^a(t)以及抽样序列 x(n)三者的关系
x(n)=xa(t)t=nT=xa(nT),−∞<n<∞
1.4.6 时域信号的插值重构
如果满足奈奎斯特抽样定理,即信号谱的最高频率小于折叠频率,则抽样后不会产生频谱混叠,可以由信号的抽样值经插值而重构原信号xa(t)。当∣Ω∣<Ωs/2时,只存在k=0一项,即有
X^a(jΩ)=T1Xa(jΩ),∣Ω∣<2Ωs 将X^a(jΩ)送人作为重构用的理想低通滤波器 H(jΩ)
H(jΩ)=⎩⎨⎧T,0,∣Ω∣<2Ωs∣Ω∣⩾2Ωs 就可以得到原信号频谱,即
Ya(jΩ)=X^a(jΩ)H(jΩ)=Xa(jΩ) 所以输出端即为原模拟信号。
信号重构的抽样内插公式:
ya(t)=m=−∞∑xa(mT)π(t−mT)/Tsin[π(t−mT)/T] xo(mT)经此公式而得到连续信号xa(t),而sin[π(t−mT)/T]/[π(t−mT)/T]称为内插函数,如图
在抽样点mT上,函数值为 1,在其余抽样点上,函数值为零,不影响其他抽样点。也就是说,xa(t) 等于各xa(mT)乘以对应的内插函数的总和。在每一个抽样点上,只有该点所对应的内插函数不为零,这使得各抽样点上信号值不变,而抽样点之间的信号则由各加权抽样函数波形的延伸叠加而成,如图
因而内插函数 h(t)的作用是在各抽样点之间起连续插值的作用。
1.4.7 正弦信号的抽样
正弦型信号 xa(t)=Asin(Ωot+φ)=Asin(2πfat+φ)的频谱在 f=fo处为δ函数,一般来说,正弦信号的抽样频率必须满足fs>2fo. 因为若取 fs=2fo,则有以下几种情况发生:
当φ=0时,一个周期抽取的两个点为x(0)=x(1)=0,相当于xα(0)和xα(π) 两个点,故不包含原信号的任何信息。
当φ=2π时,则有x(0)=A,x(1)=−A,此时从x(n)可以恢复xa(t)。
当φ为已知,且 0<φ<2π时,恢复的不是原信号,但经过变换后,可得到原信号。
当φ为未知数时,抽样后不能恢复出原信号xa(t)。所以,至少要取 fs>2fo,避免产生不确定性。
处理周期性正弦序列时,应该注意:
对抽样后的离散周期性的正弦序列作截断时,其截断长度必须为序列周期的整数倍,才不会产生频域的泄露。
离散正弦序列不宜补零后作频谱分析,否则会产生频域泄露。
考虑到作DFT时,当要求数据个数为N=2P时,正弦信号一个周期中最好抽取4个点。
满足∣ω2−ω1∣=2πk,k为整数,的两个正弦序列是相同的,有以下两个结论:
两个不同频率的模拟正弦型信号,如果用同一抽样频率fs对其抽样,得到的序列可能是相同的序列,我们没法从序列中区分出它们分别来源于哪一个模拟正弦型信号。
同一个模拟正弦型信号,如果用两个不同的抽样频率抽样后,所得到的序列仍可能是相同的,我们没法确定其原抽样频率。