diskrete Cosinus Transformation DCTII

[Marian G. (houwdy)]

Hallo Forengemeinde.

ich hab ein Problem an dem ich jetzt schon ewig hänge und einfach nicht 
weiterkomme.

Und zwar geht es um eine Hausübung.
Aufgabenstellung ist eine diskrete Cosinustansformation zu schreiben.
zuerst war das ganze in C gefragt und jetzt mit der C++ syntax 
allerdings nicht objektorientiert.

Soweit so gut. Ich hab das ganze programmiert mit Codeblocks compiliert 
aber ich bekomme andere Werte. sie sind zwar ähnlich von der Struktur, 
aber eben falsch...

könnte mir bitte jemand kurz die zwei Funktionen vergleichen und mir auf 
die Sprünge helfen? ich hab die beiden Projekte noch einmal komplett im 
Anhang beigefügt für den Fall das es Fragen gibt.

hier der funktionierende Code

 #include "dct.h"
#include <stdlib.h>
#include <stdio.h>
#include <math.h>

// Function to compute the discrete cosine transform DCT-II of a
// dataseries and store the result in another dataseries of equal size
void dct_ii (dataseries p, dataseries q) {
  // Copy dataseries size to time saver, check for equal sizes
  int n = p.size;
  if (q.size != n) {
    fprintf (stderr, "dct_ii: dataseries sizes don't match, exiting\n");
    exit (-1);
  }
  // Perform actual DCT-II
  const double pi = 3.1415926535897932;
  for (int i=0; i<n; ++i) {
    q.value[i] = 0.0;
    for (int j=0; j<n; ++j) {
      q.value[i] = q.value[i] + cos (pi / n * i * (j + 0.5)) * p.value[j];
    }
  }
  return;
}

und hier der mit falschen Ergebnissen.

 #include <iostream>
#include <vector>
#include <cstdio>
#include <math.h>

using namespace std;

void dct(vector<double> &eingabe, vector<double> &ausgabe) // Berechnung der Cosinustransformation
{
  int n = eingabe.size();
  if (ausgabe.size() != n) {
    cout<<"dct_ii: dataseries sizes don't match, exiting"<<endl;
    return;
  }
  // Perform actual DCT-II
  const double pi = 3.1415926535897932;
  for (int i=0; i<n; ++i) {
    ausgabe[i] = 0.0;
    for (int j=0; j<n; ++j) {
      ausgabe[i] = ausgabe[i] + cos (pi / n * i * (j + 0.5)) * eingabe[j];

    }
  }
  return;
}

Danke im Vorraus!
Marian

[Marian G. (houwdy)]

keiner eine Idee?