Fourierreihen und Fouriertransformation - Fachhochschule ...

Fourierreihen und Fouriertransformation 62
Ein möglicher zweiter Matlab-Code mit Figurenoutput:
X=v:0.01:v+2*pi;
Y=[ones(size(X')) cos(X') cos(2*X') cos(3*X') cos(4*X')
sin(X') sin(2*X') sin(3*X')]*FK;
hold on
plot(x,y,'ko','MarkerFaceColor','k')
plot(X,Y','LineWidth',1.5)
grid on
title(['Funktion f(t) mit ', num2str(n) , ' Abtastpunkten']);
pause
close
Abbildung 38: Interpolierendes Fourierpolynom
des diskret abgetasteten Rechtecksignals.
Und schliesslich das Amplitudenspektrum , d. h. das Zeichnen der Amplitudenwerte
2 2
= α + β
als senkrechte Balken an den Frequenzstellen k ω0
.
% Amplitudenspektrum
FK=FK';
a0=FK(1);
a=FK(2:n/2+1);
b=[FK(n/2+2:n) 0];
ampl=[a0 sqrt(a.^2+b.^2)];
hold on
for k=0:n/2
plot([k k],[0 ampl(k+1)],'k','LineWidth',2)
end
plot([-1 n/2+1],[0 0],'k','LineWidth',1.5)
grid on
title(['Amplitudenspektrum von f(t)=' ,func ,' bei
',num2str(n), ' Abtastpunkten']);
axis([-1 n/2+1 -0.1*max(ampl) 1.1*max(ampl)]);
Ak
k
k