matthiasm@104: f = 50:800;
matthiasm@104: fs = 44100;
matthiasm@104: t = (0:(2*fs))/fs;
matthiasm@104: nFreq = length(f);
matthiasm@104: 
matthiasm@104: for iFreq = 1:nFreq
matthiasm@104:     f(iFreq)
matthiasm@104:     x = zeros(size(t));
matthiasm@104:     for iAdd = 1:100
matthiasm@104:         x = x + cos(2*pi*f(iFreq)*t*iAdd)*0.6^(iAdd-1);
matthiasm@104:         if f(iFreq)*iAdd > fs/2
matthiasm@104:             break
matthiasm@104:         end
matthiasm@104:     end
matthiasm@104:     x = x / max(abs(x));
matthiasm@104:     wavwrite(x, fs, sprintf('/Users/matthiasm/data/pyin/simplewavs/%iHz.wav', f(iFreq)));
matthiasm@104: end