%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%PETIT TUBE DE KUNDT - CALCUL DU COEFFICIENT D'ABSORBTION 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% micro 1 près  du HP en première configuration, puis près échantillon après inversion
% micro 2 près échantillon en première configuration, puis près HP après inversion
[nom,chemin]=uigetfile;
toto = [chemin,nom];
load(toto);
%load kundtessai1.mat; %récupération de la fonctionde transfert H1 mesurée
h1=donnee;
h1=h1;
fq=0:Fe/Ne:Fe/Ne*(Ne-1)/2.56;
[nom,chemin]=uigetfile;
toto2 = [chemin,nom];
load(toto2);
%load kundtessai2.mat; %récupération de la fonctionde transfert H1 mesurée après inversion des micros
h2=donnee;
h2=1./h2;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%fonction de transfert moyennée après échange de micros
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

h=(abs(h1).*abs(h2)).^0.5.*exp(j.*(angle(h1)+angle(h2))/2); 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%calcul des paramètres d'absorption et de réflexion
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Hr=h(1:max(size(fq)),1)' ;% fonction de transfert moyennée sur bande de fréquence d'analyse
e=input('épaisseur de l"échantillon en m (exemple : 40e-3) :      ');
%e=40e-3 %épaisseur échantillon
d=input('distance échantillon - micro le plus éloigné en m :       ');
%d = 73e-3 + 25e-3 - e %distance échantillon - micro le plus éloigné
s=input('distance entre les microphones en m (25 e-3 pour le tube du LISMMA) :       ');
%s=25e-3 ;%distance inter-micros
c=342 ;%célérité du son
k=2*pi.*fq/c ;% nombre d'onde
Gr=exp(2*j*d.*k).*(Hr-exp(j.*k*s))./(exp(-j.*k*s)-Hr); %coefficient de réflexion complexe
gr=abs(Gr);
alphar=1-gr.^2; %coefficient d'absorption
Zr=(1+Gr)./(1-Gr); %impédance de paroi


%%%%%%%%%%%%%%%%%%%%%
% TRACES
%%%%%%%%%%%%%%%%%%%%%
figure(1);
%hold on;
plot(fq,alphar);
title('Tube à ondes stationnaires - méthode des deux microphones');
xlabel('fréquence (Hz)');
ylabel('coefficient d"absorption');
axis([0 max(fq) 0 1]);
%igure(2);
%plot(fq,abs(h1),fq,abs(h2));