数字信号处理知识点总结.docx

数字信号处理辅导一、离散时间信号和系统的时域分析(一)离散时间信号(1)基本概念信号：信号传递信息的函数也是独立变量的函数，这个变量可以是时间、空 间位置等。连续信号：在某个时间区间，除有限间断点外所有瞬时均有确定值。模拟信号：是连续信号的特例。时间和幅度均连续。离散信号：时间上不连续，幅度连续。常见离散信号-数字信号：幅度量化，时间和幅度均不连续。(2)基本序列(课本第710页)1, n>00, n<01, n>00, n<0u(n)= <f 1 n = 01)单位脉冲序列 3()= '2)单位阶跃序列0,X.a u(n)f 1 0<n< N-l3)矩形序列Rn()= '一一 4)实指数序列八 0, n<0,n>Nx()=5)正弦序列= Asin(g + e)6)复指数序列(3)周期序列1)定义：对于序列(")，若存在正整数N使%(") = %(" +N),-8八<8则称85)为周期序列，记为元()，N为其周期。注意正弦周期序列周期性的判定(课本第10页)2)周期序列的表示方法：a.主值区间表示法b.模N表示法3)周期延拓设x5)为N点非周期序列，以周期序列L对作x()无限次移位相加，即可得到 周期序列元,即00元二e x(n-iL)/=co当 L N N 时，x(n) = x(n)RN ()当 L v N 时，x() w x(n)RN (n)(4)序列的分解序列共朝对称分解定理：对于任意给定的整数M,任何序列15)都可以分解成关于c = A1/2共辗对称的序列乙和共辄反对称的序列乙之和，即时域序列离散傅里叶级数变换（DFST）5（八）1Rn5)N3(k)声加e n %5)NB(k-ni)cos(2;rm/ N)/?n()N6(k-m) + 6(k-N + m) / 2/ N)Rn ()-jNdk-/n)-3(k-N + m) / 24.频域采样定理设序列的傅里叶变换为X（&）,在区间。24）内对X（,g）进行N点等间 隔采样（采样间隔为2/N）得到序列X/）,且X（口对应的IDFT为/5），则、N（"）= Z M，2+rN）r二-00这是因为，在频域内对x（/。）等间隔采样，导致时域序列周期延拓，并 且在区间。2幻采样得到的序列X（Z）的IDFT是原序列以N为周期进行周期延 拓后的主值序列。若序列的长度为M,那么只有当频域采样点数时一,才 有赤（） = %（），此时才能由频域采样序列X（Q恢复X（e%。（三）连续信号傅里叶变换（CFT）、序列傅里叶变换（SFT）、离 散傅里叶级数变换（DFST）、离散傅里叶变换（DFT）的关系,八抽样4 （ M训也）叫 ACFTSFtJ周期延拓：,、卷积周期延拓 誓石 刈）加） * _a /（）|,取主值,sftT dfstTdft.» 1周暨，,八抽样4 （ M训也）叫 ACFTSFtJ周期延拓：,、卷积周期延拓 誓石 刈）加） * _a /（）|,取主值,sftT dfstTdft.» 1周暨，人（网A文网二人的Q=2/4取主值各个变量对应关系:%:0N 1k数字角频率0：02127ik/N数字频率/：01k/N模拟频率帆/N模拟角频率Q: 02兀九个fsco = C1TS, cos =。工=2兀于丁 , co = 271kl N 编者按：为什么要有DFT?我们从外界接收到的信号都是连续信号，但是在现代人类都用计算机对信号 进行处理，而计算机只能识别离散的值，所以需要对接收到的连续信号进行采样 截短得到离散的序列。但是，一个域的离散必然导致另一个域的周期延拓，当对 时域的连续信号进行采样时，其频谱必然进行周期延拓，所以序列的傅里叶变换 是连续周期的，这样计算机就没法对其频谱进行分析。这时，对时域信号进行周 期延拓，又会使其频谱离散化。经过两个域的分别离散化和周期延拓，这时得到 的就是DFST的对应关系。那么，分别对两个域取主值，就可得到适合计算机处 理的时域和频域序列。DFT就应运而生。(一家之言，仅供参考)(四)卷积的计算L循环卷积与线性卷积(有限长序列的卷积)设有限长序列()的长度为N, h(n)的长度为M,它们线性卷积结果为, 长度为4=N + M-1；循环卷积结果为”()，长度为L。则两类卷积有如下对 应关系：(设N2M)(1)当L = N时ycM =ycM =yl (n) + y,(" + N),0 < n < Af - 2y()， M -l<n< N-1(2)当时 o纥5)=切5)(3)当时K()二K()二>/() + %( + £), 0<n<L-L-l y ，-L<n<L-2.重叠保留法和重叠相加法(无限长序列得卷积)(1)重叠保留法基本思路：将两个序列中长度较长或无限长的序列均匀分段，计算各个有 限长的子序列与另一短序列的线性卷积，最后将结果重叠相加起来输出。(重叠 的是卷积结果)设有限长序列h(n)的长度为M, x()为无限长序列，计算步骤：1)将x()均匀分段，每段长度为Ncok=0鼠()，kN<n<(k + l)N-lXk 5) = x(n)RN(n-kN) = <| 0,else2)计算每段子序列与短序列的线性卷积设x(«) = xkn-kN),即计算%(八)与h(ri)的线性卷积y()3)将各子序列线性卷积的结果移位后相加得总输出8令"()=义(一攵N)，贝1J y()= Z”(力)k=0(2)重叠保留法基本思路：将两个序列中长度较长或无限长的序列在时间上有重叠地分段, 计算各个有限长的子序列与另一短序列的线性卷积，最后保留每段结果中间N 个点，相加输出。(重叠的是较长的序列)设有限长序列力5)的长度为M, xn为无限长序列，计算步骤：1)将5)有重叠地分段(每一段由kN向前重叠M-1个点)，每 段长度为N+M-1kN M + 14n4(k + l)N l0,else2)计算每段子序列与短序列的线性卷积设乂5) = xk(n + kN-M +1),即计算%()与h(n)的线性卷积 叉5)，乂()的长度为N+2M-2,将前M-1个点去掉，后M-1 个点去掉，保留中间N个点得力()3)将各子序列线性卷积的结果移位后相加得总输出00即、() = £”()k=0说明：重叠保留和相加法必须掌握，公式可以不必记忆，明白其算法思想，会计 算即可。而且计算时注意三步走(写在卷子上)，否则答案正确也没分(与数学 归纳法一样，有固定格式)。(五)用DFT进行频谱分析的误差1泄漏现象产生原因：用DFT进行分析时，隐含对序列在时域加窗截断，使得信号的原有 频率的能量向其他频率上泄漏减少方法：(1)加大窗长，增加实际DFT计算的点数；(2)变换时域所加窗函数的形式.栅栏现象产生原因：DFT只计算3 = 2/N,左=0,1,2,川-1的频谱减少方法：在序列末尾加零以增加DFT的点数2 .混叠现象产生原因：序列截断以及米样频率不完全满足米样定理减少方法：以较高的采样频率对信号进行采样，之后序列通过数字低通滤波器, 降低采样频率后再进行DFT分析.DFT的分辨率：参数选择的一般原则：a.若已知信号的最高频率防止混叠，选定采样频率£22ymaxb.根据频率分辨率F,确定所需DFT的长度N =工/尸c.和N确定以后，即可确定相应模拟信号的时间长度7； = £/N = NT,这里T是 采样周期。(六)离散时间信号的抽取和内插.离散时间信号的整数倍抽取时域：yri) = x(Dn)i Z)-l兀 k频域：y(")=Zx(Jk)Dh。整数倍抽取将导致数字频谱的展宽.离散时间信号的整数倍内插朗” ,、吊7)， = 0，±/，±2/,时域：口() = ' I0,else频域：V(&) = X(eW)序列相邻采样点之间插零将导致数字频谱压缩说明：即使抽取和内插的公式记不住，也要学会画图分析其过程各种域和各种变换的关系总结如图1一5所示,取取OFSIDFS向 期 越 拓采样数字频域期散周期时域工取主值周期模拟时域傅氏交换模拟频域高散时蔡傅氏变换XgDillIDTFT佰 区 同离散画期领域x()= x(ri) + x(ri)- 00 c < ooC-Cz/(")= g x(n) - x* (M - n)并且xe (n) = x(n) + x* (M - ri)(4)序列的运算1)基本运算运导性质描述序列相乘y()= % ()2 ()y =ax(n)序列相加y(n) = X()+ %2()序列翻转y() = x(-力)(将x()以纵轴为对称轴翻转)尺度变换y()= x(mn)(序列x()每隔m-1点取一点形成的序列)用单位脉冲 序列表示Q0x(")= Z x(i)b(一，)/=002)线性卷积：将序列()以y轴为中心做翻转，然后做m点移位，最后与()对应点相乘求和翻转、移位、相乘、求和8定义式:y(n) = Z F (加)/(. 一 m)- x (九)* Z()"7=-CO线性卷积的计算：A、图解B、解析法C、不进位乘法(必须掌握)3)单位复指数序列求和(必须掌握)N-1£/所= n=()N-1£/所= n=()_ e-论Ne-j(oN/2 jcoN/2 _ jcoNH)jsNi2 dNQ _jsN/2 ) / (2 j)e-j(0/ej0)/2 - /2)e-jco/2 (ej(t)/2 - e-j(0f2 )/(2，)-X/V-D/2 Csin(oN / 2)sin(69/2)如果=2兀k/N,那么根据洛比达法则有sin(3 N / 2) = Nb 也=0)(或 NS(N)(Z = N)sin(69/2)可以结合作业题3. 22进行练习(5)序列的功率和能量能量：E=£|M)72=-001 N功率：P=lim-|x(h)|22N + 16(6)相关函数一一与随机信号的定义运算相同(二)离散时间系统.系统性质(1)线性性质定义：设系统的输入分别为15)和82()，输出分别为y和()，即yi(n) = Txl(n), y2(n) = Tx2(ri)统的输对于任意给定的常数。、b,下式成立y() = Taxx (n) + bx2 (n)= ayi(ri) + by2(ri)则该系统服从线性叠加原理，为线性系统，否则为非线性系统。判定系统的线性性质时，直接用定义(2)时不变性质统的如果系统对输入信号的运算关系在整个运算过程中不随时间变化，则称该系 统是时不变系统。即对任意给定的整数i,若下式成立：yn -i) = Tx(n -/)则称该系统为时不变系统，否则为时变系统。判定系统的时不变性质时，直接用定义(3)系统的因果性定义：如果系统n时刻的输出序列只取决于n时刻及以前的输入序列，而与n时 刻以后的输入序列无关，则称该系统具有因果性质，即系统是因果系统，否则是 非因果系统。离散时间LTI系统具有因果性的充要条件是：系统的单位脉冲响应力伽)满足h(n) = 0, 。(4)系统的稳定性定义：对任意有界的输入，系统的输出都有界，则该系统是稳定的，否则是不稳 定的。离散时间LTI系统具有因果性的充要条件是：系统的单位脉冲响应5)满足绝00对可和,即/=O0(5)对离散时间LTI系统的描述(1)时域：差分方程1 2) Z域：系统函数”(z)2.信号过系统y(n) = h(n) * x(n)用线性卷积的相关知识计算，信号系统学的基本性质可以套用二、离散时间信号和系统的频域分析(一)离散时间信号1.序列傅里叶变换(Sequence Fourier Transform)(即本书中的离散时间信号 的傅里叶变换)(1)定义8SFT：= SFTx(n) =x(n)e'/7y-oo< < ooISFT： x()= ISFTX(ejeo) =Xd,s < < oo说明：1、物理意义：序列傅里叶变换本质上是序列的一种分解，它将一般序列分解为 无穷多个数字角频率-肛划中的复指数序列。称X(e，0)为序列(八)的频谱，其 模I X(") |称为幅频特性，其幅角argX(")=伙。)称为相频特性。2、尽管序列伽)是离散时间信号，但它的序列傅里叶变换对数字角频率而言 却是连续函数，因此，序列15)的傅里叶变换是连续的。003、X("3+2%)= z= X(")/?=-<»由上式可知，序列傅里叶变换X(e，。)是以2为周期的周期函数，其原因正是由 于*对。而言以2万为周期，即数字角频率相差2万的所有单位复指数序列等 价。因此，对的所有单位复指数序列只有一个周期。对于离散时间信 号，由于的周期性，使得口 = 0或2万的整数倍都表示信号的直流分量，而的奇数倍表示信号的最高频率。（2）性质名称性质描述线性性质SFTaxx () + 法2() =(n) + /?*5FTx2(/i)时移性质SFTx(n - m) = /8n频移性质SFTejn.x(n) =共朝对称性质SFTxR(n) = Xe(ej6),SFTjXl(n) = X")S尸取()=ReX(0), SFTxo(n) = j!mX(en线性卷积性质SFTx(n) * y(n) = SFTx(n)SFTy(n)帕斯瓦尔定理001二-oo乙兀相乘性质SFTx(n)y(n)=X(ej0)Y(ej(c°-0)d0序列乘以nSFT/?.x(n) = jdX(ejC0)/ dco（3）基本序列的傅里叶变换序列傅里叶变换gl)112痛(Rn5)"j/(NT)/2 sin侬)/ sin(-)anu(n) (| a < 1)(l-f )T*3 （2万/怒为有理数）245（69一 4）cos。源兀/ 4为有理数）7lS(69 COq ) + 3(69 + CO )sin0n（2兀/ 4为有理数）一）万（69 例）- 3（69+4）u(n)（12一加尸+痛（2. Z变换（不熟悉的复习信号系统相关内容，或本书2.3相关内容）（1）定义8ZT： X(z) = ZTx(n) = xnz-n Rx_ <| z |< Rx+n=-<x)IZT： x(«) = /ZTX(z) = (f) XCz)zdz R <z<R +2兀 j Jc(2)性质课本49页表(3)收敛域与基本序列Z变换课本45页表、表3.离散时间信号Z变换与SFT的关系Z变换是由SFT推广得到的，反过来,如果某序列的Z变换的收敛域包括z = * , 则也可以通过ZT求得序列的SFT。即00X(z)&= 2>5)/& = 乂(0) n=-<x>上式表明，SFT正是序列的ZT在z = /。的值(二)离散时间系统.系统函数的收敛域与系统因果性和稳定性当且仅当系统函数H(z)的收敛域为小于单位圆的某个圆的园外时，系统是因 果稳定的。1 .系统函数的零极点分布与系统因果性和稳定性若系统是因果稳定的，则H(z)的极点必定在单位圆内。2 .系统函数的零极点分布对系统频率响应特性的影响1、对极点而言：当单位圆上的点转到某个极点附近时，|"(，3)|在这附近出 现峰值。极点越靠近单位圆，振幅特性的峰值越大，当极点出现在单位圆上时， 振幅特性将出现无穷大，系统不稳定。2、对零点而言：当单位圆上的点转到某个零点附近时，| "()0)1在这附近出 现谷点。当零点出现在单位圆上时，振幅特性为零。零点可以位于单位圆外，不 影响稳定性。两个概念1、最小相位系统：系统H(z)的全部零极点都在单位圆内，某点在单位圆上逆时 针旋转一周时，系统的相位变化最小。2、最大相位系统：H(z)的全部零点在单位圆外，系统的相位变化最大。说明：处于坐标原点的零极点不影响系统的幅频响应；利用零极点分析系统的幅频响应, 仅对低阶系统有效。(三)离散时间信号与模拟(连续)时间信号.时域关系设连续时间信号4，离散时间信号则(/) = %" t=nT1 .频域关系8X(小儿皿=下 E XJ(Q-mQJ/ 77：=-00在时域对信号抽样，其频域的特征就是频谱以采样频率。，为周期进行周期延拓。一个域的离散必然导致另一个域的周期延拓一个域的周期延拓必然导致另一个域的离散对应变量的关系：口-单位：rad。-单位：Hza)= QT由于Q W Qs.,所以0= 2三、离散傅里叶变换(DFT)(一)离散傅里叶级数变换(DFST)说明：周期序列不满足绝对可和的条件，不适用于序列傅里叶变换的定义式, 但是它可以展开成离散傅里叶级数(Discrete Fourier Series, DFS),利用离散傅 里叶级数可以得到周期序列的离散傅里叶变换表示式。L定义NTDFST：又() = £35)叫/，00<%<871=()IDFST：IDFST：1 N-1元5)=不2又伏)跖欣.2.24一/nk.1nk注：1、周期单位复指数序列nV丁卜n周期单位复指数序列对n、k而言都是以N为周期的，即=Wnk 9_00</1<00=卬丁 ,8<n,k<s阅成+N)=w，_8<z<002、周期为N的周期序列无()可以分解成N个周期复指数序列的和，这些周期复指数序列的数字角频率为里（2 =。】2,双-1）周，它们的幅度和相位由离散 N傅里叶级数任2决定。2.基本周期序列的离散傅里叶级数变换3.周期序列的离散傅里叶变换时域序列离散傅里叶级数变换（DFST）$5)11N$(k)jnin e nN$(k - m)cos(2mzi / N)N6(k-m) + (k + m)/2sinQ/rmn / N)-jNS(k - m)- 3* + m) / 227r2兀xd）=$£x*）/-千 k）IV k =oo1V可类比信号系统中周期信号的傅里叶变换，具体推导过程见课本76页。（二）离散傅里叶变换（DFT）L定义N-1DFT： X(Z) = £x()叫;，0WZWN-172=01 N-1IDFT： x(n) = YX(k-nk ,Q<n<N-l要点:（1） DFT没有实际的物理含义，但是可以理解为SFT的等间隔采样，即X（Q = X（e加）| 27r ,0<k<N-lco-kN（2）变换区间：O,N-1,有限长N点（3）变换结果：与序列长度N有关，当N足够大时，|X（左）|的包络趋近于|X（0）曲线（4）频谱分析的意义:|X（左）|表示以=（2»/N）左频点的幅度谱线，如果x（）是模拟信号的采样，采样间隔为T, co = QT = 2兀f IT ,则k与相应的模拟频率的关系为：cok=-k = 27ifkT即人二-七。对模拟频率域而言，N点DFT意味着频域采样间隔为一匚Hz。所 a NTNT以用DFT进行谱分析时，称尸=,为频率分辨率。而NT表示时域采样的区间 NT长度（即观察时间或记录长度NT）,显然为了提高分辨率就必须是记录长度 足够大。（5） DFT的隐含周期性1 ） DFT是SFT的等间隔采样，而X（e，。）以2%为周期；2）阅=叱广'必的周期性3）时域抽样，频域周期延拓；频域采样，时域周期延拓2 .DFT的主要性质3.基本序列的离散傅里叶变换性质时域（l（几）、y（ji）频域(X(k)、Y(k)线性性质咐（）+ 也（）aX,(k) + bX2(k)时域循环移位性质x(n + m)NRN(n)W/"X(左)频域循环移位性质Wlx(n)X(k + l)NRN(k)时域循环卷积%（）区工2（）X1伏冰2伏）频域循环卷积% ()/()拉心）”收）复共轨序列的DFTx(n)*X*(N-k)共辗对称性“X£k)jX,(k)XrXq,(k)jx1(n)X°p(k)帕斯瓦尔定理Nli N-