基于卷积同态系统解决混响问题.pdfVIP

New新Te技ch术nology 基于卷积同态系统解决混响问题 口丁海峰 袁展光 海军蚌埠士官学校 一 、 引言 n2=24，MATLAB程序如下 ： 在通话环境声吸收缺乏的情况下进行语音通信时，混响 a=0．4；al=0．5；a2=O．7；N=64；n0=8；n1=16；n2=24； 会导致接受的语音信号含糊不清 ，混响严重时甚至无法正常 forn=l：nO ％生成信号 xl(n) 通话 1【]，因此 ，如何减少混响对声音接收系统的影响是一个 xl(n)=power(a，n-1)；end 非常重要环节。本文利用了卷积 同态系统 对混响信号进行 forn=(n0+1)：N 了有效的恢复 。 x1(n)=O；end 二 、信号恢复的方法 x2(1)=l；x2(nl+1)=a1；x2(n2+1)=a2；％生成信号 x2(n) 假设一个信号可以表示成某基本波形经不同延 时后得 x=conv(xl，x2)；xx(1)=log(x(1))； ％生成复倒谱 xx(n) 到的若干个互相交叠的复现波形之和，我们将这类信号的失 真叫做混响。它的数学模型可以表示为 ：首先我们给出序列 forn=l：nl一1 ％生成窗函数 L(n) x(n)，它是 由有用信号 x，(n)经过不同延时并乘以不同的系数 L(n)=l；end 之后叠加得到 的波形 ： forn=n1：N x()=X (，z)+∑ax．( 一n)：xl(，z)木 ：() L(n)=0；end yy=xx．*L； ％线性滤波 ： ()=6(，z)+∑a6(一n) y(1)=exp(yy(1))； ％复倒谱 的逆运算 设 xj(n)：10， l，…。-1l ， n=n，…，N一1 l ee=sum(abs(y—x1)) ％误差函数 2 ()：6(n)+alfi(n—n1)+a28(n—n2) 程序仿真运行结果如下 ： 则其复倒谱为当n0时 ， ㈤ = 一 ㈣ k=O n 当n：0时 ， ()=lg[x(O] 当 n0时 ， f，z=0 我们可 以用一个线性非频变滤波器将这两个变量分离出 来 。由于频域和时域之间存在着对称关系，上述两个分量变 化的速率就体现在时域上。所 以，想要得到 XIn1分量 ，可 以用一个只让复倒谱短时分量通过 的滤波器进行滤波 ，该滤 波器可 以设计成一个短时通 窗 ： ()：J0【， o’l’…， n n … N 一1 ， 相应 的时域表达式为 ： ()= ，() () 。 当y()经过 的卷积 同态系统特征系统 D+的逆系统 D．1 后，我们理论上将恢复 (，z)=X(，z)。因为 )是无限时宽 的，乘 以，f，z1后相 当于进行 了截断操作 ，这也使y(一)的长度 为 N。事实上 ，对于长度为 N的输入序列 ，只需要 N个采样 图2 信号 y(n)与信号 X1(n)之间误差函数 即可确定。