Advertisement
Guest User

Untitled

a guest
Nov 23rd, 2017
68
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.45 KB | None | 0 0
  1. class DFT : public sinusoida {
  2. double sumRe, sumIm;
  3. double *y, *sumMod;
  4.  
  5. public:
  6. DFT() :sinusoida() {
  7. sumMod = new double[N];
  8.  
  9. }
  10.  
  11. void generateDFT() {
  12. for (int i = 0; i < N; i++) {
  13. for (int j = 0; j < N; j++) {
  14. sumRe += x[i] * cos(2 * M_PI * j * i);
  15. sumIm += y[i] * sin(2 * M_PI * j * i);
  16. }
  17.  
  18. sumMod[i] +=(sqrt (pow(sumRe, 2) + pow(sumIm, 2)));
  19. time[i] = i * (fs / N);
  20. sumRe = 0;
  21. sumIm = 0;
  22. }
  23.  
  24.  
  25. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement