নিষ্পাপ পদ্ধতি ব্যবহার করে বিচ্ছিন্ন ফুরিয়ার ট্রান্সফর্ম কম্পিউট করার জন্য C++ প্রোগ্রাম

বিচ্ছিন্ন ফুরিয়ার ট্রান্সফর্মে (DFT), একটি সসীম তালিকা একটি ফাংশনের সমান ব্যবধানের নমুনাগুলিকে জটিল সাইনোসয়েডগুলির একটি সীমাবদ্ধ সংমিশ্রণের সহগগুলির তালিকায় রূপান্তরিত করা হয়। তারা তাদের ফ্রিকোয়েন্সি অনুসারে অর্ডার করে, যার একই নমুনা মান রয়েছে, নমুনাকৃত ফাংশনকে তার আসল ডোমেন (প্রায়শই সময় বা একটি লাইন বরাবর অবস্থান) থেকে ফ্রিকোয়েন্সি ডোমেনে রূপান্তর করতে।

অ্যালগরিদম

Begin
   Take a variable M and initialize it to some integer
   Declare an array function[M]
   For i = 0 to M-1 do
      function[i] = (((a * (double) i) + (b * (double) i)) - c)
   Done
   Declare function sine[M]
   Declare function cosine[M]
   for i =0 to M-1 do
      cosine[i] = cos((2 * i * k * PI) / M)
      sine[i] = sin((2 * i * k * PI) / M)
   Done
   Declare DFT_Coeff dft_value[k]
   for j = 0 to k-1 do
      for i = 0 to M-1 do
         dft_value.real += function[i] * cosine[i]
         dft_value.img += function[i] * sine[i]
      Done
   Done
   Print the value
End

উদাহরণ কোড

#include<iostream>
#include<math.h>
using namespace std;
#define PI 3.14159265
class DFT_Coeff {
   public:
   double real, img;
   DFT_Coeff() {
      real = 0.0;
      img = 0.0;
   }
};
int main(int argc, char **argv) {
   int M= 10;
   cout << "Enter the coefficient of simple linear function:\n";
   cout << "ax + by = c\n";
   double a, b, c;
   cin >> a >> b >> c;
   double function[M];
   for (int i = 0; i < M; i++) {
      function[i] = (((a * (double) i) + (b * (double) i)) - c);
      //System.out.print( " "+function[i] + " ");
   }
   cout << "Enter the max K value: ";
   int k;
   cin >> k;
   double cosine[M];
   double sine[M];
   for (int i = 0; i < M; i++) {
      cosine[i] = cos((2 * i * k * PI) / M);
      sine[i] = sin((2 * i * k * PI) / M);
   }
   DFT_Coeff dft_value[k];
   cout << "The coefficients are: ";
   for (int j = 0; j < k; j++) {
      for (int i = 0; i < M; i++) {
         dft_value[j].real += function[i] * cosine[i];
         dft_value[j].img += function[i] * sine[i];
      }
      cout << "(" << dft_value[j].real << ") - " << "(" << dft_value[j].img <<" i)\n";
   }
}

আউটপুট

Enter the coefficient of simple linear function:
ax + by = c
4 5 6
Enter the max K value: 10
The coefficients are:
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)
(345) - (-1.64772e-05 i)