Subversion Repositories dblclockfft

% Read the file

fid = fopen('fft_tb.dbl','r');

ftlen = fread(fid, 1, 'int32');

raw = fread(fid, [2 inf], 'double');

fclose(fid);

% Convert the raw doubles into complex values

datc = raw(1,:)+j*raw(2,:);

% Reshape the matrix into one line per FFT

ndat = reshape(datc, ftlen*2, floor(length(datc)/(ftlen*2)));

truth  = ndat((ftlen+1):(2*ftlen), :);

output = ndat(1:ftlen,:);

% Create a time axis, for use in plotting if desired

tm = 0:(ftlen-1);

% Now, the data from the test is ready for inspection

k = 1; plottst; k