Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- class DFT : public sinusoida {
- double sumRe, sumIm;
- double *y, *sumMod;
- public:
- DFT() :sinusoida() {
- sumMod = new double[N];
- }
- void generateDFT() {
- for (int i = 0; i < N; i++) {
- for (int j = 0; j < N; j++) {
- sumRe += x[i] * cos(2 * M_PI * j * i);
- sumIm += y[i] * sin(2 * M_PI * j * i);
- }
- sumMod[i] +=(sqrt (pow(sumRe, 2) + pow(sumIm, 2)));
- time[i] = i * (fs / N);
- sumRe = 0;
- sumIm = 0;
- }
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement