Fourier-Reihe

Fourier-Reihe

Fourier_Reihe.m — HTML, 625 bytes

Dateiinhalt

%
% Fourier-Zerlegung eines Rechtecksignals mit 0.2 Hz
% 
t = [0:1/100:10];       % 100 Hz Abtastrate für 10 Sekunden
t1 = [3.75:1/100:13.75];
d = 0:1/0.2:10;         % 0.2 Hz: Frequenz des Rechtecksignals
%
% Berechnung des Rechtecksignals (Signal Processing Toolbox)
%
y1 = 2*pulstran(t1,d,'rectpuls',2.5)-1;
plot(t,y1,'m')
hold on
%
% Berechnung der Fourier-Reihe
%
w_0 = 0.2*2*pi;
y = 4/pi*sin(w_0*t);
for i=3:2:511
    y=y+4/pi/i*sin(i*w_0*t);
end
plot(t,y)
grid on
xlabel('t [s]')
ylabel('f(t)')
title('f(t) = 4/\pi *[sin(\omega_0t) + 1/3sin(3\omega_0t) + 1/5sin(5\omega_0t) + ... ]')
hold off