Fourierreihen und Fouriertransformation - Fachhochschule ...

Fourierreihen und Fouriertransformation 72
end
a(n/2)=real(fy(n/2+1))/n;
an=a'
% Bezeichnungen Literatur Index mit 0 beginnend
for i=1:n/2-1
b(i)=-2/n*imag(fy(i+1));
end
bn=b'
% Plot des interpolierten Signals
% t-Intervall ist [0,2N]=[0,n]
% Abtastzeitpunkte tk=k/(2N) , k = 0,..., 2N-1
t=-1:0.05:n+1; % Bereich zum plotten
yf=a0;
for i=1:n/2-1
yf=yf+a(i)*cos(i*2*pi/n*t)+b(i)*sin(i*2*pi/n*t);
end
yf=yf+a(n/2)*cos(pi*t);
plot(t,yf,'LineWidth',1.5);
pause
% Falls Periode der Abtastzeitpunkte = 2*pi
% t- Intervall ist [0,2pi]
% Abtastzeitpunkte tk=k*pi/N , k = 0,..., 2N-1
t=[0:7]*2*pi/n;
subplot(3,1,2)
plot(t,y,'o','MarkerFaceColor','k');
title(['Periode = 2pi']);
axis([-pi/2 5/2*pi min(y)-1 max(y)+1])
grid on
% Plot des interpolierten Signals
t=-1:0.05:2*pi+1; % Bereich zum plotten
yf=a0;
for i=1:n/2-1
yf=yf+a(i)*cos(i*t)+b(i)*sin(i*t);
end
yf=yf+a(n/2)*cos(n/2*t);
hold on
plot(t,yf,'LineWidth',1.5);
pause
close
Output des obigen Matlab-Programms:
>> dfi3_Dreieck
a0 =
4
an =