L'intégration numérique à l'aide de la Règle de Simpson sur des données discrètes

Je suis à la recherche pour l'intégration numérique avec matlab. Je sais qu'il y a un trapz fonction dans matlab, mais la précision n'est pas assez bon. Par la recherche en ligne, j'ai trouvé il y a un quad fonction de là, il semble accepter l'expression symbolique comme entrée. Mes données sont discrètes et unidimensionnel. Est-ce une façon d'utiliser quad sur mes données? Merci.

  • quad() est une mise en œuvre de la règle de Simpson, qui est à l'école secondaire en mathématiques. Est-il rien de vous empêche de mise en œuvre de la règle de Simpson-vous?
  • s'agit-il seulement d'utiliser les simpson est la règle, je pense que c'est adaptative Simpson en quadrature. Mais de toute façon, je ne sais pas pourquoi, le quad est plus rapide que ma mise en œuvre. J'ai beaucoup de données à intégrer et à la recherche d'un moyen plus rapide.
  • Il a été quelques années depuis que mes méthodes numériques sûr, mais autant que je me souvienne, la seule différence entre la règle de Simpson et l'adaptation du genre, c'est que Adaptative de la règle de Simpson s'applique à une variable d'espace entre les points d'échantillonnage. Puisque vous n'êtes pas l'intégration d'une expression symbolique vous ne pouvez pas varier l'intervalle d'échantillonnage - vous avez déjà vos données et vous ne pouvez pas interpoler des points de plus entre elle. Donc régulier de la règle de Simpson sera aussi bon que vous pouvez obtenir.
  • Concernant la vitesse: Avez-vous envisagé d'utiliser l'exécution en parallèle des outils disponibles dans les nouvelles versions de MATLAB? La boucle for parallèle, parfor, est mort, facile à utiliser et diffuser votre charge CPU le plus grand nombre de cœurs (ou nœuds d'exécution) que possible.
InformationsquelleAutor user1285419 | 2012-03-26