数字信号处理2-1
第十一讲 数字信号处理-DIT—FFT
r 0
r 0
r 0
2 x1 (r )(WN ) rk j 2 2 N j
r 0
2 k x2 (r )(WN ) rk WN
2 WN e
e
2 N /2
WN / 2
3. 求出子序列的DFT
X(k )
N / 2 1 r 0 rk x1 (r )WN / 2 rk k x2 (r )WN / 2WN r 0 N / 21
5. 求出后半部的表示式
rk r WN / 2 WN (/k2 N / 2)
k X(k ) X1 (k ) WN X 2 (k )
N / 21 N / 21 N r ( N / 2 k ) rk X 1 ( k ) x1 (r )WN / 2 x1 (r )WN / 2 X 1 (k ) 2 r 0 r 0 N / 21 N / 21 N r ( N / 2 k ) rk X 2 ( k ) x2 (r )WN / 2 x2 (r )WN / 2 X 2 (k ) 2 r 0 r 0
看出:后半部的k值所对应的X1(k),X2(k)则完全重复了前半部分的 k值所对应的X1(k),X2(k)的值。
又 W
(N / 2 k ) N
W
N /2 N
W W
k N
k N
N k 后半部分: X (k ) X 1 (k ) WN X 2 (k ) 2
6. 结论2
频域中的N个点频率成分为:
作图要素: (1)左边两路为输入 (2)右边两路为输出 (3)中间以一个小圆表示加、 减运算(右上路为相加输出、 右下路为相减输出)
数字信号处理第2章
Z变换与拉氏变换的关系:
这一关系实际上是通过 到了Z平面。
若将Z平面用极坐标表示
标表示
,代入
将S平面的函数映射
,S平面用直角坐 ,得:
上述关系表明: z 的模 r 仅与 s 的实部 相对应, z 的幅角 则仅与 s 的虚部 对应。
映射关系:
Z变换与拉氏变换的关系
0 0,2 (S平面实轴映射到Z平面的正实轴)
解:
,求它的傅立叶变换。
其幅度谱和相位谱分别为:
典型例题
❖ 例2 已知序列的傅立叶变换如下,求它的反变换。
解:
显然序列 h(n)不是绝对可和的,而是平方可和 的 ,但其依然存在傅立叶变换。 Parseval定理
典型例题
❖ 例3 证明复指数序列 x(n) e j0n 的傅立叶变换为:
证:根据序列的傅立叶反变换定义,利用冲击函 数 的性质,有:
即序列绝对可和
某的有 立些序些叶既列序变不,列换满若虽依足引然然绝入不存对频满在可 域足。和的以见的冲上后条击条例件函件。也数,不但满满,足足其平平傅方方立可可叶和和变条,换件其傅
也存在。如
、某些周期序列,见后例。
序列傅立叶变换的定义
5.常用序列的傅立叶变换
序列
(n)
傅立叶变换
1
1
典型例题
❖ 例1 已知
A形k(式k=求0,X取1(…:z),N)B,(此z) A( z )
时
为了方bi 便z i通常利用
i0
N
1 ai z i
X(z)/z的
i 1
若序列为因果序列,且N≥M,当X(z)的N个极点都是单
极点时,可以展开成以下的部分分式的形式:
则其逆Z变换为:
数字信号处理(方勇)第二章习题答案
2-1 试求如下序列的傅里叶变换: (1))()(01n n n x -=δ (2))1(21)()1(21)(2--++=n n n n x δδδ (3)),2()(3+=n u a n x n10<<a(4))4()3()(4--+=n u n u n x(5)∑∞=-⎪⎭⎫⎝⎛=05)3(41)(k nk n n x δ(6)()6cos ,14()0,n n x n π⎧-≤≤=⎨⎩其他解: (1) 010()()j n j j nn X e n n ee ωωωδ∞--=-∞=-=∑(2) 2211()()122j j nj j n X e x n e e e ωωωω∞--=-∞==+-∑ωsin 1j +=(3) 2232()(2)1j j nj nn j nj n n a e X e a u n ea eaeωωωωω-∞∞---=-∞=-=+==-∑∑, 10<<a(4) []4()(3)(4)j j nn X e u n u n eωω∞-=-∞=+--∑∑-=-=33n nj e ω∑∑==-+=313n n j n nj e eωωωωωωωj j j j j e e e e e --+--=--111134=⎪⎭⎫ ⎝⎛⎪⎭⎫⎝⎛=----ωωωω21sin 27sin 1137j j j e ee(5) 3350011()(3)44n kj jn j k n k k X e n k e e ωωωδ∞∞+∞--=-∞==⎛⎫⎛⎫=-= ⎪ ⎪⎝⎭⎝⎭∑∑∑∑∞+=--⎪⎭⎫⎝⎛-=⎪⎭⎫⎝⎛=033411141k j kj e e ωω(6) 44336441()cos 32j j j jn jn n n X e nee e e ππωωωπ---=-=-⎛⎫==+ ⎪⎝⎭∑∑994()()4()()3333001122j j n j j n n n e e e e ππππωωωω--++===+∑∑ ()9()9334()4()33()()3311112211j j j j j j e e e e e e ππωωππωωππωω-+-+-+⎡⎤⎡⎤--⎢⎥⎢⎥=+⎢⎥⎢⎥++⎢⎥⎢⎥⎣⎦⎣⎦2-2 设信号}1,2,3,2,1{)(---=n x ,它的傅里叶变换为)(ωj e X ,试计算(1)0()j X e (2)()j X ed πωπω-⎰(3)2()j X e d πωπω-⎰。
山东大学 DSP数字信号处理PPT 第二章z变换 习题讲解
1 1 z2
X z
4
1
1 4
z
2
1
5 4
z 1
3 8
z
2
解:对X z的分子和分母进行因式分解,得
1 1 z2
X z
4
1
1 4
z
2
1
5 4
z 1
3 8
z
2
1
1 2
z 1
1
1 2
z 1
1
1 4
z 2
1
1 2
z 1
1
3 4
z 1
1 1 z1
2
1
1 2
jz
1
1
1 2
2-13 研究一个输入为x(n)和输出为 y(n)的 时域线性离散移不变系统,已知它满足
y(n 1) 10 y(n) y(n 1) x(n) 3
并已知系统是稳定的。试求其单位抽样 响应。
y(n 1) 10 y(n) y(n 1) x(n) 3
解:对差分方程两边取z变换
z1Y (z) 10 Y (z) zY (z) X (z) 3
在围线c外有单阶极点 z 1/ 4,
且分母阶次高于分子阶次二阶以上
x(n)
Re
s
F
(
z) z 1 /
4
z
1/
4
(
z 2)zn1 z 1/4
z 1 /
4
7 4
1 4
n 1
7
4n
x(n) 8 (n) 7 4n u(n 1)
j Im[z]
C
1/ 4
0
Re[z]
③部分分式法
X (z) z
jz
数字信号处理 第二章习题
1 为因果序列,故收敛域为: z 2
8
(2) (n n0 ) n0 0
解:
X ( z)
n
x(n) z n
n
(n n0 ) z n
X ( z) z
n0
1 n n0 (n n0 ) 0 other
1 n0 z
z 0.5 左边序列 0.5 z 2 双边序列 右边序列 z 2
16
采用围线积分法求解:
3 2 X ( z) 1 1 0.5 z 1 2 z 1 3(1 2 z 1 ) 2(1 0.5 z 1 ) 5 7 z 1 1 1 (1 0.5 z )(1 2 z ) (1 0.5 z 1 )(1 2 z 1 )
z1 1, z2 2
X(z)的收敛域为
左边序列 z 1 1 z 2 双边序列 z 2 右边序列
24
F ( z) X ( z) z
n 1
z ( z 3) ( z 3) n 1 z zn ( z 1)( z 2) ( z 1)( z 2)
z 2
21
当收敛域为: z 2 0.5
1 n n 1 x(n) 3( ) u (n) 2 u (n 1) 2
22
收敛域为: z 2
右边序列
n 0 ,围线c内有2个1阶极点
x(n) Re s[( z 0.5) F ( z), 0.5] Re s[( z 2) F ( z), 2] ( z 0.5) 5z 7 zn ( z 0.5)( z 2) ( z 2)
双边序列
n 0 ,围线c内有1个1阶极点
数字信号处理习题解答
数字信号处理习题(xítí)解答第1-2章:1. 判断下列(xiàliè)信号是否为周期信号,若是,确定其周期。
若不是,说明(shuōmíng)理由(1)f1(t) = sin2t + cos3t(2)f2(t) = cos2t + sinπt2、判断下列序列是否为周期(zhōuqī)信号,若是,确定其周期。
若不是(bùshi),说明理由(1)f1(k) = sin(3πk/4) + cos(0.5πk)(2)f2(k) = sin(2k)(3)若正弦序列x(n)=cos(3πn /13)是周期的, 则周期是N=3、判断下列信号是否为周期信号,若是,确定其周期; 若不是,说明理由(1)f(k) = sin(πk/4) + cos(0.5πk)(2)f2(k) = sin(3πk/4) + cos(0.5πk)解1、解β1 = π/4 rad,β2 = 0.5π rad 由于2π/ β1 = 8 N1 =8,N2 = 4,故f(k) 为周期序列,其周期为N1和N2的最小公倍数8。
(2)β1 = 3π/4 rad,β2 = 0.5π rad由于2π/ β1 = 8/3 N1 =8, N2 = 4,故f1(k) 为周期序列,其周期为N1和N2的最小公倍数8。
4、画出下列函数的波形(1).(2).解5、画出下列函数的波形x(n)=3δ(n+3)+δ(n+1)-3δ(n-1)+2δ(n-2)6. 离散线性时不变系统单位阶跃响应,则单位响应=?7、已知信号(xìnhào),则奈奎斯特取样(qǔyàng)频率为( 200 )Hz。
8、在已知信号(xìnhào)的最高频率为100Hz(即谱分析范围(fànwéi))时,为了避免频率(pínlǜ)混叠现象,采样频率最少要200 Hz:9. 若信号的最高频率为20KHz,则对该信号取样,为使频谱不混叠,最低取样频率是40KHz10、连续信号:用采样频率采样,写出所得到的信号序列x(n)表达式,求出该序列x(n) 的最小周期解:,11、连续信号:用采样频率100s f Hz = 采样,写出所得到的信号序列x(n)表达式,求出该序列x(n) 的最小周期长度。
(完整word版)数字信号处理答案第二章
第二章2.1 判断下列序列是否是周期序列。
若是,请确定它的最小周期.(1)x (n )=Acos(685ππ+n ) (2)x (n)=)8(π-ne j(3)x (n)=Asin(343ππ+n ) 解 (1)对照正弦型序列的一般公式x (n )=Acos (ϕω+n ),得出=ω85π。
因此5162=ωπ是有理数,所以是周期序列。
最小周期等于N=)5(16516取k k =。
(2)对照复指数序列的一般公式x (n )=exp[ωσj +]n,得出81=ω。
因此πωπ162=是无理数,所以不是周期序列。
(3)对照正弦型序列的一般公式x (n)=Acos(ϕω+n ),又x (n)=Asin (343ππ+n )=Acos (-2π343ππ-n )=Acos(6143-n π),得出=ω43π.因此382=ωπ是有理数,所以是周期序列。
最小周期等于N=)3(838取k k =2.2在图2.2中,x (n )和h(n)分别是线性非移变系统的输入和单位取样响应。
计算并列的x (n )和h (n)的线性卷积以得到系统的输出y(n ),并画出y(n)的图形。
(a)1111(b)(c)111110 0-1-1-1-1-1-1-1-1222222 33333444………nnn nnnx(n)x(n)x(n)h(n)h(n)h(n)21u(n)u(n)u(n)a n ===22解 利用线性卷积公式y(n )=∑∞-∞=-k k n h k x )()(按照折叠、移位、相乘、相加、的作图方法,计算y(n)的每一个取样值。
(a ) y (0)=x (O)h (0)=1y (l )=x (O )h(1)+x (1)h (O)=3y (n)=x(O)h (n )+x (1)h(n-1)+x(2)h (n —2)=4,n ≥2 (b) x(n )=2δ(n )-δ(n-1)h(n)=-δ(n)+2δ(n —1)+ δ(n —2)y(n )=-2δ(n)+5δ(n —1)= δ(n-3) (c ) y (n )=∑∞-∞=--k kn k n u k u a)()(=∑∞-∞=-k kn a=aa n --+111u (n )2。
(完整word版)数字信号处理第二章习题解答
数字信号处理第2章习题解答2.1 今对三个正弦信号1()cos(2)a x t t π=,2()cos(6)a x t t π=-,3()cos(10)a x t t π=进行理想采样,采样频率为8s πΩ=,求这三个序列输出序列,比较其结果。
画出1()a x t 、2()a x t 、3()a x t 的波形及采样点位置并解释频谱混淆现象。
解:采样周期为2184T ππ== 三个正弦信号采样得到的离散信号分别表示如下:1()cos(2)cos()42a n x n n ππ=⋅=2()cos(6)cos()42a n x n n ππ=-⋅=-3()cos(10)cos()42a n x n n ππ=⋅=输出序列只有一个角频率2π,其中1()a x n 和3()a x n 采样序列完全相同,2()a x n 和1()a x n 、3()a x n 采样序列正好反相。
三个正弦信号波形及采样点位置图示如下:tx a 1(t )tx a 2(t )tx a 3(t )三个正弦信号的频率分别为1Hz 、3Hz 和5Hz ,而采样频率为4Hz ,采样频率大于第一个正弦信号频率的两倍,但是小于后两个正弦信号频率的两倍,因而由第一个信号的采样能够正确恢复模拟信号,而后两个信号的采样不能准确原始的模拟信号,产生频谱混叠现象。
2.3 给定一连续带限信号()a x t 其频谱当f B >时,()a X f 。
求以下信号的最低采样频率。
(1)2()a x t (2)(2)a x t (3)()cos(7)a x t Bt π解:设()a x t 的傅里叶变换为()a X j Ω(1)2()a x t 的傅里叶变换为22()[()]Ba a BX j X j d ππωωω-⋅Ω-⎰因为22,22B B B B πωππωπ-≤≤-≤Ω-≤ 所以44B B ππ-≤Ω≤即2()a x t 带限于2B ,最低采样频率为4B 。
数字信号处理课后习题答案
(修正:此题有错,
(3)系统的单位脉冲响应 而改变,是两个复序列信号之和)
(4)
(修正: 随上小题答案
(修正:此图错误,乘系数应该为 0.5,输出端 y(n)应该在两个延迟器 D 之间)
1-25 线性移不变离散时间系统的差分方程为
(1)求系统函数 ; (2)画出系统的一种模拟框图; (3)求使系统稳定的 A 的取值范围。 解:(1)
(2)
(3)
解:(1)
(2)
(3)
1-7 若采样信号 m(t)的采样频率 fs=1500Hz,下列信号经 m(t)采样后哪些信号不 失真? (1) (2) (3) 解:
(1)
采样不失真
(2)
采样不失真
(3)
,
采样失真
1-8 已知
,采样信号 的采样周期为 。
(1) 的截止模拟角频率 是多少?
(2)将 进行 A/D 采样后, 如何?
(3)最小阻带衰减 5-4
由分式(5.39)根据 A 计算 ,如下: 由表 5.1 根据过度带宽度 计算窗口:
单位脉冲响应如下:
单位脉冲响应如下:
其中 为凯泽窗。 5-5 答:减小窗口的长度 N,则滤波器的过度带增加,但最小阻带衰减保持不变。 5-6:图 5.30 中的滤波器包括了三类理想滤波器,包括了低通,带通和高通,其响应的单位
(1)
,
(2)
1-18 若当 时
;时
(1)
,其中
(2) 证明:
,收敛域
,其中 N 为整数。试证明: ,
(1) 令 其中
,则 ,
(2)
,
1-19 一系统的系统方程及初时条件分别如下: ,
(1)试求零输入响应 ,零状态响应 ,全响应 ; (2)画出系统的模拟框图 解: (1)零输入响应
数字信号处理-第2章第1讲 离散时间信号和离散时间系统
当a>1时 当-1<a<0时 当a< -1时
2.2 常用序列
5、正弦序列
x(n) Asin(n )
x(n) xa (t) tnT Asin(nT ) T / fs 2 f / fs 单位rad, 单位rad / s
6、复指数序列
一阶后向差分: y(n) y(n) y(n 1) 二阶后向差分: 2 y(n) y(n) y(n 1)
y(n) 2 y(n 1) y(n 2) 用延时算子:Dy(n) y(n 1) y(n) y(n) Dy(n) (1 D) y(n) 1 D 2 y(n) y(n) y(n 1) (1 D) y(n) (1 D)Dy(n) (1 D)2 y(n)
卷积和
卷积和的定义
1. 交换律 2. 结合律
y(n) x(k)h(n k) x(n) h(n) k
y(n) h(n)x(n k) h(n) x(n) k
y(n) [x(n) h1(n)]*h2(n)
[x(n) h2(n)]*h1(n) x(n) [h1(n)*h2(n)]
线性非移变系统稳定的充要条件是满足绝对可 和的条件:
S h(n) n
证明:
(1)充分性
当 x(n) M得
y(n) h(k)x(n k) h(k) x(n k)
k
k
M h(k) 得证 k
(2)必要性
x(n) e( j)n
数字频率又叫归一化频率
x(n) en cos(n) jen sin(n)
《信号、系统与数字信号处理》第二章 连续时间信号与系统的频域分析
0 21
/4
/2
(b)相位图
图2.1-2例2.1-2的频谱图
二、指数形式的傅里叶级数
利用欧拉公式将三角形式的傅里叶级数,表示为 复指数形式的傅氏级数
其中
f t F n1 e jn1t
n
F n1
1 T
t0 T t0
f t e jn1tdt
F n1 是复常数,通常简写为 Fn 。
21t
5
4
2
sin
1t
1 2
sin
31t
解:将 f t 整理为标准形式
f
(t)
1
2cos 1t来自4cos 21t
5
4
1 2
cos
31t
2
1
2
cos
1t
4
cos
21t
4
1 2
cos
31t
2
振幅谱与相位谱如图2-1所示。
cn
2
1
1
1/2
0 1 21 31
(a) 振幅图
n
/4
31
第二章 连续时间信号与系统的频域分析 ——Fourier变换
2. 1 周期信号的傅里叶级数分析 2. 2 非周期信号的频谱--傅里叶变换 2. 3 傅里叶变换的性质及定理 2. 4 系统的频域分析方法 2. 5 无失真传输系统与滤波
LTI系统分析的一个基本任务,是求解系统对任意 激励信号的响应,基本方法是将信号分解为多个基本信 号元。
一、三角形式傅里叶级数
周期信号: f t f t nT
其中
T
是信号的最小重复时间间隔,f1
1 是信号的基波频率。 T
若 f t 满足狄里赫利条件,则 f t 可以展开为三角形
数字信号处理——第2章 离散时间傅里叶变换与Z变换
• 总结:
①序列ZT的收敛域以极点为边界(包含0 和 ②收敛域内不含任何极点,可以包含0 ③相同的零极点可能对应不同的收敛域,即: 不同的序列可能有相同的ZT ④收敛域汇总:右外、左内、双环、有限长z平面
)
常见典型序列z变换
序列 Z变换 收敛域
z a
z b
注意:只有z变换和它的收敛域两者在一起才和序列相对应。 其它序列见P54: 表2-1 几种序列的z变换
2.3
z反变换
Z反变换: 从X(z)中还原出原序列x(n)
X ( z ) ZT [ x ( n)]
n
x (n) z n
实质:求X(z)幂级数展开式
Z反变换的求解方法: 留数定理法
部分分式法
长除法
1. 留数定理法
根据复变函数理论,可以推导出
x ( n)
1 2 j
X ( z ) z n 1dz
1 1 3z 1
n
z 2
2 n u ( n)
z 3
3
n
n
u (n 1)
x n 2 u n 3 u n 1
3. 幂级数法(长除法)
如果序列的ZT能表示成幂级数的形式,则序列x(n) 是幂 级数 说明: ①这种方法只对某些特殊的ZT有效。 ②如果ZT为有理函数,可用长除法将X(z)展开成幂级 数。 若为右边序列(特例:因果序列),将X(z)展开成负幂 级数; 若为左边序列(特例:反因果序列),将X(z)展开成正 幂级数; 中
z z 1 1 X z 1 z 2 z 3 1 2z 1 3 z 1
1 ZT [a u (n)] z a 1 1 az 1 n ZT [a u (n 1)] z a 1 1 az
数字信处理课后习题matlab作业
数字信号处理MATLAB习题数字信号处理MATLAB 习题M1-1 已知1()cos(6)g t t π=,2()cos(14)g t t π=,3()cos(26)g t t π=,以抽样频率10sam f Hz =对上述三个信号进行抽样。
在同一张图上画出1()g t ,2()g t 和3()g t 及抽样点,对所得结果进行讨论。
解:从以上两幅图中均可看出,三个余弦函数的周期虽然不同,但它们抽样后相应抽样点所对应的值都相同。
那么这样还原回原先的函数就变成相同的,实际上是不一样的。
这是抽样频率太小的原因,我们应该增大抽样频率才能真实还原。
如下图:f=50Hz程序代码f=10;t=-0.2:0.001:0.2;g1=cos(6.*pi.*t);g2=cos(14.*pi.*t);g3=cos(26.*pi.*t);k=-0.2:1/f:0.2;h1=cos(6.*pi.*k);h2=cos(14.*pi.*k);h3=cos(26.*pi.*k);% subplot(3,1,1);% plot(k,h1,'r.',t,g1,'r');% xlabel('t');% ylabel('g1(t)');% subplot(3,1,2);% plot(k,h2,'g.',t,g2,'g');% xlabel('t');% ylabel('g2(t)');% subplot(3,1,3);% plot(k,h3,'b.',t,g3,'b');% xlabel('t');% ylabel('g3(t)');plot(t,g1,'r',t,g2,'g',t,g3,'b',k,h1,'r.',k,h2,'g.',k,h3,'b.')xlabel('t');ylabel('g(t)');legend('g1(t)','g2(t)','g3(t)');M2-1 利用DFT的性质,编写一MATLAB程序,计算下列序列的循环卷积。
《数字信号处理》朱冰莲二单元答案
2-1(4)0.5[()(10)]n u n u n --解:()101101091901(0.5)0.5[0.5[()(10)]]0.5||010.50.5n n nn z z u n u n z z z z z z ---=----===>--∑Z=0为九阶极点,Z=0.5为一阶极点;零点为()2100.50,1,2,9k j k z ek π==2-1(5)2-2(8)(8)21(1)(1)(1)(21)11111()()[()(1)]2[(1)(21)]1(1)()[]2()[]11111()[]11Nnnnnn n n N n N N N N N N X z x n zn u n u n N zNz n u n N u n N zd z z z d z z z N z dz z z z dz z z d z z z dz z z ∞∞∞----=-∞=-∞=+=-∞-+-+--+-+------==---+-------=--+---------=--+--+∑∑∑∑22(1)22212221112222222221()2()[]1112(1)2()(1)1(21)2(1)2(1)222222(21)(1)2(1)N N N NN N N N N NN N N N N N N N N z z d z z N z z dz z z z N z Nz N z z z z N z Nz z z N z Nz Nz Nz Nz N N N z z z z z z z z z ----------+-++++++-+----+---++-=+--+-+--+++--+-++=--+==-222122121(1)(1)(1)N N N N N z z z z z z ---+-=--通过验证知,上式分子中包含2)1(-z ,故在z=1处,零极点抵消。
故收敛域为 |z|>0,零点:无 极点:z=0或者:或者2-4()()()()112111111311125212121(1)2()(()2)();21(2)0.5()(2())(--1);21(3)()()()+2(--1);2112,12112112()n n n n n n z X z z z z z z x n u n z x n u n z x n u n u n X z z z z z z x n u -----------==+-+-->=-<=-=-==+<<----∴=-当时,为右序列,当时,为左序列,当0.5<<2时,为双边序列,1()22(--1)()2(--1)n n n u n u n u n +-=-- 2-5(1)+1(1)()12(1)(2)12|z|<2()()2(1)n X Z Z Z Z Z Z Z x n u n u n -==+----=----因1< 2-5(7)07()()()(1)6(4)(6)nn n z x n n n n n δδδδδδ-→-→=+-+-+-()(n )12-81111[*()]*(*)1,1(1)()(1)()*()(1)()1*(*),1(1)n n Z x n X z jX z j z x n j u n x n j u n jX z z j z-++-=+>-+∴=+∴=--∴=>--由(z)=2-12(1)直接法:111110.9911(),0.110(10.1)(10.1)10.1110()(0.1)()(10)(1)1(),0.110.1()(0.1)(1)1[()()],11n n n X z z z z z z x n u n u n Y z z zy n u n Z x n y n z z-----==-<<----∴=+--=<-∴=---=<-复卷积公式法: 110.990.99(),0.110(10.1)(10.1)(0.1)(10.1)1(),0.110.10.111[()()]()()210.99112(0.1)(10.1)0.110.992(0.1)(10.1)(0.1)c c c zX z z z z z z zY z z z z z Z x n y n X v Y dv j v v zv vdvz z j v v v vzdvj v v z v πππ--==<<----==<--==<------⎰⎰⎰收敛域120.1101max(0.1,10)10,0.1,100.990.99[()()]Re []Re [](0.1)(10.1)(0.1)(0.1)(10.1)(0.1)9.910.01(100.1)(1)11[()()],11v v zz v c v v zz zZ x n y n s s v v z v v v z v z z z z z z z z Z x n y n z z ==-<<∴==∴=+------=+=<----∴=<-为积分围线包围了两个一阶极点: 2-13(3)0000223()*() n 0()*()0n 0()*()n :()|z|>a;()()max(,0)||1()*()()*212()n n n n n n cn cx n y n x n y n x n y n a parsval az X Z z a Y Z z a V x n y n X V Yjavv j v a ππ∞=-∞∞=-∞∞=-∞-∞=-∞=<=>===-=<=-∑∑∑∑⎰⎰0()直接法当(V/1)*(1/V)dV=(1/v)dV 当n >=0时,只有一0000000002220222112v=a1=(()2()()n 1(()2()()(()n n cn n n cn n n n av d av v v a v jv a dz v a a av d av v v a v jv a dz v a d a v v dz v a ππ--------=---+-⎰⎰v=a v=av=a0个二次极点(1/v)dV (1/v))|=当n <0时,有两极点v=a,v=0(1/v)dV (1/v))|)|=02-14(3)22-22-214(1),()()11(2)()2(0)2(3)|()|2|()|54()(4)[()]()||2|()|280j n jw jw n jw jw n X e x n X e x X e dw x n dX e DTFT nx n jdwdX e dw nx n dwππππππππππππ∞=-∞-∞-=∞∞-=∞-=========∑⎰∑⎰∑⎰2-15=0==1()()()*()|2=()()11=x(0)y(0)=[()][()()]22n jw jw m m jw jw jw X e Y e dw x n y n x m y m X e dw X e Y e dw πππππππππ-∞-∞--=-⎰∑⎰⎰证明:因为两者都是因果的序列2-21(1)对于线性移不变系统,若输入序列为jwn e n x =)(,则系统输出为)13112()()()()(*)()()(--=====-∞-∞=-∞-∞=-∑∑jw jwnjw jwnm jwmjwnm m n jw e ee H eem h eem h n h n x n y(2)系统的频率响应为)](arg[11311)(jw e H j jwjw jwe ee e H =-+=-- 1111111z ()11/3(11/3)221()1||()11/311/311/331()=2()()()3n Y z z z H z z X z z z z h n u n n δ-------+--+====-+>---⋅-解:(1)对差分方程进行变换: (1-1/3z )Y(z)=(1+1/3z )X(z)所以(3)当系统输入为正弦序列是,则输出为同频的正弦序列,其幅度受频率响应加权,相位为输入信号相位与系统相位响应之和。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
n
x(n) a u (n)
n
18
逆Z变换
用长除法求 X ( z) [(1 az)(1 az 1 )]1 , | a 1 || z || a |的逆Z变换 例:
解: ∵收敛域| a 1 || z || a |为环域, ∴x(n)必为双边序列。
X ( z) 1 (1 az )(1 az )
x ( n)
2j
1
X ( z) z
n 1
dz
c
式中C为收敛域中的一条逆时针环绕原点的闭合曲线。
9
逆Z变换
若被积函数 X ( z ) z n 1是有理分式,一般采用留数定理来计 算围线积分 。根据留数定理, x(n) 等于围线C内全部极 点留数之和,即:
x(n) Re s[ X ( z ) z
收敛 不定 发散
对于不同的序列 x(n) ,可求得相应的收敛域。
6
Z变换的收敛域
收敛域内不包含任何极点,在极点处,X(z)为无穷大, Z变换不收敛。 有限长序列的收敛域为整个Z平面, 可能除开z=0, z=。 右边有限长序列: X(z)=x(1)z-1+x(2)z-2+···· |z|>0 左边有限长序列: X(z)=x(-1)z1+x(-2)z2+···· |z|< 如果是右边序列,并且|z|=位于收敛域内,那么, |z|>也位于收敛域内。
s k
[( z zi )
s
X ( z) z
]
z zi
k 1,, s
15
逆Z变换
例: 已知 X ( z)
z
2
( z 2)( z 0.5)
z 2
,求X(z)的原序列。
解: 将X(z)变为X(z)/z的形式并化为部分分式
X ( z) z z ( z 2)( z 0.5) A1 z2 A2 z 0.5
4
Z变换的定义
例3:求序列 x (n)= 解: ( z ) ( 1 ) X 3
n n
(1/3)|n| 的Z变换。
j Im[ z ]
Z平面
n
z
n
n
n
( )
n
1
1 n 3
z
n
(1) z 3
n n 0
0
( 1 z ) ( 31z ) 3
A0 A( z ) X (0 Re s[ i 1 ai z z aN
i 1
,0]
X ( z) 1 若序列为因果序列,且N≥M,当X(z)的N个极点都是单 Ak (1 z k z ) X ( z ) z z Re s[ , zk ] z 极点时,可以展开成以下的部分分式的形式:
n 1 n 0
1 3
Re[ z]
3
收敛域
|z|<3时,第一项收敛于
z z 3
,对应于左边序列。
1 3
|z|>1/3时,第二项收敛于
z z
,对应于右边序列。
z z
1 3
当
1 3
| z | 3
时:
X ( z)
z z 3
8z 3 ( z 3)( z 1 ) 3
ST
n
x( n ) z
双边Z变换 单边Z变换
n
拉氏变换与Z变换: X ( z ) z e X ( s )
2
Z变换的定义
例1:求序列 x (n)= an u(n) 的Z变换。
j Im[ z ]
解:
X ( z)
n
a u k
10
逆Z变换
如果 z k 为单阶极点,按留数定理:
Re s[ X ( z) z
k
n 1
, zk ] ( z zk ) X ( z ) z
k
n 1 z zk
如果 z k m阶极点,则其留数为: 为
z zk (m 1)! dz 在具体利用留数定理进行围线积分计算时,应根据被 积函数的特点及n值灵活选用公式来计算,可使问题 得以简化。 例如,在n小于某一值时,被积函数在围线内部z=0处 可能具有高阶极点,这时采用围线外部的极点进行计 算将方便得多。
j Im[ z ] j Im[ z ]
n
x ( n) z
0
n
x ( n) e
n n 0
j
j Im[ z ]
0
Re[ z ]
收敛域
0
Re[ z ]
收敛域
0
Re[ z ]
收敛域 右边序列的收敛域 左边序列的收敛域
双边序列的收敛域
8
逆Z变换
逆Z变换 从给定的Z变换表达式(包括收敛域)求原序列的过程称为 逆z变换。其实质是求X(z)的幂级数展开式各项的系数。 逆Z变换的三种基本方法 围线积分法 部分分式展开法 长除法(幂级数展开法) 围线积分法
X ( z ) x ( n) z
n 0 n
x ( n) e
n n 0
j
越大收敛越快。
所以,收敛域在圆外。
7
Z变换的收敛域
如果是左边序列,并且|z|=位于收敛域内,那么, 0<|z|< 的全部 z 值也位于收敛域内。
X ( z)
所以,收敛域在圆内。 如果是双边序列,收敛域由圆环组成。
x(n) Re s[ 1 n z
1
j Im[ z ]
围线C
z )( z a a a(a ( za a )( z ) )
|n|
z1 , a ]
n
1
a
n
1 a
2
n |z|=|a| 0
0
a
1/a
Re[ z ]
在收敛域内作包围原定的围线C a 收敛域 当 x(n) n n 0 时,只有一个单阶极点z=a, 1 a 2 1 a n 其围线积分为: 1 z , a] x(n) Re s[ n0 2
Re s[ X ( z ) z
n 1
, zk ]
1
d
m 1
[( z zk ) X ( z ) z m 1
m
n 1
]
11
逆Z变换
例1:已知某序列的Z变换为: X ( z) (1 az 1 ) 1 求原序列x(n)
z a
解:
x ( n)
1 2j 1 2j
(1 az
c
1
)
1
z
n 1
dz
za
c
1
z dz
n
由于收敛域为 z a ,可知该序列必定是因果序列。 并且当 n 0时,z=0处不是极点,被积函数仅有单阶 极点a,在收敛域内取围线C包含极点a,可求得:
12
逆Z变换
例2: 已知序列的Z变换为:
X ( z ) [(1 az )(1 az )] , | a || z || a |
2 3 3
1
1 1 a
2
[
a za
1 1 az
]
对右边序列 ∴右边序列为: 对应右边序列。用长除法
1 2
n az n z a z a aa 将其展开成z的负幂级数, 1 az 1 将其展开成z的正幂级数, x ( n) n0 x( ) a na n0 2 z将分母多项式按降幂排列 2 将分母多项式按升幂排列 1 - az 1 a 1 2 1 a aa z
x(0) x(1) z
1
x(2) z
2
若把X(z)展开成z-1的幂级数之和,则该级数的各系数 就是序列 x(n) 的值。 在具体进行长除法时,要根据收敛域先确定序列是左 边序列还是右边序列。对于左边序列Z变换为z的正幂 级数,分子分母多项式应按升幂排列展开;对于右边 序列,Z变换为z的负幂级数,分子分母应按降幂排列 进行展开。
第二章 离散时间信号与系统的 变换域分析
序列的 Z变换
序列的傅里叶变换
离散时间系统变换域分析
希尔伯特(Hilbert)变换
1
§1 序列的Z变换
Z变换的定义 抽样信号 x( n) x( nT )
ˆ x(t ) x(t ) T (t ) x(t ) (t nT )
n
a
n
Z平面
0
为保证收敛,则
a z
1
a
Re[ z ]
收敛域
z | z | | a |
或 | z || a |
X ( z) 1 1 ( )
a z
za
若 a = 1, 则
u ( n)
z z 1
| z | 1
3
Z变换的定义
例2:求序列 x(n)= -an u(-n-1)的Z变换。 解: ( z ) a nu (n 1) z n X
5
零点:0,极点:3,1/3
Z变换的收敛域
Z变换的收敛域 对于任意给定的序列 x(n) ,使其Z变换收敛的所有 z值的集合称为 X (z ) 的收敛域。 其收敛的充要条件是满足绝对可和条件,即:
n
x ( n) z
n
根据级数收敛的阿贝尔定理
lim
n
an
n
1 1 1
a( z a)( z a ) 1 a
13
逆Z变换
部分分式展开法 说明:1、 X(z)较简单时可按算术展开求各系数 B( z ) 用部分分式展开法求反Z变换,X ( z ) 通常为有理分式。 Ak(k=0,1…,N) 。 A( z ) 1、单极点 X(z)较复杂时可按留数定理求各系数 M 2、 i bi z Ak(k=0,1…,N),此时为了方便通常利用 X(z)/z的 B( z ) X ( z) b i 0N 形式求取: X ( z) N