Verfasst am: 26.07.2010, 18:04
Titel: Frequency Domain LMS-Algorithmus
Hallo Zusammen,
ich muss ein Matlab-Code erstellen für ein Frequency Domain LMS-Algorithmus. Das ist ein adaptive Filter, nur im Frequenzbereich implementiert.
Ich habe alle Schritte vor mir, nur mit dem Syntax komme ich nicht klar Ich bitte um Hilfe!! Es ist dringend!
vielleicht solltest du dann dein Problem etwas konkretisieren oder sagen was du möchtest, damit dir geholfen werden kann.
Ich habe mich vor einiger Zeit auch mit dem LMS-Algorithmus vor dem Hintergrund der Systemidentifikation beschäftigt, vielleicht kann ich dir weiterhelfen...
danke erstmal für deine schnelle Antwort.
So ich möchte ein Block LMS Algorithmus implementieren, dh:
- die Eingangsfolge x[n] wird in Blöcke der Länge C eingeteilt, wobei C die DFT-Länge ist. Um von einen Block zum anderen zu kommen, machen wir eine Verschiebung um L Werte.
-das Verfahren (LMS) wird dann für jeden Block im Frequenzbereich wiederholt. Sinn der Sache, ist es dass im Frequenzbereich eine Faltung einer einfachen elementweise Multiplikation entspricht!! (Rechenaufwand wird weniger )
So wird gemacht:
1.Filterung mit dem Overlap-Save-Verfahren im Frequenzbereich:
X[k]=FFT( x[k])
Y[k]=W[k] .* X[k]
y[k]=IFFT(y[k])
Fehlervektor: e[k]=d[k]-y[k]
E[k]=FFT(e[k])
2.Adaption des Filters im Frequenzbereich
-Leistungschätzung und Schrittweite
-Aufdatierung des Koeffizientenvektors
wobei
k=Blockindex, dh
x[1] ist en Vektor, der die ersten C Werte enthält
x[2] ist en Vektor, der um L Werte im Vergleich zu x[1] verschoben wurde und enthält ebenso C Werte
so allgemein: x[k]=[x[kL+L-C],........,x[kL+L-1]]
w[k] ist dann der für diesen Block entsprechende Gewichtsvektor. Er hat die Länge N, deswegen wird er um C-N Nullen auf die Länge C ausgedehnt: w[k]=[w1,w2,...,wN,0,0,....0]
Die Blockweise Verarbeitung bringt es mit sich, dass auch der Filterausgang jeweils in Blöcken der Länge L vorliegt.
Ich kann weiter machen um Einzelheiten zu klären
Ich hoffe das ist erstmal ein Anfang um mich weiter zu helfen.
Ich kann auch gerne Fragen zum Algorithmus beantworten.
was hast du denn bis jetzt geschrieben? Du solltest den Code den du bis jetzt hast schon mal posten, damit man darüber diskutieren kann... Ich denke nicht das jemand hier für Dich programmieren wird...
ein und schau dir das Beispiel an. Ich glaube hier im Forum hat jemand auch mal eine Art Tutorial zur FFT geschrieben.
Deine "x"-Matrix kannst du mit einer for Schleife erstellen. Dies sollte anhand eines Minimalbeispiels keine größeren Probleme darstellen. Beispiele:
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.
Ich bitte um Hilfe!! Es ist dringend!
)
RSS