5ndft_vhdl/radix2.m

function output = radix2(input, wn)
output = 0*input;
for i = 1:size(input,1)/2

    s1 = [input(i,:); input(i+(size(input,1)/2),:)*wn(i)];
    output(i,:) = s1(1,:)+s1(2,:);
    output(i+(size(input,1)/2),:) = s1(1,:)-s1(2,:);
end