% 
%  Periodisches Signal aus Grundton und 3 Obertönen 
% 
function signalTon()
   % Zeitintervall von 0 bis 2 s, Sample-Rate 1/40000 s
   t = 0 : 1/40000 : 2;
   f   = 220;      % Grundfrequenz 220 Hz
   A   = 1;        % 1 als Maximalamplitude

   % Array y mit Funktionswerten an den t-Stützpunkten
   y = A * ( sin(2*pi*f*t)   + 1/3 * sin(2*pi*3*f*t) + ... 
       1/5 * sin(2*pi*5*f*t) + 1/7 * sin(2*pi*7*f*t) );

   % Ausgabe der Funktion als 2D-Plot
   plot( t, y );
   axis( [0 0.02 -1.2 1.2] ) % fester Ausgabebereich

