版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
《数字信号处理》辅导
一、离散时间信号和系统的时域分析
(一)离散时间信号
(1)基本概念
信号:信号传递信息的函数也是独立变量的函数,这个变量可以是时间、空
间位置等。
连续信号:在某个时间区间,除有限间断点外所有瞬时均有确定值。
模拟信号:是连续信号的特例。时间和幅度均连续。
离散信号:时间上不连续,幅度连续。常见离散信号一一序列。
数字信号:幅度量化,时间和幅度均不连续。
(2)基本序列(课本第7——10页)
1,n-Q1,72>0
1)单位脉冲序列<?(«)=<2)单位阶跃序列〃(〃)=<
0,〃。00,72<0
04〃4N-1
3)矩形序列&4)实指数序列anu(n)
n<0.n>N
5)正弦序列》(〃)=弦sin(g〃+6)6)复指数序列x(n)=<Jconcern
(3)周期序列
1)定义:对于序列x(〃),若存在正整数N使X(〃)=X(”+N),-8<〃<OO
则称x(〃)为周期序列,记为文(“),N为其周期。
注意正弦周期序列周期性的判定(课本第10页)
2)周期序列的表示方法:
a.主值区间表示法
b.模N表示法
3)周期延拓
设x(〃)为N点非周期序列,以周期序列L对作x(")无限次移位相加,即可得到
周期序列£(〃),即
00
x(n)-Zx(n-iL)
/=-co
当LNN时,x(〃)=£(〃)&,(〃)当L<N时,x(n)*x(n)Nx(n)
(4)序列的分解
序列共也对称分解定理:对于任意给定的整数M,任何序列x(“)都可以分解成
关于c=M/2共轨对称的序列七(〃)和共加反对称的序列x0(n)之和,即
x(n)=xe(n)+%〃(〃),一oov〃<oo
并且
以〃)=尹⑺+””一〃)]W〃)=/(〃T(M-〃)]
(4)序列的运算
1)基本运算
运舁性质描述
序列相乘y(〃)=xx(n)x2(«)y(ri)=ax(n)
序列相加y(n)=xl(n)+x2(n)
序列翻转),(〃)=x(—〃)(将x(〃)以纵轴为对称轴翻转)
尺度变换y(〃)=x(mn)(序列x(〃)每隔m-1点取一点形成的序列)
用单位脉冲
M〃)=
序列表示/=—00
2)线性卷积:
将序列x(〃)以y轴为中心做翻转,然后做m点移位,最后与x(〃)对应点相
乘求和——翻转、移位、相乘、求和
00
定义式:y(〃)=E芭(施)/(〃-M=■(")*—(〃)
,n=-ao
线性卷积的计算:A、图解
B、解析法
C、不进位乘法(必须掌握)
3)单位复指数序列求和(必须掌握)
N—1[_~j(oN(e"oN/2_-joN/2)j°N/2(ej/N/2_-j(oN/2)/(2))
£日刖=ee(e
J<a/2/2/2
'l-e-_e-加2(e*w2_e-网2)-e->(e>-e->)/(2/)
〃=0
„->(N-I)/2sin(oN/2)
一C
sin(tw/2)
如果①=2成/N,那么根据洛比达法则有
四3N/2)=Nb(())(&=0)(或N6(N)*=TV))
sin(6y/2)
可以结合作业题3.22进行练习
(5)序列的功率和能量
能量:£=^|x(n)|2
n=-oo
1N
功率:P=lim——£|x(»)|2
"2N+1,占
(6)相关函数一一与随机信号的定义运算相同
(二)离散时间系统
1.系统性质
(1)线性性质
定义:设系统的输入分别为王(〃)和工2(〃),输出分别为X(〃)和为(〃),即
%(〃)=71一(〃)],y2(n)=T[x2(n)]
统的输对于任意给定的常数。、b,下式成立
y(n)=T[ax[(n)+bx2(n)]=ayi(n)+by2(n)
则该系统服从线性叠加原理,为线性系统,否则为非线性系统。
判定系统的线性性质时,直接用定义
(2)时不变性质
统的如果系统对输入信号的运算关系在整个运算过程中不随时间变化,则称该系
统是时不变系统。即对任意给定的整数i,若下式成立:
y(n-i)^T[x(n-i)]
则称该系统为时不变系统,否则为时变系统。
判定系统的时不变性质时,直接用定义
(3)系统的因果性
定义:如果系统n时刻的输出序列只取决于n时刻及以前的输入序列,而与n时
刻以后的输入序列无关,则称该系统具有因果性质,即系统是因果系统,否则是
非因果系统。
离散时间LTI系统具有因果性的充要条件是:系统的单位脉冲响应〃(〃)满足
h(n)=0,H<0
(4)系统的稳定性
定义:对任意有界的输入,系统的输出都有界,则该系统是稳定的,否则是不稳
定的。
离散时间LTI系统具有因果性的充要条件是:系统的单位脉冲响应出〃)满足绝
对可和,即力⑺K00
/=-00
(5)对离散时间LTI系统的描述
(1)时域:差分方程
(2)Z域:系统函数H(z)
2.信号过系统
y(ri)=h(n)*x(n)
用线性卷积的相关知识计算,信号系统学的基本性质可以套用
二、离散时间信号和系统的频域分析
(一)离散时间信号
1.序列傅里叶变换(SequenceFourierTransform)(即本书中的离散时间信号
的傅里叶变换)
(1)定义
SFT:X(eJ(u)=SFT[x(n)]=Jx(〃)"""",一oo<<y<oo
"=-8
ISFT:x(〃)=ISFT[X(ejM)]^-^-^X(ejM)ejmidco,-co<n<oo
说明:
1、物理意义:序列傅里叶变换本质上是序列的一种分解,它将一般序列分解为
无穷多个数字角频率[-肛乃]中的复指数序列。称X(e*")为序列x(〃)的频谱,其
模IX(*)|称为幅频特性,其幅角arg[X(e"")]=伙⑼称为相频特性。
2、尽管序列x(〃)是离散时间信号,但它的序列傅里叶变换对数字角频率。而言
却是连续函数,因此,序列x(")的傅里叶变换是连续的。
+Mn
3、X(e"i)=力x(n)e-^=X(e初)
n=-ao
由上式可知,序列傅里叶变换X(e,。)是以2〃为周期的周期函数,其原因正是由
于一.对。而言以2%为周期,即数字角频率相差2%的所有单位复指数序列等
价。因此,对-8</<8的所有单位复指数序列只有一个周期。对于离散时间信
号,由于的周期性,使得口=0或2〃的整数倍都表示信号的直流分量,而万的奇
数倍表示信号的最高频率。
(2)性质
名称性质描述
线性性质SFT[axx(〃)+bx2(n)]=a»SFT[x](〃)]+b^SFT[x2(n)]
时移性质SFT[x(n-m)]=”加"・SFT[x(〃)]
频移性质SFT[ej,^^x(n)]=
共也对称性质SbK(〃)]=XWSFT[jx,(〃)]=
SFT[xe(W)]=Re[X(j)],SFn1(〃)]=jIm[X(e>)]
线性卷积性质SFT[x(n)*y(")]=SF7Ix(n)].SF7Iy(“)]
帕斯瓦尔定理
«=-<x>乙?I
相乘性质SF7Ix(〃)y(”)]=(J"
序列乘以nSfT[〃・尤(〃)]=j[dX(e"°)/dco]
(3)基本序列的傅里叶变换
序列傅里叶变换
3(九)1
12病(①)
%(〃)e-j3(NT)/2sin(等)/sin(?)
anu(ri)(|a\<1)(l-ae-j<aY'
(24/为为有理数)2君(g一%)
cos刃()〃(2乃//为有理数)7T信(69—例)+$(/+4)]
singnQi/例)为有理数)一力r[S(<y-%)-5(切+%)]
u(n)(1_6-肚尸+痛(0)
2.Z变换(不熟悉的复习信号系统相关内容,或本书2.3相关内容)
(1)定义
1
ZT:X(z)==ZT[x(n)]=£x(n)Z-'Rx_<|z|<Rx+
n=-<x)
IZT:x(〃)=/ZT[X(z)]=—!-$X(z)z"〃zR
2兀产
(2)性质——课本49页表2.3.3
(3)收敛域与基本序列Z变换——课本45页表2.3.1、表2.3.2
3.离散时间信号Z变换与SFT的关系
Z变换是由SFT推广得到的,反过来,如果某序列的Z变换的收敛域包括z=e%
则也可以通过ZT求得序列的SFTo即
X(z)L/=£宜〃)""""=乂(〃3)
rt=-00
上式表明,SFT正是序列的ZT在z=e""的值
(二)离散时间系统
1.系统函数的收敛域与系统因果性和稳定性
当且仅当系统函数H(z)的收敛域为小于单位圆的某个圆的园外时,系统是因
果稳定的。
2.系统函数的零极点分布与系统因果性和稳定性
若系统是因果稳定的,则H(z)的极点必定在单位圆内。
3.系统函数的零极点分布对系统频率响应特性的影响
1、对极点而言:当单位圆上的点转到某个极点附近时,|"(/。)|在这附近出
现峰值。极点越靠近单位圆,振幅特性的峰值越大,当极点出现在单位圆上时,
振幅特性将出现无穷大,系统不稳定。
2、对零点而言:当单位圆上的点转到某个零点附近时,|"(,力|在这附近出
现谷点。当零点出现在单位圆上时,振幅特性为零。零点可以位于单位圆外,不
影响稳定性。
两个概念一一
1、最小相位系统:系统H(z)的全部零极点都在单位圆内,某点在单位圆上逆时
针旋转一周时,系统的相位变化最小。
2、最大相位系统:H(z)的全部零点在单位圆外,系统的相位变化最大。
说明:处于坐标原点的零极点不影响系统的幅频响应:利用零极点分析系统的幅频响应,
仅对低阶系统有效。
(三)离散时间信号与模拟(连续)时间信号
1.时域关系
设连续时间信号五⑺,离散时间信号尤(〃),则
2.频域关系
18
乂(小儿皿ZX“"Q—,0)]
/7n=-00
在时域对信号抽样,其频域的特征就是频谱以采样频率为周期进行周期延拓。
一个域的离散必然导致另一个域的周期延拓
一个域的周期延拓必然导致另一个域的离散
对应变量的关系:0-单位:radQ-单位:Hz
co=dT
由于。《Q,,所以(91rax=QT=24
三、离散傅里叶变换(DFT)
(一)离散傅里叶级数变换(DFST)
说明:周期序列不满足绝对可和的条件,不适用于序列傅里叶变换的定义式,
但是它可以展开成离散傅里叶级数(DiscreteFourierSeries,DFS),利用离散傅
里叶级数可以得到周期序列的离散傅里叶变换表示式。
L定义
N-t
DFST:元(Q=Z双〃),一8<攵<8
z?=0
1NT
IDFST:'(〃)=—Z又(口卬消
N〃=()
注:1、周期单位复指数序列W*=/如,跖,*=ak
周期单位复指数序列对n、k而言都是以N为周期的,即
+N}k
W^"=%。-00<〃,左<8
卬产)=阅大,一00<〃,%<00
W产=w『,-g<n,k<8
2、周期为N的周期序列以〃)可以分解成N个周期复指数序列的和,这些周期复
指数序列的数字角频率为笆(%=0,1,2,…,N-l)周,它们的幅度和相位由离散
N
傅里叶级数外决定。
N
2.基本周期序列的离散傅里叶级数变换
时域序列离散傅里叶级数变换(DFST)
$(〃)1
1N$(k)
声,““N8{k-m)
eN
cos(2幻%〃/N)N[3(k-m)+3(k+m)]/2
sin(2^mn/N)-jN政k-m)-S(k+m)]/2
3.周期序列的离散傅里叶变换
X(e"")=?之也烦。-手)
/VJt=—oo
可类比信号系统中周期信号的傅里叶变换,具体推导过程见课本76页。
(二)离散傅里叶变换(DFT)
1.定义
NT
DFT:X(k)=^x(n)W^,0<k<N-l
n=0
1NT
IDFT:x(〃)=—£X(A)照成,0W〃4N—1
要点:
(1)DFT没有实际的物理含义,但是可以理解为SFT的等间隔采样,即
X(Q=X(e"°)|2.,0(kWN—1
(
0=Nk
(2)变换区间:[O,N-1],有限长N点
(3)变换结果:与序列长度N有关,当N足够大时,|X(Q|的包络趋近于
曲线
(4)频谱分析的意义:
|X(Q|表示?=(2万/N乂频点的幅度谱线,如果x(〃)是模拟信号的采样,采样间
2rr
隔为T,(o=QT=2兀门T,则k与相应的模拟频率的关系为:co=—k=2^fT
kNk
即人=_七一。对模拟频率域而言,N点DFT意味着频域采样间隔为一1-"z。所
kNTNT
以用DFT进行谱分析时,称尸=」-为频率分辨率。而NT表示时域采样的区间
NT
长度(即观察时间或记录长度NT),显然为了提高分辨率就必须是记录长度
足够大。
(5)DFT的隐含周期性
1)DFT是SFT的等间隔采样,而X(e%以2%为周期;
2)胱=叱尸曲的周期性
3)时域抽样,频域周期延拓;频域采样,时域周期延拓
2.DFT的主要性质
性质时域(%(〃)、y(n))频域(X(k)、Y(k))
线性性质axl(n)+bx2(n)aXt(k)+bX2(k)
11m
时域循环移位性质x((n+m))NRN(n)WNx(k)
频域循环移位性质X((Z+)))N%(幻
时域循环卷积Xj(n)0x2(«)xgx.k)
频域循环卷积尤1(〃)/(〃),Xi(k)BX式k)
复共辗序列的DFTx(〃)*X»(N-k)
共辗对称性%(〃)XR(Q
%,(〃)jX人k)
XR(〃)Xf)
A(«)X〃(k)
帕斯瓦尔定理N-liN—l
〃=oNk=o
3.基本序列的离散傅里叶变换
时域序列离散傅里叶级数变换(DFST)
3(九)1
%(〃)N33
j—mnN次k-ni)
eNRJn)
cos(2/rmn/N)RN(n)N[3(k-m)+^k-N+m)]/2
上一
sin(2乃〃2〃/N)RN(n)-jN[3(k-m)—6(N+m)]/2
4.频域采样定理
设序歹Ux(〃)的傅里叶变换为X(e"〃),在区间[0,2万)内对X(e"")进行N点等间
隔采样(采样间隔为2万/N)得到序列X(Z),且X(幻对应的IDFT为无w(〃),则
00
XN(〃)=Zx(n+rN)
r=—oo
这是因为,在频域内对X(eW)等间隔采样,导致时域序列x(〃)周期延拓,并
且在区间[0,2万)采样得到的序列X(幻的IDFT是原序列以N为周期进行周期延
拓后的主值序列。若序列的长度为M,那么只有当频域采样点数时,才
有4(〃)=双〃),此时才能由频域采样序列X(A)恢复X3")。
(三)连续信号傅里叶变换(CFT)、序列傅里叶变换(SFT)、离
散傅里叶级数变换(DFST)、离散傅里叶变换(DFT)的关系
截短,雕延拓驷5
抽样
—►也)侧-----xN(n)
V,八押期延拓而
人(阳一
Q=2乃/%取主值
各个变量对应关系:
左:0〜N-lk
数字角频率。:0〜2万21Z/N
数字频率F:O~1k/N
模拟频率f-.O~fskf」N
0rrl^
模拟角频率。:0〜—X
CD=C1TS,0)$=。>1=2兀于*,0=2nkiN
编者按:为什么要有DFT?
我们从外界接收到的信号都是连续信号,但是在现代人类都用计算机对信号
进行处理,而计算机只能识别离散的值,所以需要对接收到的连续信号进行采样
截短得到离散的序列。但是,一个域的离散必然导致另一个域的周期延拓,当对
时域的连续信号进行采样时,其频谱必然进行周期延拓,所以序列的傅里叶变换
是连续周期的,这样计算机就没法对其频谱进行分析。这时,对时域信号进行周
期延拓,又会使其频谱离散化。经过两个域的分别离散化和周期延拓,这时得到
的就是DFST的对应关系。那么,分别对两个域取主值,就可得到适合计算机处
理的时域和频域序列。DFT就应运而生。(一家之言,仅供参考)
(四)卷积的计算
1.循环卷积与线性卷积(有限长序列的卷积)
设有限长序列x(")的长度为N,〃(〃)的长度为M,它们线性卷积结果为力(〃),
长度为4=N+M-1;循环卷积结果为”(〃),长度为L。则两类卷积有如下对
应关系:(设)
(1)当L=N时
国(")+M(〃+N),0<n<M-2
I%(〃),M-l<n<N-l
(2)当LNL.时
(3)当时
/、\y,(n)+y,(n+L),0<n<Lg-L-l
M(〃),Lx-L<n<L-\
2.重叠保留法和重叠相加法(无限长序列得卷积)
(1)重叠保留法
基本思路:将两个序列中长度较长或无限长的序列均匀分段,计算各个有
限长的子序列与另一短序列的线性卷积,最后将结果重叠相加起来输出。(重叠
的是卷积结果)
设有限长序列〃(〃)的长度为M,x(〃)为无限长序列,
计算步骤:1)将x(〃)均匀分段,每段长度为N
x(")=Zx«(〃)
k=0
x(〃),kN<n<(k+\)N-l
Xk(n)=X(〃)RN5—kN)=\
0,else
2)计算每段子序列与短序列的线性卷积
设匕(〃)=队,即计算■(〃)与人(〃)的线性卷积:(〃)
3)将各子序列线性卷积的结果移位后相加得总输出
令%(〃)=叉(〃一ZN),则加)=£”(〃)
k=0
(2)重叠保留法
基本思路:将两个序列中长度较长或无限长的序列在时间上有重叠地分段,
计算各个有限长的子序列与另一短序列的线性卷积,最后保留每段结果中间N
个点,相加输出。(重叠的是较长的序列)
设有限长序列〃(〃)的长度为M,x(n)为无限长序列,
计算步骤:1)将x(〃)有重叠地分段(每一段由kN向前重叠M-1个点),每
段长度为N+M-1
x(〃),kN-M+l<n<(k+l)N-l
%(〃)=',
n(),else
2)计算每段子序列与短序列的线性卷积
设x'k(n)=xk(n+kN-M+1),即计算x'k(n)与h(n)的线性卷积
叉(〃),义(〃)的长度为N+2M-2,将前M-1个点去掉,后M-1
个点去掉,保留中间N个点得以(〃)
3)将各子序列线性卷积的结果移位后相加得总输出
即>(〃)=£”(〃)
k=0
说明:重叠保留和相加法必须掌握,公式可以不必记忆,明白其算法思想,会计
算即可。而且计算时注意三步走(写在卷子上),否则答案正确也没分(与数学
归纳法一样,有固定格式)。
(五)用DFT进行频谱分析的误差
1.泄漏现象
★生原因:用DFT进行分析时,隐含对序列在时域加窗截断,使得信号的原有
频率的能量向其他频率上泄漏
减少方法:(1)加大窗长,增加实际DFT计算的点数;
(2)变换时域所加窗函数的形式
2.栅栏现象
产生原因:DFT只计算。=2mt/N«=0,l,2,…,N—1的频谱
减少方法:在序列末尾加零以增加DFT的点数
3.混叠现象
产生原因:序列截断以及米样频率不完全满足采样定理
减少方法:以较高的采样频率对信号进行采样,之后序列通过数字低通滤波器,
降低采样频率后再进行DFT分析
4.DFT的分辨率:参数选择的一般原则:
a.若已知信号的最高频率防止混叠,选定采样频率X>2fmm
b.根据频率分辨率F,确定所需DFT的长度N=£//
c.和N确定以后,即可确定相应模拟信号的时间长度〃=£/N=NT,这里T是
采样周期。
(六)离散时间信号的抽取和内插
1.离散时间信号的整数倍抽取
时域:y(〃)=x(Dn)
iD-l.a)-2
频域:y(,3)=7Tzx(Jk)
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2026陕西西安西京初级中学教师招聘备考题库含答案详解(轻巧夺冠)
- 2026福建龙岩学院附属中学招聘心理健康教师的1人备考题库含答案详解(典型题)
- 2026年山东师范大学第二附属中学第二批公开招聘人员备考题库(4名)附答案详解(达标题)
- 2026河北秦皇岛市市直医疗卫生单位第二批招聘工作人员36人备考题库附答案详解
- 2026年马鞍山市和县文化旅游体育局度校园招聘备考题库带答案详解(完整版)
- 2026浙江丽水莲都电信业务项目部招聘智家工程师6人备考题库含答案详解(培优)
- 2026浙江宁波前湾新区招聘事业编制教师(第四批)24人备考题库及答案详解(名师系列)
- 2026青海海南州贵南县文化馆演职人员招聘4人备考题库及1套完整答案详解
- 2026北京顺义区卫生健康人才服务中心第一批招聘编外工作人员1人备考题库含答案详解(典型题)
- 2026陕西西安市未央区辛家庙社区卫生服务中心招聘4人备考题库及答案详解(各地真题)
- 2026年急性缺血性卒中患者早期管理-美国心脏协会指南
- 2026年台州临海市粮食收储有限公司招聘工作人员5人建设笔试备考题库及答案解析
- DB42∕T 1528-2026 湖北省农业用水定额
- 2026年两会政府工作报告精神宣讲《聚焦两会关注民生》
- 船舶代理公司考核制度
- 2025年建信期货招聘笔试真题及答案
- 2025年北京纪委监委公开遴选公务员笔试试题及答案解析
- GMP计算机系统用户权限管理操作规程
- 2025年7月陕西省普通高中学业水平合格性考试语文试题及答案
- 2025河南编导考试真题及答案
- DB44 07∕T 70-2021 地理标志产品 新会陈皮
评论
0/150
提交评论