hifi.slovanet.sk

[miero]

dustin, co tak na zaciatok uplne vypustit virtualne zvukovky a urobit to pre mataa. alebo je to narocnejsie na vykon a ty to chces do meracej stanice?

[opa]

Počítá se s variantou, že DUT bude mít přenos menší než 1 ? Potom by musela být možnost kalibrační dělič (potenciometr) přesunout do kalibrační cesty místo DUT, aby simuloval zeslabení v signálové cestě.
Nebo by za (před) DUT musel být předzesilovač s malým zkreslením, který dotáhne přenos na hodnotu větší nebo rovnou 1.

A ještě jeden drobný problém. Řekli jsme si, že ta korekce zahrnuje zkreslení vstupu i výstupu zvukovky najednou. Pokud ovšem zkreslení vstupu bude nezanedbatelné a DUT nebude mít zcela rovnou fázovou charakteristiku, tak se část korekce příslušející vstupu nepotká ve fázi se zkreslením, které vzniká až v té vstupní části a THD poleze nahoru.
Když to přeženu, tak pasivní RC obvod nulovým zkreslením, který nějak zatočí s fází, bude vykazovat zhoršené THD.

[dustin]

[PMA]

[dustin]

[dustin]

[PMA]

Protoze z jednoho PC je tam smycka. Si to vyzkousej, bez traf, optiky nebo symetriky to nejde. Chce to zkusit, delal jsem to mockrat. Proste z teorie na zem.

[dustin]

[dustin]

[opa]

No tak koukám, že je to rozsáhlejší projekt, než jsem myslel a že je patrně na všechno pamatováno. Takže omezím své kibicování a budu držet palce, ať se daří.
S programováním bohužel pomoct neumím, tak se budu držet svého analogu.

[PMA]

[miero]

dustin, trochu som vylepsil spectrum_analyser z playrec-u. uz sa na to da pozerat :D
- https://github.com/miero/playrec-examples/commit/39983a056cf602cde412fede36f6db9568b79b05

skusim sa pohrat s dialogmi, a pridat tam nejake tlacitka start/stop/screenshot a povedzme vypis harmonickych

PMA, ak sa najde niekto ochotny platit (resp. sponzorovat) volne dostupny projekt, tak parada :-)

[dustin]

[dustin]

[dustin]

Miero, ty jsi drak. Parádní práce, díky moc. Otestujeme to tedy na tomhle skriptu?

[miero]

skusam i to mataa, ale tam je to "na rozjazd" az prilis pokrocile a integrovane

napriklad je tam i prototyp funkcie[1] merajucej/pocitajucej GedLee metriku[2] (z merania sumu odvodi prenosovu funkciu a z nej vypocita to cislo Gm), ale tu som nerozbehol ani na loopbacku ... chce to najprv nastudovat manual a riadne nakalibrovat

takze, early verziu najradsej s playrec ... je to dynamicke, co to sa najviac podoba klasickym programom :-)

EDIT: tak uz mam namerane i to Gm na loopbacku (pozor nenakalibrovane), ale chce to i premyslat a nie len copy&paste :-D

[1] https://github.com/mbrennwa/mataa/blob/master/mataa_tools/mataa_measure_GedLee.m
[2] http://www.gedlee.com/Papers/THD_.pdf

[opa]

[dustin]

Miero, jak prosím chápeš interpolaci v https://github.com/mbrennwa/mataa/blob/master/mataa_tools/mataa_measure_GedLee.m ?

[dustin]

Se spolehlivou funkcí interp1 souvisí https://www.mathworks.com/matlabcentral/answers/421846-interpolation-of-noisy-data - hodnoty musí být distinct, tak tam přidávají maličký šum.

[miero]

este som nasiel toto pre MATLAB:
- https://www.mathworks.com/help/signal/ref/tfestimate.html

[dustin]

Nevím, jestli by nám vrátila rovnou parametry harmonických, to asi ne. To by měl být gain vs. frekvence. Ale třeba se bude taky hodit, je dobré o tom všem vědět, abychom příp. neprogramovali kolo

[miero]

tie harmonicke by sme nasledne dostali z tej tf. no bolo by napr. zaujimave zobrazit ako sa meni tf s amplitudou. u lomenej ciary je to nazornejsie. a mohlo by to ist aj lepsie interpolovat pre rozne urovne.

[PMA]

Možná jsem úplně mimo, pracujete jen s reálnou částí FFT, nebo s reálnou + imaginární?

[miero]

Na realne hodnoty zo vstupu zvukovej karty pouzivame diskretnu fourierovu transformaciu, takze imaginarna cast (zaporne frekvencie) je zhodna s realnou castou. Avsak v kazdnom bine je komplexna hodnota, ktorej absolutna hodnota urcuje amplitudu toho binu.

[dustin]

Miero, prosím tě, co dělám blbě:

https://github.com/miero/nonlinear-compensation/commit/4e410537af7c332c805e968d46f5c8af14a1d8f4#diff-f6e45fa8399189abdde597f5a23c60b2R18 - teď generuje finální fft jenom z 240 vzorků, tj. zubaté grafy a nízká čísla amplitud - IMO protože hodně široký bin

Když zvednu na 5000 period tj. 200000 vzorků, parádní grafy, ale zase se to zasekne na smyčce https://github.com/miero/nonlinear-compensation/commit/4e410537af7c332c805e968d46f5c8af14a1d8f4#diff-f6e45fa8399189abdde597f5a23c60b2R51 , protože cykluje přes všechny výstupy fft.

Co bych tedy prosím měl volat, abych to hodilo původní grafy a výpis harmonických s dobrým rozlišení? Pracovně jsem 10 period změnil na 5000 a hledám ve 100 výsledcích, jako bylo předtím. To šlape OK.

Máme nějakou funkci, přes kterou bychom mohl vygenerovat tu aproximaci harmonických?

Našel jsi v octave nějaký způsob, jak run (F5) s parametry konkrétního skriptu, aby si to pamatovalo? Zatím musím pokaždé do skriptu přidávat konkrétní wav místo wavPath = arg_list{1} (to mi hodí chybu, nebo někdy má octave další parametr (např. --force-gui pro nejnovější octave z flatpacku) a ten je pak v arg_listu taky).

Měli bychom se asi domluvit, jak ty skripty tedy budeme spouštět. Osvědčil se mi include nějakého neverzovaného ignorovaného config souboru s nastavením základních parametrů, přičemž bychom verzovali nějaký example.

Předem díky moc za rady, nechci ti to rozvrtat, ale potřebujeme to sdílitelné s různým nastavením a zároveň spustitelné v GUI octave kvůli debugu.

[dustin]

Ještě prosba - netušíš, proč jsou amplitudy v peaks https://github.com/miero/nonlinear-compensation/commit/4e410537af7c332c805e968d46f5c8af14a1d8f4#diff-f6e45fa8399189abdde597f5a23c60b2R56 zaokrouhlená na celá čísla? Nikde jsem to nenašel, sorted2 je float, peaks vrací pouze celá čísla, i když tam u amplitudy žádné zaokrouhlení není:

[miero]

jaj, som tu pisal[1], ze si nemerguj ten moj posledni getHarmonics() ... ze som to rozbil a nefunguje to :-(

s tymi celymi cislami je matlab/octave hrozny ... ked sa ti tam raz dostanu do nejakeho pola pri definici, tak uz to konverti vsetko na int ... som to ladil dlho a zda sa stale nedoladil ... jednoduchsie je mat vsetko double a ked treba tak ceil/floor/round ... popripade kazdy int pri vypocte a ulozeni do pola skonvertit do double

tebe asi prave vadi toto, zmen si to docasne asi na round:
- maxfreq = int32(maxfreq);

alebo pouzi staru verziu, ktora vacsinou funguje... i ked niekedy zamlci harmonicke

mam rozpracovane, ze sa nebude vsetko triedit, ale len par binov okolo ocakavanej harmonickej

premyslal som tiez nad tymi parametrami a napada mi akurat extrahovat zo spustitelnych skriptov funkcie

[1] http://hifi.slovanet.sk/bb/viewtopic.php?p=317571#317571

[dustin]

Promiň, z dalších příspěvků jsem nevím proč pochopil, že už to běží dobře. No problema, resetnu HEAD o ten merge zpět, nikam jsem to nepushoval.

Na int můžeme ty zaokrouhlovací funkce, čistý int asi zatím nikde nepotřebujeme.

Tou konfigurací jsem myslel obyčejný include nějakého souboru s nastavenými proměnnými, pokud by existoval - tím by se přeskočily všechny ty testy na argvars, něco takového.

[miero]

dustin, aha aka prasarna funguje z prikazoveho riadku Octave :-D

[dustin]

Ale není to F5