% verze 1.0, Jan Chvojka
clear, close all
fname = 'violin-g-4-non-vibrato.wav';
[x,fs] = audioread(fname); 
x = x(:,1);

% vibrato
b0 = 0;  bd = 1;
max_zpozdeni = 0.0005;
fmd = 4; D = round(max_zpozdeni*fs);

% flanger
b0 = 0.5; bd = 0.5;
max_zpozdeni = 0.005;
fmd = 0.1; D = round(max_zpozdeni*fs);


y = zeros(size(x));
for n = 1:length(x)   
    d = round(D*cos(2*pi*fmd*n/fs));
    if (n-d)>0 && (n-d)<=length(x)
        y(n) = b0*x(n) + bd*x(n-d);
    end   
end

soundsc(x,fs);
pause(3.5);
soundsc(y,fs);