Code covered by the BSD License  

Highlights from
MatLab Solutions: "Introduction to Digital Signal Processing: A Computer Laboratory Textbook".

image thumbnail

MatLab Solutions: "Introduction to Digital Signal Processing: A Computer Laboratory Textbook".

by

 

29 Oct 2012 (Updated )

These files are the MatLab solutions of exercises contained in the above DSP lab textbook.

radix4fft(x)
function X = radix4fft(x)
% Radix-4, Decimation-In-Time (DIT), FFT Computation function.
% Attention: Length of input signal x should be divisable by 4 and >=4.
% Attention: Only the first iteration of the algorithm is implemented here.

% Make sure x is row vector for subsequent processing:
x = x(:).';
N = length(x);

% Generate sequences e[n], f[n], g[n] and h[n].
e = zeros(1,N/4);
f  = zeros(1,N/4);
g = zeros(1,N/4);
h = zeros(1,N/4);

ind1 = 1:4:N-3;
ind2 = 2:4:N-2;
ind3 = 3:4:N-1;
ind4 = 4:4:N;

e = x(ind1); % N/4 point signal made of samples: x[0], x[4], x[8],   ..., x[N-4]. 
f = x(ind2);  % N/4 point signal made of samples: x[1], x[5], x[9],   ..., x[N-3].
g = x(ind3); % N/4 point signal made of samples: x[2], x[6], x[10], ..., x[N-2]. 
h = x(ind4); % N/4 point signal made of samples: x[3], x[7], x[11], ..., x[N-1].

% Compute their N/4-point DFT's:
E1 = fft(e);
F1 = fft(f);
G1 = fft(g);
H1 = fft(h);

% Create the periodic extensions of total size N of these 4 DFT's:
E = [E1 E1 E1 E1];
F = [F1 F1  F1 F1];
G = [G1 G1 G1 G1];
H = [H1 H1 H1 H1];

% Create the phasor vector:
 k = 0:N-1;
 W = exp(-1i*2*pi*k/N);

% Apply the formula.
X = E + W.*F + (W.^2).*G + (W.^3).*H;

Contact us