Hangi IIR filtreleri bir Gauss filtresine yakındır?


15

Bu yüzden son zamanlarda Bessel filtrelerinin, diğer yaygın türlerle birlikte listelenmesine rağmen, gerçekten farklı bir "sınıf" a ait tuhaf bir şey olduğu konusunda şaşkına döndüm ve bu konuda daha fazla bilgi edinmeye çalışıyorum.

Dikdörtgen büyüklük yanıtı, ideal bant etki alanı tepkisini temsil eder, çünkü geçiş bandı sıfırdır ve durdurma bandı sonsuz zayıflamaya sahiptir. Diğer yandan, Gauss büyüklüğü tepkisi, dürtü tepkisinde ve adım tepkisinde herhangi bir aşma meydana gelmemesi için ideal zaman alanı tepkisini temsil eder. Uygulamada elde edilen yanıtların çoğu, bu ideal olanlara yaklaşımlardır kaynak

Yani bir brickwall filtresi , samimi bir işleve sahip bir evrişimdir ve bu frekans etki alanı özelliklerine sahiptir:

  • Düz geçiş bandı
  • Sıfır durdurma bandı
  • Sonsuz hareket hızı / geçiş bandı yok

Her iki yöndeki sonsuz kuyruklar nedeniyle nedensel değildir ve gerçekleştirilemez. Bu IIR filtreleri ile yaklaşık olarak düzenlenir ve sipariş arttıkça yaklaşım artar:

  • Butterworth (maksimum düz geçiş bandı)
  • Chebyshev (durdurma bandı veya geçiş bandı dalgalanması ile maksimum kalkış hızı )
  • Eliptik (durdurma bandı ve geçiş bandı dalgalanması ile maksimum kalkış hızı )
  • Legendre (monotonik geçiş bandı ile maksimum kalkış hızı)

Sipariş arttıkça brickwall yanıtına yaklaşan 4 filtre türü


Gauss filtre bir Gauss fonksiyonu ile konvolüsyon olduğunu ve bu zaman alanı özelliklere sahiptir:

  • Sıfır aşma
  • Minimum yükselme ve düşme süresi
  • Minimum grup gecikmesi

Sinc fonksiyonu ile aynı nedenlerle gerçekleştirilemez ve bu IIR filtreleriyle sipariş arttıkça daha yakın bir şekilde tahmin edilebilir:

  • 1 ve 2'ye göre Bessel (maksimum düz grup gecikmesi)

e12(πω)2

Bessel, sipariş arttıkça Gauss yanıtına yaklaştığını ve Gauss için tahmin ettiğini filtreledi

Yani sorularım:

Her şey yolunda mı? Öyleyse, Gaussian'a yaklaşan başka IIR filtreleri var mı? Ne için optimize edilmişler? Belki de aşımı en aza indiren biri?

"IIR Gaussian" için arama yaparsanız birkaç şey bulabilirsiniz (Deriche? Van Vliet?), Ama gerçekten bir Bessel ile aynı olup olmadıklarını veya başka bir mülk için optimize olup olmadıklarını bilmiyorum.


aslında, "IIR" dediğimde gerçekten "fiziksel olarak gerçekleştirilebilir analog filtreler" demek istediğimi düşünüyorum?
endolith

Yanıtlar:


5

Deriche ve van Vliet filtreleri sezgiselleri vardır. Her iki durumda da, RMS farkını veya bir Gaussian'dan gelen filtrenin dürtü yanıtının maksimum farkını en aza indirmek için kutupların ve sıfırların yerlerini seçerler.

Her iki filtre de nedensel-nedensel çiftlerdir. Bu yüzden faz hatası veya grup gecikmesi olmadığını düşünüyorum, ancak bunları verilerde olduğu kadar ileriye doğru da çalıştırabilmeniz gerekir. Bu onları görüntü işlemede popüler kılar, ancak belki de başka yerlerde uygulanabilirliklerini kısıtlar.

Sezgisel oldukları, onları değiştiren kağıtların zenginliği ile kanıtlanmıştır. Örneğin, bir google araması (Deriche gazetesine bağlantı ararken) , Deriche türevi filtrenin tam 0 DC yanıtı olmadığı sorununu çözmeye çalışan bir arama çıktı . Sınır koşullarının doğru bir şekilde başlatılmasıyla ilgili bazı ilginç konular da vardır .

Aşağıdaki genel bakışı iyi bir kaynak olarak buldum: Dave Hale, Özyinelemeli Gaussian ters lters , Colorado Dalga Madenleri Okulu Maden Merkezi raporu CWP-546.


3

Sanırım ayrık zamanlı IIR filtreleri için mevcut analitik çözümlerin güzel bir özetini verdiniz. Ama aynı zamanda ideal frekans seçici filtre özelliklerine yaklaşan filtre listesine Bessel filtreleri de eklerdim. Büyüklük yanıtı, aynı sıradaki diğer filtre türleri kadar keskin bir geçiş göstermez, ancak bu, geçiş bandındaki neredeyse doğrusal bir faz için ödediğiniz fiyattır. Bu nedenle Bessel filtresi, frekans seçici büyüklük yanıtı ile iyi bir faz yanıtı arasında bir uzlaşmadır.

