Verfasst am: 15.05.2009, 13:07
Titel: Plotten einer Fourierreihe
Hallo ich habe folgendes Problem. Ich habe ein Rechtecksignal abgetastet und einen Vektor der Länge N erhalten. Dann habe ich diesen mit der komplexen Fouriermatrix multipliziert, um den Fouriertransformierten Vektor zu erhalten. jetzt weiß ich ja das a_0=2*c_0 ist und a_k=2*Re(c_k) und b_k=-2*Im(c_k). Das habe ich auch noch hinbekommen! jetzt wollte ich daraus mir die Fourierreihe in abhängigkeit von N berechnen und anschließend Plotten! Bloss da scheitere ich dran schon seit mehren Stunden. Könnte mir einer weiter helfen, wie ich die Fourierreihe aufstelle? Oder vll auch mal gucken ob irgendwas falsch ist in meinen Programm?!
Code:
clear
%Es Wir ein Rechtecksiganl eingelesen wobei die Länge abgefragt wird
K=input('Geben Sie die halbe Länge einer Periode des Rechtecksignals ein: ');
%Eingabe des Vektors
N=input('Geben Sie die Länge des abtastenden Vektors ein: ');
%Vektorelemente auf Null setzen
a=linspace(0,0,N);
%Einlesen des Vektors
a(1)=0.5;
for k=2:N/2
a(k)=1;
end;
a(N/2+1)=0;
for k=N/2+2:N
a(k)=-1;
end;
a(N+1)=-0.5;
%Plotten des Abtastvektors und des Rechtecksignals
plot([00 K K 2*K 2*K],[011-1-10],0:2*K/N:2*K,a,'r.'),axis([-K/42*K+K/4-1.51.5]);
%Transponieren des Vektors a
b=zeros(N,1);
b=a';
%Erzeugung der komplexen Fouriermatrix
W=zeros(N);
l=1;
while l<N+2 for k=1:N+1
W(k,l)=exp(i*(k*l*2*pi)/N);
end;
l=l+1;
end;
%Korrektur der Matrix
W(1,:)=1;
W(:,1)=1;
%Fouriermatrix-Vektor Multiplikation
c=zeros(N,1);
for k=1:N+1
c(k)=0;
for l=1:N+1
c(k)=c(k)+W(k,l)*a(l);
end;
end;
%Überführung in die reellen Zahlen
d=zeros(N+1,1);
e=zeros(N,1);
f=zeros(N,1);
d(1)=2*c(1);
for k=1:N+1
e(k)=2*real(c(k));
f(k)=-2*imag(c(k));
end;[\code]
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.
RSS