DIY - faktiskt.se

[lilltroll]

I vårat specialfall när vi söker en FIR kompensator är polynomet b=konstant. Vi ska bara invertera frekvensgången hos h(k) och hitta en stabil AR modell. Alles

[lilltroll]

I vårat specialfall när vi söker en FIR kompensator är polynomet b=konstant. Vi ska bara invertera frekvensgången hos h(k) och hitta en stabil AR modell. Alles

Nu börjar jag minnas något långt långt bak.


Om vi söker:


och r(n) är h(n) diskreta autokorrelationsekvens

så löses det med *=komplexkonjugat


Den första matrisen är ett Toeplitz system av linjära ekvationer.

Vektorn a innehåller lösningen till problemet.
Det funna lösningen bli en minfas variant eftersom alla poler måste ligga innanför enhetscirkeln för att bibehålla stabilitet.
Det hela heter Levinson-Durbin recursion och finns färdig i MATLAB som a = levinson(r,n)

Nu har vi bra verktyg för att gå vidare !!!

Den kräver bara O(n^2) flops och blir då snabbare än / för stora n

Kom gärna med alternativa metoder att testa med!
Var det ungerfär i den här riktningen du tänkte Johan S, eller har du något annat spännande i påsen???

[lilltroll]

Det är ju ovan som är AR Yule-Walkers lösningsmetod
Yule-Walker ekvationerna löses m h a Levinson-Durbin rekursion

Då så, inte konstigt att Yule-Walker är bättre än Hilbert transform, samt att Yule-Walker gick mycket snabbare än invfreqz - skönt när bitarna faller på plats!

Nu har jag ialla fall fattat någonting

Nu ska jag ge mig själv bakläxa att läsa Kapitel 12.3 Parametric methods for power spectrum estimation, Digital Signal Processing av Proakis & Manolakis. En mycket bra bok.

[lilltroll]

Det här var riktigt svårt, inte MATLAB delen utan DSP delen. Något geni hade definerat FIR-coef. baklänges i DSP programmet. Varför göra så ? Dessutom varnar inte DSP när den inte hinner med, den gör saken halvklar och börjar på nästa. Mycket förvirrande.

Men tillslut - genom att ha använt supermetoden minimization of sum |B-H*A|^2*Wt har jag med enbart 256 tappar i ett FIR filter + ett digitalt IIR filter som delningsfilter åstakommit detta som nu kör i realtid på DSP:n

SCANSPEAK 8545 + digital minfas EQ + digitalt delningsfilter
Impulsvar


Tonkurva + dist

!!! OBS, okalibrerat!!!

Jag är riktigt imponerad av resultatet. Med 4 ggr längre FIR filter (1024) tappar så löser man det sista! Man ser det lättast i impulsvaret att det inte är perfekt utan där finns lite skräp. Nästa steg är att falta in ett allpass filter som ger konstant grupplöptid runt delningsfrekvensen för hela kalaset. Man ska dock vara försiktig här, så det inte föringer otrevligt i någon annan riktning. Lyckas man så får man dock en perfekt loobning onaxis tillsammans med diskanten.

Dags att undersöka vad som händer under 200 Hz.

PS. Kaffekoppen, det lilla du fick höra med digitala delningsfilter så gick tyvärr kompensatorn baklänges Nästan -retande DS. Nu fungerar det bättre!

[Kaffekoppen]

PS. Kaffekoppen, det lilla du fick höra med digitala delningsfilter så gick tyvärr kompensatorn baklänges Nästan -retande DS. Nu fungerar det bättre!

Perfekt.. då kan man ju be Martin komma dit och testköra sin superdistade L18 SEAS midd Kan vi kolla hur högt den kan användas utan att 3:e och 5:e disten kommer...

Ps... ville mest se hur man gjorde sakerna och det var rätt imponerande det.. ds

[JohanS]

Snyggt jobbat! Har du några mätningar på hur det ser ut off axis?

[lilltroll]

Jag tyckte det lät lite klent a basen, så jag körde gamla Tombstone på DSP´s utgång. Mycket riktigt så tappar den några dB under 100 Hz, vilket är mer än FIR-ripplet som man inte kommer ifrån utan att öka antalet tappar.
När man kör minsta kvadratmetoder, så optimeras det på linskala så optimeringen blir inte straffad mycket för att den strular till det lite i basen misstänker jag utan att veta säkert. Nu kan man vikta optimeringen, det kan vara ett sätt. Dessutom råkade jag använda en udda längd. Det blev 257 istället för 256.
Det är bara 1 nollställe eller i mitt udda fall, 1 nollställepar som ska flyttas lite i kompensatorn.

Får labba lite till.

Ska kolla hur det ser ut off axis, det brukar vara ett avslöjande test. Kanske får jag inte kompensera för uppbrytningarna så hårt i praktiken, för det straffar sig i alla andra riktningar så pass att det inte är värt det fina utseendet onaxis.

[lilltroll]

Det här är inte lätt Det är matematiskt svårt!

Men i alla fall. Nu har jag förbättrat alla metoderna, samt lagt till en metod som bygger på Levinson-Durbin med modifierad autokorrelation (Lilltroll special). Den är snabb och bra för mellan långa filter. Tackar Johan S för inspirationen!
Hilbert konvergerar mot rätt lösning ju fler tappar man använder.

Impulsvar


Här har jag beställt en digital EQ inom intervallet [200 15k] Hz
Alla metoder klarar det inom 1 dB. Den med kortast filter har det svårast = helt naturligt. Det trixiga är att lösningen får inte skena utanför området. Under 200 Hz så är den blå kurvan den korrekta - Hilbert transformen är alltid perfekt i basen. Med Levinson-Durbin blir det en liten bashöjning, och med de andra 2 "mini..." blir det en puckel runt 100 Hz.
Anledningen till att min Hilbert fungerade dåligt tidigare beror på mycket skit bakom spakarna. Det här är får svårt för att jag ska skriva rätt från början. Men nu är det bättre. Ny Trolldist kommer !

[lilltroll]

Jag som trodde jag skulle bygga högtalare i den här tråden. Det är nästan bara DIGSIG nu, men i alla fall.

Lite Lön för mödan.

Liten trevlig låda på 16 liter + lite dämpull och en port avstämd på 44 Hz (Scanspeak 8545). Nu ser man på den icketrunkerade grundtonen att basen inte faller. För att spana in det noga så kommer mätningar inne i lådan.

Tonkurva


Ca 30 grader offaxis
Tonkurva


Det blir bättre o och bättre för varje försök

[lilltroll]

Hmm - 16 liter var lite lite - men med LEGO nemas problemas:
Nedan i en 21 liters låda.
portens resonansfrekvens har nu sjunkit till 37 Hz

Tonkurva


Mätning är gjord inne i lådan, nya Trolldist 7.52 klarar av att räkna om det. Vid 150 Hz har vi ett tryckminima från en stående våg i lådan vid mikrofonen. Men under 150 Hz så är nog metoden ganska säker.

Behöver ett lite längre rör - använder det längsta 3" rör jag har som port. HIFIkit på måndag!

[i]

Verkar oerhört seriöst det här !
(Inte för att jag kan skilja på det här och Fatilarkalkylen... )

[lilltroll]

Jag har lyckats få till realtidsfiltrering i MATLAB under WINDOWS !!!
(Det enda som behövs är DAQ toolboxen)
MATLAB filtrerar nu i realtid på CPU:n !!!

Och det fungerar i praktiken. Inga hack mellan paketen.

Källan är en extern CD-spelare som går in på ljudkortet. Utgången går sedan direkt till förstärkaren.

TOKCOOOLT

Man kan t o m ändra filterparametrarna i realtid !

Jag trodde aldrig jag skulle kunna få det här att fungera under Windows !!!

OK, vart hittar man modifierade drivrutiner som låser upp alla kanaler till SB-kort. Har haft en sådan till SB-Live som fungerade bra.

[lilltroll]

Vadå? Om man dedikerar ett extra USB ljudkort som delningsfilter och ger processen hög prio så borde gemeneman kunna använda detta hemma på sin dator som kör i bakgrunden bara !?
Sedan kan man fortsätta att använda sin dator precis som vanligt. Man kopplar digital-ut från det vanliga ljudkortet till digital-in på "delningsfilterkortet".

[lilltroll]

Nu kan man välja på att köra ARMA modeller dvs filtrering med IIR filter i realtid eller MA modeller och då kör den FFT baserad FIR-filtrering i realtid

Men nu ska jag nog sova - om det går

[lilltroll]

Nu har jag bytt till CA18RLY för att se om jag gillar mellanregistret bättre. Seriekopplade med en ca 0.75 mH spole och fick en förbluffande fin tonkurva. Då spar man på det dynamiska området som man har tillgängligt i kompensatorn, dessutom fås en ev. linjariserande verkan på den elektriska sidan en bit upp i frekvens.

Tonkurva + dist

OBS okalibrerad, men 2.83 V RMS

Det måste vara lätt att bygga passiva filter till denna !!!

Den har inga vassa distpeakar under 4 kHz, är det månne bra för ljudupplevelsen!?
Jag trodde först att 18W8545 skulle låta bättre, men när en elgitarr väl träffar en peak så låter det riktigt illa. Kanske är det här bättre.

[lilltroll]

På med lite digitala leksaker !

Delningsfrekvens 2500.00 Hz
(Så säker är nog inte klockan på ett Soundblaster kort)

On axis, tonkurva och dist


30 grader off axis, tonkurva och dist


Ser riktigt lovande ut. Nivån peakar aldrig högre off axis än on axis vilket förebådar goda tider!? (Mikrofonen är nog lite längre bak, då medelnivån har sjunkit)

[lilltroll]

Här kommer det första digitala 2-vägs försöket.

CA18RLY+ScanSpeak 9300.

Mätte upp det på 50 cm onaxis, skapade digital EQ.

Sedan kontrollmättes systemet på 1 meters avstånd, dvs i en helt annan punkt. I detta försök användes 4:e ord. digitala Linkwitz-Riyley filter med en gränsfrekvens på 2.5 kHz. Dessa har den egenskapen att de summerar till 1 omm de två bidragen ligger i fas. För alla andra faslägen minskar nivån.



Vad kan man då se från detta.

1: Jag har inte kompenserat för fasen än. Detta leder till att huvudloben riktar sig ner mot golvet. Man ser några dB's minskning on axis p g a detta. Det löser sig efter att ett lämpligt allpassfilter har faltas in i EQ:n.

2: Kompensatorn för diskanten byggdes på 1 enda mätpunkt - utan någon som helt smoothing av kurvan. Baffeldiffraktionen komenseras nu i 1 punkt, 50 cm framför diskanten (on axis). I så gott som alla andra punkter är det nu sämre. Det finns många olika lösningar på detta.
Exempel kommer någon annan dag.

3: Disten är lägre för diskanten än för CA18RLY, kanske tål diskanten att delas lägre än 2500 Hz!

Nu ska jag lyssna lite på resultatet från omgång 1.

Jag är ändå imponerad att allt fungerar, även om det är en bra bit kvar till mål

Vad bör jag seriekoppla diskanten med för att minska disten på den elektriska sidan

[lilltroll]

Tycker nog att CA18RLY låter helt ok även på höga ljudtrycksnivåer - så länge den inte får kämpa med basen - inte dess styrka. Skulle nog inte bygga en tvåvägare med denna.

Det är en aha upplevelse att kunna ändra delningsfrekvensen i realtid. Delning på 1500 låter inte bra från SCANSPEAK 9300 - Delning på 8 kHz LÅTER INTE BRA från CA18RLY. Kom ihåg att tonkurvan fortfarande är rak i lyssningpunkten oberoende av delningen m h a allt det digitala joxet.

Måste dock fila på allpassfiltret så loobning blir bra

Gonatt butiken

[JohanS]

Lite funderingar:

Elementet kanske gör sig bäst i en 2.5-vägare, så att det får lite hjälp när baffelstödet inte vill ställa upp längre?

Hur tänkte du gå till väga rent praktiskt med allpassfiltret, det känns spontant krångligare att manipulera med fasen än med tonkurvan?

Riktigt kul att det börjar likna ett komplett delningsfilter redan, kan du även se hur disten ökar när det börjar låta sämre vid olika delningar?

[Kaffekoppen]

Hej Johan! Ser också ca 18rly som rysligt intressant kandidat till MTM/MMT

[lilltroll]

Så där ja, nu har jag äntligen fått fart på multikanal ASIO-ljudkortet. Dessutom har antivikningsfiltrena kommit till A/D D/A omvandlarna.

[lilltroll]

Lite funderingar:

Elementet kanske gör sig bäst i en 2.5-vägare, så att det får lite hjälp när baffelstödet inte vill ställa upp längre?

Hur tänkte du gå till väga rent praktiskt med allpassfiltret, det känns spontant krångligare att manipulera med fasen än med tonkurvan?

Riktigt kul att det börjar likna ett komplett delningsfilter redan, kan du även se hur disten ökar när det börjar låta sämre vid olika delningar?

Allpassfilter delen är lätt att göra matematiskt, svårare kan vara att veta hur man vill ha det.

Om du faltar ett allpassfilters impulsvar med sig själv baklänges, (motsvarar akf) vad få du då? Jo en impuls. Så man kan i den enklaste modellen falta digital EQ med hela det kompenserade systemets impulsvar baklänges.

I högtalarfallet så begränsar man dock kompensatorn i frekvensdomänen, och bivillkoret är att slippa långa förringningar.
Men det går bra att bara rita fasen i frekvensdomänen (med magnituden 1) och sedan beräkna impulsvaret i tidsdomänen för det.

Lättast är nog att beräkna grupplöptiden för systemet samt korrigera den till ett konstant värde runt delningsfrekvensen. Vilket motsvarar att det akustiska centret hamnar i samma punkt sett onaxis för det frekvens intervallet.

[Haakan_W]

väldigt intressant att kunna justera delningsfrekvens i realtid....

samt att kunna mäta disten i samma veva....du är helt inne på rätt spår här lilltrollet!

snart måste vi labba med monstret

[Kaffekoppen]

hoppas lilltroll får tid att titta hit igen...

[Haakan_W]

hoppas lilltroll får tid att titta hit igen...

ja detta är superintressant!

[lilltroll]

Tid kommer det alltid mera av

[Svante]

Tid kommer det alltid mera av

Fast ibland känns det som om man bor granne med ett svart hål...

[Zalve]

Tycker nog att CA18RLY låter helt ok även på höga ljudtrycksnivåer - så länge den inte får kämpa med basen - inte dess styrka. Skulle nog inte bygga en tvåvägare med denna.

Gonatt butiken

Tycker den passar sig utmärkt för en 2vägare, om man tänker sig lite aktivt basstöd...

[lilltroll]

Tid kommer det alltid mera av

Fast ibland känns det som om man bor granne med ett svart hål...

Betänk att det finns en gräns en bit ifrån det svarta hålets centrum, och passerar man den gränsen så finns det ingen väg tillbaka, i alla fall inte i samma form som man kom in

PS. Vad heter gränsen

[Naqref]

Händelsehorisonten.