数字信号处理_第三章
x( L 1) x( L 2) y (0)c x(0) y (1) x(1) x(0) x( L 1) c y (2)c = x(2) x(1) x(0) x( L 1) x( L 2) x( L 3) y ( L 1)c
DFTx2 (n) X 2 (k )
二、循环移位性质
1、序列的循环移位(圆周移位)定义: 一个有限长序列 x(n) 的圆周移位定义为
y(n) xn mN RN n
(1) 先将x(n)作 周 期 ~ x延 (n) 拓 xnN
~ n mN (2) 延 拓 后 再 进 x (n 行 m移 ) x位
1 e
e
k j 38
sin(k / 2) sin(k / 8)
15 j
0k 7
2kn 16
(2)N 16 时 X (k ) x(n) W
n 0 N 1 nk N
R4 (n)e
n 0
e
n 0
3
kn j2 16
1 e
4k j2 16 k j2 16
~ 周期序列 x (n) 是有限长序列x(n)的周期延拓。
x (n) 0 n N 1 ~ x(n) 其它 0
或
x(n) ~ x (n) RN (n)
x (n) 的主值序列。 有限长序列x(n)是周期序列 ~
二、DFT的隐含周期性
如:
0
x(n)
n N-1
~ x ( n)
yc (n) h(n) L x(n) [ h(m) x((n m))L ]RL (n)
m 0 L 1
其中L为循环卷积区间长度,L≥max(N,M)
表示方法:
L 或
*
循环卷积的计算方法
2、循环卷积的计算方法——矩阵相乘
x(n) 序列: {x(0), x(1), x(2) …, x(L-1)} x(n)的循环倒相序列: 令n = 0, m = 0, 1, …, L-1,x((n-m))L形 成的序列为
...
0 N-1
...
n
定义从n=0 到(N-1)的第一个周期为主值序列或区间。
二、DFT的隐含周期性
~ 3、频域周期序列 X (k ) 与有限长序列X(k)的关系
~ X ( k ) X k N ~ X ( k ) X ( k ) RN ( k )
~ 周期序列 X (k ) 是有限长序列X(k)的周期延拓。 ~ 有限长序列X(k)是周期序列 X (k ) 的主值序列。
时域圆周移位的性质
3、时域圆周移位的性质
设有限长序列为 x(n), 0 n N 1, 若xm (n) xn mN RN n
km 则有X m (k ) DFT[ xm (n)] WN X (k )
即:有限长序列的圆周 移位只引入一个相移 W
mk N
e
x(n) x(0), x(1), x(2), x(3) 1,1,1,1
解: h(n)与x(n)的4点循环卷积矩阵形式为
yc (0) 1 y (1) 2 c yc (2) 3 yc (3) 4
4 3 2 1 10 1 10 1 4 3 2 1 4 1 10 3 2 1 1 10
(3) 最 后 取 主 y( 值 n) 序 xn 列 mN RN n
循环移位
x ( n)
0
N-1
~ x (n) x((n))N
n
周期延拓
0
N-1
n
循环移位
循环移位
2、循环移位的含义 (1)主值区间: n=0~N-1; (2)当某序列值从此区间一端移出时,与它相同 的序列值又从此区间的另一端移进来; (3)如果把x(n)首尾排列(n=0~N-1)在一个N等分 的圆周上,序列的移位就相当于x(n) 在圆上旋转, 故又称作圆周移位。当围着圆周观察几圈时,看到 的就是周期序列 ~ x ( n)
nN
n1
或n模N等于n1
例 : (1) n 25, N 9,
(2) n 4, N 9,
259 7
49
5
二、DFT的隐含周期性
~ 2、有限长序列x(n)和周期序列 x (n) 的关系
~ x ( n)
m
x(n m N) xnN
j
2 mk N
, 对信号的幅度没有影响 。
频域圆周移位的性质
4、频域圆周移位的性质(调制特性)
IDFT[ X ((k l ))N RN (k )] W x(n)
nl N
2nl 1 DFT x(n) cos [ X ((k l )) N X ((k l )) N ]RN (k ) N 2
DT FT: X (e jw ) x(n)e jwn
n 0 N 1
DFT: X (k ) x(n)e
n 0
N 1
j
2 kn N
kn x(n)WN n 0
N 1
X (k ) X (e )
j
2 k N
DFT与序列傅里叶变换的关系
X(k)也可以看作序列x(n)的傅里叶变换X(ejω)在 区间[0, 2π]上的N点等间隔采样,其采样间 隔为ωN=2π/N
连续: 不适合 计算机 处理
正 : X (e jw )
n
jwn x ( n ) e
反 : x(n)
1 2
X (e
jw
)e
jwn
dw
由DTFT到DFT
离散时间、离散频率的傅立叶变换(DFT) 由上述分析可知,对DTFT,要想在频域上 离散化,那么在时域上必须作周期延拓。 对长度为M的有限长序列x(n),以N为周期 延拓(N≥M)。
x((0))L , x((1))L , x((2))L , , x((L 1)) L x(0), x( L 1), x( L 2), , x(1)
循环卷积的计算方法
令n = 1, m = 0, 1, …, L-1,x((n-m))L形成的序
列为
x((1))L , x((0)) L , x((1)) L , , x((L 2)) L x(1), x(0), x( L 1), , x(2)
(4)可约性: WiN WN
in
n
例1、已知 x(n) R4 (n),分别求8和16点DFT ( 1 ) N 8时 解:
X (k ) x(n)W
n 0 N 1 nk N
R4 (n)e
n 0
7
j
2kn 8
e
n 0
3
kn j 28
1 e
4k j 28 k j 28
该序列相当于x(n)的循环倒相序列向右循环移1 位。 再令n = 2, m = 0, 1, …, L-1,此时 得到的序列又是上面的序列向右循环移1位。依 次类推,当n和m均从0变化到L-1时,得到一个
关于x((n-m))L的矩阵如下:
循环卷积的计算方法
循环卷积矩阵:
பைடு நூலகம்
循环卷积的计算方法
yc (n) x(n) L h(n)
二、DFT的隐含周期性
这里的周期延拓仅看作数学处理方法,或 者说借助时域周期延拓实现有限长序列频 谱的离散化。 在DFT中,有限长序列都是作为周期序列 的一个周期来表示的,总是隐含周期性。
三、DFT与序列傅里叶变换、Z变换的关系
若x(n)是一个有限长序列,长度为N
N 1 n 0
Z变换:
循环卷积的计算方法
h(n)与x(n)的8点循环卷积矩阵形式为
yc (0) 1 y (1) c 2 yc (2) 3 y (3) c 4 y (4) 0 c yc (5) 0 y (6) 0 c 0 yc (7) 0 1 2 3 4 0 0 0 0 0 1 2 3 4 0 0 0 0 0 1 2 3 4 0 0 0 0 0 1 2 3 4 4 0 0 0 0 1 2 3 3 4 0 0 0 0 1 2 2 1 1 1 3 3 4 1 6 0 1 10 0 0 9 0 0 7 0 0 4 1 0 0
注意:离散傅里叶变换(DFT)只对有限长 序列作周期延拓或周期序列成立。
……
……
3.1 离散傅里叶变换的定义及物理意义
一、 DFT的定义
设x(n)是一个长度为M的有限长序列,则定义x(n) 的N点离散傅里叶变换为
N 1 n 0 2 kn N N 1 n 0
正变换: X (k ) DFT[ x(n)]N x(n)e
引言
DFT要解决两个问题:
一是频谱的离散化; 二是算法的快速计算(FFT)。 这两个问题都是为了使计算机能够实时处 理信号。
引言
离散
连续
周期性
非周期性
时域周期化→频域离散化 时域离散化→频域周期化
序列的傅里叶变换
离散时间、连续频率的傅立叶变换(序列 的傅立叶变换)
x(n) -1 0 1 2 t 时域离散、非周期 频域连续、周期
j
kn x(n)WN
反变换: x(n) IDFT[ X (k )]N
《数字信号处理》 完整加精版
采用抽象算法表达:由软件程序虚拟实现。 在采用硬件电路实现时,由于不需要考虑 物理环境对信号的影响,可以在设计中尽可
能采用低功耗高密度集成。
数字系统的特点
信号采用数字序列表达后,对模拟信号难以 进行的很多处理能够方便地实现,例如: 对信号的乘法调制和各种编码调制、信号的时 间顺序处理、信号的时间压缩/扩张、复杂标准 信号的产生…
时间变量与对应的函数值采用两个相等长度的序列 (一维向量)表示。 两个序列可以进行直接数值设臵:
例:n=[0 1 2 3 4 5 6 7];
x=[1 2 4 6 5 3 1 0];
数字信号的MATLAB表达
坐标区间设臵: n=[n1:n2] 只取整数,设定起点和终点;
信号函数设臵:其序列长度由n序列限定; x=3*n x=exp(j*(pi/8)*n)
设臵好坐标序列t和信号序列x后,可以采 用下列作图语句画出连续时间信号图形: plot(t,x) 该语句通过将离散的信号点之间用直线连 接得到连续图形。
模拟信号的作图表达
例:MATLAB程序
t=[0:0.1:10];x1=[zeros(1,30) ones(1,40) zeros(1,31)]; x2=2-0.3*t;x3=exp(j*(pi/8)*t);x4=exp(-0.2*t).*cos(2*pi*t);
欠采样导致的问题
s N
若原始频谱与镜像频谱混叠,产生混叠失真,则
信号不可恢复!
采样定理
待采样信号必须为带限信号
X 0
M
采样频率应大于信号最高频率的2倍
2 s 2M N Ts
Nyquist 频率
重建滤波器(低通)截止频率应满足:
数字信号处理第三版第3章.ppt
x2 (n) N•DFT X 2 (k )
y(n) x1(n) x2 (n) Y (k) DFT[ y(n)]
1 N 1
N l0
X1(l) X 2 ((n L))N RN (n)
1 N1 N l0
X 2 (l) X1((n L))N RN (n)
第3章 离散傅里叶变换(DFT)
• DFT与Z变换和DTFT关系图解说明
z e WNk
j 2 k
e N
j
2 k
N
2 k
N
第3章 离散傅里叶变换(DFT)
• DFT与Z变换和DTFT关系举例说明
第3章 离散傅里叶变换(DFT)
• DFT的隐含周期性
N 1
1768年3月21日傅里叶生于法国荣纳省欧塞尔。其父亲 是裁缝,且很早就父母双亡,小时候在天主教受的教育。 毕业后在军队中教授数学。
1795年他到巴黎高等师范教书。 1798年随拿破仑东征,任下埃及的总督。 1801年,远征军失败后回到法国,任伊泽尔省长官。 1822年当选为科学院秘书,发表《热的分析理论》一文。在文中首次提出 并证明了将周期函数展开为正弦级数的原理,从而奠定了傅里叶级数(FS)与傅 里叶变换(FT)的理论基础。二者后被统称为傅里叶分析(FA)。 为了使FA应用于工程实际,人们又提出了离散傅里叶变换(DFT),但因计 算量太大而在较长时间内并未得到广泛应用,直到1965年美国Coo1y和Tukey两 人提出快速傅里叶变换(FFT)之后,FA才真正从理论走向实践,成为大家爱不 释手的一种数学工具。 1830年5月16日病逝于巴黎。
,求它的N点DFT。
第3章 离散傅里叶变换(DFT)
数字信号处理答案(第三版)清华大学
数字信号处理教程课后习题答案目录第一章离散时间信号与系统第二章Z变换第三章离散傅立叶变换第四章快速傅立叶变换第五章数字滤波器的基本结构第六章无限长单位冲激响应(IIR)数字滤波器的设计方法第七章有限长单位冲激响应(FIR)数字滤波器的设计方法第八章数字信号处理中有限字长效应第一章 离散时间信号与系统1 .直接计算下面两个序列的卷积和)n (h *)n (x )n (y =请用公式表示。
分析:①注意卷积和公式中求和式中是哑变量m ( n 看作参量), 结果)(n y 中变量是 n ,; )()()()()(∑∑∞-∞=∞-∞=-=-=m m m n x m h m n h m x n y ②分为四步 (1)翻褶( -m ),(2)移位( n ),(3)相乘,; )( )( 4n y n n y n 值的,如此可求得所有值的)相加,求得一个(③ 围的不同的不同时间段上求和范一定要注意某些题中在 n00 , 01()0 , ,()0,n n n a n N h n n n n x n n n β-⎧≤≤-=⎨⎩⎧≤⎪=⎨<⎪⎩其他如此题所示,因而要分段求解。
)(5.0)(,)1(2 )()4()(5.0)(,)2( )()3()()(,)( )()2()()(,)( )()1(3435n u n h n u n x n R n h n n x n R n h n R n x n R n h n n x n n n =--==-=====δδ2 .已知线性移不变系统的输入为)n (x ,系统的单位抽样响应 为)n (h ,试求系统的输出)n (y ,并画图。
分析:①如果是因果序列)(n y 可表示成)(n y ={)0(y ,)1(y ,)2(y ……},例如小题(2)为)(n y ={1,2,3,3,2,1} ;②)()(*)( , )()(*)(m n x n x m n n x n x n -=-=δδ ;③卷积和求解时,n 的分段处理。
数字信号处理知识点整理Chapter3.
第三章 自适应数字滤波器3.1 引言滤波器的设计都是符合准则的最佳滤波器。
维纳滤波器参数固定,适用于平稳随机信号的最佳滤波;自适应滤波器参数可以自动地按照某种准则调整到最佳。
本章主要涉及自适应横向滤波器.....、自适应格型滤波器........、最小二乘自适应滤波器..........。
3.2 自适应横向滤波器自适应...线性组合....器.和自适应....FIR ...滤波器...是自适应信号......处理的基础.....。
3.2.1 自适应线性组合器和自适应FIR 滤波器自适应滤波器的矩阵表示式 滤波器输出:()()()1N m y n w m x n m -==-∑n 用j 表示,自适应滤波器的矩阵形式为T T j jj y ==X W W X 式中1212,,,,,,,TTN N w w w x x x ⎡⎤⎡⎤==⎣⎦⎣⎦W X误差信号表示为T T j j j j jj j e d y d d =-=-=-X W W X 与维纳滤波相同,先考虑最小均方误差准则:()2222T T j j j j dx xx E e E d y E e ⎡⎤⎡⎤⎡⎤=-=-+⎣⎦⎣⎦⎢⎥⎣⎦R W W R W2j E e ⎡⎤⎣⎦称为性能函数....,将其对每个权系数求微分,形成一个与权系数相同的列向量: 2221222,,,Tj j jj xx dx N E e E e E e w w w ⎡⎤⎡⎤⎡⎤⎡⎤∂∂∂⎣⎦⎣⎦⎣⎦⎢⎥∇==-∂∂∂⎢⎥⎣⎦R W R令梯度为零,可得最佳权系数此时最小均方误差为:22*min T j j dx E e E d ⎡⎤⎡⎤=-⎣⎦⎣⎦W R 要求2minj Ee ⎡⎤⎣⎦和最佳权系数*W ,先求自相关矩阵xx R 和互相关矩阵dx R 。
3.2.2 性能函数表示式及几何意义3.2.3 最陡下降法3.2.1给出了要求2minj Ee ⎡⎤⎣⎦和最佳权系数*W 的理论求解方法,但实际很难应用。
数字信号处理教程课后习题及答案
6.试判断:
是否是线性系统?并判断(2),(3)是否是移不变系统?
分析:利用定义来证明线性:满足可加性和比例性, T [a1 x1 (n ) + a 2 x2 (n )] = a1T [ x1 (n )] + a2T [ x2 (n )] 移不变性:输入与输出的移位应相同 T[x(n-m)]=y(n-m)。
,
(2)x(n) = R3(n)
,
(3)x(n) = δ (n − 2) ,
(4)x(n) = 2n u(−n − 1) ,
h(n) = R5(n) h(n) = R4 (n) h(n) = 0.5n R3(n) h(n) = 0.5n u(n)
分析:
①如果是因果序列 y (n ) 可表示成 y (n ) ={ y (0) , y(1) , y(2) ……},例如小题(2)为
y1 (1) = ay1 (0) + x1 (1) = 0 y1 (2) = ay1 (1) + x1 (2) = 0
┇
8
y1(n) = ay1(n − 1) + x1(n) = 0 ∴ y1 (n) = 0 , n ≥ 0 ii) 向 n < 0 处递推,将原方程加以变换
y1(n + 1) = ay1(n) + x1(n + 1)
结果 y (n ) 中变量是 n ,
∞
∞
∑ ∑ y (n ) =
x ( m )h (n − m ) =
h(m)x(n − m) ;
m = −∞
m = −∞
②分为四步 (1)翻褶( -m ),(2)移位( n ),(3)相乘,
数字信号处理高西全课后答案ppt
详细描述
线性时不变系统是指系统的输入和输出之间存在线性关系,并且系统的特性不随时间变化而变化。这种系统的行为可以用线性常系数微分方程来描述,同时它的输出不依赖于输入的时间函数,只依赖于输入的初始状态。
线性时不变系统
VS
频域分析可以揭示信号的频率成分和频率域中的每个成分与原始信号之间的关系。通过在频域中对信号进行分析和处理,可以实现信号的滤波、去噪、压缩和恢复等功能。
频域分析在信号处理、图像处理、通信系统等领域得到广泛应用。例如,在图像处理中,频域分析可以用于图像滤波、边缘检测等任务;在通信系统中,频域分析可用于调制解调、频谱分析等。
详细描述
04
第四章 傅里叶变换与频域分析
傅里叶变换的定义
傅里叶变换是一种将时间域信号转换到频域的方法,通过将信号分解成一系列不同频率的正弦和余弦函数的线性组合。
傅里叶变换的性质
傅里叶变换具有一些重要性质,包括线性、对称性、可逆性、Parseval等式等。这变换的定义与性质
离散时间信号
定义
如果信号仅在离散时间点上有定义,则该信号称为离散时间信号。
例子
数字音频、图像数据等。
数学表示方法
通常使用序列形式来表示,例如y[n] = sin(n)。
01
03
02
连续时间信号的数学表示方法
离散时间信号的数学表示方法
其他表示方法
信号的数学表示方法
03
第三章 系统分析基础
总结词
快速傅里叶变换(FFT)算法的基本思想
根据算法实现方式的不同,可以分为按时间抽取(DIT)和按频率抽取(DFT)两种FFT算法。
经典:数字信号处理2-Z变换
| x(n) || z |n
n
n
收敛域在哪个平面:Z平面/复3 数平面 (z是个
收敛域形状有哪些:全平面,圆内4,圆外,圆环 对应于什么序列:有限长,左边5,右边,双边 11 a<|z|<bn)
正、逆Z变换:逆变换
定义和对比
正正::X (Xz)( z) x(nx)(znn),z n , Rx Rz x Rxz Rx n n
z (0, ] z [0, )
n1 >- , n2
z (0, )
n1 >- , n2 = n1 0, n2 =
z (Rx, ) z (Rx, ]
n1 =- , n2 +
z (0, Rx )
n1 =- , n2 0 n1 =- , n2 =
z [0, Rx ) z (Rx, Rx )
X(z)极点
0 +∞
0, +∞ Rx-, +∞ Rx0, Rx+
Rx+ Rx-, Rx+
5
正、逆Z变换:收敛域(有限长)
xn
xn
xn
n1 0 n2
0 z
Im(z)
n1
n2
0 z
Im(z)
n1
n2
0 z
Im(z)
Re(z)
Re(z)
Re(z)
6
正、逆Z变换:收敛域(右边)
Im(z)
X
(
z)
2 (1 z1
)
+
(1
1 0.5z
1 )
正、逆Z变换:逆变换
展开法例题(P57 例3.2.3)
1 X ( z) (1 z 1)(1 0.5z 1) , | z | 1
数字信号处理离散傅里叶变换DFT(第三章)
0kN 1
比较上面二式可得关系式
X (k )X (z )
2 j k N z e
,
0kN 1 0kN 1
( 3 .1 .3 ) ( 3 .1 .4 )
X (k )X ( ej)
2 , k N
离散傅里叶变换的定义
图 3.1.1
X(k)与X(e jω)的关系
离散傅里叶变换的定义
(kmN )n X(k mN) x(n)WN n0 kn x(n)W N X (k) n0 N1 N1
同理可证明(3.1.2)式中
x(n+mN)=x(n)
离散傅里叶变换的定义
实际上, 任何周期为N的周期序列
x(n )
都可以看作长度
为N的有限长序列x(n)的周期延拓序列, 而x(n)则是
循环移位性质 1. 序列的循环移位 设x(n)为有限长序列, 长度为N, 则x(n)的循环移位定义为 y(n)=x((n+m))NRN(N) (3.2.2)
X ( k ) x ( n )W
n0
7
kn 8
e
n0
3
j
2 kn 8
e
3 j k 8
设变换区间N=16, 则
16
sin( k ) 2 , k 0,1, , 7 sin( k ) 8
3 j 2 kn 16
kn X (k ) x(n)W16 e n 0 n 0
N称为DFT变换区间长度,一般N≥M, 通常称(3.1.1)式和(3.1.2) 式为离散傅里叶变换对。
离散傅里叶变换的定义
下面证明IDFT[X(k)]的唯一性。 把(3.1.1)式代入(3.1.2)式有
数字信号处理(第四版)第三章-下ppt
Nyquist frequency:
Also folding frequency, cut off frequency
Nyquist rate:
Over sampling
Under sampling Critical sampling
Example 3.18
14
Digital Signal Processing
9
Digital Signal Processing
© 2013 Jimin Liang
Discrete-Time Signals in Frequency Domain 3.8 Sampling theorem Form 1 of Gp(jOmega)
10
Digital Signal Processing
© 2013 Jimin Liang
Discrete-Time Signals in Frequency Domain 3.8 Sampling theorem
G1( j)
G1 p ( j)
10
(a)
6
0
G 2( j)
Hr ( j)
6
(d )
0.1
20 c 6 0 6 c
P( j)
Ga ( j)
1
1
…
m
0 (a)
m
…
T
0
T
2T
3T
(b)
G p ( j )
1T
T
...
0
T
2T
T mmΒιβλιοθήκη 13Fig.3.15
Digital Signal Processing
数字信号处理DFT
DFT[x(n)]
1 N
X1(k)
X 2 (k )
3.3 频率域采样
1、对于有限长(N点)序列x(n)
Z[x(n)]=X (Z)
取单位圆上的Z变换
X (e j ) FT[x(n)]
会引起x(n)周期延拓
频域一周内等间隔N点采样
~
~
X (k) DFS[x(n)]
X (k) DFT[x(n)]
00
tt
T TpTp / N
XXXaa(a((jjfjff)))
fs 2 fc
T (2 fc )1
N 1
X ( jf ) T xa (nT )e j2 fnT n0
令f kF
00
F fs2/fNsc
N 1
j 2 kn
X ( jkF) T xa (nT )e N
n0
N 1
j 2 kn
时域信号 频域信号
连续的 非周期的
周期的 离散的
三 离散时间、连续频率的序列傅里叶变换
x(nT) T
X e j 或 X (e jT )
---
-T 0 T 2T t
0
时域信号 频域信号
离散的 周期的 非周期的 连续的
---
s
2 T
四 离散时间、离散频率的离散傅里叶变换
x(nT)=x(n)
Tp
1 F
1
W k4 4
1 W4k
4 0
k 0 k 1, 2,3
x(n)的8点DFT为
X (k)
7
x(n)WNkn
3
j 2 kn
e8
n0
n0
x(n)的幅频特性特性曲线及其4点、8点、16点DFT.