Bir Gauss filtresini IIR filtrelerle yaklaşıklaştırmak için, bahsettiğiniz Bessel filtresinin dışında herhangi bir analitik çözüm bilmiyorum. Ancak Bessel filtresinin bir Gauss filtresine yakın olması gerektiğini unutmayın, bu yüzden böyle bir filtreye yaklaşmada gerçekten ne kadar iyi olduğundan emin değilim. Bu amaç için gerçekten bir IIR filtresi istiyorsanız, Gauss filtresinin sayısal yaklaşımını tercih etmenizi öneririm. Bunu nasıl yapacağınız için birkaç seçenek vardır.

Frekans alanında Gauss filtresine yaklaşmayı deneyebilirsiniz. Sorun, istenen faz yanıtı ile ilgili bazı kararlar vermeniz gerektiğidir. Minimum faz yanıtı ile saf büyüklükte bir yaklaşım, büyük olasılıkla çok zayıf zaman alanı özellikleri ile sonuçlanacaktır. İstenen doğrusal bir faz belirtirseniz, karmaşık bir yaklaşım problemi elde edersiniz (çünkü karmaşık frekans cevabını büyüklük ve faz ile yaklaşık olarak tahmin edersiniz). Böyle bir yaklaşım probleminin çözülmesi oldukça zor olsa da, literatürde yöntemler mevcuttur.

Daha basit ve muhtemelen daha iyi bir yaklaşım, zaman alanındaki Gauss filtresine yaklaşmaktır. Prony'nin yöntemi iyi bir başlangıç ​​noktası olacaktır.

Lütfen bunların sadece konu hakkındaki düşüncelerim olduğunu unutmayın. Ben bir IIR Gauss filtresi kendim tasarlamaya çalışmadım. Buna karşı çok iyi nedenler olmadıkça aslında FIR uygulaması için giderdim.

DÜZENLEME: Bessel filtresinin bir Gaussian'a yaklaşıp yaklaşmadığı sorusuyla ilgili birkaç açıklama daha. Bessel filtresinin bir Gauss filtresine yaklaşırken en aza indirdiği anlamlı bir hata kriteri bilmiyorum. Yine de bunu öğrenmekten memnuniyet duyarım. İnsanlar bir Bessel filtresinin dürtü yanıtının Gaussian'a benzediğini veya frekans cevabının bir Gaussian'a benzediğini iddia edebilir, ancak Bessel'in herhangi bir anlamda bir Gauss filtresine yaklaştığını ve yaklaşık hataların gittiğini gösteren bir kanıt görmedim. filtre sırası arttıkça sıfır olur. Bir Gauss'a diğer standart filtrelerden (Butterworth, Chebyshev, vb.) Daha benzer olduğunu inkar etmiyorum, ancak bu soru için önemli değil.

Octave'de (işlev besself) tasarlanan Bessel filtrelerinin (sipariş 5, 10, 15, 20) dört dürtü tepkisine bakınız. Gördüğünüz gibi, kuyruktaki zil sesi artan filtre sırası ile azalmıyor ve bu filtrelerin Gauss'a nasıl yaklaştığını ve eğer öyleyse hangi iyimserlik ölçütüne göre göremiyorum. Ancak, kimse beni bu konuda aydınlatabilirse, çok mutlu olurum.

resim açıklamasını buraya girin


Bu gibi yerlerden alıyorum: "Bessel-Thomson filtrelerinin dürtü yanıtı, filtre sırası arttıkça Gauss'a yöneliyor " robots.ox.ac.uk/~sjrob/Teaching/SP/l3.pdf " analog Bessel filtresi bir Gauss filtresine yaklaşık değerdir ve filtre sırası arttıkça yaklaşık değer artar. " dsprelated.com/showmessage/130958/1.php
Endolit

... ve bir Gauss'un Fourier dönüşümü bir Gausscu olduğu için, diğerleri gibi bir tuğla duvar tepkisine yaklaştığını söylemenin doğru olduğunu düşünmüyorum.
endolith

Bir Gauss'un Fourier dönüşümü bir Gauss'dur, bundan kuşkusuz. Ama bildiğim kadarıyla bir Gaussian'a anlamlı bir şekilde yaklaşmayan Bessel filtrelerinden bahsediyoruz. Daha fazla bilgi eklemek için cevabımı düzenledim.
Matt L.

Bununla ilgili soruya biraz daha ayrıntı ekledik. Grafiklerinizi saplar yerine çizgi eğrileri olarak ve aynı Y ekseniyle yeniden yapabilir misiniz? Sırayla şekil değiştiriyor gibi görünüyor ama söylemek zor.
endolith

1
Ben 60 siparişe kadar bessel filtreleri hesaplamaya çalıştı ve gaussians tarafından dürtü yanıtlarını uydurma ve benim kod ham olmasına rağmen, onlar "gahoş" azalan ve sipariş azalan hata ile gaussian yaklaşıyor gibi görünüyor. 10. sıra: imgur.com/1qNsHeg 60. sıra: imgur.com/BgmFzZp sipariş arttıkça hata: imgur.com/cpHDDJs hesaplamalar yanlış olabilir. seninkini nasıl hesapladın
endolith
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.