harj34.m

% Analysoidaan signaalin taajuus käyttämällä positiivisesta
% negatiiviseksi transition ajanhetkiä

x = load( 'data34.txt' );
hold off
plot( x(:,1), x(:,2), 'b-' )

% Ensimmäinen transitio
for i=1:size(x,1)-1
    if x(i,2) > 0 && x(i+1,2) < 0
        t1 = x(i,1);
        istart = i;
        break;
    end
end
    
% Toinen transitio
for i=istart+1:size(x,1)-1
    if x(i,2) > 0 && x(i+1,2) < 0
        t2 = x(i,1);
        break;
    end
end

% Merkitään triggeripisteet kuvaan
hold on
plot([t1 t1], [-3 3], 'r-')
plot([t2 t2], [-3 3], 'r-') 

T = t2-t1;
freq = 1/T