test.cpp

很好的数学函数库

// Test.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include "math.h"
#include "stdlib.h"
#include "time.h"
#include "memory.h"

#define PI (3.1415926535897932384650288419716039987510)
#define N (16)
int X[N];

int main(int argc, char* argv[])
{
  srand((unsigned)time(NULL));
  int i;
  for(i=0; i<N; i++){
    X[i]=rand();//i+1;
  }

  float x[N];
  memset(x, 0, sizeof(x));
  int k;
  int n;
  for(k=0; k<N; k++){
    for(n=0; n<N; n++){
      x[k]+=X[n]*(cos(2*PI*n*k/N));
    }
    //x[k]/=N;
  }

  float y[N];
  memset(y, 0, sizeof(y));
  for(n=0; n<N; n++){
    for(k=0; k<N; k++){
      y[n]+=x[k]*(cos(2*PI*n*k/N)+sin(2*PI*n*k/N));
    }
    y[n]/=N;
  }

/*
  float f[N];
  memset(f, 0, sizeof(f));
  for(k=0; k<N; k++){
    f[k]=X[k]/y[k];
  }
*/

  return 0;
}