时域表示
计算如下有限长序列的
L1范数、
L2范数和
L∞范数:
(a)
{x1[n]}={4.50−2.68−0.143.912.62−0.43−4.813.21−0.55}
(b)
{x2[n]}={0.922.343.371.90−2.59−0.753.483.33}
解:
∥x∥=(−∞∑∞∣x[n]∣p)1/p
∥x∥1=−∞∑∞∣x[n]∣∥x∥2=(−∞∑∞∣x[n]∣2)1/2
∥x∥∞=∣x∣max
(a)
∥x1∥1=22.85∥x1∥2=9.1396∥x1∥∞=4.81
(b)
∥x2∥1=18.68∥x2∥2=7.1944∥x2∥∞=3.48
请参考这篇文章离散时间信号的时域表示
序列的运算
基本运算
已知序列
x[n]y[n]w[n]={2,0,−1,6,−3,2,0},={8,2,−7,−3,0,1,1},={3,6,−1,2,6,6,1},−3−5−2≤n≤3≤n≤1≤n≤4
上述序列在给定区间以外的样本值都为零。生成以下序列:
(a)c[n]=x[n+3](b)d[n]=y[n−2](c)e[n]=x[−n](d)u[n]=x[n−3]+y[n+3](e)v[n]=y[n−3]⋅w[n+2](f)s[n]=y[n+4]−w[n−3](g)r[n]=3.9w[n]
解:
(a)
序列
x[n]左移
3个单位,则
c[n]={2,0,−1,6,−3,2,0},−6≤n≤0
(b)
序列
y[n]右移
2个单位,则
(b)d[n]={8,2,−7,−3,0,1,1},−3≤n≤3
(c)
序列
x[n]反褶,则
e[n]={0,2,−3,6,−1,0,2},−3≤n≤3
(d)
x[n]右移
3个单位得到
x[n−3]={2,0,−1,6,−3,2,0},0≤n≤6
y[n]左移
3个单位得到
y[n+3]={8,2,−7,−3,0,1,1},−8≤n≤−2
则
u[n]={8,2,−7,−3,0,1,1,0,2,0,−1,6,−3,2,0}−8≤n≤6
(e)
y[n]右移
3个单位得到
y[n−3]={8,2,−7,−3,0,1,1},−2≤n≤4
w[n]左移
2个单位得到
w[n+2]={3,6,−1,2,6,6,1},−4≤n≤2
则
v[n]={0,0,−8,4,−42,−18,0,0,0}−4≤n≤4
(f)
y[n]左移
4个单位得到
y[n+4]={8,2,−7,−3,0,1,1},−9≤n≤−3
w[n]右移
3个单位得到
w[n−3]={3,6,−1,2,6,6,1},1≤n≤7
则
s[n]={8,2,−7,−3,0,1,1,0,0,0,−3,−6,1,−2,−6,−6,−1},−9≤n≤7
(g)
r[n]={11.7,23.4,−3.9,7.8,23.4,23.4,3.9},−2≤n≤4
卷积运算
证明一个长度为
M的序列与一个长度为
N的序列进行卷积,可得到一个长度为
(M+N−1)的序列。
参考这篇文章中卷积后的长度部分
设
x[n]、
y[n]、
w[n]分别表示长度为
N、M和
L的三个序列,每个序列的第一个样本都出现在
n=0处,序列
x[n]∗y[n]∗w[n]的长度是多少?
解:
借用上题的结论,两序列卷积后的长度为两序列长度之和减一。
令
h[n]=x[n]∗y[n],则序列
h[n]的长度为
H=M+N−1,所以序列
h[n]∗w[n]的长度为
H+L−1=M+N+L−2,所以序列
x[n]∗y[n]∗w[n]的长度为
M+N+L−2
求下面序列与其自身的卷积
x[n]={1,−1,1},−1≤n≤1
解:
所以
x[n]∗x[n]={1,−2,↑3,−2,1}−2≤n≤2
可以观察到
x[n]它是左右对称的,事实上序列自己与自己卷积得到的序列都是左右对称的。
不理解上述计算卷积算法的,请参考用多项式乘法快速计算卷积。
设
y[n]=x1[n]∗x2[n]且
v[n]=x1[n−N1]∗x2[n−N2],试用
y[n]来表示
v[n]。
解:
首先从数学的角度对公式进行推导
y[n]=x1[n]∗x2[n]=m=−∞∑∞x1[m]x2[n−m]
v[n]=x1[n−N1]∗x2[n−N2]=m=−∞∑∞x1[m−N1]x2[n−m−N2]
令
k=m−N1得到
v[n]=k=−∞∑∞x1[k]x2[n−(N1+k)−N2]=k=−∞∑∞x1[k]x2[(n−N1−N2)−k]=y[n−N1−N2]
重点是怎么理解这一个结果,该结果表明,如果序列
x1[n]时移
N1,
x2[n]时移
N2,那么卷积后的序列
y[n]时移
N1+N2。
设
x[n]={1}是定义在
0≤n≤N−1长度为
N的序列,而
h[n]={1}是定义在
0≤n≤N−3范围内长度为
N−2的序列,不进行卷积运算,求
y[n]=x[n]∗h[n]最大样本的值的位置。
解:
由于序列
x[n]、h[n]的每一项的数值都是
1,所以
x[n]与
h[n]各项对应的乘积也为
1,所以要使
y[n]最大,就得包含更多的项相加,由于
h[n]的长度小于
x[n],所以最大包含
h[n]的长度
N−2项,所以由
y[n]=x[n]∗h[n]=m=−∞∑∞h[m]x[n−m]
所以有
0≤m≤N−30≤n−m≤N−1
由于要保证要取到所有的
h[n]才能保证
y[n]最大,所以对于任意
0≤m≤N−3存在
0≤n−m≤N−1
即
mmax≤n≤mmin+N−1⇒N−3≤n≤N−1
所以
N−3,N−2,N−1处
y[n]有最大值为
N−2。
圆周运算
考虑序列
g[n]={−3,0,4,9,2,0,−2,5},−4≤n≤3
(a)
求
g[n]向右圆周平移5个样本周期得到的序列
h[n]
(b)
求
g[n]向左圆周平移4个样本周期得到的序列
w[n]
解:
(a)
h[n]={9,2,0,−2,5,−3,0,4}−4≤n≤3
(b)
w[n]={5,−3,0,4,9,2,0,−2}−4≤n≤3
序列的分类
基于对称性的分类
证明实数值序列
x[n]的平均功率
Px是其偶部和奇部平均功率
Pxev和
Pxod之和。
解:
Px=2K+11n=−K∑K∣x[n]∣2=2K+11n=−K∑K(xev[n]+xod[n])2=2K+11(n=−K∑K(xev[n])2+n=−K∑K(xod[n])2+2n=−K∑Kxev[n]xod[n])=Pxev+Pxod+2K+12n=−K∑K((x[n]+x[−n])(x[n]−x[−n]))=Pxev+Pxod+2K+12(n=−K∑Kx[n]2−n=−K∑Kx[−n]2)=Pxev+Pxod
上式用到了
n=−K∑Kx[n]2−n=−K∑Kx[−n]2=0
求序列
x[n]={−1+j3,2−j7,4−j5,3+j5,−2−j},−2≤n≤2的共轭对称部分和共轭反对称部分。
解:
x[−n]={−2−j,3+j5,4−j5,2−j7,−1+j3},−2≤n≤2
x∗[−n]={−2+j,3−j5,4+j5,2+j7,−1−j3},−2≤n≤2
xcs[n]=(x[n]+x∗[−n])/2={−1.5+j2,2.5−j6,4,2.5+j6,−1.5−j2}
xca[n]=(x[n]−x∗[−n])/2={0.5+j,−0.5−j,−j5,0.5−j,−0.5+j}
证明当所有
n≥0时,因果序列
x[n]可以从它的偶部
xev[n]中完全恢复出来,而当所有
n>0才可从它的奇部
xod[n]中恢复出来。
解:
由于序列
x[n]为因果序列,则
x[−n]=0,n>0。由于
xev[n]=21(x[n]+x[−n])
当
n=0时,
x[0]=xev[0]
当
n>0时,
x[n]=2xev[n]
所以
x[n]=⎩⎪⎨⎪⎧2xev[n],xev[0],0,n>0n=0n<0
而
xod[n]=21(x[n]−x[−n])
当
n=0时,
x[0]=0
当
n>0时,
x[n]=2xod[n]
所以
x[n]={2xod[n],0,n>0n≤0
因果复序列
y[n]能够从它的共轭对称部分
ycs[n]中完全恢复出来吗?能从它的反共轭对称部分
yca[n]中恢复出来吗?证明你的结论。
解:
因为
y[n]为因果序列,所以
y∗[−n]=0,n>0,由
ycs[n]=21(y[n]+y∗[−n])
当
n=0时,
yre[0]=ycs[0]
当
n>0时,
y[n]=2ycs[n]
由于在
n=0时,只能恢复出实部,所以由共轭对称部分不能完全恢复出
y[n]。
同理,共轭反对称部分在
n=0只能恢复出
yim[0],也不能完全恢复出
y[n]。
考虑以各自共轭对称和共轭反对称之和的形式表示的两个复值序列
h[n]和
g[n],即
h[n]=hcs[n]+hca[n]且
g[n]=gcs[n]+gca[n]。试确定下列序列是共轭对称还是共轭反对称的。
(a)
y1[n]=hcs[n]∗gcs[n]
(b)
y2[n]=hca[n]∗gcs[n]
(c)
y3[n]=hca[n]∗gca[n]
解:
(a)
y1[n]y1∗[−n]=m=−∞∑∞hcs[m]gcs[n−m]=(m=−∞∑∞hcs[−m]gcs[−(n−m)])∗=m=−∞∑∞hcs∗[−m]gcs∗[−(n−m)]=m=−∞∑∞hcs[m]gcs[n−m]=y1[n]
(b)
y2[n]=−y2∗[−n]
(c)
y3[n]=y3∗[−n]
能量信号与功率信号
计算长度为
N的序列
x[n]=sin(2πkn/N)的能量,
0≤n≤N−1。
解:
ε2=n=0∑N−1sin2(2πkn/N)=n=0∑N−121−21cos(4πkn/N)=2N−21n=0∑N−1cos(4πkn/N)
令
C=n=0∑N−1cos(4πkn/N),S=n=0∑N−1sin(4πkn/N)
则
C+jS=n=0∑N−1e−jN4πkn=1−e−jN4πk(1−e−jN4πkN)=0
这说明
C=0,所以
ε2=2N
计算下列序列的能量:
(a)
xa[n]=Aαnμ[n],∣α∣<1,
(b)
xb[n]=n21μ[n−1]
解:
(a)
εa2=n=0∑∞A2(α2)n=1−α2A2
(b)
εb2=n=1∑∞n41=90π4
周期序列和非周期序列
求下列周期序列的基本周期:
(a)
x~a[n]=ej0.25πn
(b)
x~b[n]=Re(ejπn/8)+Im(ejπn/5)
解:
(a)
0.25π2π=8所以周期
N=8。
(b)
π/82π=16所以
N1=16
π/52π=10所以
N2=10
所以周期为二者的最小公倍数
N=LCM(N1,N2)=LCM(16,10)=80
其他分类
证明序列
x[n]=n(−1)n+1μ[n]不是绝对可和的。
解:
n=−∞∑∞∣x[n]∣=n=1∑∞n1=∞
证明下面序列是绝对可和的。
(a)
x1[n]=αnμ[n−1]
(b)
x2[n]=nαnμ[n−1]
其中
∣α∣<1
解:
(a)
n=1∑∞∣α∣n=1−∣α∣∣α∣<∞
(b)
n=1∑∞∣nα∣n=n=1∑∞n∣α∣n=∣α∣+2∣α∣2+3∣α∣3+...+n∣α∣n+...=(∣α∣+∣α∣2+∣α∣3+...+∣α∣n+...)+(∣α∣2+∣α∣3+...+∣α∣n+...)+(∣α∣3+...+∣α∣n+...)+...=1−∣α∣∣α∣+1−∣α∣∣α∣2+1−∣α∣∣α∣3+...+1−∣α∣∣α∣n+...=1−∣α∣1n=1∑∞∣α∣n=(1−∣α∣)2∣α∣<∞
下面序列中,哪些序列是有界序列?
(a)
x[n]=Aαn,其中
A和
α是复数,且
∣α∣<1
(b)
v[n]=(1−n21)μ[n−1]
解:
(a)
∣x[n]∣=∣A∣∣α∣n
对于
n<0存在
n使得
∣x[n]∣=∞,所以
x[n]不是有界序列
(b)
∣v[n]∣=∣1−n21∣,n≥1
则
∣v[n]∣<1
所以序列
v[n]是有界序列。
典型序列及其表示
使用单位阶跃序列
μ[n]表示序列
x[n]=1,−∞<n<∞。
解:
x[n]=μ[n]+μ[−n−1]
单位抽样序列
δ[n]与单位阶跃序列
μ[n]的关系。
解:
δ[n]=μ[n+1]−μ[n]
μ[n]=m=0∑∞δ[n−m]=m=−∞∑nδ[m]
试用单位阶跃序列
μ[n]表示序列
x[n]={2,1,n≥2n<2
解:
x[n]=μ[n−2]+μ[−n+3]
求下列卷积和的闭式表示:
(a)
αnμ[n]∗μ[n]
(b)
nαnμ[n]∗μ[n]
解:
(a)
αnμ[n]∗μ[n]=m=0∑nαm=1−α1−αn+1,n≥0
(b)
nαnμ[n]∗μ[n]=m=0∑nmαm,n≥0
信号的相关
求下面序列的自相关序列,并证明它们均为偶序列,自相关序列最大值的位置在哪里?
(a)
x1[n]=αnμ[n],∣α∣<1
(b)
x2[n]={1,0≤n≤N−10,其它
解:
(a)
ra[l]=n=−∞∑∞x1[n]x1[n−l]=n=−∞∑∞αnμ[n]αn−lμ[n−l]
当
l≥0时,则
ra[l]=n=l∑∞α−l(α2)n=α−l1−α2α2l=1−α2αl
当
l<0时,则
ra[l]=n=0∑∞α−l(α2)n=α−l1−α21=1−α2α−l
所以
ra[l]=⎩⎪⎨⎪⎧1−α2αl,1−α2α−l,l≥0l<0=1−α2α∣l∣
当
l=0时,
ra[l]最大,最大为
1−α21
ra[−l]=1−α2α∣l∣=ra[l]
所以
ra[l]为偶序列。
(b)
解:
rb[l]=n=−∞∑∞xb[n]xb[n−l]
如果
l+N−1<0,即
l<−N+1则
rb[l]=0
如果
0≤l+N−1<N−1,即
−N+1≤l<0,则
rb[l]=n=0∑l+N−11=l+N
如果
0≤l≤N−1,则
rb[l]=n=l∑N−11=N−l
如果
l>N−1,则
rb[l]=0
所以
rb[l]=⎩⎪⎨⎪⎧N−∣l∣,0,∣l∣≤N−1∣l∣>N−1
当
l=0时,
rb[l]有最大值为
rb[0]=N。
而
rb[−l]=rb[l]
所以序列
rb[l]为偶序列。
计算下列各周期序列的自相关序列以及它们的周期:
(a)
x~1[n]=cos(πn/M),其中
M为正整数
(b)
x~2[n]=n模6
(c)
x~3[n]=(−1)n
解:
(a)
π/M2π=2M
所以序列
x~1[n]的周期为
2M,则
r~[l]=2M1n=0∑2M−1cos(Mπn)cos(Mπ(n−l))=2M1n=0∑2M−1cos(Mπn){cos(Mπn)cos(Mπl)+sin(Mπn)sin(Mπl)})=2M1cos(Mπl)n=0∑2M−1cos2(Mπn)+4M1sin(Mπl)n=0∑2M−1sin(M2πn)=2M1cos(Mπl)n=0∑2M−1cos2(Mπn)=2M1cos(Mπl)n=0∑2M−121(1+cos(M2πn))
考虑
n=0∑2M−121(1+cos(M2πn))并令
N=2M,则
n=0∑N−121(1+cos(N4πn))=2N+21n=0∑N−1cos(N4πn)
令
C=n=0∑N−1cos(N4πn),S=n=0∑N−1sin(N4πn),则
C+jS=n=0∑N−1ej4πn/N=0
这表示
C=0,所以
n=0∑2M−121(1+cos(M2πn))=2N=M
所以
r~[l]=21cos(Mπl)
r~[l]的周期为
2M。
(b)
x~2[n]={0,1,2,3,4,5},0≤n≤5
r~[l]=61n=0∑5x~[n]x~[n−l]
所以
r~[0]r~[1]...=61(x~[0]x~[0]+x~[1]x~[1]+x~[2]x~[2]+x~[3]x~[3]+x~[4]x~[4]+x~[5]x~[5])=655=61(x~[0]x~[5−1]+x~[1]x~[0]+x~[2]x~[1]+x~[3]x~[2]+x~[4]x~[3]+x~[5]x~[4])=640
(c)
x~3[n]={1,−1},0≤n≤1
r~[l]=21n=0∑1x~[n]x~[n−l],0≤l≤1
r~[0]=1,r~[1]=−1