数字信号处理程佩青第三版课件第八章数字信号处理中的有限字长效

合集下载

数字信号处理

数字信号处理
2《Digital Signal Processing—A Computer Approach Third Edition》 Mitra
3《Digital Signal Processing》A.V.Oppenheim 4…….
4
第一章 数字信号处理概述
1.1 数字信号处理技术 1.2 数字信号与连续时间信号的关系 1.3 数字信号处理的分析方法 1.4 A/D、D/A原理 1.5 模拟信号的数字滤波
12
1.4 A/D、D/A原理
1.4.1 A/D原理与抽样定理
模拟信号的抽样 抽样信号的频谱 无失真抽样条件 前置预滤波器的作用 A/D变换的指标
.4.2 D/A原理和重构定理
重构定理 一种D/A变换器原理
13
1.4.1 A/D原理与抽样定理
A/D 将模拟信号转变为数字信号
s
Ya (
j)

FT

ya (t) X a ( j)G(
ya (t) xa (t)
j)

Xa(
j) (*)
X a ( j)
19
讨论
1、(*)式成立的条件:
s 2m
s
1
T
k
Xa(
j
jks )
Xˆ a ( j) s
当m s / 2
Xˆ a ( j)
18
m s / 2
时信号的提取
xˆa (t)
G( j)
Xˆ a ( j)
ya (t)
G(
j)

T , 0,


1 2
s


1 2
s

(完整word版)《数字信号处理》课程教学大纲

(完整word版)《数字信号处理》课程教学大纲

课程编号15102308《数字信号处理》教学大纲Digital Signal Processing一、课程基本信息二、本课程的性质、目的和任务《数字信号处理》课程是信息工程本科专业必修课,它是在学生学完了高等数学、概率论、线性代数、复变函数、信号与系统等课程后,进一步为学习专业知识打基础的课程。

本课程将通过讲课、练习使学生建立“数字信号处理”的基本概念,掌握数字信号处理基本分析方法和分析工具,为从事通信、信息或信号处理等方面的研究工作打下基础。

三、教学基本要求1、通过对本课程的教学,使学生系统地掌握数字信号处理的基本原理和基本分析方法,能建立基本的数字信号处理模型。

2、要求学生学会运用数字信号处理的两个主要工具:快速傅立叶变换(FFT)与数字滤波器,为后续数字技术方面课程的学习打下理论基础。

3、学生应具有初步的算法分析和运用MA TLAB编程的能力。

四、本课程与其他课程的联系与分工本课程的基础课程为《高等数学》、《概率论》、《线性代数》、《复变函数》、《信号与系统》等课程,同时又为《图像处理与模式识别》等课程的学习打下基础。

五、教学方法与手段教师讲授和学生自学相结合,讲练结合,采用多媒体教学手段为主,重点难点辅以板书。

六、考核方式与成绩评定办法本课程采用平时作业、期末考试综合评定的方法。

其中平时作业成绩占40%,期末考试成绩占60%。

七、使用教材及参考书目【使用教材】吴镇扬编,《数字信号处理》,高等教育出版社,2004年9月第一版。

【参考书目】1、姚天任,江太辉编,《数字信号处理》(第二版),华中科技大学出版社,2000年版。

2、程佩青著,《数字信号处理教程》(第二版),清华大学出版社出版,2001年版。

3、丁玉美,高西全编著,《数字信号处理》,西安电子科技大学出版社,2001年版。

4、胡广书编,《数字信号处理——理论、算法与实现》,清华大学出版社,2004年版。

5、Alan V. Oppenheim, Ronald W. Schafer,《Digital Signal Processing》,Prentice-Hall Inc, 1975.八、课程结构和学时分配九、教学内容绪论(1学时)【教学目标】1. 了解:什么是数字信号处理,与传统的模拟技术相比存在哪些特点。

《数字信号处理教程》程佩青(第三版)清华大学出版社课后答案

《数字信号处理教程》程佩青(第三版)清华大学出版社课后答案

结果 y (n ) 中变量是 n ,


