语音记录包括由墙壁反射引起的回波。使用真正的凯斯特鲁姆来过滤掉。

在录音中，一个人说Matlab®这个词。加载数据和采样率， $$F_{\mathrm{S.}}=7.4.18.H\mathrm{Z.}$$。

加载MTLB.％听到，键入soundsc（MTLB，FS）

通过添加到记录延迟的信号副本来模拟回声 $$\mathrm{\delta .}$$样品并通过已知因素衰减 $$\mathrm{\alpha .}$$： $$y（\mathrm{N.}）=X（\mathrm{N.}）+\mathrm{\alpha .}X（\mathrm{N.}-\mathrm{\delta .}）$$。指定0.23秒的时间滞后和衰减因子为0.5。

timelag = 0.23;delta = round（fs * timelag）;alpha = 0.5;ORIG = [MTLB; Zeros（Delta，1）];echo = [零（delta，1）; mtlb] * alpha;mtecho = orig +回声;

绘制原始，回声和所产生的信号。

t =（0：长度（mtecho）-1）/ fs;子图（2,1,1）图（t，[over echo]）图例（'原来的'那'回声'）子图（2,1,2）绘图（T，MTecho）传奇（'全部的'）xlabel（'时间''）

％听到，键入soundsc（mtecho，fs）

计算信号的真实倒谱系。绘制薄型克拉姆并注释其最大值。综注在回声开始到达时具有尖峰。

C = RCEPS（MTECHO）;[px，locs] = findpeaks（c，'临界点'，0.2，'minpeakdistance'，0.2）;CLF图（T，C，T（LOC），PX，'o'）xlabel（'时间''）

通过输出的IIR系统过滤信号来取消回声， $$\mathrm{W.}$$，obeys. $$\mathrm{W.}（\mathrm{N.}）+\mathrm{\alpha .}\mathrm{W.}（\mathrm{N.}-\mathrm{\delta .}）=y（\mathrm{N.}）$$。绘制过滤的信号并将其与原件进行比较。

dl = locs（2）-1;mtnew =滤波器（1，[1，1，DL-1），Mtecho）;子图（2,1,1）绘图（t，orig）图例（'原来的'）子图（2,1,2）绘图（T，MTNew）图例（'过滤'）xlabel（'时间''）

％听到，键入soundsc（mtnew，fs）