Witam.
Projektuję plugin VST jako przykład filtra cyfrowego.
Początkowo muszę sygnał podzielić na 3 odrębne pasma (LP, BP, HP) , przy czym pasmo jest zmiennym parametrem programu
Do obliczenia współczynników filtra zastosowałem wzory z Zielińskiego:
Np. Dla LP:
(*) h(n)=2*F0(sin(Ω0n)/Ω0n) -k<n<k
h(0)=2*F0
F0-częstotliwość graniczna unorm., Ω0-pulsacja unormowana
i przemnożyłem przez okno Kaisera.
Na wykresie amplitudowo-częstotliwościowym dostałem bardzo dobre wyniki -90 dB osiągam mniej więcej 20 % za zadanym F0 (filtr LP).
Ale wykres fazowy wygląda jak z sejsmografu w Chile w 1960r.
Napisali w "Zielińskim" że aby uzyskać liniową fazę potrzeba symetrii odpowiedzi impulsowej. Z założenia (*) jest symetryczna, okno Kaissera również. Tylko czy symetria względem środka (k) np: h(k-1)=h(k+1) ?? czy taka h(k)=h(k+1) ??
I czy to w ogóle prawda że symetria h(n) jest gwarantem linowej fazy ??
Z góry dzięki za pomoc.
[/u]
Projektuję plugin VST jako przykład filtra cyfrowego.
Początkowo muszę sygnał podzielić na 3 odrębne pasma (LP, BP, HP) , przy czym pasmo jest zmiennym parametrem programu
Do obliczenia współczynników filtra zastosowałem wzory z Zielińskiego:
Np. Dla LP:
(*) h(n)=2*F0(sin(Ω0n)/Ω0n) -k<n<k
h(0)=2*F0
F0-częstotliwość graniczna unorm., Ω0-pulsacja unormowana
i przemnożyłem przez okno Kaisera.
Na wykresie amplitudowo-częstotliwościowym dostałem bardzo dobre wyniki -90 dB osiągam mniej więcej 20 % za zadanym F0 (filtr LP).
Ale wykres fazowy wygląda jak z sejsmografu w Chile w 1960r.
Napisali w "Zielińskim" że aby uzyskać liniową fazę potrzeba symetrii odpowiedzi impulsowej. Z założenia (*) jest symetryczna, okno Kaissera również. Tylko czy symetria względem środka (k) np: h(k-1)=h(k+1) ?? czy taka h(k)=h(k+1) ??
I czy to w ogóle prawda że symetria h(n) jest gwarantem linowej fazy ??
Z góry dzięki za pomoc.
[/u]