Discussion :

[Kichott]

Bonjour,

Je rencontre un petit problème avec la fonction, pwelch, cette fonction détermine la dsp, créer une figure et plot le résultat.

Mon problème est que je suis entrain entrain une petite interface ac guide.

Pour les autres graphiques il n'y a pas de problèmes, mais là je ne sais pas comment faire pour que ma fonction pwelch s'affiche dans "l'axes" que je veux...

Quelqu'un a une idée ?

Merci
Jonathan

[Caro-Line]

Comme a priori on ne peut définir l'axe de destination dans la fonction PWELCH, il doit suffire de faire en sorte que l'axe de destination soit mis en axe courant avant l'appel à cette fonction.
Tout simplement :

Code :

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

1
2
axes(handle de l'axe souhaité);
pwelch(...);

[Kichott]

C'est ce que j'ai fais mais ça ne marchais pas.

Enfin ca marche maintenant que je test chez moi avec une version assez récente de Matlab

Je vais réessayer mais j'ai l'impression qu'avec la version du taf (7.0..) ca ne marchait pas. Je confirmerai...

[Kichott]

Bon malheuresement je confirme

Code :

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

1
2
axes(handle de l'axe souhaité);
pwelch(...);

Ceci ne fonctionne pas avec la version R14 a mon boulot...
Quelqu'un aurait eventuellement une autre solution ? (A part demander une nouvelle version de Matlab )

[Caro-Line]

Ceci ne fonctionne pas avec la version R14 a mon boulot...

Etrange...
Quel est le comportement exactement ?
Le graphe se fait obligatoirement dans un nouvel objet axes ou alors il se met au hasard ?

Je n'ai pas cette fonction sur ma version mais est-elle éditable ?
Peut-être y a-t-il quelques lignes à modifier dedans pour améliorer le comportement.

Hum, ça doit venir de là :
Version 6.3 (R14SP2) Signal Processing Toolbox Software
=> Spectral Analysis Changes

Pas sur que tu puisses faire quelque chose (à part donc changer de version)

[duf42]

Bonjour,

Pourquoi ne pas récupérer les valeurs de PWELCH et les tracer ensuite à l'aide d'un PLOT :

Code :

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

1
2
3
4
5
6
 
% Calcul avec PWELCH
[ydata, xdata] = pwelch(...);
 
% Tracé des données
plot(axe_handles, xdata, ydata);

A tester...

Bonne journée,
Duf

[Kichott]

Etrange...
Quel est le comportement exactement ?
Le graphe se fait obligatoirement dans un nouvel objet axes ou alors il se met au hasard ?

Je n'ai pas cette fonction sur ma version mais est-elle éditable ?
Peut-être y a-t-il quelques lignes à modifier dedans pour améliorer le comportement.

Hum, ça doit venir de là :
Version 6.3 (R14SP2) Signal Processing Toolbox Software
=> Spectral Analysis Changes

Pas sur que tu puisses faire quelque chose (à part donc changer de version)

C'est assez bizarre, il se met dans un axes que j'ai dans mon gui mais juste le titre la courbe ne s'affiche pas.

Duf je crois avoir déjà essayer ça mais je réessayerai Lundi...

[Jerome Briot]

Duf je crois avoir déjà essayer ça mais je réessayerai Lundi...

Cela marche sous MATLAB 6.5.1 car dans PWELCH il y a :

Code :

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

1
2
3
4
5
6
7
8
9
if nargout==0, % Plot when no output arguments are specified      
   yscale = 'db';
   titlestring = 'Welch PSD Estimate';
   psdplot(Pxx,w,units,yscale,titlestring);    
elseif nargout == 1,
   varargout = {Pxx};
elseif nargout == 2,
	varargout = {Pxx,w};
end

Donc si tu mets des arguments de sortie, ke tracé n'apparait pas.

Le problème vient de la fonction PSDPLOT qui contient la ligne suivante :

Code :

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

newplot;

Tu fais donc une copie de PSDPLOT que tu enregistres sous MYPSDPLOT.M

Et dans MYPSDPLOT.M, tu supprimes la ligne ci-dessus.

Ensuite tu fais comme ceci pour tracé dans l'axe courant :

Code :

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

1
2
3
4
5
Fs = 1000;   t = 0:1/Fs:.296;
x = cos(2*pi*t*200)+randn(size(t));
[Pxx,W] = pwelch(x,[],[],[],Fs,'twosided');
 
mypsdplot(Pxx,W,'Hz','db','Welch PSD Estimate');

ou pour tracer dans un axe en particulier :

Code :

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

1
2
3
4
5
6
Fs = 1000;   t = 0:1/Fs:.296;
x = cos(2*pi*t*200)+randn(size(t));
[Pxx,W] = pwelch(x,[],[],[],Fs,'twosided');
 
axes(...)
mypsdplot(Pxx,W,'Hz','db','Welch PSD Estimate');

[Kichott]

Code :

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

1
2
3
4
5
% Calcul avec PWELCH
[ydata, xdata] = pwelch(...);
 
% Tracé des données
plot(axe_handles, xdata, ydata);

En effet cela trace quelque chose mais ce n'est pas le bon résultat.

Moi dans pwelch, j'ai ça :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
function varargout = pwelch(x,varargin)
 
error(nargchk(1,7,nargin));
error(nargoutchk(0,3,nargout));
 
esttype = 'psd';
if nargin > 1 & isstr(varargin{end}) & strcmpi(varargin{end},'ms'),
    esttype = 'ms';
    varargin(end)=[];
end
 
% Possible outputs are:
%       Plot
%       Pxx
%       Pxx, freq
[varargout{1:nargout}] = welch(x,esttype,varargin{:});

Je ne pense pas que j'ai les données x et y en sortie.

[Kichott]

Code :

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

1
2
3
4
5
6
Fs = 1000;   t = 0:1/Fs:.296;
x = cos(2*pi*t*200)+randn(size(t));
[Pxx,W] = pwelch(x,[],[],[],Fs,'twosided');
 
axes(...)
mypsdplot(Pxx,W,'Hz','db','Welch PSD Estimate');

Dut j'ai aussi essayé ça, il y a bien un trace qui s'affiche mais ce n'est pas le bon résultat non plus.

[Kichott]

Bon comme je trouvais pas de solution, j'ai finalement j'réussi a changer de version au boulot, donc je n'ai plus ce problème.