Funktion som beskriver Fletcher-Munsonkurvorna?

av **Naqref** » 2004-03-24 15:44

Undrar om det finns några matematiskt beskrivna överföringsfunktioner för dessa vackra kurvor?

Någon som har koll på det? Kan ju vara olika beroende på vilken normalnivå man ska beskriva kurvorna på.

av **Wolfie** » 2004-03-24 15:53

Fletcher-muns... vad för något?

Beskriv gârn va de är té för...

av **Naqref** » 2004-03-24 16:00

ah!

Det är dessa:
Bild

Som visar hur högt vi upplever ljud med olika toner subjektivt.

av **Wolfie** » 2004-03-24 16:12

Bra fråga... Det borde ju rimligtvis finnas någon nerd som beskrivit funktionen matematiskt...

Annars så får man väl göra det själv?

Skicka in amplituderna vid de olika frekvenserna i Matlab

Läs sedan hjälpen för system identification (hjälpfilerna finns inkluderade i system identification toolbox) så hittar du säkert något lämpligt sätt att få fram en överföringsfunktion. En ARX-modell borde räcka då frekvenssvaret inte är sådär våldsamt komplicerat.

Förutom att det varierar i frekvens beroende på amplitud... Darn...

Annars så kan du bygga ett adaptivt rekursivt FIR-filter som kopplas som en "adaptive noise canceller" (finns färdiga under DSP-toolboxen). Det fina med en sådan är att den så gott den kan sätter sina vikter till att efterlikna överföringsfunktionen till det som den ska försöka "tracka".
Jag kan vid tillfälle se om jag hittar någon som jag byggt under tiden jag pluggade på KAU (borde ligga sparade på något gammal hårddisk någonstans).

av **Naqref** » 2004-03-24 17:49

Tackar! Kollar lite i hjälpen för system identification!

Behöver bara ha en kurva åt gången nämligen.

av **Svante** » 2004-03-24 18:10

Nagref:
Nyfiken... Vad vill du göra?

av **Naqref** » 2004-03-24 18:36

Tro det eller ej så mekar jag ihop lite ljudfiler till LTS nya hemsida som man kan hämta ner för att göra tester (ser ut som ett sammanträffande och det var faktiskt helt okorrelerat med det som försigår i tråden härintill!

) och tänkte då att man skulle kanske kunna ha några signaler korrigerade för örats känslighet så att man kanske mer rimligt kan bedöma subjektivt avvikelser. Om det nu blir så bra vet jag inte men det är värt ett försök!

av **Svante** » 2004-03-24 19:28

Aha, en loudnessknapp i datorn! Det ser vi fram emot. Perfekt blir det kanske inte, örat är lite olikt på sinusar och sammansatta ljud var det nån som sa till mig nån gång. Men lite nästankompensation kan ju vara bättre än ingen alls. Och förmodligen finns en del att lära under resans gång.

av **Wolfie** » 2004-03-24 21:40

Naqref: För att göra det så enkelt som möjligt för dig...

kör kommandot

>> fdatool

och se vad som händer

Du kan här grafiskt sätta ut det frekvenssvaret som du vill ha hos ett filter och vips så får du automagiskt filterkoefficienterna!
Går att välja lite olika typer av filter, beroende på tycke och smak... Men lämpligt är väl FIR filter som är faslinjära? :wink:

av **Wolfie** » 2004-03-24 21:59

Detta verktyg är nästan helt obegränsat över vad som går att göra
i filterväg, men det kan vara lite bökigt i början.. Så ta god tid på dig och gå igenom exempel och hjälpen så kommer du nog en bra bit ganska så snabbt.

Har för mig att det ska gå att bygga multiband filter i nästan obegränsade former vilket kan vara vad du söker...

Om inte så finns ju filter() funktionen som genererar ett "analogt" filter. filt() är den digitala motsvarigheten. Där går oxo att att i en vektor beskriva frekvenssvaret mha olika graders amplitud vid olika frekvenser...

Tips!! Gå in på www.ee.kau.se och leta reda på lite färdiga matlabfilter m.m.

Finns en del "labinfo" under kurserna signalbehandling, DSP, Reglerteknik mfl...

av **Naqref** » 2004-03-24 22:56

låter som väldigt intressanta saker! Ska undersökas imorgon! Kan ju användas vid framtida filterberäkningar oxå...

av **Wolfie** » 2004-03-25 21:56

Kollade lite gamla labbar igårkväll...

lek lite med kommandot remez()

>>help remez %gäller nog här

Det är ett fiiiint verktyg för att tillverka faslinjära FIR-filter

"Parks-McClellan optimal equiripple FIR filter design.
B=REMEZ(N,F,A) returns a length N+1 linear phase (real, symmetric coefficients) FIR filter which has the best approximation to the desired frequency response described by F and A in the minimax sense."

Ett annat rätt så användbart kommando...

>>freqz

" [H,W] = FREQZ(B,A,N) returns the N-point complex frequency response
vector H and the N-point frequency vector W in radians/sample of
the filter:
jw -jw -jmw
jw B(e) b(1) + b(2)e + .... + b(m+1)e
H(e) = ---- = ------------------------------------
jw -jw -jnw
A(e) a(1) + a(2)e + .... + a(n+1)e
given numerator and denominator coefficients in vectors B and A. The
frequency response is evaluated at N points equally spaced around the
upper half of the unit circle. If N isn't specified, it defaults to
512."

Sen så finns ju färdiga funktioner för Butterworth eller Chebyshev filter.. butter() cheby1, cheby2 där man anger gränsfrekvens och ordning så får man en färdig överföringsfunktion... (simulerar analoga filter).

fir1 och fir2 är inte heller dumma....

Sen så finns ju självklart funktioner som laplace, ilaplace, ztrans och fourier för att göra livet enkelt...

Vill man veta rötterna hos ett jobbigt polynom...
skriv roots(x) och vips...

Sade jag någon gång att jag gillar Matlab :?: