clear; min = 10; max = 20; rate = 1/(max*4*1000); result = zeros(0, 2^16); for n = min:max noise = sin(2*pi*n*1000*(0:rate:2)); result = [result noise]; end wavwrite(result, max*4*1000, 16, '10khz_20khz_1sec_sample.wav');