Neden evrişim gereklidir ya da evrişimin arkasındaki felsefe nedir?


15

Dijital görüntü restorasyonu alanında çalışıyorum. Evrişim ile ilgili her şeyi okudum, bir LTI sistemi için dürtü tepkisini bilirsek , çıktısını sadece girdi ve dürtü yanıtı arasında evrişim kullanarak bulabiliriz.

Birisi bana arkasındaki ana matematiksel felsefenin ne olduğunu söyleyebilir mi? Deneyiminiz bana sadece internette sörf yapmaktan daha fazlasını anlatacak.



3
Bu soruyu reddediyorum çünkü (veya küçük varyasyonları) tekrar tekrar bu sitede pichenettes'in yorum durumları olarak soruldu ve cevaplandı. Bunun yerine bu sitede "internet surfed" olmalıdır.
Dilip Sarwate

Yanıtlar:


14

Evrişim Düşüncesi

Konuyla ilgili en sevdiğim sergim Brad Osgood'un Fourier Dönüşümü üzerine verdiği derslerde . Evrişim tartışması 36: 00'da başlar, ancak tüm dersin izlemeye değer ek bir bağlamı vardır.

Temel fikir, Fourier Dönüşümü gibi bir şeyi her zaman doğrudan tanımla çalışmak yerine tanımladığınızda, hesaplamaları basitleştiren daha üst düzey özellikler türetmenin yararlı olmasıdır. Örneğin, böyle bir özellik, iki fonksiyonun toplamının dönüşümünün, dönüşümlerin toplamına eşit olması, yani

F{f+g}=F{f}+F{g}.

Bu, bilinmeyen bir dönüşüme sahip bir işleviniz varsa ve bilinen dönüştürmelere sahip işlevlerin toplamı olarak ayrılabilirse, temel olarak cevabı ücretsiz olarak alırsınız.

Şimdi, iki dönüşümün toplamı için bir kimliğimiz olduğundan, iki dönüşümün çarpımı için kimliğin ne olduğunu sormak doğal bir sorudur, yani

F{f}F{g}= ?.

Cevabı hesapladığınızda, evrişimin ortaya çıktığı ortaya çıkıyor. Tüm türetme videoda verilir ve sorunuz çoğunlukla kavramsal olduğundan, burada tekrar özetlemeyeceğim.

