% dft.m
% Berechnung und grafische Darstellung der DFT 
% 26.3.2002
% Rolf Hofstetter, Thomas Borer



clear all												%alle Variablen löschen
close all												%alle Figuren-Windows schliessen
Ta=0.0001;												%Abtastperiode
t=0:Ta:0.03;											%Zeitvektor
fa=1/Ta;													%Abtastfrequenz
w0=2*pi*100;											%Grund-Kreisfrequenz
s=2*sin(4*w0*t)+5*sin(7*w0*t);					%Signal mit gleicher Länge wie der Zeitvektor t
N=200;													%Fensterlänge bzw. Anzahl Abtastwerte im Zeitfenster
x=s(1:N);												%Signal innerhalb des gewählten Zeitfensters
subplot(2,1,1);plot(t(1:N),x);grid; 			%Darstellung des Signals in der oberen Figur
fax=abs(fft(x));										%Bestimmen der DFT-Absolutwerte von x
f=0:(fa/N):fa-(fa/N);								%Frequenz-Vektor
subplot(2,1,2);stem(f(1:50),fax(1:50));grid; %Darstellung der DFT-Absolutwerte in der unteren Figur


