Obsah fóra hifi.slovanet.sk hifi.slovanet.sk


 
 FAQFAQ   HľadaťHľadať   Zoznam užívateľovZoznam užívateľov   Užívateľské skupinyUžívateľské skupiny   RegistráciaRegistrácia 
 NastaveniaNastavenia   Súkromné správySúkromné správy   PrihláseniePrihlásenie 
https://www.tnt-audio.com/gif/smallstereo_white.gif

Impulsní odezva z frekvenčního průběhu

 
Pridať novú tému   Zaslať odpoveď    Obsah fóra hifi.slovanet.sk -> Dokumentácia - How to ...
Zobraziť predchádzajúcu tému :: Zobraziť nasledujúcu tému  
Autor Správa
opa
Hifi inventar


Založený: 24 február 2007
Príspevky: 10768
Bydlisko: Praha

PríspevokZaslal: Ut september 22, 2020 13:07:56    Predmet: Impulsní odezva z frekvenčního průběhu Odpovedať s citátom

Mám nějaký frekvenční průběh, definovaný stovkou bodů 50-20000 Hz.
Fázi nemám, ale mohu klidně předpokládat, že je všude nulová (vzniklá chyba mi nevadí).
Potřeboval bych z toho vypočítat impulsní odezvu, nejlépe přímo v programu REW.
Načíst a zobrazit průběh lze, ale nevím, jak dál.
Ale i jakákoliv jiná cesta, třeba i v jiném programu, by mi pomohla.
Uměl by prosím někdo poradit ?

_________________
Cui bono ?
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu
miero
Hifi inventar


Založený: 08 september 2010
Príspevky: 9926
Bydlisko: Praha

PríspevokZaslal: Ut september 22, 2020 13:29:24    Predmet: Odpovedať s citátom

navod pre matlab, resp. octave:

citácia:
If you only have real, positive frequency data... then that's all you got. Imaginary parts are all zero, so that is fine... you have data with zero phase.

You will still need to create a symmetrical dataset (i.e., "FR_data"). In order to create your "FR_data" dataset, you will take your current dataset, say "FR_positive_frequency_only" and do this:
FR_data = zeros(1,1000);
FR_data(2:501) = FR_positive_frequency_only;
FR_data(502:1000) = fliplr(FR_postive_frequency_only(1:end-1));

These amplitudes correspond to frequencies according to:
f_shift = ifftshift(f);

Your impulse response dataset is then just:
IR_data = ifft(FR_data); % No need to use ifftshift when constructing FR_data like above

Here the amplitude at 0 frequency (i.e., FR_data(1)) is just set to zero.

Note: if your "FR_positive_frequency_only" dataset is a column matrix, then use flipud instead of fliplr.

Zdroj: https://www.mathworks.com/matlabcentral/answers/38316-obtain-the-impulse-response-from-a-frequency-response#
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu
BV
Hifi inventar


Založený: 06 október 2006
Príspevky: 5216
Bydlisko: Trnava

PríspevokZaslal: Ut september 22, 2020 15:20:33    Predmet: Odpovedať s citátom

A v čom je problém? Naimportovať kmitočtový priebeh do REW (.txt. .csv, .frd..) a zobraziť potom impulznú dozvu. A tú exportovať ako wav,. txt.
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
opa
Hifi inventar


Založený: 24 február 2007
Príspevky: 10768
Bydlisko: Praha

PríspevokZaslal: Ut september 22, 2020 17:32:46    Predmet: Odpovedať s citátom

BV: To jsem si myslel taky, ale nic. Dokonce jsem pro jistotu doťukal ručně třetí sloupeček s fází (všude nuly), ani to nepomohlo.
Přikládám testovací soubor, jenom se musí přejmenovat na TXT.txt. Jinak to nešlo dát jako příloha.
Normálně zobrazí frekvenční průběh i tu fázi, ale impulsní odezva je prázdná. Při pokusu o export napíše, že nejsou data. Prostě žádný výpočet neproběhnul a nevím, jak k tomu program přinutit.

_________________
Cui bono ?
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu
godz
Hifi expert


Založený: 05 máj 2007
Príspevky: 1001
Bydlisko: Trenčín

PríspevokZaslal: Ut september 22, 2020 18:12:07    Predmet: Odpovedať s citátom

S REWom moc nepracujem ale na sto percent to funguje v HOLMimpulse. A tuším sa to dá urobiť aj v Xoveri.
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu
BV
Hifi inventar


Založený: 06 október 2006
Príspevky: 5216
Bydlisko: Trnava

PríspevokZaslal: Ut september 22, 2020 18:45:23    Predmet: Odpovedať s citátom

Je to "záhada", asi nejaký bug.. Ale prehnal som to cez HolmImpulse, načítal txt a vyexportoval impulzku ako wav. Je to v prílohe s príponou .pdf, treba zmeniť na wav.
Na obrázku je ako to vyzerá v HolmImpulse.


Naposledy upravil BV dňa Ut september 22, 2020 18:50:11, celkom upravené 2 krát.
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
PMA
Pavel Macura


Založený: 18 január 2007
Príspevky: 19729

PríspevokZaslal: Ut september 22, 2020 18:45:39    Predmet: Odpovedať s citátom

100 bodu, není to málo, milý Watsone?
_________________
Dispre2-JFET preamp, Z-30_2020 power amp, CNO-T25 repro
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
PMA
Pavel Macura


