Viadeo Twitter Google Bookmarks ! Facebook Digg del.icio.us MySpace Yahoo MyWeb Blinklist Netvouz Reddit Simpy StumbleUpon Bookmarks Windows Live Favorites 
Logo Documentation Qt ·  Page d'accueil  ·  Toutes les classes  ·  Toutes les fonctions  ·  Vues d'ensemble  · 

fastfouriertransformer.cpp Example File

demos/mobile/guitartuner/src/fastfouriertransformer.cpp
 /****************************************************************************
 **
 ** Copyright (C) 2011 Nokia Corporation and/or its subsidiary(-ies).
 ** All rights reserved.
 ** Contact: Nokia Corporation (qt-info@nokia.com)
 **
 ** This file is part of the QtDeclarative module of the Qt Toolkit.
 **
 ** $QT_BEGIN_LICENSE:BSD$
 ** You may use this file under the terms of the BSD license as follows:
 **
 ** "Redistribution and use in source and binary forms, with or without
 ** modification, are permitted provided that the following conditions are
 ** met:
 **   * Redistributions of source code must retain the above copyright
 **     notice, this list of conditions and the following disclaimer.
 **   * Redistributions in binary form must reproduce the above copyright
 **     notice, this list of conditions and the following disclaimer in
 **     the documentation and/or other materials provided with the
 **     distribution.
 **   * Neither the name of Nokia Corporation and its Subsidiary(-ies) nor
 **     the names of its contributors may be used to endorse or promote
 **     products derived from this software without specific prior written
 **     permission.
 **
 ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
 ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
 ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
 ** A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
 ** OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
 ** SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
 ** LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
 ** DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
 ** THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
 ** (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
 ** OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE."
 ** $QT_END_LICENSE$
 **
 ****************************************************************************/

 #include "fastfouriertransformer.h"
 #include "math.h"

 #define STIN  inline
 #define __STATIC

 #include "fftpack.c"

 // called by __ogg_fdrffti
 __STATIC void drfti1(int n, float *wa, int *ifac);
 void __ogg_fdrffti(int n, float *wsave, int *ifac);
 void __ogg_fdcosqi(int n, float *wsave, int *ifac);
 // called by drftf1
 STIN void dradf2(int ido,int l1,float *cc,float *ch,float *wa1);
 // called by drftf1
 STIN void dradf4(int ido,int l1,float *cc,float *ch,float *wa1,
         float *wa2,float *wa3);
 // called by drftf1
 STIN void dradfg(int ido,int ip,int l1,int idl1,float *cc,float *c1,
         float *c2,float *ch,float *ch2,float *wa);
 // called by drftf1
 STIN void drftf1(int n,float *c,float *ch,float *wa,int *ifac);
 void __ogg_fdrfftf(int n,float *r,float *wsave,int *ifac);
 STIN void dcsqf1(int n,float *x,float *w,float *xh,int *ifac);
 void __ogg_fdcosqf(int n,float *x,float *wsave,int *ifac);
 STIN void dradb2(int ido,int l1,float *cc,float *ch,float *wa1);
 STIN void dradb3(int ido,int l1,float *cc,float *ch,float *wa1,
                           float *wa2);
 STIN void dradb4(int ido,int l1,float *cc,float *ch,float *wa1,
         float *wa2,float *wa3);
 STIN void dradbg(int ido,int ip,int l1,int idl1,float *cc,float *c1,
             float *c2,float *ch,float *ch2,float *wa);
 STIN void drftb1(int n, float *c, float *ch, float *wa, int *ifac);
 void __ogg_fdrfftb(int n, float *r, float *wsave, int *ifac);
 STIN void dcsqb1(int n,float *x,float *w,float *xh,int *ifac);
 void __ogg_fdcosqb(int n,float *x,float *wsave,int *ifac);

 FastFourierTransformer::FastFourierTransformer(QObject *parent) :
         QObject(parent),
         m_waveFloat(0),
         m_workingArray(0),
         m_ifac(0),
         m_last_n(-1)
 {
 }

 FastFourierTransformer::~FastFourierTransformer()
 {
     if (m_waveFloat != 0) {
         delete [] m_waveFloat;
     }
     if (m_workingArray != 0) {
        delete [] m_workingArray;
     }
     if (m_ifac != 0) {
         delete [] m_ifac;
     }
 }

 /**
   * Prepares the arrays to be of length n.
   */
 void FastFourierTransformer::reserve(int n)
 {
     Q_ASSERT(n>0);
     if (m_waveFloat != 0) {
         delete [] m_waveFloat;
     }
     if (m_workingArray != 0) {
        delete [] m_workingArray;
     }
     if (m_ifac != 0) {
         delete [] m_ifac;
     }
     m_workingArray = new float[2*n+15];
     m_waveFloat = new float[n];
     m_ifac = new int[n];
     __ogg_fdrffti(n, m_workingArray, m_ifac);
     m_last_n = n;
 }

 /**
   * Calculates the Fast Fourier Transformation (FFT).
   */
 void FastFourierTransformer::calculateFFT(QList<qint16> wave)
 {
     const int n = wave.size();
     if (m_last_n != n) {
         reserve(n);
     }
     for (int i = 0; i < n; i++) {
         m_waveFloat[i] = (float) wave.at(i);
     }

     __ogg_fdrfftf(n, m_waveFloat, m_workingArray, m_ifac);

 }

 /**
   * Returns the index which corresponds to the maximum density
   * of the FFT.
   */
 int FastFourierTransformer::getMaximumDensityIndex()
 {
     const int halfN = m_last_n / 2;
     float maxDensity = 0;
     int maxDensityIndex = 0;
     float densitySquared = 0.f;
     for (int k = 1; k < halfN; k++) {
         // Here, we calculate the frequency k/N.
         // k=1, the wave oscillation time is N, and the frequency
         //      is 1/sample.
         // k=2, the wave oscillation time is N/2, and the frequency
         //      is 2/sample.
         // k=3, the wave oscillation time is N/3, and the frequency
         //      is 3/sample.
         // Note, that the documentation is for Fortran, so indexes in the
         // documentation does not match.
         // The sine and cosine coefficients are obtained thus as follows:
         const float cosCoefficient = qAbs(m_waveFloat[2*k-1]);
         const float sinCoefficient = qAbs(m_waveFloat[2*k]);

         densitySquared = sinCoefficient*sinCoefficient + cosCoefficient*cosCoefficient;
         if (densitySquared > maxDensity) {
             maxDensity = densitySquared;
             maxDensityIndex = k;
         }
     }

     if (m_cutOffForDensitySquared < maxDensity) {
         return maxDensityIndex;
     }
     else {
         return -1;
     }
 }

 /**
   * Sets the cutoff density.
   */
 void FastFourierTransformer::setCutOffForDensity(float cutoff)
 {
     m_cutOffForDensitySquared = cutoff*cutoff;
 }
