clear all; close all; Fs=8000; N=0.5; f1=440;f2=440; cut1=5; cut2=0; x=0:1:N*Fs-1; t=x/Fs; y1=cos(2*pi*t*f1); y2=cos(2*pi*t*f2); myx=(1:length(1))/length(1)*Fs; plot(myx,abs(fft(y1))); myOut=[]; for idx=1:5; myOut=[ myOut y1(1+cut1:4000-cut2) y2(1+cut1:4000-cut2)]; end; plot(myOut(3900:4100),'-o'); sound(myOut, Fs);