∑ ∑ y (n ) =
x ( m )h (n − m ) =
h(m)x(n − m) ;
m = −∞
m = −∞
②分为四步 (1)翻褶( -m ),(2)移位( n ),(3)相乘,
(4)相加,求得一个 n 的 y(n) 值 ,如此可求得所有 n 值的 y(n) ;
10
T [ax1(n)+ bx2 (n)] =
n

[ax1
(n
)
+
bx2
(n
)]
m = −∞
T[ax1(n) + bx2(n)] = ay1(n) + by2(n)
∴ 系统是线性系统
解:(2) y(n) =
[x(n )] 2
y1(n)
= T [x1(n)] = [x1(n)] 2
y2 (n) = T [x2 (n)] = [x2 (n)] 2
(3) y(n) = δ (n − 2) * 0.5n R3(n) = 0.5n−2 R3(n − 2) (4) x(n) = 2n u(−n −1) h(n) = 0.5n u(n)
当n ≥ 0 当n ≤ −1
∑ y(n) = −1 0.5n−m 2m = 1 ⋅ 2−n
m = −∞
3
y(n) = ∑n 0.5n−m 2m = 4 ⋅ 2n
+ 1)

x1 (n
+ 1)]
=
−a n
综上 i) , ii) 可知: y1 (n) = −a nu(−n − 1)
(b) 设 x(n) = δ (n − 1)
i)向 n > 0 处递推 ,

高西全-丁玉美-数字信号处理课件(第三版)

高西全-丁玉美-数字信号处理课件(第三版)

出版信息
出版社:清华大学出版社 出版时间:2019年 作者:王志强、李志刚、张志强 内容简介:本书主要介绍数字信号处理的基本概念、原理和方法,以及其在通信、雷达、 图像处理等领域的应用。
主要内容
数字信号处理的 基本概念和原理
数字信号处理的 应用领域
数字信号处理的 算法和实现
数字信号处理的 发展趋势和挑战
感谢观看
汇报人:PPT
信号处理在音频处理中的应用
添加标题
添加标题
信号处理在图像处理中的应用
添加标题
添加标题
信号处理在雷达系统中的应用
04
学习资源
习题答案
教材配套习题 答案
教师提供的习 题答案
网络资源: 如CSDN、 GitHub等
同学之间的互 助解答
教学PPT
课件形式:图文并茂,动 画演示,互动问答等
课件内容:数字信号处理 基础知识、应用案例、实 验操作等
课件特点:简洁明了,逻 辑清晰,易于理解
课件下载:提供课件下载 链接,方便学生课后复习
和预习
学习笔记
教材:数字信号处理课件(第三 版)
学习资料:教材、课件、实验指 导书、习题集等
添加标题
添加标题
添加标题
添加标题
课程内容:数字信号处理基础知 识、数字信号处理算法、数字信 号处理应用等
学习工具:MATLAB、Python 等编程工具,数字信号处理软件 等
实验指导书
实验目的:掌握数字信号处理的基本概念和原理 实验内容:包括信号的采样、量化、编码、传输、解码等 实验步骤:详细描述每个实验的步骤和注意事项 实验结果:对实验结果进行分析和讨论,提出改进意见
05
使用指南

数字信号处理-程佩青第三版课件

数字信号处理-程佩青第三版课件

xa(t) 0
xa(nT)
t
2T
0
t
T
这里 n 取整数。对于不同的 n 值,xa(nT) 是 一个有序的数字序列,该数字序列就是离散时间信 号。注意,这里的n取整数,非整数时无定义,另 外,在数值上它等于信号的采样值,即
x(n) xa (nT ), n
离散时间信号的表示方法:公式表示法、图形 表示法、集合符号表示法,如
称该系统是因果系统。 因果系统是指输出的变化不领
先于输入的变化的系统。
对于线性时不变系统,具有因果性的充要条件是 系统的单位取样响应满足:

稳定系统
稳定系统是指对于每个有界输入x(n),都产生有 界输出y(n)的系统。即如果|x(n)|≤M(M为正常数), 有|y(n)|<+∞,则该系统被称为稳定系统。
x(n) ...1,2,3,7,8,9,...
二、常用序列
1. 单位抽样序列(n)
(t) 1/
0 t
(n)
1
0
n
(t)
(1)
t
0
2. 单位阶跃序列u(n)
u(n) 0
u(t)
1

