% Frequency response of "echo" channel in mus.m

f=linspace(500,550,501);

H = 1 + 0.5*exp(j*2*pi*0.4*f);
aH = abs(H);
pH = angle(H)/pi*180;

subplot(211)
plot(f, aH);
