% Eecoute de battements
%-------------------------------------------------------------

clear all
Fe = 22050;

%-------------------------------------------------------------
% deux notes proches de l'unisson
%-------------------------------------------------------------
h = 1:3;
f1 = 330; f2 = 332; 
fr1 = h*f1; fr2 = h*f2;
a = [1000 1 .1];
p = 0*a; T = 2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe);
disp('unisson desaccorde');
soundsc(s,Fe);
pause(T+1);

%-------------------------------------------------------------
% deux notes proches d'une quinte, avec puis sans h2 et h3
%-------------------------------------------------------------
h = 1:4;
f1 = 330; f2 = 330*1.5+1; 
fr1 = h*f1; fr2 = h*f2;
a = [100 50 20 2];
p = 0*a; T = 2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe);
disp('une quinte descaccordee');
soundsc(s,Fe);
pause(T+1);

h = 1:4;
f1 = 330; f2 = 330*1.5+1; 
fr1 = h*f1; fr2 = h*f2;
a = [100 0 0 2];
p = 0*a; T = 2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe);
disp('idem sans les harmoniques 2 et 3');
soundsc(s,Fe);
pause(T+1);


%-------------------------------------------------------------
% tierces majeures
%-------------------------------------------------------------

f1 = 220*3/4;
h = 1:5;
a = [100 30 20 15 10];
p = 0*a; T = 2;

% tierce majeure Zarlino
f2 = f1*5/4; 
fr1 = h*f1; fr2 = h*f2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe);
disp('tierce Zarlino');
soundsc(s,Fe);
pause(T+1);

% tierce majeure temperee
f2 = f1*2^(4/12); 
fr1 = h*f1; fr2 = h*f2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe);
disp('tierce temperee');
soundsc(s,Fe);
pause(T+1);

% tierce majeure Pythagore
f2 = f1*81/64; 
fr1 = h*f1; fr2 = h*f2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe);
disp('tierce Pythagore');
soundsc(s,Fe);
pause(T+1);

%-------------------------------------------------------------
% accords majeurs
%-------------------------------------------------------------

f1 = 220;
h = 1:5;
a = [100 30 20 15 10];
p = 0*a; T = 2;

% tierce majeure Zarlino
f2 = f1*5/4; 
fr1 = h*f1; fr2 = h*f2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe) ...
    + synthad(a,1.5*fr1,p,T,Fe);
disp('accord majeur Zarlino');
soundsc(s,Fe);
pause(T+1);

% tierce majeure temperee
f2 = f1*2^(4/12); 
fr1 = h*f1; fr2 = h*f2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe) ...
    + synthad(a,1.5*fr1,p,T,Fe);
disp('accord majeur tempere');
soundsc(s,Fe);
pause(T+1);

% tierce majeure Pythagore
f2 = f1*81/64; 
fr1 = h*f1; fr2 = h*f2;
s = synthad(a,fr1,p,T,Fe) + synthad(a,fr2,p,T,Fe) ...
    + synthad(a,1.5*fr1,p,T,Fe);
disp('accord majeur Pythagore');
soundsc(s,Fe);
pause(T+1);
