hwdp

close all
clc

my_path = what;
my_file=dir([my_path.path '/sounds/']);
pathname_pliki=([my_path.path '/sounds/']);

L = (length(my_file)-2);

%F = zeros(1:(length(my_file)-2));
%J = zeros(1:(length(my_file)-2));
%S = zeros(1:(length(my_file)-2));


for i=1:(length(my_file)-2)
    filename=my_file(i+2).name;
    name = [pathname_pliki filename];
    [s,fs]= audioread(name);

    [Jitt, Shim, Fn0] = jittshim(s, fs);

    if i < 108
        F(i) = Fn0;
        J(i) = Jitt;
        S(i) = Shim;
    else
       Fx(i) = Fn0;
       Jx(i) = Jitt;
       Sx(i) = Shim;
    end
end


figure(1)
plot(J(:), S(:),'.');
hold on;
plot(Jx(:), Sx(:), '*');
xlabel('jitter[%]');
ylabel('shimmer[%]');
axis([0 10 0 10]);
grid on;
hold off;


figure(2)
plot3(J(:), S(:), F(:), '.');
hold on;
plot3(Jx(:), Sx(:), Fx(:), '*');
grid on;
xlabel('jitter[%]');
ylabel('shimmer[%]');
zlabel('freqs [Hz]');
hold off


symbols = ['*', 's'];