Bu şekilde evrişime yaklaşmanın anlamı, bunun Laplace Dönüşümünün (Fourier Tranform'un özel bir durum olduğu) lineer sabit katsayılı sıradan diferansiyel denklemleri (LCCODE) cebirsel denklemlere dönüştürmesinin temel bir parçası olmasıdır. LCCODE'nin analitik olarak izlenebilir hale getirilmesi için böyle bir dönüşümün mevcut olması, sinyal işlemede incelenmesinin nedeninin büyük bir parçasıdır. Örneğin, Oppenheim ve Schafer'ı alıntılamak için :

Matematiksel olarak karakterize edilmeleri nispeten kolay oldukları ve yararlı sinyal işleme işlevleri gerçekleştirmek üzere tasarlanabildikleri için, doğrusal kaydırma-değişmez sistemlerin sınıfı kapsamlı bir şekilde incelenecektir.

Sorunun bir cevabı, LTI sistemlerini analiz etmek ve / veya sentezlemek için dönüşüm yöntemlerini kullanıyorsanız, er ya da geç, kıvrım ortaya çıkacaktır (örtülü ya da açıkça). Evrişim girişine yönelik bu yaklaşımın diferansiyel denklemler bağlamında çok standart olduğunu unutmayın. Örneğin, Arthur Mattuck'ın bu MIT dersine bakın . Çoğu sunum ya evrişim integralini yorum yapmadan sunar, daha sonra özelliklerini türetir (yani bir şapkadan çıkarır) ya da integralin garip formu hakkında bilgi verir, ters çevirme ve sürükleme, zaman tersine çevirme vb. .

Prof. Osgood'un yaklaşımını sevmemin sebebi, benim görüşüme göre matematikçilerin ilk başta bu fikre nasıl geldiğine dair derin bir kavrayış sağlamanın yanı sıra tsouris'lerden de kaçınmasıdır. Ve alıntı yapıyorum:

Dedim ki, "F ve G'yi zaman alanında birleştirmenin bir yolu var mı, böylece frekans alanında spektrumlar çoğalıyor, Fourier dönüşümleri çoğalıyor mu?" Ve cevap, evet, bu karmaşık integral tarafından var. Çok açık değil. Sabah yataktan kalkmaz ve bunu yazmazsınız ve bunun bu sorunu çözeceğini düşünürsünüz. Nasıl elde ederiz? Diyelim ki, sorun çözüldü, ne olması gerektiğini görün ve sonra zaferi ilan etmenin zamanı geldiğini fark etmeliyiz. Ve zafer ilan etme zamanı.

Şimdi, iğrenç bir matematikçi olarak, izlerinizi kapatıyorsunuz ve "Şey, ben sadece bu fonksiyonla iki fonksiyonun evrimi tanımlayacağım" diyorsunuz.

LTI Sistemleri

Çoğu DSP metninde, evrişim genellikle farklı bir yoldan getirilir (bu, dönüştürme yöntemlerine herhangi bir göndermeden kaçınır). Ölçekli ve kaydırılmış birim impulsların toplamı olarak rastgele bir giriş sinyali ifade ederek ,x(n)

(1)x(n)=k=x(k)δ(nk),

nerede

(2)δ(n)={0,n01,n=0,

doğrusal zamanla değişmeyen sistemlerin tanımlayıcı özellikleri , dürtü yanıtı içeren bir evrişim toplamına yol açar . Bir LTI operatör tarafından tanımlanan sistem ise L olarak ifade edilir y ( n ) = L [ x ( n ) ] , daha sonra repective özellikleri uygulayarak, yani doğrusalh(n)=L[ δ(n) ]Ly(n)=L[ x(n) ]

(3)L[ ax1(n)+bx2(n) ]Transform of the sum of scaled inputs=aL[ x1(n) ]+bL[ x2(n) ]Sum of scaled transforms,

ve zaman / vardiya değişmezliği

(4)L[ x(n) ]=y(n) impliesL[ x(nk) ]=y(nk),

sistem şu şekilde yeniden yazılabilir

y(n)=L[k=x(k)δ(nk)]Tranform of the sum of scaled inputs=k=x(k)L[δ(nk)]Sum of scaled transforms=k=x(k)h(nk).Convolution with the impulse response

Bu, evrişimi sunmanın çok standart bir yoludur ve bunun için mükemmel zarif ve kullanışlı bir yöntemdir. Benzer türevler Oppenheim ve Schafer , Proakis ve Manolakis , Rabiner ve Gold'da bulunabilir ve eminim diğerleri. [Standart tanıtımları daha ileri gider] Bazı derin kavrayış onun mükemmel cevap Dilip verilir burada .

(1)

(fg)(n)f convolved with g=k=f(k)g(nk),

(1)xδnn+0n×1

Bu perspektiften, evrişim özünde bir delta işlevi fikri ile ilişkilidir (yani, delta işlevinin kimlik öğesi olarak sahip olduğu bir ikili işlemdir). Evrişimle olan ilişkisini dikkate almadan bile, sinyalin tanımı büyük ölçüde delta işlevi fikrine bağlıdır. Öyleyse soru daha sonra, delta işlevi için fikri nereden edindik? Anlayabildiğim kadarıyla, en azından Fourier'in Analitik Isı Teorisi üzerine, dolaylı olarak göründüğü kadar geri gidiyor . Daha fazla bilgi için kaynaklardan biri Alejandro Domínguez'ın Kökeni ve Konvolüsyon Tarihi üzerine yazdığı makaledir.

Şimdi, bunlar doğrusal sistem teorisi bağlamında fikre ana yaklaşımlardan ikisi. Biri analitik içgörü, diğeri sayısal çözümü destekler. Her ikisinin de evrişimin öneminin tam bir resmi için yararlı olduğunu düşünüyorum. Bununla birlikte, ayrık durumda, doğrusal sistemleri tamamen ihmal ederek, evrişimin çok daha eski bir fikir olduğu duygusu vardır.

Polinom Çarpımı

Ayrık katlama olduğu fikrine biri iyi sunum sadece polinom çarpma içinde kıtabının tarafından verilen bu ders etrafında 5:46 başlayarak. Bu perspektiften, fikir, konumsal sayı sistemlerinin (sayıları dolaylı olarak polinomlar olarak temsil eden) tanıtılmasına kadar uzanır. Z-dönüşümü sinyalleri z cinsinden polinomlar olarak temsil ettiğinden, bu bağlamda da kıvrım ortaya çıkacaktır - Z-dönüşümü karmaşık analizlere başvurmadan resmi olarak bir gecikme operatörü olarak ve / veya Laplace'ın özel bir durumu olarak tanımlansa bile Dönüştür .


çok değerli rehberiniz için teşekkürler efendim, takip etmem için bana doğru yolu gösterdiniz. Bu yardımın bana başkaları için iyi bir insan olmanın nasıl başladığını öğretti .... :)
Mayank Tiwari

Bu büyük tesadüf, sorulan davadaki evrişimi yapmanız gerektiğini nasıl açıklıyor? Her alanda, argümanları zaman alanına geri döndürdüğünüzde evrişime dönüşen bir işlem olduğuna inanıyorum. Yanıt almak için zaman alanında muiltiplication yapmamız gerekebilir mi? Neden zaman taramaları yerine frekansları çarpalım?
Val

1
OP'nin LTI sistemlerine göre dürtülerin rolü hakkında bir soru sorduğunu göz önünde bulundurarak, konvülsiyonun nereden geldiğiyle ilgili bir soruyu motive etmek için bunu örnek olarak kullanarak soruyu okudum - bir dürtü hesaplamadaki rolü değil kendi başına cevap. Sorduğun şey bu mu?
datageist

1
Evrişime ihtiyacımız olduğunu söylemek, Fourier çarpımıyla özdeş olduğunu söylemek, Fourier çarpımına neden ihtiyacımız olduğunu bilmediğimizde bana saçma geliyor. Dürtü yanıtı verildiğinde, bu, herhangi bir kara büyü yerine fourier temelinde zaman alanı ve evrişim anlamına gelir. Bu soruya yapılan göndermenin bunu netleştirebileceğini düşünmüyorum. Her durumda, genel, temel (yani felsefi) sorulara “yerelleştirilmiş cevaplar” vermek iyi değildir. Soru-Cevap, gelecekteki ziyaretçiler için yararlı olmalıdır.
Val

Val'in yukarıdaki yorumu doğru. Fourier dönüşümleri icat edilmeden / keşfedilmeden önce doğrusal sistemlerin nasıl çalıştığını merak ediyorum. Duygusal olmayan cansız bir nesne nasıl bu kadar karmaşık bir formül keşfetti?
Dilip Sarwate

6

xt. Doğrusal sistemler, x (t) giriş işlevi için yanıt işlevini y (t) elde etmek için bu giriş girişlerinin her biri için çıkışları toplayabileceğiniz kadar havalıdır. Çıkış darbesinin y (t = 10), h (0) * x (10) 'a katkıda bulunan acil giriş x (10)' a bağlı olduğunu biliyorsunuz. Ancak, önceki darbe, x (9) çıkışına x (9) * h (1) katkısı ve daha önceki giriş değerlerinden katkılar da vardır. Önceki girdilerden katkı ekledikçe, bir defalık argüman azalırken diğeri artar. Bu MAC y içine tüm katkıları (10) = h (0) * x (10) + h (1) * x (9) + h (2) * x (8) + ..., bir kıvrım olup.

Y (t), h (t) ve x (t) fonksiyonlarını vektör olarak düşünebilirsiniz. Matrisler lineer cebirdeki operatörlerdir. Girdi vektörü (bir sayı dizisi) alır ve çıktı vektörü (başka bir sayı dizisi) üretir. Bu durumda, y vektörü x olan bir evrişim matrisi ürünüdür,

y=[y0y1y2]=[h000h1h00h2h1h0][x0x1x2]=Hx

Şimdi, evrişim bir Toeplitz matrisi olduğu için, Fourier özüne sahiptir ve bu nedenle evrişim operatörü (lineer operatörler matrislerle temsil edilir, ancak matris de temele dayanır) fourier alanında güzel bir diyagonal matristir,

Y=[Y0Y1Y2]=[λ0000λ1000λ2][X0X1X2]=diag(H)X

Unutmayın, çok daha fazla sıfır ve dolayısıyla çok daha basit hesaplama. Bu sonuç "evrişim teoremi" olarak bilinir ve ilk yanıtın yanıtlandığı gibi, fourier alanında çok daha basittir. Ancak, bu, "evrişim teoreminin", dört katlı temelin ve doğrusal operatörlerin her yerde evrişim ihtiyacından ziyade felsefe.

y[currentTime]=k1x[time1]+k2x(time2)+by[time1]


Bir kaç not: Sürekli açık vaka için bu açıklamayı nasıl genişletirsiniz (açıkçası ayrık zamanlı davadan önce gelir)? Ayrıca, hızlı evrişim için Fourier dönüşümü tabanlı yöntemler kullanan birçok gerçek zamanlı uygulama vardır. Gerçek zamanlı uygulamalar için çıktıların her zaman birer birer hesaplandığını söylemek doğru değildir.
Jason R

Bununla birlikte, evrişim matrisinin Toeplitz yapısının Fourier temelinde çapraz bir temsili kabul ettiğini ima eden güzel bir iş.
Jason R

Evet, gerçek zamanlı olarak fourier transfrom kullanıyor olabilirsiniz. DSP uzmanı olmaktan çok uzakım. Sadece kıt uygulama ve DSPGuide okumaktan aldığım fikri ifade ettim. Her neyse, fourier'in evrim felsefesi ile ilgisi olmadığını vurgulamak istiyorum. Rahatsız edici olduğu için fourier ile ilgili tartışmanın tümünü kaldırmam gerekebilir. Konvolüsyon zaman alanında doğaldır ve Fourier ne kadar havalı olursa olsun, herhangi bir Fourier olmadan gereklidir.
Val

f(x)dxf(x)dx=limdx0(f(x)dx)
Val

@JasonR Sürekli ayarda, Toeplitz matrisini vardiya-değişmez bir operatörle değiştirirsiniz. Daha sonra Fourier temel fonksiyonlarının bu operatörü köşegenleştirdiğini gösterebilirsiniz.
lp251

4

Önceki cevaplar gerçekten iyi olmasına rağmen, rakamlar nedeniyle görselleştirmeyi daha kolay hale getirdiğim evrişim hakkındaki görüşümü eklemek istiyorum.

Belirli bir giriş sinyali için bir sistemin çıkış sinyalinin belirlenebileceği herhangi bir yöntem olup olmadığı merak edilir. Evrişim, sistemin doğrusal ve zamanla değişmeyen (LTI) olması şartıyla, bu sorunun cevabıdır.

Rastgele bir sinyalimiz olduğunu varsayın . Daha sonra, , aşağıdaki akıl yürütme yoluyla kaydırılmış birim impulsların ölçekli bir toplamına ayrıştırılabilir. M birimini, örnekleri tarafından olarak değiştirilmiş bir birim impuls ile çarpın . Yana 0'a eşit olduğu yerde de hariç bu olur ve çok kez tüm değerler , 0 ile eşit değildir ve 1 ile eşittir . Sonuçta elde edilen dizi, değerinde değerine eşit bir impulsa sahip olacaktır.s[n]s[n]s[n]mδ[nm]δ[nm]n=ms[n]nmnmn=ms[m]. Bu işlem aşağıdaki Şekil'de açıkça gösterilmiştir.

resim açıklamasını buraya girin

Gibi bu matematiksel yazılabilir farklı bir gecikme ile aynı prosedür tekrar edilmesi verir

s[n]δ[nm]=s[m]δ[nm]
m
s[n]δ[nm]=s[m]δ[nm]

Değeri , bu anda ekstre edilir. Bu nedenle, bu çarpma tüm olası gecikmelerde tekrarlanırsa ve üretilen tüm sinyaller birlikte toplanırsa, sonuç dizisinin kendisi olur. s[m]<m<s[n]

s[n]=+s[2]δ[n+2]+s[1]δ[n+1]+s[0]δ[n]+s[1]δ[n1]+s[2]δ[n2]+=m=s[m]δ[nm]

Özetle, yukarıdaki denklem, , her birim dürtü nin genliği olduğu ölçekli birim darbelerinin bir toplamı olarak yazılabileceğini belirtir . Böyle bir toplamın bir örneği aşağıdaki Şekil'de gösterilmiştir.s[n]δ[nm]s[m]

resim açıklamasını buraya girin

Darbe yanıtı olan bir LTI sistemine girdi olarak verildiğinde ne olacağını düşünün .h[n]

resim açıklamasını buraya girin

Bu giriş-çıkış sırasına yol açar.

resim açıklamasını buraya girin

Yukarıdaki prosedür sırasında, dürtü yanıtı olan bir LTI sistemine giriş için çıkışını tanımlayan ünlü evrişim denklemini çalıştık .r[n]s[n]h[n]

Evrişim çok mantıklı ve basit bir süreçtir, ancak birçok DSP öğrencisinin açıklama şekli nedeniyle kafa karıştırıcı bulabilir. Geleneksel bir yöntemi ve daha sezgisel bir yaklaşımı anlatacağız.

Geleneksel Yöntem


Evrişim denklemini tanımladıktan sonra çoğu kitap, aşağıdaki adımlar aracılığıyla uygulanmasını önerir. Her bireysel zaman kayması ,n

[Flip] Denklemi olarak düzenleyerek, dürtü yanıtını değişkeninin bir fonksiyonu olarak düşünün , elde etmek için yaklaşık çevirin .r[n]=m=s[m]h[m+n]mh[m]m=0h[m]

[Shift] elde etmek için zaman kayması için , kaydırma ile pozitif sağdan birim ve negatif için kalan . h[m+n]nh[m]nnn

[Çarpma] ürün sekansı elde etmek için sekansını sekansı ile çarpın . s[m]h[m+n]s[m]h[m+n]

[Toplama] zamanındaki evrişim çıktısını elde etmek için yukarıdaki ürün dizisinin tüm değerlerini toplayın .n

[Tekrarla] Yukarıdaki adımları her olası değeri için tekrarlayın .n

İki sinyal arasındaki evrişim örneği ve , her bir için sonucunun gösterildiği aşağıdaki Şekil'de gösterilmiştir .s[n]=[211]h[n]=[112]r[n]n

Yukarıdaki sinyal gösteriminde bir değişikliğe dikkat edin. Gerçek sinyaller ve zaman indeksi bir fonksiyonudur, ancak evrişim denklemi bu sinyallerin her ikisini de zaman endeksi . Öte yandan , noktası ile çarpılmadan önceki zaman kaymasını temsil etmek için kullanılır . çıkışı , uygulanan kaydırma olan zaman indeksinin bir fonksiyonudur .s[n]h[n]nmnh[m]s[m]r[n]nh[m]

resim açıklamasını buraya girin

Ardından, bir sinyalin çevrilmesinin gerekli olmadığı daha sezgisel yönteme dönüyoruz.

Sezgisel Yöntem


Evrişimi anlamak için başka bir yöntem daha var. Aslında, bu kıvrım denklemi, yani türetilmesi üzerine kuruludur, çıkış bulmak olarak r[n] s[n]=[2

r[n] = +s[2]h[n+2] +s[1]h[n+1] +s[0]h[n] + s[1]h[n1] + s[2]h[n2] +
Yukarıdaki Şekil ile aynı örneği çözelim, burada ve . Bu, aşağıdaki Tablo'da gösterilmiştir.sa [ n ] = [ - 1s[n]=[2 11]h[n]=[112]

resim açıklamasını buraya girin

Böyle bir yöntem aşağıdaki Şekil'de gösterilmiştir. Uygulama açısından bakıldığında, her iki yöntem arasında bir fark yoktur.

resim açıklamasını buraya girin

Özetle, evrişim bize bir LTI sisteminin belirli bir girdiye tepki olarak nasıl davrandığını anlatır ve yukarıdaki sezgisel yöntem sayesinde evrişimin zaman alanında çarpma olduğunu (ve sinyalin çevrilmesi gerekli olmadığını) söyleyebiliriz. bu kez etki alanı çarpımı bellek içerir. Saygısız gelir ve hangi frekans etki olur, sen benim kitaptan bir örneği bölümüne indirebilirsiniz ayrıca daha derin düzeyde anlamak için burada .

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.