(complexe) Cela dépend de la façon de nommer le coefficient de l'ondelette de morlet, de la valeur de réglage appropriée et du matériau, j'ai donc essayé de l'organiser autant que possible.

Je ne peux pas vraiment l'organiser, mais ...

1. Je ne sais pas si les variables et expressions de l'expression d'ondelettes sont légèrement différentes et pointent exactement vers la même chose.

ref. D'après le matériel de matlab (probablement) Teolis (1998) Traitement du signal informatique avec ondelettes. P66

formule

morlet formula orig

Nom de variable

γb: variance nommée (largeur de bande). Identique à fb dans le matériau matlab. Cela semble être une dispersion temporelle. Mais qu'est-ce que la dispersion temporelle ici? Elle s'appelle γc: fréquence centrale. C'est presque la même. Je comprends que c'est la valeur fixée par factor = np.arange (8,30) comme mne. j: unité imaginaire

ref. stack overflow

formule

w (t, f 0 </ sub>) = A ・ exp (−t 2 </ sup> / 2 * σ 2 </ sup>) ・ exp (2 * i * π * f 0 </ sub> * t) = A ・ exp (2 iπ f 0 </ sub> t- (t 2 </ sup> / 2σ 2 </ sup> >))

variable

f 0 </ sub>: fréquence centrale (γc ci-dessus, fc dans la boîte à outils d'ondelettes matlab) Il est décrit comme σ: σ = m / 2πf 0 </ sub>. Comme décrit dans [Tallon-Baudry (1999)] 3, la signification de l'écart type (sigmat) par rapport au temps. Par rapport à l'équation de Teolis (1998) ci-dessus, Cela correspond en supposant que γb = 2σ 2 </ sup>. m: constant. Compromis entre la résolution temporelle et fréquentielle. Je l'ai mis à 7. A = 1 / (σ (2π) 0,5 </ sup>) = 1 / (σ√2π): La partie avant e (exp) est remplacée par une variable. Introduite sous le nom de facteur de normalisation Oui.

réf. [Tallon-Baudry (1999)] 3, p154 Encadré 2 Description

formule

w (t, f) = A. exp (- t 2 </ sup> / 2σ t </ sub> 2 </ sup>). exp (2i π ft): C'est Identique au débordement de pile.

variable

f: Identique à f 0 </ sub> ci-dessus. Fréquence centrale σ t </ sub>: Identique à σ ci-dessus. Décrit comme l'écart type du domaine temporel.

** A = 1 / (σt √π) 1/2 **: Il est écrit comme un facteur de normalisation comme [stack overflow] 2. La formule est un peu différente. [Tallon-Baudry (1997)] 11 Il a la même description, on ne sait pas si c'est intentionnel ou faux + copie.

ref. matlab wavelet toolbox

formule

ψ (x) = 1 / √πfb * ** exp (x 2 </ sup> / fb) ** * exp (2iπfcx): le copier-coller n'a pas fonctionné, alors copiez-le. Comparé aux deux ci-dessus, ** exp (x 2 </ sup> / fb) ** n'a pas de ** moins **. Cette page est la seule, donc j'ai l'impression d'avoir fait une erreur. .. D'autres correspondent à Teolis (1998), stack overflow

variable

x: identique à t ci-dessus fc: Identique à f0, f, γc ci-dessus. Fréquence centrale fb: paramètre de bande passante nommé, identique à γb dans Teolis (1998).

ref. matlab help

formule

wavelet formula C'est aussi un peu différent. Le dénominateur et la molécule sont inversés dans la partie avant e (A: facteur de normalisation qui apparaît ci-dessus). La différence entre le dénominateur et la molécule n'est que cette page, donc c'est probablement faux. D'autres correspondent à Teolis (1998), stack overflow

variable

Le nom de fc et fb est le même que [matlab wavelet toolbox] 4

réf. [Échange de pile] 7 wavelet cmor fb fc ou quelque chose comme ça va frapper le dessus.

formule

wavelet formula another Le moins de e -j2pifct </ sup> semble inutile.

Résumé de 1

Jugement complet (complexe) La formule d'ondes de Morlet est Teolis (1998) Computational Signal Proceesing with Wavelets. P66, Stack overflow semble correct. Teolis (1998) est un texte. La formule est vraiment légèrement différente selon le matériau, et j'ai l'impression qu'elle est en fait pleine d'erreurs.

2. Quel est un bon paramètre pour les ondelettes. (Principalement avec mne python à l'esprit)

[sortie sur le terrain] 9(analyse temps-fréquence, section d'ondelettes)

cfg.width. Il détermine la ** largeur ** des ondelettes en nombre de cycles. Dans l'échantillon

cfg.width      = 7; 

C'était. Aussi, The spectral bandwidth at a given frequency F is equal to F/width*2 the wavelet duration is equal to width/F/pi F est la fréquence, qui semble être la même que la fréquence centrale (f, f 0 </ sub>, γc, fc) ci-dessus.

--Je ne sais tout simplement pas ce que signifient la bande passante spectrale et la durée des ondelettes. --Il semble que ce n'est pas une relation si simple que la bande passante spectrale corresponde à σ f </ sub> et que la durée de la forme d'onde corresponde à σ t </ sub>. Il existe de nombreux termes et c'est déroutant.

Tallon-Baudry(1999) ** σ f </ sub> = 1 / (2πσ t </ sub>) ** (Équation 2.1) La relation est montrée.

--σ t </ sub> est décrit comme un écart type (domaine temporel) .Il semble être le même que σ t </ sub> et σ qui sont apparus plusieurs fois. --σ f </ sub> est également décrit comme un écart type (domaine fréquentiel).

** f / σ f </ sub> (rapport constant: supérieur à cinq en pratique) ** (équation 2.2).

--f semble être la même que la fréquence centrale (F dans le champ trip, f 0 </ sub>, γc, fc qui apparaît ci-dessus).

Roach & Mathalon (2008) ** σ t </ sub> = 1 / (2πσ f </ sub>) ** (Équation 2.3) Il en est de même si l'équation (2.1) de [Tallon-Baudry (1999)] 3 est modifiée.

Description des variables

σ t </ sub>: Décrit comme SD temporel. σ f </ sub>: bande passante spectrale autour de sa fréquence centrale, f0, qui a un SD de σ f </ sub>, on peut donc dire que c'est SD de fréquence. L'ondelette de Morlet est l'enveloppe gaussienne ( Etant donné que le mot «enveloppe normalement distribuée» apparaît fréquemment, SD peut être facile à utiliser comme un index sous forme d'ondelette. Juste au bas de la même page spectral bandwidth around any given center frequency is mσf. Il y a aussi. m est multiplié. Pour en déduire spectral bandwith = mσf , and has σf in it. J'imagine que ça veut dire, je suis troublé par la description qui confond le texte.

** c = f 0 </ sub> / σ f </ sub> ** (Équation 2.4) Et ** c (rapport constant) ** sont définis ** (probablement la même chose que cfg.width) **. Cela indique également la même chose, seule la variable c est affectée à l'équation 2.2. On pense que oui.

--Il est recommandé dans le texte que c soit un nombre d'environ 5,6,7, de sorte que σ t </ sub> contienne des ondes avec un ou plusieurs cycles ("constan"). , c, est généralement recommandé pour être supérieur à 5, et est souvent mis à des valeurs de 6 ou 7, ce qui correspond à un σ t </ sub> qui, par conséquent, au moins un cycle sinusoïdal complet pour une fréquence particulière ") .- -Dans [stack overflow] 2, la formule ** σ = m / (2 π f 0 </ sub>) ** apparaît, et (σ = σ t </ sub> (Peut-être), qui correspond à l'équation (2.5) ci-dessous, il semble donc que m et le rapport constant c dans cette équation pointent vers la même chose.

Un peu à part: Explication que σ t </ sub> contient des ondes avec un ou plusieurs cycles lorsque c est d'environ 7 (peut être faux)


Lorsqu'il est transformé sous la forme de σ t </ sub> = en utilisant les équations (2.3) et (2.4)

** σ t </ sub> = c / (2πf) = c / (2π) * 1 / f ** (équation 2.5) est obtenu.

** Variable inconnue m **

[Roach & Mathalon (2008)] 10 l'explique comme un facteur m. ** fenêtre temporaire pour toute fréquence donnée (f 0 </ sub>) est m σ t </ sub> ** (équation 2.6) Je ne sais pas à quoi se réfère la fenêtre temporelle. Cependant, j'estime que la fenêtre temporelle est la même que la durée de l'ondelette qui apparaît dans [field trip] 9 Transformons l'expression.

À partir des équations (2.6) et (2.5) fenêtre temporelle = mc / (2π) * 1 / f = mc / (2π) / f. Je répète la description dans la section sortie sur le terrain ci-dessus,

** durée des ondelettes ** est égale à ** largeur / F / pi ** (équation 2.7) Puisque F = f, en supposant c = largeur et m = 2, fenêtre temporelle = durée de l'ondelette Tient. ** fenêtre temporaire ** (par Roach & Mathalon (2008)) ** = durée des ondelettes ** (par sortie sur le terrain) semble bonne.

[Roach et Mathalon (2008)] 10 spectral bandwidth around any given center frequency is mσf Est-ce la même chose que la bande passante spectrale dans [field trip] 9? À partir de Fieldtrip: bande passante spectrale = F / largeur * 2

Si la largeur précédente = c est correcte, alors c = f 0 </ sub> / σ f </ sub> (équation 2.4) Parce qu'il y a.

Roach & Mathalon (2008) bande passante spectrale = mσ f </ sub> = mf 0 </ sub> / c = F / largeur * m (équation de parcours) Ainsi, on peut dire que m = 2 indique la même chose.

Vous pouvez également dire que fieldtrip cfg.width = c (rapport constant).

Il y a trop de termes, c'est déroutant. Au fait, il n'y a pas de description explicite correspondant à ce "m" lors de la visite sur le terrain. Est-il traité de manière appropriée en interne?

MNE python

Référence: si vous jetez un œil au script morlet que tfr_morlet utilise en interne

Il y a un morlet def dans mne / time_frequency / tfr.py dans le script

tfr.py


sigma_t = this_n_cycles / (2.0 * np.pi * f)

Par conséquent, à partir de l'équation (2.5), je pense que this_n_cycles (= n_cycles) = c (rapport constant).

Résumé de 2

--cfg.width of fieldtrip et n_cycles of mne python, c (constant ratio. Equation 2.4) in the document, m in [stack overflow] 2 semble indiquer la même chose. Comprendre la signification des paramètres " Je ne peux pas, mais j'aimerais lire le rapport constant et essayer de le comprendre autant que possible.Notez que m dans [Roach & Mathalon (2008)] 10 est différent. --σ t </ sub> est temporaire SD, qui est différente de la durée de l'ondelette (probablement = fenêtre temporaire) .Wavelet duration = mσ t </ sub> ([Roach & Mathalon (2008)] [ À partir de 10]

Recommended Posts

(complexe) Cela dépend de la façon de nommer le coefficient de l'ondelette de morlet, de la valeur de réglage appropriée et du matériau, j'ai donc essayé de l'organiser autant que possible.
L'arbre.plot_tree de scikit-learn était très simple et pratique, j'ai donc essayé de résumer comment l'utiliser facilement.
J'ai essayé de faciliter la modification du paramètre du proxy authentifié sur Jupyter
J'ai essayé de récupérer les données de l'ordinateur portable en le démarrant sur Ubuntu
Je n'ai pas compris le redimensionnement de TensorFlow, alors je l'ai résumé visuellement.
J'ai essayé de comprendre comment utiliser les pandas et la colinéarité multiple en utilisant l'ensemble de données Affaires comme thème.
[Rails] v1.0 est sorti sur google-cloud-vision de gem, j'ai donc essayé de le soutenir
Il fait froid, j'ai donc essayé de permettre d'allumer / d'éteindre automatiquement le chauffage AC avec Raspberry Pi!
[Linux] [C / C ++] Comment obtenir la valeur d'adresse de retour d'une fonction et le nom de fonction de l'appelant
Comment exécuter Python sous Windows sans polluer autant que possible l'environnement (édition Scoop)
Je ne suis pas sûr de la différence entre les modules, les packages et les bibliothèques, alors j'ai essayé de les organiser.
J'ai résumé comment changer les paramètres de démarrage de GRUB et GRUB2
J'ai essayé de comprendre la normalisation spectrale et la décomposition de valeurs singulières qui contribuent à la stabilité du GAN.
J'ai essayé de pousser les documents Sphinx vers BitBucket et de les refléter automatiquement sur le serveur Web