Discussion :

[titi_du_34]

Bonjour à tous,

J'explique tout d'abord ce que j'ai fait :

Génération d'un sinus glissant de 150 Hz à 1 kHz par pas de 100 Hz (soit 9 sinus à 9 fréquences différentes générés) via un HP puis récupération de la réponse à ce sinus avec deux microphones via audacity.
Exportation de la réponse sous matlab.

Ce que je veux faire :

1) la réponse de mes deux micros me donnant deux vecteurs de 402633 pour mes 9 fréquences, je souhaite tout d'abord transformer ce vecteur en une matrice représentant mes signaux à mes 9 fréquences différentes (sachant que la longueur du premier sinus est un peu plus courte que celles des autres signaux qui ont la même durée, car j'ai coupé le début du premier signal qui faisait un peu de bruit)

2) Une fois cela fait, j'aimerai trouver le moyen d'appliquer un filtre passe bande de butterwoerth glissant qui va ainsi me traiter la réponse de tous mes sinus (afin d'enlever toutes distorsion du signal)

3) Enfin le filtre passe bande atténuant toujours les fréquences filtrées (cf images jointes), je souhaite ne conserver que les oscillations à peu près constantes (un exemple d'image est donné avec un sinus pur à 1000 Hz)

Je joins mon travail en attaché.

Encore merci à ceux qui accepteront de jeter un oeil à ma demande

[titi_du_34]

Personne pour m'aider :'( ??

Bon j'ai quand même résolu le problème du découpage de mon signal global x1 qui a pour dimension 1587600x1 en une matrice Mic_1 de dimension 88200 points pour 18 fréquences à traiter.
Si quelqu'un peut m'aider pour le filtrage de mes 18 signaux à mes 18 fréquences différentes : voici mon code :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
Q=5; n=4; ftype='bandpass'; fe=44100;
f1=[]; f2=[]; y1=[]; A=[]; B=[];
 
for f=150:50:1000
f1=[f1 f.*(1-(1/(2*Q)))];
f2=[f2 f.*(1+(1/(2*Q)))];
Wn=[f1 f2]/(fe/2);
[A,B]=butter(n,Wn,ftype);
y1=filter(A,B,Mic_1).*hamming(M)
 
end

Encore merci

[Invité]

Bonjour,

Code :

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

1
2
f1=[f1 f.*(1-(1/(2*Q)))];
f2=[f2 f.*(1+(1/(2*Q)))];

Ici tu ajoutes au fur et à mesure les valeurs des itérations dans un tableau. Plutôt comme ceci :

Code :

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

1
2
f1 = f.*(1-(1/(2*Q)));
f2 = f.*(1+(1/(2*Q)));

Attention aux noms :

Code :

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

1
2
[B,A]=butter(n,Wn,ftype);
y1=filter(B,A,Mic_1).*hamming(M)

tu les as inversé dans les deux lignes, ce qui au final donne un résultat juste, mais autant ne pas s’emmêler les pinceaux si jamais le code vient à évoluer.

[titi_du_34]

Un grand merci à vous Jérôme pour votre aide , du coup ça marche