现代信号处理_第13讲.pptVIP

第七章 同态信号处理 ; 狭义线性系统概念可推广到广义线性系统范畴 若干个元素x1，x2，?构成一个集合 G，元素间定义一种运算 ?。G 被称为群需满足4个条件： 封闭性：x1?x2?G 结合性：x1?(x2 ?x3)= (x1?x2)?x3; 群G 构成了一个矢量空间 ;;;;; 乘法同态系统：输入、输出矢量空间矢量间的运算是乘法运算(?、?为乘法运算，?、?为指数运算) 根据相乘信号的形式，分为： 实数乘法同态系统 复数乘法同态系统 ;;;;;;(2)、复数乘积同态系统; 第三个子系统T-1?[?]为对复数进行指数运算： ＝ln[|y(n)|]+jarg[y(n)] 则有： y(n)=exp[ ]=exp{ln[|y(n)|]+jarg[y(n)]} =exp{ln[|y(n)|]}exp{jarg[y(n)]} =|y(n)|exp{jarg[y(n)]} 对于ln[x1(n)x2(n)]，有： ln[x1(n)x2(n)]=ln[|x1(n)x2(n)|]+jarg[x1(n)x2(n)] 其中ln[|x1(n)x2(n)|]=ln[|x1(n)|]+ln[|x2(n)|] 而arg[x1(n)x2(n)]却不一定等于arg[x1(n)]+arg[x2(n)] 因此对复数乘积同态系统，需注意复对数的多值性和解析性的问题 这个问题将下面的内容中一并讨论; 信号和噪声的关系除相加、相乘外，可以为卷积： 语音信号是声带源和声道冲激响应的卷积 地震波是地震源波形和地壳冲激响应的卷积 处理这类信号，使用卷积同态系统 卷积同态系统：输入、输出矢量空间中矢量间的运算是卷积运算(?、?为卷积运算);;;;;;;;;; 对卷积同态系统第一个子系统的时间序列x(n) 先进行傅里叶变换，再取对数，然后求反傅里叶变换 或先进行Z变换，再取对数，然后求反Z变换 第一个子系统的输出 称为x(n)的复倒谱 ; 推论： X(z)的对数 不是唯一的(复对数的多值性问题) ; 对 求反Z变换，得到x(n)的复倒谱 在 的收敛域内， 是z的解析函数 若 是稳定、因果的， 的收敛域是某个圆外且包括单位圆 ? 在单位圆上是解析的 现在问题就变为要求 是?的连续函数;其中A为常数 N1、N2为单位圆内、单位圆外的零点数目 N3、N4为单位圆内、单位圆外的极点数目 ai、bi、ci、di的模均小于1; 第二项ln zM的反Z变换为： ; 结论：复倒谱具有以下性质 1．复倒谱的长度总是无限长的 2．复倒谱??幅度至少按1/|n|的速度衰减 3．最小相位序列在单位圆外无零点和极点， 即：N2＝N4＝0，因此其复倒谱是因果序列;最小相位序列的复倒谱是因果序列： ; 的傅里叶变换： ;;②直接计算法 先求x(n)的N点离散傅里叶变换X(k)，k=0～N-1，然后求X(k)的复对数 ，最后利用离散傅里叶反变换得到 采用N点离散傅里叶变换和离散傅里叶反变换，算出的 是理论值 以N为周期延拓的结果 由于 是无限长序列，以N为周期进行延拓必然会造成混叠 好在 的幅度至少按1/|n|的速度衰减，只要N取得足够大，这种混叠造成的失真还是可以忽略的 当N不够大时，可以在x(n)序列后面补零后进行计算，从而使算出的 更接近理论值 实际应用中可使用快速傅里叶变换;;;;;;;