Discussion :

[xMaudinette]

Bonjour

J'ai créé un sinus sweep dans Scilab, afin de faire de la simulation d'enceinte. Seulement, je voudrais récupérer l'enveloppe du signal ou sinon les valeurs max de la matrice que j'ai créé.
Le problème c'est que max(A) ne permet de récupérer qu'une seule valeur, or j'en ai besoin de plusieurs.
Bien sur, si je connaissais la période du signal ca serait bien mais la particularité du sweep c'est que justement sa fréquence augmente en fonction du temps.
Le sinus dure 1min, avec une fréquence d'échantillonnage à 96k donc la matrice est assez grande (juste un peu, 5,7M CA VA)

J'ai essayé différentes méthodes, j'ai créer un petit algo 'enveloppe', le problème, après enregistrement du sweep à travers mon enceinte, il y a du bruit, et c'est plutot compliqué dans les basses fréquences.

Des suggestions ?


Voici comment je crée mon sweep:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
fe=96000;
T=60;
w1=20;
w2=20000;
te=1/fe;
n=0:(T*fe);
tailleFic=(fe*T)+4800;
 
K=(T*w1)/log(w2/w1);
L=T/log(w2/w1);
 
sweep=sin(K*(exp(n/(fe*L))-1));
f=(K/L)*exp(n/(fe*L));

[Vincent PETIT]

Salut,
Tu peux nous montrer une image du signal ? Je n'ai pas Scilab là où je suis.

[xMaudinette]

Salut,
Tu peux nous montrer une image du signal ? Je n'ai pas Scilab là où je suis.

Le signal n'est pas affiché en entier, ca devient un gros bloc rouge à la fin vu que la période est vraiment petite.
Et les abscisses sont en samples

[Vincent PETIT]

Ok mais c'est un chirp ça !

C'est une pure modulation de fréquence et du coup son enveloppe est rectangulaire, c'est juste les min et max que tu connais déjà, -1 et +1.


En revanche si le signal eu été celui ci, donc modulé en amplitude et en fréquence alors l'enveloppe est le signal en pointillé qui passe par tout les maximums et minimums.
Cette enveloppe est obtenue grâce à la transformée de Hilbert (de mémoire mais je suis quasi sur)





Exemple d'enveloppe rectangulaire d'un chirp.

[xMaudinette]

Je sais bien que son enveloppe est rectangulaire, c'est parce que je ne l'ai pas encore passé à travers l'enceinte, pour l'instant je ne fait que tester.
Je veux trouver les indices des maximums pour pouvoir récupérer l'enveloppe du signal passé dans l'enceinte.

J'ai essayer la fonction de Hilbert, mais quand il y a du bruit c'est vraiment pas top :/

[Vincent PETIT]

Ah ok !

[...]Je veux trouver les indices des maximums pour pouvoir récupérer l'enveloppe du signal passé dans l'enceinte.[...] afin de faire de la simulation d'enceinte. [...]

De microphone plutôt !
Parce que depuis le début, je pensais que tu voulais synthétiser un signal pour l'envoyer sur des enceintes alors que tu simules une réception. C'est pour ça que j'avais l'impression que tu connaissais déjà la réponse à ta question.
Bref !

J'ai essayer la fonction de Hilbert, mais quand il y a du bruit c'est vraiment pas top :/

C'est vraiment pas top à ce point ?

C'est la méthode la plus classique pourtant. Tu as une capture d'écran ?

[xMaudinette]

Non non, de la simulation d'enceinte! Enfin plus précisement : de la simulation d'enceinte+micro+acoustique de la salle!
Je passe un sweep dans une enceinte, j'enregistre avec un micro oui, mais j'obtiens la réponse fréquentielle de l'enceinte et non du micro.
C'est un système de boite noire en fait : je connais l'entrée (le sinus sweep), et la sortie mais ce que je veux c'est la réponse fréquentielle!

Et si je connaissais l'indice des max du sweep, je peux récupérer ma réponse fréquentielle!
Après enregistrement du sweep dans mon enceinte j'obtient un fichier de la forme :


Avec Hilbert, je peux obtenir la RF :

Seulement il y a trop de bruit dans les BF j'ai besoin que ça soit plus précis, ya 15 dB de bruit!