0001 N=2048; original fft size 0002 pfbTaps=4 0003 M = pfbTaps.*N; % Window length (we'll use 0004 0005 n = [0:M-1]; % time indices for sinusoid and FFT 0006 wxT = 2*pi/(4); % Sinusoid frequency in rad/sample (1/4 sampling rate) 0007 A = 1; % Sinusoid amplitude 0008 phix = 0; % Sinusoid phase 0009 0010 0011 x = randn(1,M)+A * exp(j*wxT*n+phix); % the complex sinusoid itself: [1,j,-1,-j,1,j,...] 0012 wf=hamming(N.*pfbTaps) 0013 plot(abs(fft(x,N)))