Chasse au Loup

Protocoles

Une liste des différents protocoles expérimentaux mis en oeuvre durant notre stage.

Trouver la fréquence du corps sujette au loup

On fait un premier enregistrement du phénomène du loup. Grace à un spectrographe, on fait une première estimation de la fréquence du corps en repérant le battement. Dans le même temps, on repère une zone vibrante de la table manuellement.

On place un micro sur cette zone. On fait un balayage de fréquence autour de la valeur repérée à l'aide d'un GBF branché à un amplificateur et un haut-parleur pour trouver la fréquence où la zone est la plus vibrante.

Cartographie de la table

On place sur le violoncelle un quadrillage léger transparent. On fixe la valeur du GBF à la fréquence du corps trouvée, et on déplace le micro sur chaque point du quadrillage pour enregistrer l'amplitude de vibration. On en tire une carte grâce au programme Octave ci-dessous :

clear

A = load("vib_carte.txt")

%----------MOYENNAGE DES 0

k=0;

%Prend la matrice a(13x8) contenue dans A(14x9)
for (i=2:13)
   for (j=2:8)
     if (A(i,j)==0)
         k=0;
       if (A(i-1,j)!=0)
         A(i,j)=A(i,j)+A(i-1,j);
        k=k+1;
       endif
  
        if (A(i,j-1)!=0)
         A(i,j)=A(i,j)+ A(i,j-1);
        k=k+1;
       endif
  
       if (A(i+1,j)!=0)
        A(i,j)=A(i,j)+A(i+1,j);
        k=k+1;
       endif
  
       if (A(i,j+1)!=0)
        A(i,j)=A(i,j)+A(i,j+1);
        k=k+1;
       endif
       
       if (A(i+1,j+1)!=0)
         A(i,j)=A(i,j)+A(i+1,j+1);
         k++;
        endif
       
        if (A(i-1,j+1)!=0)
         A(i,j)=A(i,j)+A(i-1,j+1);
         k++;
        endif
       
        if (A(i+1,j-1)!=0)
         A(i,j)=A(i,j)+A(i+1,j-1);
         k++;
        endif
       
        if (A(i-1,j-1)!=0)
         A(i,j)=A(i,j)+A(i-1,j-1);
         k++;
        endif
       
       if (k!=0)
         A(i,j)=A(i,j)./k;
       endif
      endif
    endfor
  endfor
 
i=0; j=0;

%Prend le coin inferieur gauche de A
A(14,1) = (A(13,1)+A(14,2)+A(13,2))/3;
%A(14,1) = (A(13,1)+A(14,2))/2;

%Prend le coin inferieur droit de A
A(14,9) = (A(13,9)+A(14,8)+A(13,8))/3;
%A(14,9) = (A(13,9)+A(14,8))/2;

%Prend le bord sup du A
for (i=1)
   for(j=2:8)
     if(A(i,j)==0)
       k=0;
        
       if (A(i,j-1)!=0)
         A(i,j)=A(i,j)+ A(i,j-1);
        k=k+1;
       endif
  
       if (A(i+1,j)!=0)
        A(i,j)=A(i,j)+A(i+1,j);
        k=k+1;
       endif
  
       if (A(i,j+1)!=0)
        A(i,j)=A(i,j)+A(i,j+1);
        k=k+1;
       endif
       
       if (A(i+1,j+1)!=0)
         A(i,j)=A(i,j)+A(i+1,j+1);
         k++;
        endif
       
       if (A(i+1,j-1)!=0)
         A(i,j)=A(i,j)+A(i+1,j-1);
         k++;
        endif
       
       if (k!=0)
         A(i,j)=A(i,j)./k;
       endif
      endif
    endfor
  endfor

i=0; j=0;

%Prend le bord inferieur de A
for (i=14)
   for(j=2:8)
     if(A(i,j)==0)
       k=0;
        
       if (A(i,j-1)!=0)
         A(i,j)=A(i,j)+ A(i,j-1);
        k=k+1;
       endif
  
       if (A(i-1,j)!=0)
        A(i,j)=A(i,j)+A(i-1,j);
        k=k+1;
       endif
  
       if (A(i,j+1)!=0)
        A(i,j)=A(i,j)+A(i,j+1);
        k=k+1;
       endif
       
       if (A(i-1,j+1)!=0)
         A(i,j)=A(i,j)+A(i-1,j+1);
         k++;
        endif
       
       if (A(i-1,j-1)!=0)
         A(i,j)=A(i,j)+A(i-1,j-1);
         k++;
        endif
       
       if (k!=0)
         A(i,j)=A(i,j)./k;
       endif
     endif
   endfor
 endfor
 
 
i=0; j=0;

%Prend le bord gauche de A
for (i=2:13)
   for(j=1)
     if(A(i,j)==0)
       k=0;
        
       if (A(i,j+1)!=0)
         A(i,j)=A(i,j)+ A(i,j+1);
        k=k+1;
       endif
  
       if (A(i+1,j)!=0)
        A(i,j)=A(i,j)+A(i+1,j);
        k=k+1;
       endif
  
       if (A(i-1,j)!=0)
        A(i,j)=A(i,j)+A(i-1,j);
        k=k+1;
       endif
       
       if (A(i+1,j+1)!=0)
         A(i,j)=A(i,j)+A(i+1,j+1);
         k++;
        endif
       
       if (A(i-1,j+1)!=0)
         A(i,j)=A(i,j)+A(i-1,j+1);
         k++;
        endif
       
       if (k!=0)
         A(i,j)=A(i,j)./k;
       endif
      endif
    endfor
  endfor

i=0; j=0;

%Prend le bord droit de A
for (i=2:13)
   for(j=9)
     if(A(i,j)==0)
       k=0;
        
       if (A(i,j-1)!=0)
         A(i,j)=A(i,j)+ A(i,j-1);
        k=k+1;
       endif
  
       if (A(i+1,j)!=0)
        A(i,j)=A(i,j)+A(i+1,j);
        k=k+1;
       endif
  
       if (A(i-1,j)!=0)
        A(i,j)=A(i,j)+A(i-1,j);
        k=k+1;
       endif
       
       if (A(i+1,j-1)!=0)
         A(i,j)=A(i,j)+A(i+1,j-1);
         k++;
        endif
       
       if (A(i-1,j-1)!=0)
         A(i,j)=A(i,j)+A(i-1,j-1);
         k++;
        endif
       
       if (k!=0)
         A(i,j)=A(i,j)./k;
       endif
      endif
    endfor
  endfor

%imagesc(A,[min(min(A)),max(max(A))])
surf(A)
shading('interp')
colorbar

Calcul de la masse et raideur du violoncelle

On place différentes masses sur un ventre de l'instrument et on effectue à chaque fois une analyse fréquencielle pour trouver la fréquence du corps. Une régression linéaire nous permet de trouver K et M (voir rapport).

Diagramme en fréquence

On place un micro devant un ventre de vibrations et on effectue analyse fréquencielle d'un signal impulsionnel avant d'enregistrer un document en trois colonnes (fréquence/dB/phase). On réitère l'expérience, puis on en sort un diagramme grâce au code suivant: