Matlab - einfache Fouriertransformation eines Sägezahns

[M & m &. (hanno85)]

Hallo!

Ich habe ein kleines Problem mit meinem Matlab Programm, es zeigt mir 
andeutungsweise einen Sägezahn, wenn ich die harmonischen Schwingen 
gemäß der Fourierformel für fallende Sägezahmsignale anwende.

hier mein Programm:

function [y,t] = fourier(ua,N,f,ap)

T=1/f;
y=0;
y2=0;
t=0:1/(f*f):T*ap;

for n = 1:N
    y = y + [(2*ua / (pi)) * (1/n) * sin(2*pi*f*n*t)/n ];
end    
y2=(-2*ua .*rem(t,T) ./T) + ua;

figure(1)
plot(t,y)
grid on
xlabel('t in s')
ylabel('A in V')
title('\it{Fouriertransformierte}','FontSize',16)
figure(2)
plot(t,y2)
grid on
xlabel('t in s')
ylabel('A in V')
title('\it{Sägezahnsignal}','FontSize',16)

[joep (Gast)]

Hab leider gerad kein Matlab da, kann mir  also dein Ergebnis nicht 
anschaun.
Aber da du ja meinst, dass er andeutungsweise einen Sägezahn anzeigt: 
hast du N groß genug gewählt?

[M & m &. (hanno85)]

Ja es ist leider kein richtiger Sägezahn, der ist an den stellen, wo er 
gerade nach oben gehen soll noch ein sinus zu erkennen und egal wie hoch 
ich die Fourier-Reihe berechne es ändert sich nichts daran : /

[Andreas S. (andreas) (Admin)]

Du hast das Gibbssche Phänomen entdeckt: 
http://de.wikipedia.org/wiki/Gibbssches_Ph%C3%A4nomen

[M & m &. (hanno85)]

ah Fehler behoben hatte ein 1/n zuviel in der Funktion! Typischer 
Anfängerfehler : /

Hier nochma die richtige Funktion, diesma hab ich die mit Vektoren 
berechnet, um die For-Schleife zu umgehen. Jedoch habe ich noch keine 
Idee, wie ich die einzelnen harmonischen Schwingungen ohne For-Schleife 
darstellen kann? Hat jmd eine Idee, am besten auch mit Vektor und 
Matrizenrechnung.

function [y,t] = fourier(ua,N,f,ap)

T=1/f;
w=2*pi*f;
n=1:N;
t=0:1/(f*f):T*ap;
y=0;
y2=0;
y3=0;

y = (2*ua / (pi)) * (1./n) * sin(w*n'*t); 

for n = 1:N
    y3 = y3 + [(2*ua / (pi)) * (1/n) * sin(2*pi*f*n*t)];
    figure(3)
    plot(t,y3)
    grid on
    hold on
    xlabel('t in s')
    ylabel('A in V')
    title('\it{harmonische Sinusschwingungen}','FontSize',16)

end  

y2=(-2*ua .*rem(t,T) ./T) + ua;
figure(1)
plot(t,y)
grid on
xlabel('t in s')
ylabel('A in V')
title('\it{Fouriertransformierte}','FontSize',16)
figure(2)
plot(t,y2)
grid on
xlabel('t in s')
ylabel('A in V')
title('\it{Sägezahnsignal}','FontSize',16)