Verfasst am: 18.05.2009, 21:33
Titel: Eine Fourier-Reihe aufstellen
Hey.
Ich hab folgendes Problem. Zu einen vorgegebenen Vektor der Länge N eines Signals,z.b. Rechtecksignal oder Sägenahnsignal, möchte ich eine reelle Fourierreihe aufstellen, also eine Funktion aus linearkombinationen trigometrischer Funktionen. Ich habe dafür folgendes im Netz gefunden(siehe Anhang). Jetzt hab ich aber Probleme das in Matlab zu verarbeiten! Ich denke es liegt irgendwo ein Fehler bei der Indizierung, weil Matlab ja nur mit positiven Integern arbeitet und, aber bei den Fourier-Reihen es ja immer bei null anfängt! Folgendes habe ich mir überlegt, aber es sieht beim Plotten nicht nach der Sägezahnfunkion aus:
Code:
clear
%%
%Es Wir ein Sägezahnsignal eingelesen wobei die Länge abgefragt wird
K=input('Geben Sie die halbe Länge einer Periode des Sägezahnsignals ein: ');
%Eingabe des Vektors
N=input('Geben Sie die Länge des abtastenden Vektors ein: ');
for k=1:N-1
b(k)=(k-1)*1/(N-2);
end;
b(N)=0.5;
plot([02*K 2*K 4*K],[0101],'g',0:2*K/(N-1):2*K,b,'ko');
A=0;
B=0;
C=0;
%Zeichnen
for k=1:N/2+1
A(k)=0;
B(k)=0;
C(1)=0;
for l=1:N
A(k)=A(k)+(1/N)*b(l)*cos((k*l*2*pi)/N);
B(k)=B(k)+(1/N)*b(l)*sin((k*l*2*pi)/N);
C(1)=(1/N)*b(l);
end;
end;
y=0;
x=linspace(0,4,1001);
for m=2:N/2
y=0;
y=C(1)+2*(y+A(m)*cos(((m-1)*2*pi*x)/(2*K))+B(m)*sin(((m-1)*2*pi*x)/(2*K)))+A(N/2)*cos((N*pi*x)/(2*K));
end;
hold on,plot(x,y);
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.