Discussion :

[Neocid]

Bonsoir,

J'ai un petit soucis avec mon programme. j'ai besoin de multiplier deux signaux sinusoïdaux mais une fois ceux-ci définit et multipliés il me fait une erreur out of memory :p

Code :

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

1
2
3
4
5
6
7
8
 
 fid = fopen ('xp.wav');   %ouverture du fichier
son = fread(fid,inf,'int32');
fclose(fid);
fe=8000;
t_porteuse=(1:length(son))/fe; % définition de la durée de la porteuse
porteuse = cos(2*pi*12800*t_porteuse); % porteuse
module = son * porteuse; % modulation

??? Error using ==> mtimes
Out of memory. Type HELP MEMORY for your options.

Une idée ?

[Jerome Briot]

Si tu veux multiplier les deux signaux éléments pas éléments, il faut faire comme ceci :

Code :

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

module = son .* porteuse; % modulation

[Neocid]

Merci dut,

j'obtiens une erreur de matrice maintenant.

comment je fais pour mettre ma "porteuse" ( voir code précédent ) dans une matrice de même taille que mon "son" ( qui est matrice 1*2 si j'ai bien compris ? )

Code :

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

1
2
3
4
 
 module = son .* porteuse;
??? Error using ==> times
Matrix dimensions must agree.

merci

[Jerome Briot]

Pourquoi n'utilises-tu pas WAVREAD pour lire le fichier wav ?

Sinon, pour connaitre la taille d'une variable, soit tu utilises la fonction WHOS, soit tu regardes dans la fenêtre Wokspace du MATLAB Desktop

[Neocid]

Je n'utilisais pas wavread parce que cela me donnait une erreur de lecture. Mais ca a l'air de marcher o_O' ( je me sens sacrément nul sur matlab ) .

Alors en faisant z = wavread('xp.wav');

j'ai les values :
z: 106150x2 double.
'son' ( fait avec fread ) : 106161x1 double
porteuse : 1x106161 double

ca devrait pas coller ? pourquoi il me fait une erreur de matrice ?

( en faisant z .* porteuse; j'obtiens encore un "Matrix dimensions must agree").

Merci.

[Jerome Briot]

Alors en faisant z = wavread('xp.wav');

j'ai les values :
z: 106150x2 double.

Si WAREAD lit correctement le fichier, le fait que z soit un vecteur Nx2 signifie que le signal est de stype stereo (2 canaux : gauche-droit).

Demande la nombre maximum d'arguments de sortie à la fonction :

Code :

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

[y,fs,nbits,opts]=wavread('xp.wav')

Et regarde les informations (surtout le nChannels) que retourne ceci :

Code :

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

opts.info

Pour jouer les contenu du fichier, il y a deux possibilités :

Code :

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

1
2
3
wavplay(y,fs);
%ou
sound(y,fs);

'son' ( fait avec fread ) : 106161x1 double
porteuse : 1x106161 double

ca devrait pas coller ? pourquoi il me fait une erreur de matrice ?

( en faisant z .* porteuse; j'obtiens encore un "Matrix dimensions must agree").

Les dimensions des vecteurs ne correspondent pas.

son : Nx1
porteuse : 1xN

Tu dois faire

Code :

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

son .* porteuse.'

[Neocid]

bonjour,

Merci pour les explications.
Malheureusement même en faisant

Code :

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

1
2
3
son .* porteuse
%ou
 porteuse .* son

il me dit toujours que ma matrice n'est pas de même taille.

Pourriez vous me renseigner sur la façon de créer mon signal sinusoïdale pur et qu'il soit contenu dans une matrice de même taille que mon 'son' svp ?

[Jerome Briot]

As-tu lu ma dernière remarque ?

Tu dois faire

Code :

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

son .* porteuse.'

[Neocid]

Oups, toutes mes excuses le ' je pensais que c'était une fin de code.

Bon en effet cela se multiplie bien et j'ai une jolie fft avec les spectres centrés sur mes fréquences de porteuse !! merci !!!

Maintenant j'obtiens une erreur lors de l'utilisation de filtres


je cherche à filtrer mon signal 'xp.wav' à la fréquence de 18200 khz.
voila mon code

Code :

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

1
2
3
4
5
6
7
 
[N, Wp] = ELLIPORD(1/fs, 18200/fs, 1, 60)
[B,A] = ELLIP(1,1,60,Wp)
Z = FILTER(B,A,z)      % z étant mon wavread('xp.wav')
 
 
??? Undefined function or method 'FILTER' for input arguments of type 'double'.

encore un soucis de matrice double.
J'ai essayer de trouver d'autre possibilité de faire des filtres ( notemment avec fir1 ) et cela me donne la même errreur

Existe t'il un moyen de filtrer un signal double ? ( ca veux pas dire grand chose en effet ) ou alors est-il possible de ne récuperer qu'un signal composé d'un seul élément ? ( en gros passer en mono le son si j'ai bien compris ? )

Merci d'avance !

[Jerome Briot]

Le nom des fonction s'écrit toujours en lettres minuscules.

Essaye ceci :

Code :

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

Z = filter(B,A,z);

[Neocid]

Damn you good !

Merci beaucoup Dut