function y = synth_complex1 % y=synth_complex1 - sum of sines global sr; % sampling rate D = 0.2; % s - duration f0 = 200; % Hz - fundamental frequency np = 10; % number of partials x = abscissa(D); y = x .* 0; for i = 1:np y = y + sin(x * i * f0 * 2 * pi); end