makewavesc.c

基于TS201 32-bit floating point DFT routine

/*______________________________________________________________________
MAKEWAVESC.C	   Test Signal Generator for Floating Point Data
______________________________________________________________________*/


#include <stdio.h>
#include <math.h>
double x[4196];  /*made global for a huge array*/
double s[4196];
double c[4196];
main()
{
  int i, n, k=0;
  char filename[25];

  double pi;
  FILE *out_file;


  printf("%c%c%c%c",27,91,50,74); /*clear screen*/
  printf("%c%c%c",27,91,72);      /*curser home*/
  printf("Interleaved Sin and Cosine Signal Generator for Optimized DFT Routine.\n");
 
  
  printf("\nEnter the output filename -------------> ");
  scanf(" %s",filename);
  out_file=fopen(filename,"w");

  printf("Enter the number of points (max 4196) -> ");
  scanf(" %d",&n);
  
  pi=3.141592654;

     
  for (i=0; i <= n-1; i++)
    {
      s[i] = (-1)*sin(2.0*pi*(double)i/(double)n);
    }
      
  for (i=0; i <= n-1; i++)
    {
      c[i] = cos(2.0*pi*(double)i/(double)n);
    }

  for (i=0; i <= n-1; i++)
    {
      x[k] = s[i];
      x[k+1] = c[i];
      k = k + 2;
    }
            /*write results*/
  printf("\n");
  for (i=0; i <= 2*n-1; i++)
    {
      fprintf(out_file,"%15.10e\n",x[i]);
      printf("%d : %15.10e\r",i,x[i]);

    }
  fclose(out_file);
  printf("\nFinished\n");
}