% trois gammes : Temperee, Pythagore et Zarlino
%----------------------------------------------

clear all; close all;

Fe = 22050; h = 1/Fe;
f0 = 110*2^(1/4); t = .3;

% gamme temperee
%---------------
fr = f0*2.^([0 2 4 5 7 9 11 12]/12);
x = [];
for k = 1:8
    x = [x, orgue(fr(k),t,Fe)];
end
x = .6*x;
sound(x,Fe);
pause(3)
%wavwrite(x,Fe,'../SONS/301gamtemp.wav');

% gamme de Pythagore
%-------------------
fr = f0*[1 9/8 81/64 4/3 3/2 27/16 243/128 2];
x = [];
for k = 1:8
    x = [x, orgue(fr(k),t,Fe)];
end
x = .6*x;
sound(x,Fe);
pause(3)
%wavwrite(x,Fe,'../SONS/302gampyth.wav');

% gamme de Zarlino
%-----------------
fr = f0*[1 9/8 5/4 4/3 3/2 5/3 15/8 2];
x = [];
for k = 1:8
    x = [x, orgue(fr(k),t,Fe)];
end
x = .6*x;
sound(x,Fe);
%wavwrite(x,Fe,'../SONS/303gamzarl.wav');