Traitement du son

Le code suivant permet charger et écouter un fichier .wav:

[s,fs]=wavread('coyote_howl.wav');
sound(s,fs);

Voici les fichiers de sons:
coyote_howl.wav
t2_hasta_la_vista.wav

Exercice 1
Charger le fichier coyote_howl.wav. Calculer sa transformée de Fourier discrète. Enlever des fréquences
a) bases (filtre passe-haut),
b) hautes (filtre passe-bas),
c) dont les coefficients sont petits en valeur absolu (compression, debruitage (dans le cas de debruitage, ajouter un bruit blanc au signal avant de commencer)).
Appliquer la transformée de Fourier discrète inverse. Afficher et écouter le résultat. Expérimenter avec des autres fichiers .wav.

Le résultat devrait sembler comme suit :

Solution ici

Traitement des images

Pour la suite on aura besoin du module SIVP qu’on peut installer en tapant dans la console : atomsInstall SIVP. Puis il faut fermer et rouvrir SciLab.

Les images à tester sont :
smiley.png
zebra.png
Joseph_Fourier.jpg

image=imread('smiley.png');//Joseph_Fourier.jpg ou zebra.png
grayscale=rgb2gray(image);//l'image est en couleur, il faut le convertir
Dgray=double(grayscale);//le type n'est pas bon, il faut le convertir
transform=fft(Dgray);
//imshow(image);
imshow(real(transform));

Exercice 2 Expérimenter avec filtres passe-bas et passe-haute et la compression pour les images.