Filtreleme ve polinom regresyon yumuşatma arasındaki farklar?


14

Klasik düşük geçişli filtreleme (bir IIR veya FIR ile) ve özellikle N'nin 1'den büyük olduğu durumlarda, yerelleştirilmiş N derece polinom regresyonu ve / veya enterpolasyon (yukarı örnekleme durumunda) ile "yumuşatma" arasındaki farklar nelerdir? ancak regresyonda kullanılan yerel nokta sayısından az.


+1 Harika bir soru, beni dövüyorsun. :-) N = 2 kullanan AFAIK, bildiğimiz doğrusal 'klasik' filtrelemeye karşılık gelir, ancak bu konuda yanlış olabilirim.
Spacey

2
iç rekonstrüksiyon vs spline enterpolasyonu: cnx.org/content/m11126/latest "spline enterpolasyonu iç enterpolasyondan daha pürüzsüzdür. Bunun nedeni, kardinal spline'ların desteğinin sinc fonksiyonundan daha kompakt olmasıdır."
endolit

Yanıtlar:


10

Hem düşük geçişli filtreleme hem de polinom regresyon düzleştirmesi bir fonksiyonun yaklaşık değerleri olarak görülebilir . Ancak, bunu yapmanın yolları farklıdır. Burada sorulması gereken kilit soru "Birini diğeri açısından yapabilir misin?" ve kısa cevap aşağıda açıklanan nedenlerden ötürü "her zaman değil" dir.

Filtreleme ile yumuşatma yaparken , frekans alanında burada , Ayrık Fourier Dönüşümü'nü (ve tersine ) belirtir . Ayrık Fourier Dönüşümü (örn. ), trigonometrik fonksiyonların toplamı olarak yaklaşık bir değerini sunar . Zaman bir düşük geçiş filtresi, düşük frekans bileşenlerinin daha az sayıda muhafaza edilir ve ani değişikliklerin düzeltilmesini. Bu , temel fonksiyonlar olarak trigonometrik fonksiyonlar kullanarak fonksiyon yaklaşımı bağlamında düşük geçişli filtrelemeyi ayarlary(n)=x(n)h(n)y=F1(F(x)F(h))FF1F(x)xhxancak, filtreleme sırasında, y (n) (filtrenin çıktısı) ve ayrıca geçmiş numunelerinin ağırlıklı toplamına (burada ağırlıklandırma tarafından belirlenen ağırlıklandırma ) bağlı olduğuna dikkat etmek için kıvrım formülünü tekrar gözden geçirmeye değer . "şekli" . (elbette IIR filtreleri için de nin geçmiş değerlerinin eklenmesi ile benzer hususlar geçerlidir )x(n)xhy(n)

Yine de bir n-dereceli polinom ile yumuşatırken , interpolantın çıktısı sadece ve (farklı) temel fonksiyonların bir karışımına ( monomiyaller olarak da adlandırılır ) bağlıdır. Bu farklı temel işlevler nelerdir? Bu bir sabiti (var ), bir çizgi ( ), bir parabol ( ) ve benzeri (bakınız Buna güzel bir örnek için). Genellikle, zaman içinde eşit uzaklıktaki örneklerle uğraşırken ve doğrulukla ilgili nedenlerle, Newton'un polinom formu kullanılır.x(n)a0x0a1xa2x2. Bunu belirtmemin nedeni, doğrusal enterpolasyon gerçekleştirirken, düşük dereceli enterpolasyon polinomunun enterpolasyon yapmak için "satırları" kullanması gibi, mevcut örneklerin lineer ağırlıklı toplamını döndüren bir filtre çekirdeği oluşturabileceğinizi görmek kolay olmasıdır. iki numune arasında. Ancak daha yüksek derecelerde, iki yaklaşım yöntemi farklı sonuçlar döndürür (temel işlevlerdeki farklılıklar nedeniyle).

Yukarıda yazdığım gibi, nin geçmiş değerlerini hesaba katmamak katı değildir. Bu ince bir nokta. Çünkü genellikle, bir polinom oluşturulurken verilen aralığın dışındaki değerler (bir sinyalin "geçmiş" ve "geleceği") dikkate alınmaz. Bununla birlikte, türevleri aralığın kenarlarına sabitleyerek bunları dahil etmek mümkündür. Ve eğer bu tekrar tekrar yapılıyorsa (örtüşmeyen bir kayar pencere gibi) etkili bir şekilde, x (n) 'nin "geçmiş örnekleri" dikkate alınacaktır. (Bu spline'ların kullandığı püf noktasıdır ve aslında bisubik enterpolasyon için bir evrişim ifadesi vardır . Bununla birlikte, lütfen yorumunun spline hakkında konuşurken farklı olduğunu unutmayın - normalizasyon hakkındaki noktayı not edin).xx(n)x