n
0
t
(n)与u(n)之间的关系
令n-k=m,有
3. 矩形序列RN(n)
N为矩形序
列的长度
R4(n)
n 012 3
4. 实指数序列
,a为实数
0<a<1
a>1
n
n
0
0
-1<a<0
a<-1
0
n0
n
a<-1或-1<a<0,序列的幅值摆动

数字信号处理教程 (第三版)程佩青 清华大学出版社dsp-ch5-1

数字信号处理教程 (第三版)程佩青  清华大学出版社dsp-ch5-1
3)存在输出到输入的反馈,递归型结构
X

二、有限阶IIR的表达式:
(其中至少有一个 ak≠0)
Y ( z) 系统函数: H ( z ) X ( z)
N
10 页
bk z k 1 ak z k
k 1
M
M
k 0 N
差分方程: y ( n ) ak y ( n k ) bk x (n k )
2、直接Ⅱ型(典范型)
13 页
只需实现N阶滤波器所需的最少的N个延时单元, 故称典范型。( N M )
X

14 页
直接Ⅰ型与直接Ⅱ型结构比较
1)直接Ⅰ型需要一个相加器,而直接Ⅱ型 需要两个相加器; 2)直接Ⅱ型需要的延迟器比直接Ⅰ型少, 因此所需的存储单元少;
3)从节约存储单元的角度来看,直接Ⅱ型

18 页
X

19 页
N 1 当M=N时,二阶因子配对方式有 ! 种 2 N 1 各二阶基本节的排列次序有 ! 种 2
X

级联型的特点:
20 页
• 调整系数 1k, 2k能单独调整滤波器的第k对零点, 而不影响其它零极点 调整系数1k , 2k 能单独调整滤波器的第k对极点, 而不影响其它零极点
方框图
流图
z
1
z
1
a
a
X
第 5 页
x(n)
z
-1
x(n-1)
x(n)
z -1
x(n-1)
支路増益
x(n)
x(n)
ax(n)
a x1(n)
a
网络 节点
ax(n)
源节点

数字信号处理第三版PPT107页

数字信号处理第三版PPT107页
数字信号处理第三版
31、园日涉以成趣,门虽设而常关。 32、鼓腹无所思。朝起暮归眠。 33、倾壶绝余沥,窥灶不见烟。
34、春秋满四泽,夏云多奇峰,秋月 扬明辉 ,冬岭 秀孤松 。 35、丈夫志四海,我愿不知老。

26、要使整个人生都过得舒适、愉快,这是不可能的,因为人类必须具备一怨环境的心情,化为上进的力量,才是成功的保证。——罗曼·罗兰

28、知之者不如好之者,好之者不如乐之者。——孔子

29、勇猛、大胆和坚定的决心能够抵得上武器的精良。——达·芬奇

30、意志是一个强壮的盲人,倚靠在明眼的跛子肩上。——叔本华
谢谢!
107

程佩青《数字信号处理教程(第三版)》课后习题答案精编版

程佩青《数字信号处理教程(第三版)》课后习题答案精编版
4
第一章 离散时间信号与系统
1 .直接计算下面两个序列的卷积和 y( n ) = x( n )* h( n )
h (n )
=
⎧an ⎨
⎩0
, 0 ≤ n ≤ N −1 , 其他n
x (n )
=
⎧⎪ β ⎨
n−n 0
⎪⎩ 0
,n0 ≤ n , n < n0
请用公式表示。
分析:
①注意卷积和公式中求和式中是哑变量 m ( n 看作参量),
y (n ) ={1,2,3,3,2,1} ;
②δ (n)* x(n) = x(n) , δ (n − m)* x(n) = x(n − m) ;
③卷积和求解时, n 的分段处理。
6
解:(1) y(n) = x(n) * h(n) = R5(n) (2) y(n) = x(n) * h(n) = {1,2,3,3,2,1}
β α
n +1
β α β =
n +1− N −n0
N−
N
α −β
y(n) = Nα n−n0 ,
(α = β )
, (α ≠ β )
如此题所示,因而要分段求解。
2 .已知线性移不变系统的输入为 x( n ) ,系统的单位抽样响应
为 h( n ) ,试求系统的输出 y( n ) ,并画图。
(1)x(n) = δ (n)
∑ ∑( ) n α m−n0 n−m = β α = β m=n0
nn β
n0
α
n β −n0
− β n0
α
β n +1 α
1

