|
hifi.slovanet.sk
Bolo zavedene manualne potvrdzovanie registracii !!!
|
|
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: 11183 Bydlisko: Praha
|
Zaslal: Ut september 22, 2020 13:07:56 Predmet: Impulsní odezva z frekvenčního průběhu |
|
|
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 |
|
|
miero Hifi inventar
Založený: 08 september 2010 Príspevky: 11533 Bydlisko: Praha
|
Zaslal: Ut september 22, 2020 13:29:24 Predmet: |
|
|
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 |
|
|
BV Hifi inventar
Založený: 06 október 2006 Príspevky: 5949 Bydlisko: Trnava
|
Zaslal: Ut september 22, 2020 15:20:33 Predmet: |
|
|
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 |
|
|
opa Hifi inventar
Založený: 24 február 2007 Príspevky: 11183 Bydlisko: Praha
|
Zaslal: Ut september 22, 2020 17:32:46 Predmet: |
|
|
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 |
|
|
godz Hifi expert
Založený: 05 máj 2007 Príspevky: 1077 Bydlisko: Trenčín
|
Zaslal: Ut september 22, 2020 18:12:07 Predmet: |
|
|
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 |
|
|
BV Hifi inventar
Založený: 06 október 2006 Príspevky: 5949 Bydlisko: Trnava
|
Zaslal: Ut september 22, 2020 18:45:23 Predmet: |
|
|
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 |
|
|
PMA Pavel Macura
Založený: 18 január 2007 Príspevky: 20466
|
Zaslal: Ut september 22, 2020 18:45:39 Predmet: |
|
|
100 bodu, není to málo, milý Watsone? _________________ mám rád hifi, ale víte, jak je to dál ..... |
|
Návrat hore |
|
|
PMA Pavel Macura
Založený: 18 január 2007 Príspevky: 20466
|
Zaslal: Ut september 22, 2020 19:11:58 Predmet: |
|
|
Vypadá to, že to REW neumí. Zkuste to ještě někdo, z přiložené FCH - korektní data i s fází. _________________ mám rád hifi, ale víte, jak je to dál ..... |
|
Návrat hore |
|
|
opa Hifi inventar
Založený: 24 február 2007 Príspevky: 11183 Bydlisko: Praha
|
Zaslal: Ut september 22, 2020 19:44:37 Predmet: |
|
|
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 |
|
|
PMA Pavel Macura
Založený: 18 január 2007 Príspevky: 20466
|
Zaslal: Ut september 22, 2020 19:47:09 Predmet: |
|
|
No já tam dal mraky bodů a nepomohlo to. _________________ mám rád hifi, ale víte, jak je to dál ..... |
|
Návrat hore |
|
|
miero Hifi inventar
Založený: 08 september 2010 Príspevky: 11533 Bydlisko: Praha
|
Zaslal: Ut september 22, 2020 20:02:54 Predmet: |
|
|
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 |
|
|
PMA Pavel Macura
Založený: 18 január 2007 Príspevky: 20466
|
Zaslal: St september 23, 2020 09:34:17 Predmet: |
|
|
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. _________________ mám rád hifi, ale víte, jak je to dál ..... |
|
Návrat hore |
|
|
PMA Pavel Macura
Založený: 18 január 2007 Príspevky: 20466
|
Zaslal: St september 23, 2020 10:10:38 Predmet: |
|
|
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. _________________ mám rád hifi, ale víte, jak je to dál .....
Naposledy upravil PMA dňa St september 23, 2020 16:35:31, celkom upravené 1 krát. |
|
Návrat hore |
|
|
PMA Pavel Macura
Založený: 18 január 2007 Príspevky: 20466
|
Zaslal: St september 23, 2020 10:30:19 Predmet: |
|
|
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. _________________ mám rád hifi, ale víte, jak je to dál ..... |
|
Návrat hore |
|
|
PMA Pavel Macura
Založený: 18 január 2007 Príspevky: 20466
|
Zaslal: St september 23, 2020 20:24:08 Predmet: |
|
|
Behringer ECM8000. Nic moc. _________________ mám rád hifi, ale víte, jak je to dál ..... |
|
Návrat hore |
|
|
opa Hifi inventar
Založený: 24 február 2007 Príspevky: 11183 Bydlisko: Praha
|
Zaslal: St september 01, 2021 20:29:04 Predmet: |
|
|
Tak jsem si počkal a REW se to mezitím naučil.
Funkce se jmenuje "Generate minimum phase"
Funguje, ale evidentně je příliš málo bodů. Nevadí, odměřil jsem to v reálu s držákem a bez držáku.
Zvlnění je +/- 0,5 dB, takže to dál nemusím řešit.
Nicméně děkuji všem za pomoc. _________________ Cui bono ? |
|
Návrat hore |
|
|
|
|
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.
|
|