% verze 1.0, Jan Chvojka clear, close all fname = 'violin-g-4-non-vibrato.wav'; [x,fs] = audioread(fname); x = x(:,1); % chorus b0 = 1; K = 50; % pocet zpozdenych kopií bd = rand(1,K)*0.5+0.5; max_zpozdeni = (rand(1,K)*30+10)*0.001; D = round(max_zpozdeni*fs); f = rand(1,K)*5; ph = rand(1,K)*2*pi; y = zeros(size(x)); for n = 1:length(x) y(n) = b0*x(n); for i = 1:length(bd) d = round(20*cos(2*pi*f(i)*n/fs+ph(i))) + D(i); if (n-d)>0 && (n-d)<=length(x) y(n) = y(n) + bd(i)*x(n-d); end end end soundsc(x,fs); pause(3.5); soundsc(y,fs);