语音录音包括由墙壁反射引起的回声。用真正的倒频谱把它过滤掉。

在录音中，一个人说了单词MATLAB®。加载数据和采样率， $$F_{\mathrm{年代}}＝7418Hz$$．

负载mtlb%听到，键入soundsc(mtlb,Fs)

通过在记录中添加一个延迟的信号副本来模拟回声 $$\Delta$$样本和衰减的已知因素 $$\alpha$$： $$y（n）＝x（n）+\alpha x（n-\Delta ）$$．指定时间延迟为0.23 s，衰减系数为0.5。

时滞= 0.23;δ=圆(Fs *时滞);α= 0.5;源自= [mtlb; 0(δ1)];回声=[0(δ1);mtlb] *α;mtEcho = origin + echo;

绘制出原始信号、回波和最终信号。

t =(0:长度(mtEcho) 1) / Fs;Subplot (2,1,1) plot(t，[origin echo]) legend(“原始”，“回声”)子图(2,1,2)图(t,mtEcho)图例(“总”)包含(“时间(s)”）

%听到，键入soundsc(mtEcho,Fs)

计算信号的真实倒谱。绘制倒谱并注释其最大值。当回声开始到达时，倒频谱有一个尖峰。

c = rceps (mtEcho);(px, loc) = findpeaks (c,“阈值”, 0.2,“MinPeakDistance”, 0.2);clf情节(t、c、t (loc), px,“o”)包含(“时间(s)”）

通过IIR系统对信号进行滤波，消除回波。 $$w$$,遵循 $$w（n）+\alpha w（n-\Delta ）＝y（n）$$．将滤波后的信号与原始信号进行比较。

dl = loc (2) 1;mtNew = filter(1，[1 zeros(1,dl-1) alpha]，mtEcho);次要情节(2,1,1)情节(t,源自)传说(“原始”) subplot(2,1,2) plot(t,mtNew) legend(“过滤”)包含(“时间(s)”）

%听到，键入soundsc(mtNew,Fs)