% Aufgabe 4.1.4: Zeitverlauf und Spektrum anzeigen
%    spektrumSign2() 
function spektrumSign2()

   % Zeitintervall von 0 bis 2 s
   % Array mit Stützpunkten im Abstand 1/40000 s
   t = 0 : 1/40000 : 2;
   f   = 220;      % Grundfrequenz 220 Hz

   % Array y mit Funktionswerten an den t-Stützpunkten
   y = sign( sin(2*pi*f*t) ) - sin(2*pi*f*t);
   
   % Ausgabe der Funktion als 2D-Plot
   subplot(2,1,1), 
   plot( t, y );
   axis( [0 1/20 -2 2] ) % fester Ausgabebereich

   % Fourier-Transformation
   c = fft( y );
   % Betrag der Frequenzamplituden  
   A = abs( c );

   % Array der Frequenzen, normiert auf 40000 Hz
   f = (0:length(c)-1).*40000/length(c);

   % 2D-Plot der Frequenzamplituden über den Frequenzen
   subplot(2,1,2)
   plot( f, A ) 
   % Frequenzbereich: 0 Hz bis 4000 Hz
   axis( [0 8000 -1 20000] )

   % Abspielen
   sound( y, 40000 );
end