Signaux sinusoïdaux et Scilab

[invite97fdb67a]

Salut moi aussi je débute avec scilab, j'ai besoin de mesurer la phase en deux signaux sinus représentant des courants de charge. Les deux signaux sont chaqun sur une voie d'un fichier wave stereo. La voie de gauche provient d'un capteur à effet hall et contient du bruit BF qui peut fausser ma mesure de phase. Je doit donc :
1 lire les deux signaux
2 filtrer les deux signaux
3 afficher les signaux
4 et au passage afficher leur spectre

Je vais poster ma progression, si ça peut t'aider toi ou un autre.
-----

[invite97fdb67a]

Voici la lecture du fichier et son affichage :

clear;
n=400; //nombre d'échantillons limité au 400 premiers du fichier
[y,Fs,bits]=wavread("C:/Cmonitor.wav",n);Fs,bits
// y contient les données des signaux
//Fs est la fréquence d'échantillonage ex.44100
//bits est le nombre de bit pour lequel le signal est codé ex. 32
subplot(2,1,1)// selectionne le graph N°1 dans une matrice de graph de 2 lignes x 1 colonne (L,C,N°)
plot2d(y(1,: )) // trace un graph 2D avec toutes les données du canal de droite
subplot(2,1,2)
plot2d(y(2,: )) // idem canal de gauche

[invite97fdb67a]

Voici l'affichage des voies et de leur FFT.
Je ne sais pas pourquoi Scilab me demande d'appuyer sur une touche pour continuer et ce à plusieurs reprises mais le résultat s'affiche. Si quelqu'un le sait je suis preneur.


clear;
n=400;
[y,Fs,bits]=wavread("C:/Cmonitor.wav",n);Fs,bits
rs=y(1,: ) //right signal
ls=y(2,: ) //left signal

subplot(2,2,1)
plot2d(rs)
subplot(2,2,2)
plot2d(ls)

rf=fft(rs); //right FFT
lf=fft(ls); //left FFT

subplot(2,2,3)
plot(abs(rf(1:length(rf)/2)))
subplot(2,2,4)
plot(abs(rf(1:length(lf)/2)))

[JPL]

Pour une nouveau problème il faut créer une nouvelle discussion et non parasiter une discussion déjà ouverte. J'ai donc séparé ton sujet.

[A voir en vidéo sur Futura]

[invite97fdb67a]

Bon ok JPL, alors j'en suis a filtrer mes deux canaux sur une même fréquence de 300Hz mais là je bloque grave.

Dans la doc scilab il n'y a pas d'exemple pour la fonction filter si tant est que ce soit la bonne fonction à utiliser.

Est-ce quelqu'un a un tuyau ou mieux un exemple ?

Merci à tous