% 
%  Signal aus MAT-File einlesen + FFT  
%    >> readSignalFFT( 'piano.mat', 0.05 ) 
%
function readSignalFFT( filename, s )
    % MAT-File einlesen 
    data = load( filename );
    t = data.t;   
    y = data.y; 
    
    % Anzahl der Daten
    len = length( t );
    tmax = t(len);   % Zeitdauer 
    SR = len/tmax;
   
    % Ausgabe der Funktion als 2D-plot
    subplot(2,1,1), 
    plot( t, y );
    axis( [0 s -1 1] ) % s sec plotten 
   
    % Fourier-Transformation
    c = fft( y );
    
    % Betrag der Amplituden  
    A = abs( c );
    
    % Array der Frequenzen, normiert auf SR Hz
    anzc = length(c);
    f = (0:anzc-1) * SR/anzc; 
    maxA = max( A ) * 1.1;
    
    % Ausgabe der Amplituden zu den Frequenzen als 2D-plot
    subplot(2,1,2)
    plot(f,A) 
    
    % Frequenzen bis max. 4000 Hz
    axis( [0 4000 -1 maxA] )
    