Založený: 18 január 2007
Príspevky: 19729

PríspevokZaslal: Ut september 22, 2020 19:11:58    Predmet: Odpovedať s citátom

Vypadá to, že to REW neumí. Zkuste to ještě někdo, z přiložené FCH - korektní data i s fází.
_________________
Dispre2-JFET preamp, Z-30_2020 power amp, CNO-T25 repro
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
opa
Hifi inventar


Založený: 24 február 2007
Príspevky: 10768
Bydlisko: Praha

PríspevokZaslal: Ut september 22, 2020 19:44:37    Predmet: Odpovedať s citátom

BV: díky, ještě s tím wav-em zápasím. V REW vychází ten tvar jinak, ale asi dělám někde chybu.

PMA: Ano, 100 bodů je málo, ale víc jich nemám. Je to kalibrační křivka mikrofonu a hledám tam odraz od držáku při kalibraci. Zdá se, že by to mohlo být 100 mm.

_________________
Cui bono ?
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu
PMA
Pavel Macura


Založený: 18 január 2007
Príspevky: 19729

PríspevokZaslal: Ut september 22, 2020 19:47:09    Predmet: Odpovedať s citátom

No já tam dal mraky bodů a nepomohlo to.
_________________
Dispre2-JFET preamp, Z-30_2020 power amp, CNO-T25 repro
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
miero
Hifi inventar


Založený: 08 september 2010
Príspevky: 9926
Bydlisko: Praha

PríspevokZaslal: Ut september 22, 2020 20:02:54    Predmet: Odpovedať s citátom

opa, bez zaruky

(EDIT) tych 99 bodov zo suboru som interpoloval na 1996 bodov (po 10Hz)

a este prilozim skript pre Octave/Matlab, keby to chcel niekto pozriet a pripadne zreprodukovat/upravit/opravit:
kód:
[a,b,c]=textread("FR_50_20k.txt","%f %f %f");
A=50:10:20000;
B=interp1(a,b,A,'spline');
l=length(A);
F=zeros(1,2*l);
F(2:l+1)=B;
F(l+2:2*l)=fliplr(B(1:end-1));
I=ifft(F);
w=fopen("imp2.txt","w");
fprintf(w,"%f\n",I(1:l));
fclose(w);
plot(1:l/10,I(1:l/10));
print imp2.pdf


Poznamka: horizontalna os ukazuje sample pri Fs 1995Hz, takze zhruba vydelit 2x pre ms.
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu
PMA
Pavel Macura


Založený: 18 január 2007
Príspevky: 19729

PríspevokZaslal: St september 23, 2020 09:34:17    Predmet: Odpovedať s citátom

Umí to VituixCad. Amplitudová charakteristika se naskenuje z grafu, přiřadí se v simulaci k měniči a program vypočte impulsní odezvu. Práce na 1 - 2 minuty.
_________________
Dispre2-JFET preamp, Z-30_2020 power amp, CNO-T25 repro
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
PMA
Pavel Macura


Založený: 18 január 2007
Príspevky: 19729

PríspevokZaslal: St september 23, 2020 10:10:38    Predmet: Odpovedať s citátom

OPA, takhle mi to vyšlo z tvých dat. Tu frekvenčku jsem jen posunul asi o 3dB dolů.

P.S.: impulsovku lze samozřejmě vyexportovat jako wav (příloha), ale IMO se z toho nic nezjistí ohledně odrazu.

_________________
Dispre2-JFET preamp, Z-30_2020 power amp, CNO-T25 repro


Naposledy upravil PMA dňa St september 23, 2020 16:35:31, celkom upravené 1 krát.
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
PMA
Pavel Macura


Založený: 18 január 2007
Príspevky: 19729

PríspevokZaslal: St september 23, 2020 10:30:19    Predmet: Odpovedať s citátom

Po velmi dlouhé době zajímavé téma, metoda se dá využít ke sledování vlivu impedance kabelu+reprobedny na temeno skokové odezvy zesilovače. Díky OPA za otevření tématu.
_________________
Dispre2-JFET preamp, Z-30_2020 power amp, CNO-T25 repro
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
PMA
Pavel Macura


Založený: 18 január 2007
Príspevky: 19729

PríspevokZaslal: St september 23, 2020 20:24:08    Predmet: Odpovedať s citátom

Behringer ECM8000. Nic moc.
_________________
Dispre2-JFET preamp, Z-30_2020 power amp, CNO-T25 repro
Návrat hore
Zobrazit informácie o autorovi Odoslať súkromnú správu Zobraziť autorové WWW stránky
Zobraziť príspevky z predchádzajúcich:   
Pridať novú tému   Zaslať odpoveď    Obsah fóra hifi.slovanet.sk -> Dokumentácia - How to ... Časy uvádzané v GMT + 1 hodina
Strana 1 z 1

 
Prejdi na:  
Nemôžete pridávať nové témy do tohto fóra.
Nemôžete odpovedať na témy v tomto fóre.
Nemôžete upravovať svoje príspevky v tomto fóre.
Nemôžete mazať svoje príspevky v tomto fóre.
Nemôžete hlasovať v tomto fóre.
Nemôžete pripojiť súbory do tohto fóra.
Nemôžete sťahovať súbory z tohto fóra.


Powered by phpBB © 2001, 2005 phpBB Group
Slovenský preklad phpBB Slovak - www.pcforum.sk
TOPlist