Matlaby Dźwięk

clear all;
close all;

[Y,FS] = audioread('glos.wav');

subplot(2,1,1);
plot(Y,'r');
title('probki');


%sound(Y,FS);

N = length(Y);
X = (0:N-1)/FS;
subplot(2,1,2);
plot(X,Y,'g');
title('czas');


time = 10 / 100; %time in ms

frame = FS * time;

allEnergy = zeros(N/frame,1);
allZeroJump = zeros(N/frame,1);


for j = 0:length(allEnergy)-1
    energy = 0;
    countZeroJumps = 0;
    for i = 1:frame
        index = i +( j* frame);
        energy = energy + (Y(index) ^ 2);
        if index + 1 <= N
            if zeroJump(Y(index),Y(index+1)) == true
                countZeroJumps = countZeroJumps + 1;
            end
        end
    end
    allEnergy(j+1) = energy;
    allZeroJump(j+1) = countZeroJumps;
end

allEnergy = allEnergy/max(allEnergy);
allZeroJump = allZeroJump/max(allZeroJump);


figure;
hold on;
plot(Y,'r');
plot((0:frame:N-1),allEnergy,'g');
plot((0:frame:N-1),allZeroJump,'b');
hold off;