β α
α β =
− n +1− n0

数字信号处理程佩青PPT

数字信号处理程佩青PPT
二阶多项式
H ( z ) A k1 1 1 1 k k z z 1 1 2 2 k k z z 2 2 A kH k ( z )
当 M N 时 , 共 有 N 2 1 节
当零点为奇数时:
有一个 2k 0
当极点为奇数时:
有一个 2k 0
H ( z ) A k1 1 1 1 k k z z 1 1 2 2 k k z z 2 2 A kH k ( z )
k 1
k 1
k 1
A为 常 数
M M 1 2 M 2
p k 和 c k 分 别 为 实 数 零 、 极 点 N N 1 2 N 2
q k , q k * 和 d k , d k * 分 别 为 复 共 轭 零 、 极 点
将共轭成对的复数组合成二阶多项式,系数即为实数。 为采用相同结构的子网络,也将两个实零点/极点组合成
N1
H(z) h(n)zn n0
N n 2 1 0 1 h (n )z n h N 2 1 z N 2 1 n N N 1 1 1 h (n )z n 2 令 n N 1 m
N n 2 1 0 1 h (n ) z n z (N 1 n ) h N 2 1 z N 2 1
h(n)偶对称,取“+”
h(n)奇对称,取“
”,且
h
N 1 2
0
N为偶数时
N1
H(z) h(n)zn
N1 2
N1
h(n)zn h(n)zn
n0
n0
nN
N1
2
2
h(n)znz(N1n)
n0
四、数字滤波器的格型结构
格型结构的优点:
1)模块化结构便于实现高速并行处理
2)m阶格型滤波器可以产生1阶到m阶的m个横向 滤波器的输出性能

数字信号处理教程程佩青笔记

数字信号处理教程程佩青笔记

数字信号处理教程程佩青笔记(原创版)目录1.教程概述2.信号处理基本概念3.数字信号处理的基本步骤4.常用数字信号处理方法5.应用案例与实践正文1.教程概述《数字信号处理教程程佩青笔记》是一本针对数字信号处理(Digital Signal Processing,简称 DSP)的教程,适用于电子工程、通信工程、计算机科学等相关专业的学生以及从事信号处理领域的研究人员和工程师。

教程以程佩青教授的课堂笔记为基础,结合实际案例和实践,系统地讲解了数字信号处理的基本概念、原理和方法。

2.信号处理基本概念信号处理是对信号进行操作、变换和处理的过程,其目的是提取有用信息、去除噪声干扰,或者将信号转换为更适合分析、传输或存储的格式。

信号可以分为模拟信号和数字信号,其中模拟信号是连续的,而数字信号是离散的。

数字信号处理就是在数字域中对信号进行操作和处理。

3.数字信号处理的基本步骤数字信号处理的基本步骤包括信号采样、量化、零阶保持、编码等。

首先,通过对连续信号进行采样,将其转换为离散信号;然后,对离散信号进行量化,即将信号的幅度转换为数字表示;接着,对量化后的信号进行零阶保持,以保持信号的连续性;最后,对信号进行编码,以便存储或传输。

4.常用数字信号处理方法常用的数字信号处理方法包括滤波、傅里叶变换、快速傅里叶变换、离散余弦变换、小波变换等。

滤波是信号处理中最基本的方法之一,可以用来去除噪声、衰减脉冲响应等;傅里叶变换和快速傅里叶变换是信号的频域分析方法,可以用来分析信号的频率成分;离散余弦变换和小波变换则是信号的时频分析方法,可以用来同时分析信号的时间和频率信息。

5.应用案例与实践数字信号处理在许多领域都有广泛的应用,如通信、音频处理、图像处理、生物医学工程等。

