Discussion :

[Krillz]

Bonjour, ce doit être tout simple mais je ne trouve pas comment faire.

Je fais une opération :

Code :

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

R=xcorr(xl,xr,'coeff');

où xl et xr sont des vecteurs ligne.

Et j'exporte dans un fichier à l'aide de save :

Code :

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

1
2
 
save('corr.dat','final','-ascii')

Néanmoins j'aimerais avoir le numéro de ligne dans mon fichier corr.dat

Code :

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

1
2
3
4
5
 
1     R(1)
2     R(2)
3     R(3)
....

Ou bien en gros, l'antécédent de chaque valeur de R, suivi de la valeur de R


Merci de votre aide!

[Loïc B.]

Pourquoi ne pas enregistrer un tableau intermédiaire avec deux colonnes et exporter ce dernier ?

[Krillz]

J'ai essayé hier mais j'ai eu une erreur de concaténation.

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
 
[x,Fs,bits] = wavread(['2000Hz40.wav']);
 xl = x(:,1);
 xr = x(:,2);
 
R=xcorr(xl,xr,'coeff');
 
 figure(2)
 stem(R)
 
 [p m]=max(R);
 phase_shift=(m-length(xl))/Fs;
 
 
x=1:1:length(R); 
 
final=[x;R]
 
save('corr2000Hz40.dat','final','-ascii')

Et l'erreur :

Code :

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

1
2
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.

[nahouto]

tu as essayé

Code :

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

final=[x,R]

au lieu de

Code :

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

final=[x;R]

?

[Krillz]

Code :

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

1
2
3
4
5
??? Error using ==> horzcat
CAT arguments dimensions are not consistent.
 
Error in ==> findphase2000onefile at 16
final=[x,R]

[nahouto]

A quoi ressemble R ?

[Krillz]

R est la corrélation (avec xcorr) de deux signaux obtenus à partir de 2 fichiers sonores (avec wavread).

Je pensais qu'il s'agissait simplement d'une fonction comme une autre, mais apparemment pas...

[nahouto]

tu pourrais donner plus de détails sur R ? Classe, Taille, exemple...

[Krillz]

Je suis désolé mais je ne saurais pas te donner plus de détails que le code que j'ai posté au dessus (c'est le code en intégralité)...

Code :

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

1
2
3
4
5
[x,Fs,bits] = wavread(['2000Hz40.wav']);
 xl = x(:,1);
 xr = x(:,2);
 
R=xcorr(xl,xr)

Plus de détails sur la fonction xcorr :
http://www.mathworks.com/access/help...ient=firefox-a

[nahouto]

Tu ne peux pas me dire ce qu'est ta variable R ? (mets un breakpoint, regarde le workspace...)

Je n'ai pas la Toolbox, donc je ne peux pas tester sur ma machine.

[Krillz]

R < 292863x1 double>

min -1
max 1

[Caro-Line]

Vu la taille de R il est étonnant que [x,R] ne fonctionne pas

Essaye ceci :

Code :

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

1
2
3
4
R=rand(292863,1);%taille que tu nous a donné
x=1:length(R);
final=[x',R];%concaténation horizontale de 2 colonnes
final(1:10,:) %pour valider

[Krillz]

Merci de vos réponses,

Caro-Line, collé tel quel dans Matlab ton code fonctionne.

Néanmoins, je ne comprends pas le "prime" dans final=[x',R] ....

[Caro-Line]

Réponse ici : Quelle est la différence entre ' et .' ?

Le plus simple pour comprendre quelque chose est encore de le tester sur de petits exemples :

Code :

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

1
2
x=[1 2 3]
x'

[Krillz]

Merci beaucoup Caro-Line !