Fourierreihen und Fouriertransformation - Fachhochschule ...

Fourierreihen und Fouriertransformation 79
fy=fy(1:kk);
kkk=min(kk,n/2);
T=1/fa;
% Berechnung der reellen Koeffizienten
a0=fy(1)/n;
for i=1:kkk-1
a(i)=2/n*real(fy(i+1));
end
for i=1:kkk-1
b(i)=-2/n*imag(fy(i+1));
end
te=0:2*pi/n:10*T; % Bereich zum plotten
subplot(2,1,1)
plot(te,x(1:length(te)))
title('Originalsignal');
axis([0 10*T -1.1*max(x) 1.1*max(x)]);
pause
hold on
% Plot des interpolierenden trig. Polynoms
yf=a0;
for i=1:kkk-1
yf=yf+a(i)*cos(i*te)+b(i)*sin(i*te);
end
subplot(2,1,2)
plot(te,yf)
title('gefiltertes Signal');
axis([0 10*T -1.1*max(x) 1.1*max(x)]);
pause
close
Abbildung 45: Die beiden Outputfenster des obigen Matlab-Codes.