[subplot] conflit plot et surf

Version imprimable

Bonjour,

Voici un code d'ouverture de fichier de signal (.gcf) dans lequel je souhaite tracer en subplot: le signal, le spectrogramme et la FFT :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
 
 
clear all
close all
 
 
RootDir=['C:\Documents\Datas\'];
 
 
%listDir = dir(RootDir)
 
FileName = '6c47e2';
 
cd([RootDir,FileName]) ;
 
rep=[RootDir,FileName] ;
ext = '*.gcf';
 
chemin = fullfile(rep,ext);
 
list = dir(chemin);
 
 
% lecture en boucle des signaux dans dossiers
for n = 1:numel(list)
 
    [list(n).name,ext(3:end)]
 
    [signal,streamID,sps,ist]=readgcffile([list(n).name]);
 
    signalRaw = signal ;
    ind = find(isnan(signal));
    signal(ind)=[] ;
 
    taille = numel(signal);
    duree = taille/sps ; % sps = sample frequency
    temps = linspace(ist,ist+(duree*(1/24/60/60)),taille) ;
 
    %% moyennage du signal = 0
    moySignal=mean(signal) ;
    signal = signal-moySignal ;
 
    %% apodisation par fenetre de hanning
    nH = 3000 ;
    pourcentage_echantillon_han = 1 ; % ! en %
    [windowHann,SigHann,maxSigx] = HannWinCell(signal,nH,pourcentage_echantillon_han);
    signal =  SigHann ;
 
 
    %% spectro
    % Compute the ifft of the signal
    L = numel(signal)
    NFFTbis = 2^nextpow2(L); % Next power of 2 from length of y
    Y = fft(signal,NFFTbis)/L;
    fbis = sps/2*linspace(0,1,NFFTbis/2+1);
 
 
    % % spectrogramme
    nfft = 128 ;
 
    [S,F,T,P] = spectrogram(signal,nfft,[],nfft,sps);
 
 
 
    figure ;
 
     subplot(3,1,1);
    plot(temps, signal,'b');
    hold on
    plot(temps,windowHann,'r')
    plot(temps,signalRaw,'g')
    legend('signal corrected', 'Hanning window','raw signal')
    %  xlim([temps(1) temps(end)]) ;
    title([list(n).name]) ;
    datetick('x','HH:MM','keeplimits')
 
    subplot(3,1,2);
    surf(T,F,10*log10(P),'edgecolor','none');
    axis tight;
    view(0,90);
 
    subplot(3,1,3);
    plot(fbis,2*abs(Y(1:NFFTbis/2+1)))
 
 
    clear signal streamID sps ist ;
end

Seulement lors de la création du spectrogramme (subplot(3,1,2)), le subplot précédent (du signal) disparait laissant celui-ci gris.. comme vous pouvez le voir:

http://img4.hostingpics.net/pics/214685MatlabPb01.png

Aucun message d'erreur, ni de warning.. je ne vois pas du tout d'où ça pourrait provenir..

Une idée?

Merci d'avance

26/11/2013, 15h56
Newenda

1 pièce(s) jointe(s)
Pour illustrer le problème, voici un .zip dans lequel test.mat contient toutes les variables nécessaires pour faire un test rapide suivant le code simplifié :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 clear all close all load test figure ; subplot(2,1,1); plot(temps, signal); subplot(2,1,2); surf(T,F,10*log10(P),'edgecolor','none');
Ça serait sympa si vous pouviez me dire si chez vous, oui ou non vous rencontrez ce problème d'affichage.

Merci d'avance

(PS: je travail avec la version R2013b)
26/11/2013, 16h15
le fab

salut

pas le temps de tester, mais ca vient peut etre de ton 'edgecolor','none' ou quelque chose comme ca
si c'est ca, tu peut regarder des fonction comme freezecolor

Merci de ta réponse.

mais malheureusement j'ai déjà essayé:
Code:

1 2 3 4 5 6 7 8 9 10 11 clear all close all load test figure ; subplot(2,1,1); plot(temps, signal); subplot(2,1,2); surf(T,F,10*log10(P));
Ça ne change rien.

mais si je change le premier plot et que j'écris par exemple :
Code:

1 2 3 4 5 6 7 8 9 10 11 clear all close all load test figure ; subplot(2,1,1); plot(rand(20000,1)); subplot(2,1,2); surf(T,F,10*log10(P),'edgecolor','none');
ça marche...

Je sais pas c'est magique non?

26/11/2013, 16h23
Adjen

Bonjour,

Testé à l'instant avec un matlab2013b sous win7 64 bits, j'ai la même chose que toi, avec ou sans edgecolor :

http://img4.hostingpics.net/pics/964165test.png
26/11/2013, 16h45
Newenda

Merci Adjen, au moins tu prouves que je ne suis pas fou !

Quelqu'un qui a une version alpha et non beta pourrait essayer pour voir? merci d'avance
26/11/2013, 16h48
le fab
cela viens peut être de la différence d'échelle entre ton signal(-1e5 5e4) et ta surface (-100 100) (ton rand lui reste entre 0 et 1)
la colormap de ta surface doit être automatiquement appliqué à ton plot

en tout cas en 2013a j'ai le même pb

une piste :
Code:

1 2 3 4 5 6 7 8 9 10 11 clear all close all load test figure ; subplot(2,1,1); plot(temps - temps(1), signal,'b-'); subplot(2,1,2); surf(T,F,10*log10(P),'edgecolor','none');
ca marche : c'est un pb d'échelle
26/11/2013, 17h06
Jerome Briot

C'est un soucis avec OpenGL

Change la propriété Renderer de l'objet Figure :

Code:

set(gcf, 'renderer', 'zbuffer')
26/11/2013, 17h20
Newenda

la commande de Dut fait tout bien marché,

Merci pour vos réponses, j'en sortais plus..

PS: c'est quand même dingue, j'ai passé ma journée à faire des recherches google sur ce problème et rien trouvé.. ça doit pas être si rare que ça quand même.