%Authors: Dylan Rumph, Jake Poteet, Brian Viel
%Title: SmokeontheWater.m
%Date 9/11/09


%sawtooth(t,.5)=triangle wave
%sawtooth(t,1)=sawtooth
%square(t,duty) = square wave, duty 0:100; 50=normal
%sin(t) = sin wave with a period
%soundsc(y,Fs) ; Fs is the sampling frequency

clear;
a=square(2*pi*220*(0:0.000125:0.5), 50);
as=sin(2*pi*233.08*(0:0.000125:0.25));
e=sawtooth(2*pi*164.81*(0:0.000125:0.5), 1);
g=sawtooth(2*pi*196*(0:0.000125:0.5), .5);
gm=sawtooth(2*pi*196*(0:0.000125:0.25), .5);
la=square(pi*220*(0:0.000125:0.5), 50);
las=sin(pi*233.08*(0:0.000125:0.25));
le=sawtooth(pi*164.81*(0:0.000125:0.5), 1);
lg=sawtooth(pi*196*(0:0.000125:0.5), .5);
lgm=sawtooth(pi*196*(0:0.000125:0.25), .5);
n = randn(1, 2000);

line1=[e, g, a, n, e, gm, as, a, n];
line2=[e, g, a, n, g, e, n, n];

line1l=[le, lg, la, n, le, lgm, las, la, n];
line2l=[le, lg, la, n, lg, le, n, n];

Echer = -.5 * (square(2* pi* 3* (0:0.000125:5))+ .5);
Expo = exp(-1*(0:0.000125:5));
Echothing = Echer .* Expo;

song=[line1,line2];
songl=[line1l, line2l];
song2=[song, songl];
songe = conv(song2, Echothing);

soundsc(song2);
soundsc(songe);