Datum: 22.04.2008 07:42
Hallo, zum besseren Verständnis der Fourier Transformation habe ich die Fouriertransformation schriftlich gerechnet. Dafür habe ich die Formel der DFT verwendet. ( http://de.wikipedia.org/wiki/Fourier-Transformation ) Allerdings klappt es noch nicht. Ich habe die Formel mit for-schleifen realisiert. Allerdings stimmt da was nicht. Vielleicht kann mir jemand einen Tip geben. Hier mein Code: time = (0.01:0.01:1) N = length(time) w=2*pi sinus = 1.*sin(w.*10.*time); for k=1:N for n=1:N for time = 0.01:0.01:1 erg = 1.*sin(w.*10.*time) end end end ergebnis = abs(erg)/N; stem(time,ergebnis) Danke Gruß Dennis
Datum: 22.04.2008 08:11
was für Code soll das sein? liegt es vieleicht daran dass du in deiner innersten for-schleife ständig "erg" überschreibst ? bei Fourierreihen integrierst du über das Produkt aus gesuchten Sinus(u. Cosin.) und Eingangssignal. das gibt dir die geuchte Intensität einer Frequenz im Einganssignal
Datum: 22.04.2008 08:31
Oh, ich programmiere mit Matlab. Ja, ich denke da hast du recht und ich überschreibe erg immer wieder. Danke für den Tip. Hast du denn noch eine Idee, wie erg nicht überschrieben wird? ;) Danke
Datum: 22.04.2008 09:28
Müsste passen.
clear j;
time = (0.01:0.01:1)
w=2*pi
sinus = 1.*sin(w.*10.*time);
N = length(sinus)
erg = zeros(N, 1);
for k=1:N
for n=1:N
erg(k) = erg(k) + ...
( sinus(n) * exp(-j*2*pi/N * (k-1) * (n-1)) ); % ???
end
end
ergebnis = abs(erg)/N;
stem(time,ergebnis)
|
Datum: 22.04.2008 09:48
Ja, coole Sache. Danke. Es klappt. Wenn du mir jetzt noch sagen kannst, wieso die Amplitude der Transformation bei 0.5 ist und nicht bei 1, wäre ich dir sehr dankbar. Danke
Datum: 22.04.2008 10:19
Weil das Spektrum konjugiert komplex ist und somit zwei "gleiche" Bin's auftreten.
Datum: 23.04.2008 14:36
Mmmmh, ok, aber so 100%ig habe ich es noch nicht verstanden. Wie meinst du das mit den zwei gleichen Bins? Liegt das an der Nyquist Frequenz? Danke
Datum: 23.04.2008 18:04
Nimm dir mal ein Blatt Papier und schreib dir für ein Testsignal mit wenigen Abtastwerten die DFT auf. Du wirst dabei feststellen, dass immer 2 Werte komplex konjugiert zueinander auftreten. Die Frequenz der beiden ist f0 und fs-f0.
Antwort schreiben
Die Angabe einer Email-Adresse ist freiwillig. Wenn Sie automatisch per Email über Antworten auf Ihren Beitrag informiert werden möchten, melden Sie sich bitte an.
Wichtige Regeln - erst lesen, dann posten!
- Suchfunktion und Betreffsuche benutzen - vielleicht gibt es schon einen ähnlichen Beitrag
- Aussagekräftigen Betreff wählen
- Im Betreff angeben um welchen Controllertyp es geht (AVR, PIC, ...)
- Groß- und Kleinschreibung verwenden
- Längeren Sourcecode nicht im Text einfügen, sondern als Dateianhang
- JPEG-Dateien (.jpg) nur für Fotos verwenden, Schaltpläne, Screenshots usw. als PNG oder GIF anhängen
Formatierung (mehr Informationen...)
- [c]C-Code[/c]
- [avrasm]AVR-Assembler-Code[/avrasm]
- [pre]vorformatierter Text (z.B. Code in anderen Sprachen)[/pre]
- [math]Formel in LaTeX-Syntax[/math]
- [[Titel]] - Link zu Artikel