|
|
| Atmos_kk |
Forum-Fortgeschrittener
|
 |
Beiträge: 93
|
|
|
|
Anmeldedatum: 23.05.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 20.07.2011, 09:18
Titel: A-Bewertungsfilter
|
 |
Moin zusammen,
ich verzweifle langsam.
Ich habe eine Akkustikmessung in dB und muss das nun mit dem richtigen Filter in dB(A) umrechnen. Direkt soll das nicht möglichsein, kann ich auch so hinnehmen aber irgendwie muss man es ja hinbekommen.
Ich denke mal, dass es mithilfe der Frequenzen irgendwie funktioniert allerdings finde ich den richtigen Filter dafür nicht.
Hat jemand ein Schlagwort für mich oder eine kleine Anleitung?
Greetz
|
|
|
|
|
|
| Maddy |
Ehrenmitglied
|
|
Beiträge: 494
|
|
|
|
Anmeldedatum: 02.10.08
|
|
|
|
Wohnort: Greifswald
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 20.07.2011, 09:55
Titel:
|
|
|
|
|
| Atmos_kk |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 93
|
|
|
|
Anmeldedatum: 23.05.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 20.07.2011, 10:21
Titel:
|
|
|
Alles klar das hab ich auch schon mal gelesen. Aber wie benutze ich das jetz? Muss ich einfach die Gleichung mit in meine Berechnung einfügen und einfach alle Frequenzen da durchjagen und anschließend mit meinen dB muliplizieren?
|
|
|
|
| DSP |
Forum-Meister
|
|
Beiträge: 2.117
|
|
|
|
Anmeldedatum: 28.02.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: R2014b
|
|
|
|
|
|
|
Verfasst am: 20.07.2011, 10:37
Titel:
|
|
Habe das noch nie gemacht, aber ich denke, dass man so vorgehen muss.
Du hast neben den dB Werten deiner Messung auch einen Frequenzvektor
|
|
|
|
| Atmos_kk |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 93
|
|
|
|
Anmeldedatum: 23.05.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 20.07.2011, 10:47
Titel:
|
|
Ja so sah auch meine Vermutung aus!
Werde dann mal loslegen und wenn es hinhaut werde ich das mal posten. Vielleicht gibts ja noch mehr Anfänger wie mich die sich darüber stundenlang den Kopf zerbrechen.
Danke an alle!
|
|
|
|
| Atmos_kk |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 93
|
|
|
|
Anmeldedatum: 23.05.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 07:38
Titel:
|
|
Moin zusammen,
ich hab da nochmal eine kleine Frage. Gehört eher zur Signalverarbeitung aber da wir hier eh schon über A-Bewertung sprechen mach ich das mal hier.
Wann bzw. wo führe ich die Filterung überhaupt durch? Sagen wir mal ich hab den Schalldruck über die Zeit gegeben dann ermittle ich daraus erstmal meinen Schalldruckpegel über die Zeit! Da ich für die a-Bewertung die Frequenzen benötige mach ich eine FFT und errechne mir meine Faktoren. Benutze ich diese nun in der FFT und führe eine ifft durch um auf dB(A) zu bekommen oder gehe ich direkt zurück zu meinem Schalldruckpegel und wende die Faktoren darauf an?
Wenn ich nichts falsch gemacht habe kommt in beiden Fällen leider nicht das Selbe heraus!
Greetz
|
|
|
|
| Maddy |
Ehrenmitglied
|
|
Beiträge: 494
|
|
|
|
Anmeldedatum: 02.10.08
|
|
|
|
Wohnort: Greifswald
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 08:54
Titel:
|
|
Die Bewertungsfunktion ist über der Frequenz gegeben und muss deshalb auch im Frequenzraum angewendet werden. Wie sollte man das auch im Ortsraum machen...
Oder verstehe ich dich falsch?
_________________
>> why
The computer did it.
|
|
|
|
| Atmos_kk |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 93
|
|
|
|
Anmeldedatum: 23.05.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 09:06
Titel:
|
|
Da bin ich mir auch nicht ganz sicher wie wir uns verstehen 
Wahrscheinlich hab ich mich doff ausgedrückt!
Also hier mal eine kurze Einteilung meines Ablaufes:
1. Einlesen der Daten
2. Aufteilen der Daten in Zeit- und Beschleunigungsvektoren
3. FFT der Beschleunigungen ( die benötige ich ja für die Formel der A-Bewertung)
4. Berechnen des Vektors Ra ( Formel für A-Bewertung )
und nun mein Problem:
Jetzt kann ich ja ohne Probleme diesen Faktor Ra
1. Anwenden im Zeitbereich auf die Beschleunigungen
2. Anwenden im Frequenzbereich auf meine Ausschläge ( oder was die y-Achse auch immer anzeigt) und anschließen müsste man da noch eine ifft machen um in den Zeitbereich zu gelangen ( Beschleunigung über Zeit!)
Nur die Ergebnisse sind nicht identisch also hab ich daraus gefolgert, dass eins von beiden falsch sein muss! Aber was und warum bekomme ich nicht raus!
Greetz
|
|
|
|
| DSP |
Forum-Meister
|
|
Beiträge: 2.117
|
|
|
|
Anmeldedatum: 28.02.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: R2014b
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 09:15
Titel:
|
|
Benötigst du denn wieder ein Zeitsignal?
Die Filterung durch Multiplikation mit der Gewichtsfunktion findet auf jeden Fall im Frequenzbereich statt. Das Ergebnis wird dann mittels ifft wieder in den Zeitbereich transformiert.
Um die Filterung im Zeitbereich durchzuführen, musst du die Gewichtsfunktion erst in den Zeitbereich transformieren und dann eine Faltung von Signal und Filter durchführen.
|
|
|
|
| Atmos_kk |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 93
|
|
|
|
Anmeldedatum: 23.05.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 09:28
Titel:
|
|
Das hört sich auf jeden Fall plausibel an. Allerdings habe ich dann direkt das nächste Problem:
Bei der A-Bewertung im Frequenzbereich fallen meine niedrigen Frequenzen vollkommen weg, was wenn man sich die Kurven ansieht völlig klar scheint. Also was soll man dann machen. Die ifft muss dann in dem Fall doch falsche Ergebnisse liefern, oder hab ich da mal wieder einen Denkfahler?
|
|
|
|
| DSP |
Forum-Meister
|
|
Beiträge: 2.117
|
|
|
|
Anmeldedatum: 28.02.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: R2014b
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 09:29
Titel:
|
|
Allerdings fällt mir jetzt noch ein, dass man so eigentlich den Faltungsatz verletzt. Das Faltungstheorem von Signal(n) und Ra(m), wobei in diesem Fall m=n ist, ergibt ein Ergebnisvektor der Länge n+m-1. So wäre es unter Benutzung der Matlabfunktion conv().
Im Frequenzbereich muss man dann eigentlich folgendes machen:
Jetzt sollten conv() und FFT Faltung identisch sein. Für conv() benötige ich aber die Gewichtsfunktion Ra im Zeitbereich, also ifft(Ra).
|
|
|
|
| Atmos_kk |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 93
|
|
|
|
Anmeldedatum: 23.05.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 09:33
Titel:
|
|
Also das übersteigt meine Kenntnisse der Signalverarbeitung sowie die mit Matlab.
Ist es also garnicht möglich den A-Bewertungsfilter einfach so zu benutzen?
|
|
|
|
| DSP |
Forum-Meister
|
|
Beiträge: 2.117
|
|
|
|
Anmeldedatum: 28.02.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: R2014b
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 09:41
Titel:
|
|
| |
 |
 |
