Discussion :

[najbet]

salut
je desire realiser un spectogramme : en fait je veux réaliser un graphique à 3 dimensions : 2 dimensions pour le spectre (Intensité et Fréquences) et une dimension pour le temps, sous matlab mais je n'y arrive pas

je dispose d'un signal que j'ai acquis d'un length de 1014300

[sango85]

La fonction spectrogram représente en X le temps, en Y la fréquence et en couleur (que l'on peut considérer comme une 3ème dimension) l'amplitude du signal... Sinon, pour de la vraie 3D, je sèche...

[nadmont]

salut
je desire realiser un spectogramme : en fait je veux réaliser un graphique à 3 dimensions : 2 dimensions pour le spectre (Intensité et Fréquences) et une dimension pour le temps, sous matlab mais je n'y arrive pas

je dispose d'un signal que j'ai acquis d'un length de 1014300

si tu as 3 séries de données F t et coeff (intensité)
tu réalises une mech(t,f,coeff)
ensuite pour avoir une vue 3D ("en volume") au lieu de couleurs pour les coeff
tu définis une vue avec 2 arguments [azimuth elevation]
exemple view([-90 80])

[najbet]

En fait c'est ce que je desire obtenir :

La fonction spectrogram représente en X le temps, en Y la fréquence et en couleur (que l'on peut considérer comme une 3ème dimension) l'amplitude du signal...

mais je ne sais pas trop coment utiliser cette fonction pour visualiser ces 3 dimention sachant que je ne dispose que de mon signal temporel

[shinobida]

En fait c'est ce que je desire obtenir :

mais je ne sais pas trop coment utiliser cette fonction pour visualiser ces 3 dimention sachant que je ne dispose que de mon signal temporel

J'ai pas la formule en tête, mais si tu veux je pourrai te la donner aprés, l'idée c'est de calculer ta trasfomré en sepctro à l'aide de ton signal (je crois que la formule c'est : de prndere ton signat x(u). h*(u-t).Exp(-j2(pi)f.u) et de l'intégré suivant ''u'' t'aura une fonction TEMPS(t)Fréquence(f) et tu pourra donc faire la représentation par par exemple mesh(...,...,...)

PS: h c'est une focntion fenetre centré en 0, et selon la largeur de cette fenetre tu aura soit une bonne résolution temporelle ou une bonne résolution fréquentielle, et le pb du spectro et du sonogram(son eq en fréquence) c'est qu'i faut faire un compromi entre ces deux résolutions

Good luck, et si t'a besoin de plus d'info , je suis la

NB: je parles du vrai spectrogramme utilisé en analyse temps-fréquentille, celui de Matlab je l'ai pas utilisé et je ne pourrai pas te dire s'il refléte le SPECTROGRAM de la TFSAP ( Time-Frequency signal analysis and processing )

[shinobida]

pour son equivalent en fréquence , le SONOGRAM voila la fonction :
intégration de ( X(v). H*(v-f).exp(-2j(pi)(v-f)t) . dv) = fonction(t,f)

good luck


PS: X est la TF de x
et H est la TF de h

[nadmont]

pour faire une représentation temps-fréquence tu as plusieurs solutions
cela dépend des outils matlab.

le plus simple, tu peux faire une analyse temps-fréquence au moyen d'une fourier-court-terme (STFT) au moyen de la "spectram" ou suivant la version "spectrogram" :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 [coeff, f, t]=specgram(signal, nfft, Fe, fenêtre, Noverlap);

** nfft, le nombre de points pour réaliser la transformation exple 1024 ;

** Fe, la fréquence d'échantillonnage ;

** fenêtre = fenêtre d'amortissement (hanning, hamming, blackman etc..), tu dois également fixer le nombre de points pour ce fenêtrage :
hannin(Nwindow) exple 256 ;

** Noverlap, le nombre de points de décalage pour ta fenêtre d'analyse sur ton signal analysé (chevauchement des fenêtres) exple Nwindow-1.

tout cela te renvoie l'intensité (coeff), la fréquence (f) et le temps (t).

Maintenant tu peux utiliser d'autres méthode au moyen de transformés en Ondelette continue (tu as de nombreuses possibilités en ondelette)

tu as d'autres méthodes en temps-fréquence conjointe, tu peux aller voir la Toolbox Time-Frequency développée par Flandrin et coll (CNRS, Rice University ) dans le cadre isis.

[ammoula1]

Bonjour,
Désolée pour le dérangement mais je veux savoir lorsque je fais

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

specgram(signal, nfft, Fe, fenêtre, Noverlap);

il suffit de passer le nom du fenêtre tout court,je veux dire par tout court,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

specgram(signal, nfft, Fe, 'hamming', Noverlap);

ou non ,
Merci pour votre aide ,