Publicité

Best Of

Actualités les plus lues

Semaine
Mois
Année
  1. « Quelque chose ne va vraiment pas avec les développeurs "modernes" », un développeur à "l'ancienne" critique la multiplication des bibliothèques 94
  2. Apercevoir la troisième dimension ou l'utilisation multithreadée d'OpenGL dans Qt, un article des Qt Quarterly traduit par Guillaume Belz 0
  3. Les développeurs ignorent-ils trop les failles découvertes dans leur code ? Prenez-vous en compte les remarques des autres ? 17
  4. Pourquoi les programmeurs sont-ils moins payés que les gestionnaires de programmes ? Manquent-ils de pouvoir de négociation ? 42
  5. Quelles nouveautés de C++11 Visual C++ doit-il rapidement intégrer ? Donnez-nous votre avis 10
  6. Adieu qmake, bienvenue qbs : Qt Building Suite, un outil déclaratif et extensible pour la compilation de projets Qt 17
  7. 2017 : un quinquennat pour une nouvelle version du C++ ? Possible, selon Herb Sutter 7
Page suivante

Le blog Digia au hasard

Logo

Créer des applications avec un style Metro avec Qt, exemples en QML et C++, un article de Digia Qt traduit par Thibaut Cuvelier

Le blog Digia est l'endroit privilégié pour la communication sur l'édition commerciale de Qt, où des réponses publiques sont apportées aux questions les plus posées au support. Lire l'article.

Communauté

Ressources

Liens utiles

Contact

  • Vous souhaitez rejoindre la rédaction ou proposer un tutoriel, une traduction, une question... ? Postez dans le forum Contribuez ou contactez-nous par MP ou par email (voir en bas de page).

Qt dans le magazine

Cette page est une traduction d'une page de la documentation de Qt, écrite par Nokia Corporation and/or its subsidiary(-ies). Les éventuels problèmes résultant d'une mauvaise traduction ne sont pas imputables à Nokia. Qt 4.7-snapshot
Copyright © 2012 Developpez LLC. Tous droits réservés Developpez LLC. Aucune reproduction, même partielle, ne peut être faite de ce site et de l'ensemble de son contenu : textes, documents et images sans l'autorisation expresse de Developpez LLC. Sinon, vous encourez selon la loi jusqu'à 3 ans de prison et jusqu'à 300 000 E de dommages et intérêts. Cette page est déposée à la SACD.
Vous avez déniché une erreur ? Un bug ? Une redirection cassée ? Ou tout autre problème, quel qu'il soit ? Ou bien vous désirez participer à ce projet de traduction ? N'hésitez pas à nous contacter ou par MP !
 
 
 
 
Partenaires

Hébergement Web