Discussion :

[jujuf1]

Bonjour à tous,

J'ai cherché sans résultats....Je cherche à calculer une densité spectrale de puissance en "unité"^2/Hz, où "unité" est l'unité du signal au lieu de dB/Hz qui semble être la méthode de calcul par défaut. Par exemple, si le signal est une accélarétion en fonction du temps ce qui est mon cas, il me faudrait avoir les valeurs en "accélération"^2/Hz et non en dB (d'ailleurs matlab prend quelle référence poour la conversion en dB ????).

Autre question: je me demande d'ailleurs comment récupérer les données concernant les abscisses et les ordonnées afin de pouvoir exporter tout ça vers un autre logiciel.


Je vous remercie par avance pour votre aide.

[aliasfdr]

matlab ne calcule que ce que tu lui demande de calculer ? donc pas de conversion sans demande de conversion !!
Tu utilise une fonction déjà faite ? peut tu nous la montrer ?
ou une fonction que t'a programmée ?

[jujuf1]

C'est une fonction toute faite. J'ai potassé l'aide et ils précisent que la fonction retourne la psd en dB/Hz.

Voici ce que j'utilise:

Soit le signal temporel AA, supposé gaussien, échantillonné à la fréquence 6000 Hz

hs=spectrum.welch('rectangular',300,50);
psd(hs,AA,'6000',fs,'NFTT',NFFT);

et ça me sort une densité specctrale en dB/Hz..., comme indiqué dans l'aide.
Il n'est nulle part fait mention de d'autres possibilités. Je ne demande pourtant pas quelquechose d'exotique.....Je ne comprend pas pourquoi Matlab n'offre pas la possibilité de garder la grandeur d'entrée...C'est pourtant courant d'avoir des DSP en (m/s²)²/Hz dans le domaine des vibrations...Au moins s'ils expliquaient en quoi consiste ces decibels.....10 log (puissance/Pref) ??, 20 log((signal dans le domaine fréquentiel)²/(reférence²)) et indiquaient la référence utilisée, je pourrais dans ce cas convertir moi même les valeurs, si toutefois on peut extraire des données d'un objet dspdata....ce que je n'ai pas réussi à faire....

Merci d'avance

[aliasfdr]

je n'ai pas trouvé de fonction psd dans ma version matlab ?
peut être que je n'ai pas la toolbox associée !
je te conseille de jeter un coup d'oeil à la fonction ou la poster pour trouver la fonction de conversion, à mon avis c'est du 20*log(grandeur)

[jujuf1]

je pense que tu as la toolbox:c'est signal processing. Cette fonction psd n'est plus détaillée, sauf en tapant help dsp: ils disent dailleurs qu'elle est obsolète. Elle est pourtant utilisée puisque j'ai utilisé ce script à partir de l'aide. Ils expliquent comment tracer des dsp dans spectrum, periodogramm etc....mais le problème, c'est que j'ai déja regardé et j'ai rien trouvé....pourtant cette possibilité doit forcément exister....ou c'est à n'y rien y comprendre...

[aliasfdr]

je n'ai pas la signal processing toolbox dans mon bureau, mais bon essaie de taper open psd dans la command window y'aura plus de details sur ta fonction peut être ..

[jujuf1]

Bah c'est super tu m'a dépanné !

Etant novice en matlab, j'avais pas eu l'idée d'ouvrir le fichier m de la fonction PSD. J'ai pu voir dedans l'opération que la fonction réalisait pour aboutir à ces fameux dB. Je vais donc pouvoir me faire une fonction perso !

Encore merci !

[aliasfdr]

content de t'avoir aidé

a+

[jujuf1]

Re bonjour.

ça fait plus d'une heure que je m'acharne à modifier le fichier psd.m en retirant le 10*log10() sans succès !!! quand je relance psd(......), ça me rebalance la meme figure (mêmes valeurs, etc....).

C'est normal ? Les fichiers appartenant à matlab sont-ils protégés ?????

Merci d'avance.

[aliasfdr]

essaie d'ouvrir un nouveau m file (copier coller )et ensuite de changer le nom de la fonction et les lignes qui t'interessent.

[Caro-Line]

Il vaut mieux en effet ne pas modifier les fichiers d'origine de MATLAB et plutôt faire un Save As et sauver sous un nom à toi (my_psd par exemple) dans un de tes répertoires de travail.
Sinon quand tu as plusieurs fichiers qui portent le même nom dans des répertoires différents tu fais :

Code :

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

which -all

Le premier qui apparait est celui qui est appelé par défaut.

[jujuf1]

merci caro pour la commande which all.

cela renvoie à deux fichier psd.m: l'un que j'ai modifié, et un autre, qui semble s'apparenter à un fichier aide.

j'avais décidé de modifier directement le fichier psd (en faisant un backup) parce que lorsque que je fais save as un autre nom, la fonction ne veut pas fonctionner, en prenant soin de modifier l'entête "...=nouveaunom" (elle doit faire appelle à d'innombrables autres fonctions très certainement). Au final, en ayant exploré chaque fichier psd, et en notant que seul celui qui m'interesse est désigné par matlab par défaut, ...., sachant que j'ai supprimé 10log10(), et que cela m'afffiche malgré tout un truc en 10log10, et bien je ne vois pas du tout le probleme, et je crois que c'est tres mal parti ......

il faudrait en fait que j'arrive à récupérer le tableau de valeur avant que matlab le plot et que je modifie moi même mes valeurs .... mais j'ai l'impression que les vecteurs plotés sont des vecteurs "temporaires"..........

Une idée ???

Merci d'avance.

[aliasfdr]

c'est normal que le save as ne marche pas car le nom de la fonction doit être le même que m file donc il faut changer le nom de ta fonction (function myfun )et ensuite l'enregistrer sous ce même nom .

Je te conseille d'utiliser des break points pour voir là où ça cloche ( le plot reste en Db )
bon courage a+

[jujuf1]

c'est ce que j'avais fait. Y'a d'autre choses qui clochent. Tant pis je laisse tomber. En tout cas merci pour tout.

[jujuf1]

Sinon j'ai trouvé ça, mais je sais pas si le script est fiable. En tout cas ça a l'air bien.

http://www.mathworks.com/matlabcentr...bjectType=FILE