yading@10: N=1024; yading@10: t = [0:1/44100:0.1-1/44100]; yading@10: x1=cos(2*pi*100*44100/1024*t); yading@10: x801 = x1(1:801); yading@10: xlength=length(x801); yading@10: fftbuffer=zeros(N, 1); yading@10: fftbuffer(1:(xlength+1)/2) = x801((xlength+1)/2:xlength); yading@10: fftbuffer(N-(xlength-1)/2+1:N) = x801(1: (xlength-1)/2); yading@10: plot(unwrap(angle(fftbuffer,1024)));