绪论:本章介绍数字信号处理课程的基本概念。
0.1信号、系统与信号处理
1.信号及其分类
信号是信息的载体,以某种函数的形式传递信息。这个函数可以是时间域、频率域或其它域,但最基础的域是时域。
分类:
周期信号/非周期信号
确定信号/随机信号
能量信号/功率信号
连续时间信号/离散时间信号/数字信号
按自变量与函数值的取值形式不同分类:
2.系统
系统定义为处理(或变换)信号的物理设备,或者说,凡是能将信号加以变换以达到人们要求的各种设备都称为系统。
3.信号处理
信号处理即是用系统对信号进行某种加工。包括:滤波、分析、变换、综合、压缩、估计、识别等等。所谓“数字信号处理”,就是用数值计算的方法,完成对信号的处理。
0.2 数字信号处理系统的基本组成
数字信号处理就是用数值计算的方法对信号进行变换和处理。不仅应用于数字化信号的处理,而且
也可应用于模拟信号的处理。以下讨论模拟信号数字化处理系统框图。
(1)前置滤波器
将输入信号x a(t)中高于某一频率(称折叠频率,等于抽样频率的一半)的分量加以滤除。
(2)A/D变换器
在A/D变换器中每隔T秒(抽样周期)取出一次x a(t)的幅度,抽样后的信号称为离散信号。在A/D 变换器中的保持电路中进一步变换为若干位码。
(3)数字信号处理器(DSP)
(4)D/A变换器
按照预定要求,在处理器中将信号序列x(n)进行加工处理得到输出信号y(n)。由一个二进制码流产生一个阶梯波形,是形成模拟信号的第一步。
(5)模拟滤波器
把阶梯波形平滑成预期的模拟信号;以滤除掉不需要的高频分量,生成所需的模拟信号y a(t)。
0.3 数字信号处理的特点
(1)灵活性。(2)高精度和高稳定性。(3)便于大规模集成。(4)对数字信号可以存储、运算、系统可以获得高性能指标。
0.4 数字信号处理基本学科分支
数字信号处理(DSP)一般有两层含义,一层是广义的理解,为数字信号处理技术——DigitalSignalProcessing,另一层是狭义的理解,为数字信号处理器——DigitalSignalProcessor。
0.5 课程内容
该课程在本科阶段主要介绍以傅里叶变换为基础的“经典”处理方法,包括:(1)离散傅里叶变换及其快速算法。(2)滤波理论(线性时不变离散时间系统,用于分离相加性组合的信号,要求信号频谱占据不同的频段)。
在研究生阶段相应课程为“现代信号处理”(AdvancedSignalProcessing)。信号对象主要是随机信号,主要内容是自适应滤波(用于分离相加性组合的信号,但频谱占据同一频段)和现代谱估计。
简答题:
1.按自变量与函数值的取值形式是否连续信号可以分成哪四种类型?
2.相对模拟信号处理,数字信号处理主要有哪些优点?
3.数字信号处理系统的基本组成有哪些?
第一章:本章概念较多,需要理解和识记的内容较多,学习时要注意。
1.1 离散时间信号 1.离散时间信号的定义
离散时间信号是指一个实数或复数的数字序列,它是整数自变量n 的函数,表示为x(n)。一般由模拟信号等间隔采样得到:()()
a
a t nT
x n x x nT n ===-∞<<∞。
时域离散信号有三种表示方法:1)用集合符号表示 2)用公式表示 3)用图形表示
2.几种基本离散时间信号(记住定义)
(1)单位采样序列
(2)单位阶跃序列
(3)矩形序列
(4)实指数序列
(5)正弦序列
ω是正弦序列数字域的频率,单位是弧度。
对连续信号中的正弦信号进行采样,可得正弦序列。设连续信号为
,它的采样值为
,因此
(重点)
这个式子具有一般性,它反映了由连续信号采样得到的离散序列,其数字频率与模拟频率的一般关系。另外需要说明的是,ω的单位为弧度,Ω的单位为弧度/秒。本书中,我们一律以ω表示数字域频率,而以Ω及f 表示模拟域频率。 例:已知采样频率F T = 1000Hz, 则序列x (n ) = cos(0.4πn) 对应的模拟频率为 ( 400π ) 弧度/s 。 说明:本题旨在理解数字频率与模拟频率之间的关系:T
F Ω
=
ω。 (6)复指数序列
复指数序列是以余弦序列为实部、正弦序列为虚部所构成的一个复数序列。 (7)周期序列(重点)
所有n 存在一个最小的正整数N ,满足:)()(N n x n x +=,则称序列)(n x 是周期序列,周期为
N 。(注意:按此定义,模拟信号是周期信号,采用后的离散信号未必是周期的)
例:正弦序列
)sin(0n ω的周期性:
当k N πω20=,k 为整数时,)sin()](sin[00n N n ωω=+,即为周期性序列。周期02ωπk
N =
,式
中,k 、N 限取整数,且k 的取值要保证N 是最小的正整数。
可分几种情况讨论如下:(1)当0/2ωπ为整数时,只要1=k ,0/2ωπ=N 就为最小正整数,即周期为0/2ωπ。(2)当0/2ωπ不是整数,而是一个有理数时,设Q P //20=ωπ,式中,P 、Q 是互为素数的整数(互为素数就是两个数没有公约数),取Q k =,则P N =,即周期为P 。(3)当0/2ωπ是无理数时,则任何k 皆不能使N 为正整数,这时,正弦序列不是周期性的。 例:X(n) = cos(0.4πn)的基本周期为( 5 )。 [说明]基本周期的定义即计算公式:k N ω
π
2=
,其中N 和k 均为整数,N 为基本周期(使得N 为
最小整数时k 取值)。本题ω = 0.4π,代入上式得到:1,
5==k N 。
3.信号运算
(1)加法:两个信号之和 由同序号的序列值逐点对应相加得到。 (2)乘法:两个信号之积 由同序号的序列值逐点对应相乘得到。
(3)移位:当
,序列右移(称为延时);当
,序列左移(称为超前)。
(4)翻转:
(5)尺度变换:或,其中M 和N 都是正整数。
当
时,序列
是通过取x(n)的每第M 个采样形成,这种运算称为下采样。对
于序列,定义如下这种运算称为上采样。
4.信号分解(重点)
任一信号x(n)可表示成单位脉冲序列的移位加权和:
简记为
1.2 时域离散系统
时域离散系统定义 []()
().x n y n T ???
→???→ []()()y n T x n =
1 线性系统(重点)
判定公式:
若1()y n =1[()]T x n ,2()y n =2[()]T x n 则1212()[()()]()()y n T ax n bx n ay n by n =+=+
2 时不变系统(重点)
判定公式:y(n)=T[x(n)] y(n-0n )=T[x(n-0n )]
例:判断下列系统是否为线性、时不变系统。(重点) (1)()()2(1)3(2)y n x n x n x n =+-+-; (2)2
()()y n x n =; 解:
(1)令:输入为0()x n n -,输出为'000'
0000()()2(1)3(2)
()()2(1)3(2)()
y n x n n x n n x n n y n n x n n x n n x n n y n =-+--+---=-+--+--=
故该系统是时不变系统。
12121212()[()()]
()()2((1)(1))3((2)(2))
y n T ax n bx n ax n bx n ax n bx n ax n bx n =+=++-+-+-+-
1111[()]()2(1)3(2)T ax n ax n ax n ax n =+-+- 2222[()]()2(1)3(2)T bx n bx n bx n bx n =+-+- 1212[()()][()][()]T ax n bx n aT x n bT x n +=+
故该系统是线性系统。
(2)2
()()y n x n = 令:输入为0()x n n -,输出为'20()()y n x n n =-,因为
2'00()()()y n n x n n y n -=-=
故系统是时不变系统。又因为
212121222
12[()()](()()) [()][()] ()()
T ax n bx n ax n bx n aT x n bT x n ax n bx n +=+≠+=+
因此系统是非线性系统。
3 线性时不变系统(LTI 或者LSI 系统)输入与输出之间关系(重点):
()[()]h n T n δ=
()()()m y n x m n m δ∞
=-∞
=
-∑
()[()()]m y n T x m n m δ∞
=-∞=-∑
y (n )=
()()m x m h n m ∞
=-∞
-∑=x (n )*h (n )
重点:线性离不变系统的输出等于输入序列和该系统的单位脉冲响应的卷积
【说明】离散时间LTI 系统的单位冲激响应h(n)为系统对单位冲激序列δ(n)的零状态响应。
单位冲激响应的概念非常重要。在时域,LTI 系统可以由其单位冲激响应h(n)唯一确定,因此,我们常常用单位冲激响应描述 LTI 系统。在这种情况下, LTI 系统的输入输出关系可以由卷积运算描述:y (n )=
()()m x m h n m ∞
=-∞
-∑=x (n )*h (n )
物理意义: 卷积和运算具有显式意义,即可以用来确定系统的输出。如果系统确定,则其单位冲激响应是唯一的。由此,可求系统对任意输入的响应。
注意:计算卷积和的关键是求和区间的确定。因此,常常需要绘制序列x(m) 和h(n-m)的图形。利用序列x(m) 和h(n-m)的图形可助我们方便地确定求和区间。 卷积的求解方法(重点):
线性卷积是一种非常重要的一种运算,对它的求解,一般我们采用作图法。线性卷积满足交换律,设两序列长度分别是N 和M ,线性卷积后序列的长度为N +M -1。 卷积的计算过程包括翻转、移位、相乘、相加四个过程。 1)将
和
用
和
表示,画出
和
这两个序列; 2)选择一个序列,并将其按时间翻转形成序列
;
3)将移位n ,得到;
4)将
和
相同m 的序列值对应相乘后,再相加。
例:设(),x n n =04n ≤≤,4()()h n R n =, ()x n 和()h n 如图1所示。求()x n 和()h n 的卷积()y n 。(重点)
n
0 1 2 3
R 4(n )
1
0 1 2 3 4
4
n
()x n
图1
解 方法一:用图解法求卷积和。
(1) 将()x n 和()h n 用()x m 和()h m 表示(图2中(a)、(b)图)。
)
(a
…
) )
(
c
)
m
)
(d
)
(g
)
(f
)
(b
(e)
图2 图解法求卷积过程
(2) 将()
h m进行反折,形成()
h m
-(图2中(c)图);将()
h m
-移位n,得到()
h n m
-(图2中(d)、(e)、
(f)图)。
(3) 将()
x m和()
h n m
-相同m的序列值相乘,再相加,得到()
y n(图2中(g)图)。
{}
()1,3,6,10,9,7,4
y n=17
n
≤≤
再讨论解析法求线性卷积。
用式
()()()
m
y n x m h n m
+∞
=-∞
=-
∑
求解上式首先要根据()
x m和()
h n m
-的非零值区间确定求和的上下限,()
x m的非零值区间为14
m
≤≤,()
h n m
-的非零值区间为03
n m
-
≤≤,或3
n m n
-≤≤,由两个非零值区间可得n的取值区间为17
n
≤≤,它们的乘积()()
x m h n m
?-的非零值区间应满足:
14
m
≤≤和3
n m n
-≤≤
因此
当1
n<、7
n>时,()0
y n=;
当13
n
≤≤时,0
(1)
()1
2
n
m
n n
y n m
=
+
=?=
∑
;
当47
n
≤≤时,
4
3
(1)(8)
()1
2
m n
n n
y n m
=-
+-
=?=
∑
。
与图解法结果一致。
y(n)用公式表示为
(1)/2
()(1)(8)/2
n n
y n n n
+
?
?
=+-
?
?
?
13
47
n
n
≤≤
≤≤
其他
方法二:当序列()
x n和()
h n的长度分别为有限长N和M时,可采用“不进位乘法”求两序列线卷积。
如图1所示:
{}
()0,1,2,3,4
x n
↑
=
,
{}
()1,1,1,1
h n
↑
=
{}
()0,1,3,6,10,9,7,4
y n
↑
=
例:两线性时不变系统级联,其单位取样响应分别为
)
(
1
n
h和)
(
2
n
h,输入为)
(n
x,求系统的输出)
(n
y。
已知:
)
(
)
(n
u
n
x=,)4
(
)
(
)
(
1
-
-
=n
n
n
h δ
δ,)
(
)
(
2
n
u
a
n
h n
=。
解:设第一个系统的输出为
)
(n
ω,则
)3
(
)2
(
)1
(
)
(
)4
(
)
(
)]
4
(
)
(
[
)
(
)
(
)
(
)
(
1
-
-
-
=
-
-
=
-
-
*
=
*
=
n
n
n
n
n
u
n
u
n
n
n
u
n
h
n
x
n
δ
δ
δ
δ
δ
δ
ω
+
+
+
因而输出为
)3
(
)2
(
)1
(
)
(
)
(
)]
3
(
)2
(
)1
(
)
(
[
)
(
)
(
)
(
3
2
1
2
-
+
-
+
-
+
=
*
-
+
-
+
-
+
=
*
=
-
-
-n
u
a
n
u
a
n
u
a
n
u
a
n
u
a
n
n
n
n
n
h
n
n
y
n
n
n
n
n
δ
δ
δ
δ
ω
4. 系统因果性和稳定性的判定(重点)
1)稳定系统:有界的输入产生的输出也有界的系统,即:若|()|
x n<∞,则|()|
y n<∞(记住!!)
线性移不变系统是稳定系统的充要条件:|()|
n
h n
∞
=-∞
<∞
∑(系统稳定的充分必要条件是系统的单位脉
冲响应绝对可和)(记住!!)
或:其系统函数H(z)的收敛域包含单位圆|z|=1(记住!!)
2)因果系统:0n 时刻的输出0()y n 只由0n 时刻之前的输入0(),x n n n ≤决定(记住!!)
线性移不变系统是因果系统的充要条件:()0,0h n n =<(记住!!)因果系统的单位脉冲响应必然是因果序列。(记住!!)
或:其系统函数H(z)的收敛域在某圆外部:即:|z|>Rx (记住!!) 3)稳定因果系统:同时满足上述两个条件的系统。 线性移不变系统是因果稳定系统的充要条件:
|()|n h n ∞
=-∞
<∞∑,()0,0h n n =<(记住!!)
或:H(z)的极点在单位圆内H(z)的收敛域满足:||,1x x z R R --><(记住!!) 例:判断线性时不变系统的因果性、稳定性,并给出依据。(重点)
(1)1
1
()()N k y n x n k N
-==
-∑;
(2)00
()()n n k n n y n x k +=-=
∑
;
解:(1)只要1N ≥,该系统就是因果系统,因为输出只与n 时刻的和n 时刻以前的输入有关。如果()x n M ≤,则()y n M ≤,因此系统是稳定系统。 (2)如果()x n M ≤,0
0()()21n n k n n y n x k n M +=-≤
≤+∑
,因此系统是稳定的。系统是非因果的,
因为输出还和x(n)的将来值有关。
注意:如果给出的是h(n),用上面要求记住的充要条件判断!
例:设某线性时不变系统的单位取样响应为)()(n u a n h n
=(a 为实数),分析系统的因果性和稳定
性。(重点) 解:讨论因果性: