Filterkoefficienter...

[Svante]

Här kommer en hitech-fråga, som ändå är lotech...

När man bygger ett andra ordningens filter digitalt så har det filterkoefficienter a0-a2 och b0-b2. Jag behöver explicita funktioner av typen a0=f(f0,Q,fs), a1=g(f0,Q,fs) osv...

Jag har de här i mina program, men jag börjar misstänka att de är fel. Jag har letat mig blå på nätet efter dem, men hittar bara hur man ska ta fram dem och inte resultatet.

Man kan ju göra det antingen med bilinjärtransform eller impulsinvarians, och egentligen vill jag ha båda.

Någon?

[MagnusÖstberg]

En fråga för lilltroll?

[Svante]

En fråga för lilltroll?

Ja, fast han kommer bara att slänga ur sig Matlabkod. Eller länkar till Matlabs manualer.

[MagnusÖstberg]

[kodapa]

Ja, fast han kommer bara att slänga ur sig Matlabkod.

Det kan väl inte bli bättre?

[Svante]

Ja, fast han kommer bara att slänga ur sig Matlabkod.

Det kan väl inte bli bättre?

Jo, Delphikod.

[E]

Varför misstänker du att de du har är fel? Är det något i MasVis som
inte står rätt till?

Mvh E*

[Svante]

Varför misstänker du att de du har är fel? Är det något i MasVis som
inte står rätt till?

Mvh E*

Nä. Jag håller på med en ljudeditor som ska kunna göra filtrering med bla 2:a ordningens filter. Ett sånt med Q=0,5 ska inte ringa, men det gör mitt.

[mike34]

Tja eftersom jag vet precis vad det är frågan om nu så måste jag lägga mig i.
SKOJAR !!!!

Lilltroll har ju lösningen så klart. Men filter har jag i kod hehe för xcore i XC . lite cryptiska metoder som lämnar dubbla returvärden. men det går väl att hacka om.

Men ringningen , tja, avrundningsfel eller overflow i variablerna ?
Mycket jag inte förstår har jag fått se sen jag satte snurr på de där filtrena.

[E]

Nä.

(Puh!)

[kodapa]

Ja, fast han kommer bara att slänga ur sig Matlabkod.

Det kan väl inte bli bättre?

Jo, Delphikod.

[-Martin-]

Svante, ring Herr Fnozkberg vettja.

[faulhund]

Här kommer en hitech-fråga, som ändå är lotech...

Ja, man gnuggar sig i ögonen, och vet inte om man läst rätt ... när självaste Svante frågar om något så grundläggande ...

Om det är LP-filter av andra ordningen är formel 7.6.5 i Fnozkberg svaret, och för högre ordningar krävs att man har partialbråksuppdelningen. (Så bra att jag har den inom räckhåll nu, sedan den varit nedpackad i 20 år!)

Man kan ju göra det antingen med bilinjärtransform eller impulsinvarians, och egentligen vill jag ha båda.

Båda samtidigt går ju inte, men du kanske menar att du vill implementera båda som alternativa val. För båda metoderna är ju fullständiga och entydiga, så i annat fall blir problemet överbestämt. En lösning för impulsinvarians skulle kanske kunna vara att (med hjälp av optimering då) approximera ett impulssvar (eventuellt givet i godtyckliga samplingpunkter i godtyckligt antal) med en ARMA-funktion av givet gradtal.

Vad har du kommit fram till?

[Svante]

Alltså, det är ju så enkelt, men jag får inte till det. Det handlar om andra ordningens filter, alltså ett polpar och ett nollställepar.

Hur går jag från Q-värden och brytfrekvenser (för pol- och nollställeparen) (vid en given samplingsfrekvens) till filterkoefficienterna a0-2 och b0-2.

Fnosken, ja, de kurserna genomsyrades av en rätt skön humor. Jag har för mig att jag bläddrade i formelsamlingen, men att jag inte hittade det jag ville, åtminstone inte explicit.

Jag är förvånad att jag inte hittar den där övergången någonstans, det måste ju vara vardagsmat för alla som bygger filter.

Faulhund: hur ser 7.6.5 ut? Jag har inte samlingen här (om det nu är saamma som jag hade).

[n3mmr]

Alltså, det är ju så enkelt, men jag får inte till det. Det handlar om andra ordningens filter, alltså ett polpar och ett nollställepar.

Hur går jag från Q-värden och brytfrekvenser (för pol- och nollställeparen) (vid en given samplingsfrekvens) till filterkoefficienterna a0-2 och b0-2.

Fnosken, ja, de kurserna genomsyrades av en rätt skön humor. Jag har för mig att jag bläddrade i formelsamlingen, men att jag inte hittade det jag ville, åtminstone inte explicit.

Jag är förvånad att jag inte hittar den där övergången någonstans, det måste ju vara vardagsmat för alla som bygger filter.

Faulhund: hur ser 7.6.5 ut? Jag har inte samlingen här (om det nu är saamma som jag hade).

Berätta gärna vem Hr Fnozkberg är, eller vad Fnosken är för kurs eller nåt!

[Tore]

Har du kollat denna?

http://www.musicdsp.org/files/Audio-EQ-Cookbook.txt

[lilltroll]

Mja, jag har en del i ANSI C kod oxå numera.

Men den första delen är "lätt", hur man trollar om från Q, fs och Gain, till a[] samt b[] i den tidskont. världen. Sen får du röja på med Z transform.

Titta mitt på denna sida
http://en.wikipedia.org/wiki/Q_factor



PS. Jag sitter mest och imp filter i asm, med omvänd pipelineordning just nu. Dvs att det sista pipelinesteget gör fetch, och det första pipelinesteget gör instruction decode. Längre bort än så kommer man inte från MATLAB!
Varför gör man så då.
Jo En branch kan då avsluta med en fetch på nästa inst. i sista steget, och i nästa operation görs då decode på denna till att börja med, utan att man förlorat en enda cykel pga branch eller missar i en branch prediction. Det medför att man kan exekvera minst 4 HW trådar parallellt utan att de någonsin krigar om L1 cache acess ... vilket lätt inses typ DS

Tiderna förändras, nu sitter Svante i något äckligt högnivå, medan jag skriver ren ASM.

[-Martin-]

Berätta gärna vem Hr Fnozkberg är, eller vad Fnosken är för kurs eller nåt!

"Fnozkberg" (eller egentligen Lars F Nozkberg om rätt ska vara rätt) var under den preudonym som Kåre Mossberg (som ledde filterkursen/kurserna på KTH) skrev sin litteratur och som användes under kurserna. Studentikååst så det står härliga till...
"Fnosk" är väl egentligen samma som (normerad) frekvens fast i "samplingsdomänen". Fnosk = 1 motsvarar väl gränsen på 1:a Nyqvist om jag inte minns fel(?)

[lilltroll]

Jo, men när Mossberg försvann inkl T3 blev S3 så fick inte FNOSKen vara kvar.
Numera motsvarar fs == 1 normerat