Discussion :

[patricx]

bonjour

S'il vous plait je créé un filtre butter

Code :

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

[B,A]=butter(20,0.6,'high');

ensuite je lis un fichier .wav

Code :

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

[elise,fs]=wavread('Elise.wav');

je créé un signal Blanc

Code :

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

sb=rand(size(elise));

je filtre le signal

Code :

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

y=filter(sb,B,A);

si j'ajoute le bruit filtré

Code :

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

c=y+elise;

j'ai cette erreur

Code :

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

1
2
3
4
5
??? Error using ==> plus
Matrix dimensions must agree.
 
Error in ==> sigBruite at 30
c=y+elise;

j'ai meme augmenté l'ordre du filtre ça ne resoud pas le problème

s'il vous plait comment y remedier

Merci

[Invité]

Bonjour,

comment y remedier

Simplement en regardant plus précisément l'emploi des fonctions que tu utilises:

y = filter(b,a,X) filters the data in vector X with the filter described by numerator coefficient vector b and denominator coefficient vector a. If a(1) is not equal to 1, filter normalizes the filter coefficients by a(1). If a(1) equals 0, filter returns an error.

[patricx]

Merci

j'ai plus d'erreur liée au taille des matrices

mais j'ai ces warnings

Code :

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

1
2
3
4
5
6
Warning: Data clipped during write to file:musique.wav 
> In wavwrite>PCM_Quantize at 287
  In wavwrite>write_wavedat at 309
  In wavwrite at 138
  In sigBruite at 36
>>

et à la fin mon son a quand meme un un arrière son

mon code

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
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
%Exercice3 Signal  bruité
function []= sigBruite()
 
%filtre passe-haut de type Butterworth en utilisant la fonction butter
[B,A]=butter(20,0.6,'high');
 
%Representation de la réponse en  frequence
[h1,w1]=freqz(B,A,512);
plot(w1/pi,abs(h1));
 
hold on;
 
%filtre passe-bas de type Butterworth en utilisant la fonction butter
[C,D]=butter(20,0.6,'low');
 
%Representation de la   réponse en  frequence
[h2,w2]=freqz(C,D,512);
plot(w2/pi,abs(h2));
 
% Lecture  du  fichier Elise.wav
[elise,fs]=wavread('Elise.wav');
 
%signal(bruit) Blanc de meme taille qu'elise
sb=rand(size(elise));
 
% Filtre du signal bruité  à  partir du  filtre passe  haut 
y=filter(B,A,sb);
 
%Ajout du signal  blanc filtré au  son  Elise 
c=y+elise;
 
% Filtre du son Elise bruité  à  partir du  filtre passe bas 
fpbs=filter(C,D,c);
 
%Ecriture du signal  filtré
wavwrite(fpbs,fs-1,'musique.wav')
%a=size(y)
%b=size(elise)

Merci

[Jerome Briot]

Une simple recherche => http://www.developpez.net/forums/d94...e/#post5312351