% create scale

clear all;

Fs = 44100;
dur = .2;
t = [0 : 1/Fs : dur-1/Fs];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% general formula for n-split equitempered octave
% f(m) = sin ( 2*pi * freq * 2^(m/n) * t)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

f0 = 400;   % fundamental frequency
n_split = 12;   % number of equidistant tones that octave is split into
n = [0:n_split];
for m = 1:length(n)
	f(m,:) = sin ( 2*pi * f0 * 2^(n(m)/n_split) * t );
end

y = [];  % prepare sound vector
for k = 1:13
	y = [y f(k,:)];
    y = wind(Fs,20,y);
end
y = .4 * y;
sound(y,Fs)

t_all = [0 : 1/Fs : dur*length(n) - 1/Fs];
plot(t_all,y)
axis([0 length(n)*dur min(y)-.1 max(y)+.1])

Published with MATLAB® 7.1