Das die niedrigen Frequenzen wegfallen, ist für mich eigentlich logisch. Im niederfrequenten Bereich hört man nichts...erst ab ca. 20 Hz...15-20kHz.
Nochmal zur Faltung...anderseits bin ich mir jetzt doch unschlüssig wegen der Erweiterung von Ra mit Nullen im Frequenzbereich. Normalerweise führt man die Faltung per FFT durch, da sie bei langen Signalen schneller ist als die Faltung im Zeitbereich. Da erweitert man aber die Impulsantwort des Filters (liegt im Zeitbereich vor) erst mit Nullen und transformiert dann in den Freq.-bereich. Da du das Filter ja aber schon im Frequenzbereich hast, könnte das erweitern mit Nullen auch falsch sein. Dann darf man das Zeisignal aber auch nicht erweitern, da H und Ra für die Multip. gleich lang sein müssen. Probier einfach mal beide Varianten.
Edit: Um deine Frage letzte Frage zu beantworten...doch! Einfach im im Frequenzbereich multiplizieren.
Zuletzt bearbeitet von DSP am 21.07.2011, 09:49, insgesamt einmal bearbeitet
|
|
|
|
| DSP |
Forum-Meister
|
|
Beiträge: 2.117
|
|
|
|
Anmeldedatum: 28.02.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: R2014b
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 09:46
Titel:
|
|
Sorry...habe noch einen Punkt vergessen.
|
|
|
|
| Atmos_kk |
Themenstarter
Forum-Fortgeschrittener
|
|
Beiträge: 93
|
|
|
|
Anmeldedatum: 23.05.11
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 21.07.2011, 09:47
Titel:
|
|
Ich sag mal so, meine Frequenzen fallen bis ungefähr 500 Hz nahezu vollständig weg und es bleibt nichts mehr übrig. Wie soll ich denn dann wieder zurück in den Zeitbereich kommen ohne Fehler zu erhalten? Das verstehe ich nicht ganz.
Oder muss ich vor der Filterung mein Signal halbieren ( hab ja eh die Spiegelung ) und dann den Filter drüber jagen?
Naja, probieren geht über studieren würde ich sagen aber das ganze ist irgendwie ein Rätsel für mich!
Greetz
|
|
|
|
|
Gehe zu Seite 1, 2, 3, 4 Weiter
|
|
Einstellungen und Berechtigungen
|
|
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
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
|  RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
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.
|
|
