Hilfe zu Fourier-Transformation

Gast #798057
Lesenswert?

Ich hab da ein Problem mit einer Aufgabe. Das Signal x(k) befindet sich 
im Dateianhang.

Bestimmen Sie die Fourier-Transformierte X(ω) von x(k).
Die Abtastfrequenz sei fs = 1 kHz. Zeichnen Sie die 
Fourier-Transformierte im Frequenzbereich 0 ≤ f ≤ 500Hz.
Angehängte Dateien:
Gast #798614
Lesenswert?

Wenn ich statt ein Dreieck ein Rechteck transformiere,dann erhalte ich 
folgende Werte:

n=0  P(n*wo)=9
n=1  P(n*wo)=0
n=2  P(n*wo)=0
n=3  P(n*wo)=9
n=4  P(n*wo)=0

[P(w)=3*cos(2*w*T)+3*cos(2*w*T)+3]

Ist dies eigentlich richtig?
Gast #798622
Lesenswert?

[P(w)=3*cos(2*w*T)+3*cos(2*w*T)+3]

Die Gleichung stimmt nicht. Hier die korrigierte Version:

[P(w)=3*2*cos(2*w*T)+3*2*cos(2*w*T)+3]

Auch die Ergebnisse sind falsch. Hier sind die richtigen Ergebnisse:

n=0  P(n*wo)=15
n=1  P(n*wo)=3
n=2  P(n*wo)=3
n=3  P(n*wo)=3
n=4  P(n*wo)=3

Mit Matlab hab ich dies mal überprüft:

abs(fft([3,3,3,3,3],6))

ans =

    15     3     3     3     3     3
#799840
Lesenswert?

Wo ist dein Problem? Hier werden lediglich einige Regeln der ZDFT 
benutzt.

x[k] nimmst du als Dreieck von -3 bis 3 an (das wir mal p[k] nennen), 
das mittels Faltung mit Diracs um 6 verschoben wird.
Diesen Ausdruck kann man ohne weiteres direkt transformieren.

Faltung wird zur Multiplikation im Frequenzbereich.

Gruß,
Stefan
Gast #799878
Lesenswert?

Heftig, ja? Was wäre, wenn du dein Skript vernünftig lesen würdest?
Du hast ein diskretes, periodisches Signal. Dessen Spektrum ist 
dementsprechend diskret und periodisch, dass heißt, dass die richtige 
Fourier-Transformation die DFS ist. Diese kannst du mittels DFT 
berechnen.
Die Länge einer Periode ist 6, das heißt:

Berechnen muss man nur X(0) .. X(3):

f=0: X(0) = .3 + .2 + .1 + 0 + .1 + .2 = .9
f=1000*2/6: X(1) = .3*exp(-j*2*pi*0*1/6) + .2*exp(-j*2*pi*1*1/6) + ... = 
.4
X(2) = 0
X(3) = .1

Prüfung:

>> fft([.3 .2 .1 0 .1 .2])

ans =

    0.9000    0.4000   -0.0000    0.1000   -0.0000    0.4000
Gast #800042
Lesenswert?

Hallo Stefan, deine angegebene Formel ist doch falsch.
X(\Omega) müsste doch so aussehen oder?

x[k] = p[k] * \sum^{\infty}_{k=-\infty} \delta[n-k \cdot 6]

X(\Omega) = P(\Omega) \cdot \frac{2\pi}{6*T} \sum^{\infty}_{k=-\infty} 
\delta(\Omega - \frac{2\pi k}{6*T})

Hallo Unit*, ich muss doch diese Formel verwenden:

X(k) = sum_{n=0}^{N-1}x(n) \cdot e^{-j2\pi\frac{nk}{N}}

-----------------------------------------------------------------------

Warum eigentlich nur X(0) .. X(3) verwenden?
Weil nach X(3) das ganze gespiegelt ist.
Gast #800049
Lesenswert?

Weil du das Spektrum im Bereich 0 .. 500 Hz darstellen musst...
f = 0Hz                <-> X(0)
f = 1000*1/6 = 166.6Hz <-> X(1)
f = 1000*2/6 = 333.3Hz <-> X(2)
f = 1000*3/6 = 500Hz   <-> X(3)

> Warum eigentlich nur X(0) .. X(3) verwenden?
> Weil nach X(3) das ganze gespiegelt ist.
Gast #800054
Lesenswert?

Hallo Stefan, deine angegebene Formel ist doch falsch.
X(w) müsste doch so aussehen oder?

[latex]
x[k] = p[k] * \sum^{\infty}_{k=-\infty} \delta[n-k \cdot 6]

X(\Omega) = P(\Omega) \cdot \frac{2\pi}{6*T} \sum^{\infty}_{k=-\infty}
\delta(\Omega - \frac{2\pi k}{6*T})
[/latex]

Hallo Unit*, ich muss doch diese Formel verwenden:

[latex]
X(k) = sum_{n=0}^{N-1}x(n) \cdot e^{-j2\pi\frac{nk}{N}}
[/latex]
-----------------------------------------------------------------------

Warum eigentlich nur X(0) .. X(3) verwenden?
Weil nach X(3) das ganze gespiegelt ist.
Gast #800244
Lesenswert?

Hallo Stefan, deine angegebene Formel ist doch falsch.
X(w) müsste doch so aussehen oder?

Hallo Unit*, ich muss doch diese Formel verwenden:
-----------------------------------------------------------------------

Warum eigentlich nur X(0) .. X(3) verwenden?
Weil nach X(3) das ganze gespiegelt ist.
Gast #800251
Lesenswert?

Hallo Stefan, deine angegebene Formel ist doch falsch.
X(w) müsste doch so aussehen oder?

Hallo Unit*, ich muss doch diese Formel verwenden:
#800328
Lesenswert?

Hi,

ja ich habe halt zeitnormiert gerechnet. Natürlich musst du in deinem 
Fall das T zu den 6 mit dazunehmen, da man es später braucht. Aber an 
der Formel ändert sich dadurch nichts. Das findest du in jeder 
ZDFT-Tabelle.

Das p[k] kannst du in diesem einfachen Fall einfach als Summe von fünf 
verschobenen Diracs mit der dementsprechenden "Höhe" betrachten. Danach 
kannst du diese porblemlos transformieren und siehst auf Anhieb den 
Cosinus.

Gruß,
Stefan
Gast #800743
Lesenswert?

Hallo Stefan,

ok ja aber warum steht in der verschobenen DIRAC Funktion ein PLUS drin 
und nicht ein MINUS. Im Beitrag "Datum: 03.03.2008 12:17" habe ich ein 
Minus verwendet. Ist dies nun falsch oder richtig?
Gast #863493
Lesenswert?

Hallo Stefan,

ich hab zu dieser Aufgabe noch eine Sache, die mir unklar ist.

Die Variable "T" ist doch der Kehrwert der Abtastfrequenz fs.

In den folgenden Rechenschritte hab ich ein weiteres Problem.
Warum ist wo = 6T ?

  
Gast #863500
Lesenswert?

Diesmal die korrekte Darstellung der Gleichungen:

  

Antwort schreiben

Bitte melde dich an, um einen Beitrag zu schreiben.

oder

Mit Google-Account einloggen

Die Registrierung ist kostenlos und dauert nur eine Minute.

Jetzt registrieren