Filtrelemenin bazı zamanlarda enterpolasyon olarak kullanılmasının nedeni, örneğin "Sinc İnterpolasyonu" durumunda, fiziksel açıdan da mantıklı olmasıdır. Bant sınırlı bir sistemin (örn. (Doğrusal) bir amplifikatör veya bir optik sistemdeki mercek ) zaman alanında ideal temsili , samimi nabızdır. Bir nabız sinyalinin frekans alanı gösterimi bir dikdörtgen "darbe" dir.. Bu nedenle, çok az varsayımla eksik bir değerin komşularına az çok yakın olmasını bekliyoruz (elbette sınırlar dahilinde). Bu, bazı n-sıralı polinomlarla (daha yüksek n için) gerçekleştirildiyse, eksik bir değerin her zaman gerçekçi olmayabilecek komşularıyla ilişkili olma yolunu "düzeltiriz" (neden ses basıncı değerleri mikrofona vurmak dalga cephesi örneğin şeklinde olacak şekilde sabitlenecek mi? Ses kaynağının her zaman doğru olmayabilecek şekilde nasıl davrandığına dair bir varsayım var. beynin işlenmesini içeren psikofizik bakış açısıyla şema (bkz. Lanczos yeniden örnekleme)x3Örneğin). Kesinlikle biri nesnel olarak eksik değerleri "tahmin etmeye" çalıştığında enterpolasyonun getirdiği kısıtlamalardan bahsediyorum.

Evrensel bir “en iyi yöntem” yoktur, karşılaştığınız enterpolasyon problemine bağlıdır.

Umarım bu yardımcı olur.

PS (İki yaklaştırma yönteminin her biri tarafından üretilen eserler de farklıdır, örneğin, Gibbs Fenomeni ve aşırı uydurma olmasına rağmen, aşırı sığdırma sorunuzun "diğer tarafında" olmasına rağmen).


+1 Mükemmel cevap. Bazı takipler: 1) Polinom uydurmada geçmiş x [n] değerlerini hesaba katmamadan bahsettiniz, ancak bu, x [n] 'nin sinüs / kosinüslerin bir toplamı olduğu hakkında söylediklerinize dayanarak tartışılan bir konu değil mi? (Geçmiş değerler hesaba katılmış olsun olmasın, bu hala geçerlidir). 2) Bu durumda 'bant sınırlı' bir şeyin fiziksel yorumu ile biraz kafam karıştı. Her şey bantla sınırlı değil mi? Yani, belirli frekansları geçecek ve diğerlerini zayıflatacak mı? Sınırsız sisteme fiziksel örnek nedir? Teşekkürler.
Spacey

1) Ne demek istediğinizi tam olarak anladığımdan emin değilim ama kıvrımdan ve polinom fittinginden çıktı elde etmek arasındaki farklara atıfta bulunuyordum. 2) Bazı durumlarda sinyaller ve sistemler aynı çerçevede ele alınır. Teorik olarak (gerçekten) beyaz gürültü ( en.wikipedia.org/wiki/White_noise ) gibi bant sınırlı olmayan ( en.wikipedia.org/wiki/… ) sinyaller vardır . Oppenheim ve Willsky tarafından Signals & Systems'da çok iyi bir tedavi mevcuttur. Burada terimi,
bandlimit-

Tamam, sorumu yeniden yazdım - sadece emin olmak için: 1) Ne kadar yüksek mertebeden polinom kullanırsak, fiziksel gerçeklikle uyuşmayabilecek noktalar arasındaki ilişkileri o kadar 'taraflı' tutarız, değil mi? (Bu durumda daha fazlası her zaman daha iyi değildir.) 2) Bant sınırlama ile ilgili olarak - Bunu neden söylediğimizi merak ediyorum, çünkü HER sistem bandı sınırlı değil, çünkü sadece belirli frekansları alıyor ve diğerlerini zayıflatıyor mu? Teşekkürler.
Spacey

Üzgünüm, bu benim dikkatimi çekti. Bu özel sorular için: 1) Gerekli değil. Verilen örnekte, monomialerin "şekli" tarafından getirilen kısıtlamalardan bahsetmiştim. 2) Sinyaller ve Sistemler çok yardımcı olacaktır. Mühendislik uygulamalarının, başka bir alanda bant dışı sınırlı sinyaller için çok iyi bir kullanımı olabilen bir matematik alt kümesi kullandığından bazı şeylerin kesin olduğu söylenir (örneğin, birbirine bağlanmış gerçekten tek tip rastgele süreç (beyaz gürültü)).
A_A

2

Güzel soru ve aydınlatıcı cevaplar. Birkaç içgörü şöyle paylaşmak istedim. Yüksek dereceli polinomların yerleştirilmesinde daha kararlı olan Legendre'nin polinom bazları (monomiyal bazların aksine) gibi dik polinom bazları da vardır. Shannon'ın enterpolasyon formülünde kullanılan iç tabanlar (aslında bir evrişim operasyonu ve dolayısıyla bir filtreleme operasyonu olarak da görülebilir) bant sınırlı Hilbert alanı için dik tabanlar olduğundan, ortogonal polinom bazlar band sınırında olmayan daha büyük bir fonksiyon sınıfına yaklaşmaya hizmet edebilir onlarla birlikte ortogonallik gücüne sahip olmak.

Polinom filtreleme (enterpolasyon değil) 1960 yılından beri Kimya literatüründe de yer almaktadır . eecs.berkeley.edu/~ee123/fa12/docs/SGFilter.pdf

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.