Ваша модель похожа на КИХ-фильтрацию, поэтому я думаю, что самый простой способ - использовать функцию filter
.
Вам необходимо определить фильтр длины d
с правильными коэффициентами:
b_echo = zeros(1,d);
b_echo(1) = 1; % original signal
b_echo(d) = k; % attenuation of the echo
y = filter(b_echo,1,x);