例如,在音频处理中,可以通过数字信号处理方法对音频信号进行降噪、均衡、混响消除等处理,以提高音质;在图像处理中,可以通过数字信号处理方法对图像信号进行滤波、边缘检测、图像增强等处理,以提高图像质量。

(完整版)数字信号处理(程佩青) 第二章 Z变换PPT文档

(完整版)数字信号处理(程佩青) 第二章 Z变换PPT文档

27
1.围线积分法(留数法)
X(z)zn-1在任意一极点zr处的留数
(1) zr 是X(z)zn-1的单(一阶)极点
R e s X z z n 1 z z r z z r X z z n 1 z z r
(2) zr 是X(z)zn-1的多重(l 阶)极点
R e s X zzn 1 z z r l 1 1 !d d z ll 1 1 z z rlX zzn 1 z z r
其中zi是X(z)的一个r阶极点 ,zk是X(z)的单极点(k=1,2……N-r),Bn是 整式部分的系数(M≥N时存在,M=N时,只有B0 项;M<N时Bn =0)。
33
各系数的求法:
(1)Bn可用长除法求得。
(2) A k 1 z k z 1X z |z z k z z k X z z z z k R e s X z z z z k
式中C是X(z)收敛域内的环绕原点的一条反时针方向的闭合
围线。
25
1.围线积分法(留数法)
比较z变换的定义式 X(z) x(n)zn
和(2.5)式
n
Xz Cnzn n
可以发现,x(n)就是罗朗级数的系数,即:
xn21j
Xzzn1dz
C
CRx,Rx
这就是用围线积分的z反变换公式:
26
1.围线积分法(留数法)
{Z:X(z)存在}=收敛区域。
注意:z变换收敛域的概念很重要,不同的 序列可能有相同的z变换表达式,但是收敛域却 不同。所以应该特别注意,只有当z变换的表达
式与收敛域都相同时,才能判定两个序列相等。
10
2. z变换的收敛域
收敛条件:按照级数理论, X(z) x(n)zn

全套电子课件:数字信号处理(第三版)

全套电子课件:数字信号处理(第三版)

5、本书的主要内容
经典的数字信号处理限于线性时不变系统理 论, 数字滤波和FFT是常用方法。
随机信号处理:基于平稳高斯随机信号 目前DSP研究热点: 时变非线性系统、非平
稳信号、 非高斯信号 处理方法的发展:自适应滤波、 离散小波 变换、 高阶矩分析、盲处理、分形、混沌
理论
课程介绍
基础理论:离散时间信号与系统(ch1)(复习和强化)
(4)可以实现多维信号处理
利用庞大的存储单元,可以存储二维的图像信号或多维的阵列信号,实现二维或 多维的滤波及谱分析等。 4G移动通信:MIMO和OFDM
缺点
(1)增加了系统的复杂性。它需要模拟接口以及比较复杂的数字系统。 (2)应用的频率范围受到限制。主要是A/D转换的采样频率的限制。 (3)系统的功率消耗比较大。数字信号处理系统中集成了几十万甚至更多的晶体管 ,而模拟信号处理系统中大量使用的是电阻、电容、电感等无源器件,随着系统的复 杂性增加这一矛盾会更加突出。
其常中用zZ为[x(复n)变]表量示,对以序其列实x(部n)为的横Z坐变标换,,虚即部为纵坐标构成的平面为z平面。
Z[ x(n)] x(n) z n n
这种变换也称为双边 Z 变换,与此相应还有单边 Z 变换,单边 Z变换只是 对单边序列(n>=0部分)进行变换的Z变换,其定义为
X ( z) x(n) z n n0
上个世纪80年代用Apple II计算机用雷米兹交替算法设计一256阶的FIR滤波 器需要20多小时。
上个世纪90年代已经可以实时地在PC机上实现音视频的编解码。
4、DSP的发展与运用(续)
DSP发展的主要表现: (1) 由 简 单 的 运 算 走 向 复 杂 的 运 算 , 目 前 几十位乘几十位的全并行乘法器可以在数 个纳秒的时间内完成一次浮点乘法运算, 这无论在运算速度上和运算精度上均为复 杂的数字信号处理算法提供了先决条件;

