Keskustelu: Ohjelmointikysymykset: C: Modified discrete cosine transform

// Modified Discrete Cosine Transform v0.1

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <string.h>

int main(void) {
   const char str[] = "abbaKabba"; // data

   int i, n, N = strlen(str)/2, mid = 0;
   double *real = calloc(sizeof(double), N), pi_n = 3.14159265359 / N, tmp_l;

   for (i = 0; i < 2*N; i++) {
      mid += str[i];
      printf("%i ", str[i]);
   } mid /= i; printf("\n\nN = %i, mid = %i\n\n", N, mid);

   // vähennetään datasta sen keskiarvo, ettei amplitudit kasva niin suuriksi

   for (i = 0; i < N; i++) { // muunnos looppi
      for (n = 0; n < 2*N; n++)
         real[i] += (str[n] - mid)*cos(pi_n * (n + 0.5 + N/2) * (i + 0.5));
      printf("\t%.3f\n", real[i]);
   }

   for (i = 0; i < 2*N; i++) { // käänteismuunnos looppi
      for (tmp_l = n = 0; n < N; n++)
         tmp_l += real[n] * cos(pi_n * (i + 0.5 + N/2) * (n + 0.5));
      tmp_l = tmp_l / N + mid;

      printf("%.3f => %X = %c\n", tmp_l, (char) round(tmp_l), (char) round(tmp_l));
   }

   free(real);
   return 0;
}