Changeset 1962 in Sophya for trunk/SophyaLib/NTools

Timestamp:

Apr 11, 2002, 6:07:44 PM (23 years ago)

Author:

lemeur

Message:

changement fft reelle en complexe, pour nombre pair d'achantillons

File:

• trunk/SophyaLib/NTools/FSAppIrrSmpl.cc (modified) (2 diffs)

trunk/SophyaLib/NTools/FSAppIrrSmpl.cc

@@ -119 +119 @@
 void FSApproximationIrregularSampling::restaureRegularlySampledSignal(int nbInterv, TVector<double>& solution) const
 {
-  if (nbInterv < 2*M_+1 || nbInterv%2 !=  0) 
-      {
-        int k; 
+  cout << " c'est celui que je corrige " << endl;
+  int k; 
+  if (nbInterv < 2*M_+1) 
+      {
         solution.ReSize(nbInterv); 
         double delta = 1./nbInterv;
@@ -132 +133 @@
     else
       {
-        int tailleTF = nbInterv/2+1;
-        TVector<complex<double> > TFf(tailleTF);
+        TVector<complex<double> > TFf(nbInterv);
+        TVector<complex<double> > bidon(nbInterv);
         TFf = complex<double>(0.,0.);
         TFf(Range(0,M_)) = coefFourier_;
-        fftIntfPtr_-> FFTBackward(TFf, solution);
+        for ( k=1; k<= M_; k++)
+          { 
+            TFf(nbInterv-k) = conj(coefFourier_(k));
+          }
+        fftIntfPtr_-> FFTBackward(TFf, bidon);
+        cout << " taille de bidon "<< bidon.Size() << endl;
+        solution.ReSize(nbInterv);
+        for (k=0; k< nbInterv; k++) solution(k) = bidon(k).real();
       }
   reshapeSignalInUsersFrame(solution);