数字信号处理教程第三版程佩青清华大学出版社dsp-ch

数字信号处理教程第三版程佩青清华大学出版社dsp-ch
IIR滤波器和FIR滤波器的设计方法是很不相同的。
X

• 对于线性相位滤波器,通常采用FIR滤波 33 页 器,其单位脉冲响应满足一定条件可以证 明其相位特性在整个频带中是严格线性的, 这是模拟滤波器无法达到的。当然.也可 以采用IIR滤波器,但必须使用全通网络对 其非线性相位特性进行相位校正,这样增 加了设计与实现的复杂性。
X

(2) 四种基本的滤波器
12 页
四种基本滤波器为低通(LP)、高通 (HP)、带通(BP)和带阻滤波器(BRF):
X

(3) 理想滤波器的数字表示
13 页
X
第 14 页
上图给出的是理想数字滤波器幅度特性; 理想滤波器是物理不可实现的;
滤波器的传输函数都是以2π为周期的; 考察频率的范围为-π~π。
阻带最小衰减: 2
2 2 0 lg H H ( ( e e j j0 s t) ) 2 0 lg H (e j s t) 2 0 lg2
其中: H(ej0) 1
当 H (e j c)2 /2 0 .7 0 7 时, 13dB
称 c 为3dB通带截止频率
X

4、表征滤波器频率响应的特征参量
• 特点是输入信号中有用的频率成分和希望滤 除的频率成分各占有不同的频带,通过一个 合适的选频滤波器达到滤波的目的。例如, 输入信号中含有干扰,如果信号和干扰的频 带互不重叠,可滤除干扰得到纯信号。
X
第 6

