xing.c

离散傅立叶变换 对48个离散的数值分析

#define N 48
#include "stdio.h"
#include<math.h>
#include "conio.h"
void main()
{
int k=1,i=0;
float x[48];
float X,q;

float XR=0;
float XI=0;
FILE *fp;
printf("input num k?\n");
scanf("%d",&k);
if((fp=fopen("in.txt","r"))==NULL){
    printf("can not open file!");
     getch();
     exit(0);
    }
    else {
    for(i=0;i<48;i++){
            fscanf(fp,"%f",&x[i]);
            printf("%f ",x[i]);
            if((i+1)%6==0)printf("\n");
            }
            getch();

    }


        for(i=0;i<=47;i++)
    {
        XR=XR+x[i]*cos(6.28*k*i/N);
        XI=XI-x[i]*sin(k*i*6.28/N);
    }
XR=XR*2/N;
XI=XI*2/N;

X=sqrt(XR*XR+XI*XI);


      if(XI==0)
      {
       printf("error\n");
        getch();
      }

else {q=atan(XR/XI);
       printf("\nOK!GO ON.\n");
        getch();
    }

printf("k=%d\n",k);
  getch();
printf("fuzhi:%f\n",X);
  getch();
printf("xiangwei:%f\n",q);
  getch();
}