main.c

来自「DSP TMS320vc5416 FFT代码」· C语言 代码 · 共 66 行

	#include "stdio.h"
	#include "math.h"

void main(void)
{
	int i,j;
	double pr[64], pi[64],fr[64],fi[64];

	for (i=0;i<64;i++)
	{
		//pr[i]=exp(-0.1*(i+0.5)); 
		pr[i] = sin(i*2) +sin(i*20);
		pi[i]=0.0;
	}

    printf("\n");

    for (i=0; i<=15; i++)
    { 
    	for (j=0; j<=3; j++)
    		printf("%e   ",pr[4*i+j]);
        printf("\n");
    }

    printf("\n");

    kfft(pr,pi,64,6,fr,fi,0,1);

    for (i=0; i<=15; i++)
    { 
    	for (j=0; j<=3; j++)
        	printf("%e   ",fr[4*i+j]);
        printf("\n");
	}
    printf("\n");
    for (i=0; i<=15; i++)
	{ 
		for (j=0; j<=3; j++)
        	printf("%e   ",fi[4*i+j]);
		printf("\n");
	}
    printf("\n");
    for (i=0; i<=15; i++)
	{ 
		for (j=0; j<=3; j++)
			printf("%e   ",pr[4*i+j]);
        printf("\n");
	}
    printf("\n");
    for (i=0; i<=15; i++)
	{ 
		for (j=0; j<=3; j++)
			printf("%e   ",pi[4*i+j]);
		printf("\n");
	}
    printf("\n");
    kfft(fr,fi,64,6,pr,pi,1,1);
    for (i=0; i<=15; i++)
	{ 
		for (j=0; j<=3; j++)
			printf("%e   ",fr[4*i+j]);
		printf("\n");
	}
	printf("\n");
}