load aiueo
Fs=11025;
clear h;
h(11025)=0.3;
h(5513)=0.5;
h(1)=1;
echoaiueoc=conv(aiueo,h);
wavplay(echoaiueoc/max(abs(echoaiueoc)))
subplot(2,1,1);plot((1:len)/Fs,aiueo);
axis([0 length(echoaiueoc)/Fs -1 1]);
subplot(2,1,2);plot((1:length(echoaiueoc))/Fs,echoaiueoc);
axis([0 length(echoaiueoc)/Fs -1 1]);