现代滤波器 :
• 如果信号和干扰的频带互相重叠,这时需要采用
另一类所谓的现代滤波器,例如维纳滤波器、卡
(ej)2 1jlnH H *((eej j))21jlnH H ((zz)1)zej
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
return 7
•定点数作加减法时结果可能会超出±1,称为 “溢出”; •乘法运算不溢出,但字长要增加一倍。 为保证字长不变,乘法后,一般要对增加的 尾数作截尾或舍入处理,带来误差。 缺点:动态范围小,有溢出;衰减比例系数不 好确定。
return 8
定点数的表示分为三种(原码、反码、补码):
设有一个(b+1)位码定点数: β0β1β2┄βb,则 ①原码表示为
另外,数字信号是通过采样和转换得到的,而转换的位数 是有限的(一般6、8、10、12、16位),所以存在量化误差 ,另外,计算机中的数的表示也总是有限的,经此表示的滤 波器的系数同样存在量化误差,在计算过程中因有限字长也 会造成误差。
return 4
量化误差主要有三种误差: ①A/D变换量化效应---信号采集时—抽样定理; ②系数的量化效应---系统函数分子分母系数的量化 表示; ③数字运算的有限字长效应--乘法运算--乘积的有效 位数比每个因子都增加,须截短或舍入。
return 5
§8.2 量化与量化误差
有限字长的二进制数表示数字系统的误差源: ①对系统中各系数的量化误差(受计算机
中存贮器的字长影响) ②对输入模拟信号的量化误差(受A/D的精
度或位数的影响) ③运算过程误差,如溢出,舍入及误差累
积等(受计算机的精度影响)
return
6
8.2.1 二进制数的表示
8.2.2 定点制的量化误差
定点制中的乘法,运算完毕后会使字长增加, 例如原来是b位字长,运算后增长到b1位,需对 尾数作量化处理使b1位字长降低到b位。
量化处理方式: 截尾:保留b位,抛弃余下的尾数; 舍入:按最接近的值取b位码。 两种处理方式产生的误差不同,另外,码
制不同,误差也不同。
return 15
N
ai Z i
H(z)
i0 N
1 bi Z i
i 1
为对于I I R这滤样波一器个形系式统,,{ 也b可}i都用为差0分时方就程是来一表个示FIR:滤波器。
N
N
y(n) ai x(n i) bi y(n i)
i0
i 1
return 2
x(n)
DF
y(n)
IIR、FIR的系统函数 网络结构形式(直接、串、并、格型结构)
正数表示:0.101 其反码为:1.010
return 10
③补码表示(正数同原码,负数则将原码中的尾 数求反加1)bx 0 i 2 i 1例: x 0.75
正数表示:0.110 取反:1.001
x 的补码:1.010
return 11
补码加法运算规律: 正负数可直接相加,符号位同样参加运算,
第八章 数字信号处理中的有限 字长效应
8.1 引言 8.2 量化与量化误差 8.3 A/D变换的量化效应 8.4 乘积误差对数字滤波器有限字长运算的影响 8.5 极限环振荡 8.6 系数量化对系数滤波器的影响
1
8.1 引言
数字滤波器的实现方法(数字计算过程):
a. 利用专用计算机(DSP系统); b.直接利用计算机和通用软件编程实现。 一个数字滤波器的系统函数一般可表示为有理函数形式:
x 0.0101 2011就是非规格化表示。
为充分利用尾数的有效位数,规格化表示为x 0.101 2010
13
return
优点: 动态范围大,一般不溢出. 缺点: 相乘、相加,都要对尾数处理作量化处理。
一般,浮点数都用较长的字长,精度较高,所 以我们讨论误差影响主要针对定点制。
return 14
“量化宽度”或“量化阶” q=2-b :代表b位字长可表 示的最小数。
一般 2-b1<<2-b, 因此正数的截尾误差为
-q≤ET≤0
return 17
2)负数
负数的三种码表示方式不同,所以误差也不同 。
原码(β0=1):
x
b1
i 2i
i1
[x]T
i
b
1
i
2
i
b1
ET x T x i 2i
1、截尾处理:
1)正数(三种码形式相同)
一个b1位的正数 x 为:
b1
x i 2i i 1
用[·]T表示截尾处理,则
b
[x]T i 2 i i 1
return 16
截尾误差
ET [x]T
x
i
b1
b
1
i
2
i
可见,ET≤0,βi全为1时,ET有最大值,
ET
b1
2
i
(2b
2b1)
ib1
b
x (1) 0 i 2i i 1
例:1.111→-0.875 , 0.010→0.25
原码的优点是乘除法运算方面,而加减法运算要增加 时间。P395
return 9
②反码表示:(正数同原码,负数则将原码中 的尾数按位求反)
b
x 0 (1 2b ) i 2i i 1
例: x 0.625
如符号位发生进位,进位的 1 丢掉。 负数以补码形式表示的原因是 将减法运算变为补码的加法运算。
return 12
(2)浮点表示 x M 2c
1 M 1
2
尾数 指数 阶数
浮点制运算: 相加 对阶
相加
归一化,并作尾数处理 相乘 : 尾数相乘, 阶码相加, 再作截尾或舍入。
尾数M的取值在[0.5, 1),只是要求规格化表示。
(1)定点表示
0 • 12 b
• 整个运算中,小数点在数码中的位置固定不变,
原则上小数点在数码中的位置是任意的,称为定点
制;e.g.六位字长:10.1001
•通常定点制总是把数限制在±1之间; 最高位为符 号位,0为正,1为负,小数点紧跟在符号位后;数 的本身只有小数部分,称为“尾数”;
•若数值较大时,可乘上一个衰减因子,保证该数在 运算中不超过1;运算后再除以该因子还原。
i b 1
0≤ET≤q
return 18
补码( 0 1)
b1
x 1 i 2i i 1 b
xT 1 i 2i i 1
b
b1
软、硬件实现
return 3
即一个输出序列是其过去 点N输出值的线性组合加上当前输 入序列与过去 点输N入序列的线性组合。 除了y与(n当) 前的输
入 有关,x(同n)时还与过去的输入和过去的输出有关,系统是 带有记忆的。
对于上面的算式,可以化成不同的计算形式,如直接计算 、分解为多个有理函数相加、分解为多个有理函数相乘等等 ,不同的计算形式也就表现出不同的计算结构,而不同的计 算结构可能会带来不同的效果,或者是实现简单,编程方便 ,或者是计算精度较高等等。
相关文档
最新文档