Discussion :

[Gazelle Arctique]

Bonjour à tous !

Je me permets de vous solliciter au sujet de mon code.
J'aimerais définir le spectre de mon signal dans différents intervalles de temps.

Voici mon signal :


Avec mon code, je suis capable de déterminer le spectre de ce signal.
Mais l'idée est vraiment de segmenter ce signal en 5 zones : [T1 T2] [T3 T4] [T5 T6] [T7 T8] [T9 T10] et de calculer le spectre correspondant à chaque zone.

L'idée a beau être claire dans ma tête, elle n'est pas claire dans mes doigts. J'ai essayé d'écrire mon code de 1001 facons, je ne suis pas encore tombée sur la bonne...

Voici une partie de mon code actuel, qui me semble le plus proche de la solution:

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
17
18
19
20
21
22
Time_RR1_16Hz = interp1(times_rr, RR_interval_ms, Time_orig, 'spline');
    Time_RR1_16Hz =  Time_RR1_16Hz - mean( Time_RR1_16Hz);
 
 
   % TRANSITION POINTS - RR
   x_RR = [1115 14400 15560 24220 24940 28740 30910 38680 39340 48010];
   y_RR = [35.03 27.06 -46.51 -40.51 25.2 30.75 -49.61 -34.13 32.19 4.636];
 
   Zone1_RR = (x_RR(1:2), y_RR(1:2));
   Zone2_RR = (x_RR(3:4), y_RR(3:4));
   Zone3_RR = (x_RR(5:6), y_RR(5:6));
   Zone4_RR = (x_RR(7:8), y_RR(7:8));
   Zone5_RR = (x_RR(9:10), y_RR(9:10));
 
   p = 256; 
   [Pxx,f] = pburg( Zone1_RR,p,1028,16); % spectrum estimate
   figure (1);
   subplot(3,1,1);
   plot(Time_orig,  Zone1_RR,'k'); xlabel('Time [s]');ylabel('RRi [ms]');
   xlim([0 3000]); ylim([-50 300]);
   subplot(3,1,[2 3]);
   plot(f,Pxx,'k'); xlim([0 0.5]); xlabel('Frequency [Hz]'); ylabel(sprintf('PSD [ms^2/Hz]'));

Seulement, j'ai une erreur de : ) } ) quand je définis les zones dans lesquelles je veux déterminer le spectre... Je dois avoir les yeux rivés dessus depuis trop longtemps.

Néanmoins, pensez-vous que je suis sur la bonne voie ?

Je vous remercie d'avance pour votre aide !

[Jerome Briot]

A tout hasard :

Code :

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

Zone1_RR = [x_RR(1:2), y_RR(1:2)];

[Gazelle Arctique]

Merci pour ta réactivité et pour ta réponse !

J'obtiens l'erreur suivante :

Error using arburg (line 56)
The length of input vector X
must at least 257.

Error in arspectra (line 76)
[a,v] = feval(method,x,p);

Error in pburg (line 88)
[Pxx,freq,msg,units,~,options,msgobj]
=
arspectra(method,x,p,varargin{:});

Error in Session11 (line 79)
   [Pxx,f] = pburg(
   Zone1_RR,p,1028,16); %
   spectrum estimate

Alors j'ai vérifié ce que Zone1_RR contient et effectivement, il y a un problème de longueur.
Zone1_RR contient 4 valeurs : 1115 14400 35,03 27,06.

L'idée serait que Zone1_RR corresponde à la portion du signal entre les points T1 et T2.

Serait-ce juste un problème d'écriture ?

Je vous remercie pour votre aide !

[Jerome Briot]

Quelque chose comme ceci ?

Code :

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

Zone1_RR = y(x_RR(1):x_RR(2));

[Gazelle Arctique]

Je viens d'essayer avec :

Code :

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

Zone1_RR =   y_RR(x_RR(1):x_RR(2));

Et j'obtiens l'erreur suivante :

"Index exceeds matrix dimensions.

Error in Session11 (line 72)
   Zone1_RR =    y_RR(x_RR(1):x_RR(2));"

Comme si j'en demandais plus que ce qu'il y a dans la matrice...

Une idée quelconque pour m'aider à avancer ?

Merci d'avance pour votre aide...

[Gazelle Arctique]

C'est bon ! je pense avoir trouvé :

Code :

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

 Zone1_RR = Time_RR1_16Hz(1115 : 14400);

Si je regarde les valeurs contenues dans Zone1_RR, il y a bien les valeurs désirées.

SAUF QUE maintenant, j'ai un problème d'affichage avec plot :

Vectors must be the same lengths

[Gazelle Arctique]

J'AI TROUVE LA SOLUTION !!!!!!

En fait, c'est tout bêtement parce que je n'avais pas sélectionné le même intervalle sur l'axe temporel !

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
   % TRANSITION POINTS - RR
   x_RR = [1115 14400 15560 24220 24940 28740 30910 38680 39340 47905];

  Zone1_y = Time_RR1_16Hz(1115 : 14400);
  Zone1_x = Time_orig (1115: 14400);
  
  
   p = 256; % AR model order
   [Pxx,f] = pburg( Zone1_y,p,1028,16); % spectrum estimate
   figure (1);
   subplot(3,1,1);
   plot(Zone1_x, Zone1_y,'k'); xlabel('Time [s]');ylabel('RRi [ms]');

ALLEZ LOUIA !

Merci pour vos pistes !