Scribd
Upload
88 views4 pages

"DSP2833x - Device.h" "Stdio.h" "Stdlib.h" "Math.h": #Include #Include #Include #Include

This C code implements a basic radix-2 FFT algorithm. It initializes global variables including input and output arrays. It calls functions to perform bit reversal ordering, one butterfly computation, and the overall FFT. The main function initializes system controls, sets the input array, calls the FFT, and enters an infinite loop.

Uploaded by

Trupti Pawar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
88 views4 pages

"DSP2833x - Device.h" "Stdio.h" "Stdlib.h" "Math.h": #Include #Include #Include #Include

This C code implements a basic radix-2 FFT algorithm. It initializes global variables including input and output arrays. It calls functions to perform bit reversal ordering, one butterfly computation, and the overall FFT. The main function initializes system controls, sets the input array, calls the FFT, and enters an infinite loop.

Uploaded by

Trupti Pawar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

EXP-11

#include "DSP2833x_Device.h"
#include "stdio.h"
#include "stdlib.h"
#include "math.h"

// external function prototypes


extern void InitSysCtrl(void);
extern void InitPieCtrl(void);
extern void InitPieVectTable(void);

void fft( int n, float x[],float y[]);


void bitreversal(int n, float x[], float y[]);

// Global Variables
int N= 4;
int M =4;
float add=0,multi=0;
int i,j,k,n;
float x[8],y[8],real[8],imag[8];

//###########################################################################
// main code

//###########################################################################

void main(void)
{
//###########################################################################
// System Initialization

//###########################################################################

InitSysCtrl(); // Basic Core Init from DSP2833x_SysCtrl.c


EALLOW;
SysCtrlRegs.WDCR= 0x00AF; // Re-enable the watchdog
EDIS; // 0x00AF to NOT disable the Watchdog, Prescaler = 64
DINT; // Disable all interrupts
InitPieCtrl(); // basic setup of PIE table; from DSP2833x_PieCtrl.c
InitPieVectTable(); // default ISR's in PIE

//###########################################################################
// FFT

// Basic one butterfly Cooley Tukey radix-2 FFT


//###########################################################################

for(i=0;i<N;i++)
{
x[i]=i+1;
y[i]=0;
}

fft(N,x,y);

while(1)
{}

void bitreversal(int n, float x[], float y[])


{

int i,j,k,n1;
j=0;
n1=N-1;
for(i=0;i<n1;i++)
{
float temp;

if(i>=j) goto end;

temp=x[j];
x[j]=x[i];
x[i]=temp;

temp=y[j];
y[j]=y[i];
y[i]=temp;

end: k=n/2;

kk: if(k>j) goto kt;


j=j-k;
k=k/2;

goto kk;

kt: j=j+k;
}

/****************************************/

void fft( int n, float x[],float y[])


{
int i,j,k;
int n1,n2,l;
float a,c,s,e;
float xt,yt;

n2=n;
for(k=0;k<M;k++)
{

n1=n2;
n2=n2/2;
e=6.283185307179586/n1;

for(j=0;j<n2;j++)
{
a = j*e;
c = cos(a);
s = -sin(a);

for(i=j;i<n;i+=n1)
{

l=i+n2;
xt=x[i]-x[l]; add++;
x[i]=x[i]+x[l]; add++;
yt=y[i]-y[l]; add++;
y[i]=y[i]+y[l]; add++;

x[l] = xt*c - yt*s; add++; multi += 2;


y[l] = xt*s + yt*c; add++; multi += 2;
}
}
}

bitreversal(n, x, y);

You might also like