%
%  SDSMT Signals Class EE312
%  Dr. Charles R. Tolle
%
%  The usefulness of the F.T. by exploring musical -- tones -- harmonics
%
%  cite equal-tempered scale
%  http://www.phy.mtu.edu/~suits/notefreqs.html
%
%  cite for clarinet tones:
%  http://www.phys.unsw.edu.au/jw/flutes.v.clarinets.html
%
%  cite for cords:
%  http://nfo.net/MFILE/
%

clear; home;
C0      =   16.35;  
Cs0_Db0 =   17.32;
D0      =   18.35;
Ds0_Eb0 =   19.45;
E0      =   20.60;
F0      =   21.83;
Fs0_Gb0 =   23.12;
G0      =   24.50;
Gs0_Ab0 =   25.96;
A0      =   27.50;
As0_Bb0 =   29.14;
B0      =   30.87;
C1      =   32.70;
Cs1_Db1 =   34.65;
D1      =   36.71;
Ds1_Eb1 =   38.89;
E1      =   41.20;
F1      =   43.65;
Fs1_Gb1 =   46.25;
G1      =   49.00;
Gs1_Ab1 =   51.91;
A1      =   55.00;
As1_Bb1 =   58.27;
B1      =   61.74;
C2      =   65.41;
Cs2_Db2 =   69.30;
D2      =   73.42;
Ds2_Eb2 =   77.78;
E2      =   82.41;
F2      =   87.31;
Fs2_Gb2 =   92.50;
G2      =   98.00;
Gs2_Ab2 =   103.83;
A2      =   110.00;
As2_Bb2 =   116.54;
B2      =   123.47;
C3      =   130.81;
Cs3_Db3 =   138.59;
D3      =   146.83;
Ds3_Eb3 =   155.56;
E3      =   164.81;
F3      =   174.61;
Fs3_Gb3 =   185.00;
G3      =   196.00;
Gs3_Ab3 =   207.65;
A3      =   220.00;
As3_Bb3 =   233.08;
B3      =   246.94;
C4      =   261.63;
Cs4_Db4 =   277.18;
D4      =   293.66;
Ds4_Eb4 =   311.13;
E4      =   329.63;
F4      =   349.23;
Fs4_Gb4 =   369.99;
G4      =   392.00;
Gs4_Ab4 =   415.30;
A4      =   440.00;
As4_Bb4 =   466.16;
B4      =   493.88;
C5      =   523.25;
Cs5_Db5 =   554.37;
D5      =   587.33;
Ds5_Eb5 =   622.25;
E5      =   659.26;
F5      =   698.46;
Fs5_Gb5 =   739.99;
G5      =   783.99;
Gs5_Ab5 =   830.61;
A5      =   880.00;
As5_Bb5 =   932.33;
B5      =   987.77;
C6      =   1046.50;
Cs6_Db6 =   1108.73;
D6      =   1174.66;
Ds6_Eb6 =   1244.51;
E6      =   1318.51;
F6      =   1396.91;
Fs6_Gb6 =   1479.98;
G6      =   1567.98;
Gs6_Ab6 =   1661.22;
A6      =   1760.00;
As6_Bb6 =   1864.66;
B6      =   1975.53;
C7      =   2093.00;
Cs7_Db7 =   2217.46;
D7      =   2349.32;
Ds7_Eb7 =   2489.02;
E7      =   2637.02;
F7      =   2793.83;
Fs7_Gb7 =   2959.96;
G7      =   3135.96;
Gs7_Ab7 =   3322.44;
A7      =   3520.00;
As7_Bb7 =   3729.31;
B7      =   3951.07;
C8      =   4186.01;
Cs8_Db8 =   4434.92;
D8      =   4698.64;
Ds8_Eb8 =   4978.03;

%flute
%h(1) = 0;
%h(2) = -18;
%h(3) = -12;
%h(4) = -18;
%h(5) = -20;
%h(6) = -35;

%Clarinet
h(1) = 0;
h(2) = -15;
h(3) = -22;
h(4) = -35;
h(5) = -35;
h(6) = -38;

%for i = 1:6,
% h(i) = 1/(1+i);
%end
%h(1) = 1;

for i = 1:6,
    mh(i) = 10^(h(i)/20);
end

% middle C --> C5
clear y;
l = 3;
t = linspace(0,l,l*8192);
disp('Pure C5 note'); 
y = sin(2*pi*C5*t);
sound(y)

pause();
disp('Carinet based C5 Harmonic note'); 
for i = 2:6,
   y = y + mh(i)*sin(2*pi*i*C5*t);
end
y = y - mean(y);
y = y ./ max(abs(y));
sound(y)
pause();

clear t;
l = 3;
t = linspace(0,l,l*8192);

% cord C5 E5 G5
clear y;
disp('Play a pure C5 E5 G5 cord');
y = sin(2*pi*C5*t) + sin(2*pi*G5*t) + sin(2*pi*E5*t);
y = y - mean(y);
y = y ./ max(abs(y));
sound(y)
pause();

% cord C5 E5 G5
clear y;
y = mh(1)*sin(2*pi*C5*t);
for i = 2:6,
   y = y + mh(i)*sin(2*pi*i*F5*t);
end

for i = 1:6,
   y = y + mh(i)*sin(2*pi*i*G5*t);
end

for i = 1:6,
   y = y + mh(i)*sin(2*pi*i*C5*t);
end

disp('Carinet based C5 E5 G5 Harmonic cord'); 
y = y - mean(y);
y = y ./ max(abs(y));
sound(y)