1.问题
上一篇中介绍了利用循环卷积计算线性卷积,但是当两个卷积项长度相差很大时,短序列需要补的0非常多,这样无助于计算量的减小,并且这种方法需要等长序列全部输入完之后才开始计算,这会造成输出有很大的延时。
但是实际中这种现象很长见,比如对音频信号进行数字滤波,滤波器抽头系数点数较少,而MIC一直在录音,需要滤波的信号可以认为是无限长的,如果等录很久后才开始进行滤波,则会造成需要存储的数据量急剧增大,同时输出的信号也会有较长的延时,这个时候怎样保证延时小的同时输出的结果跟对长信号滤波的结果相同呢
滤波器系数h(n), h ( n ) ,长度为N N,输入信号x(n) x ( n ),长度为无限长,输出y(n) y ( n )即为它们的卷积
y(n)=∑Mm=0x(n−m)h(m) y ( n ) = ∑ m = 0 M x ( n − m ) h ( m )
那怎么利用频域计算呢,两种方法,重叠相加(overlap-add)和重叠保留(overlap-save)
重叠相加法(OLA)
将x(n) x ( n )分段,每段长为M M,保证M M接近N N即可,然后将xk(n) x k ( n )补零延长到L=M+N−1 L = M + N − 1,计算L L点FFT F F T得到XK(K) X K ( K )
将h(n) h ( n )补零延长至L=M+N−1 L = M + N − 1,计算L L点FFT F F T得到H(K)。 H ( K ) 。
计算yk(k)=XK(K)∗H(K) y k ( k ) =