MIT_18.03_微分方程_Convolution_卷积_Notes
Convolution
引
已知兩個拉普拉斯變換結果F(s),G(s)F(s),G(s)F(s),G(s),
F(s)=∫0∞e?stf(t)dtG(s)=∫0∞e?stg(t)dtF(s) = \int_{0}^{\infty}e^{-st}f(t)dt\\ G(s) = \int_{0}^{\infty}e^{-st}g(t)dt F(s)=∫0∞?e?stf(t)dtG(s)=∫0∞?e?stg(t)dt
它們的乘積F(s)?G(s)F(s)\cdot G(s)F(s)?G(s)為何?
我們知道拉普拉斯變換就是連續版的冪級數
于是我們可以用冪級數的乘積類比出卷積的公式
F(x)=∑0∞anxnG(x)=∑0∞bnxnF(x)?G(x)=∑n=0∞(∑m=0nambn?m)xnF(x) = \sum_0^\infty a_nx^n\\ G(x) = \sum_0^\infty b_nx^n\\ F(x)\cdot G(x) = \sum_{n=0}^\infty(\sum_{m=0}^na_mb_{n-m})x^n F(x)=0∑∞?an?xnG(x)=0∑∞?bn?xnF(x)?G(x)=n=0∑∞?(m=0∑n?am?bn?m?)xn
冪級數的乘積的系數乘為柯西乘積,思考:
xnx^nxn項可以由1,xn1,x^n1,xnx相乘得到 可以由x,xn?1x,x^{n-1}x,xn?1相乘,x2,xn?2x^2,x^{n-2}x2,xn?2相乘,…,xm,xn?mx^m,x^{n-m}xm,xn?m相乘,…,xn,1x^n,1xn,1相乘得到,所以就得出xnx^nxn項的系數為柯西乘積
于是我們把這個式子變連續,即可得:
F(s)?G(s)=∫0∞(∫0tf(u)g(t?u)du)e?stdt=∫0∞(f?g)e?stdtF(s)\cdot G(s) = \int_{0}^{\infty}(\int_0^tf(u)g(t-u)du)e^{-st}dt = \int_{0}^{\infty}(f*g)e^{-st}dt F(s)?G(s)=∫0∞?(∫0t?f(u)g(t?u)du)e?stdt=∫0∞?(f?g)e?stdt
其中
∫0tf(u)g(t?u)du\int_0^tf(u)g(t-u)du ∫0t?f(u)g(t?u)du
稱為fff和ggg的卷積,寫作f?gf*gf?g
文章目錄
- Convolution
- 引
- @[toc]
- 證明
- 性質
- 理解卷積
- 引
- @[toc]
- 證明
- 性質
- 理解卷積
證明
用二重積分換元證明
F(s)?G(s)=∫0∞e?suf(u)du?∫0∞e?svg(v)dvF(s)\cdot G(s) = \int_{0}^{\infty}e^{-su}f(u)du\cdot \int_{0}^{\infty}e^{-sv}g(v)dv F(s)?G(s)=∫0∞?e?suf(u)du?∫0∞?e?svg(v)dv
根據小富比尼定理(Little?Fubini’s?theorem\text{Little Fubini's theorem}Little?Fubini’s?theorem)
=∫0∞∫0∞e?s(u+v)f(u)g(v)dudv=\int_{0}^{\infty}\int_{0}^{\infty}e^{-s(u+v)}f(u)g(v)dudv =∫0∞?∫0∞?e?s(u+v)f(u)g(v)dudv
令t=u+vt = u+vt=u+v, 換元
{u=uv=t?u\left\{\begin{matrix} u = u\\ v = t-u \end{matrix}\right. {u=uv=t?u?
雅可比行列式 JacobianJacobianJacobian
JT=∣10?11∣=1J_T = \begin{vmatrix} 1 & 0\\ -1 & 1 \end{vmatrix} = 1 JT?=∣∣∣∣?1?1?01?∣∣∣∣?=1
原積分區域施以逆變換T?1=(1011)T^{-1} = \begin{pmatrix} 1 & 0\\ 1 & 1 \end{pmatrix}T?1=(11?01?), 即基向量[10]\begin{bmatrix} 1\\0 \end{bmatrix}[10?] 變換為 [11]\begin{bmatrix} 1\\1 \end{bmatrix}[11?]
即得到
=∫0∞∫0te?stf(u)g(t?u)JTdudt=∫0∞e?st∫0tf(u)g(t?u)dudt=∫0∞e?st(f?g)dt=\int_{0}^{\infty}\int_{0}^{t}e^{-st}f(u)g(t-u)J_T\ dudt\\ =\int_{0}^{\infty}e^{-st}\int_{0}^{t}f(u)g(t-u)dudt\\ =\int_{0}^{\infty}e^{-st}(f*g)dt =∫0∞?∫0t?e?stf(u)g(t?u)JT??dudt=∫0∞?e?st∫0t?f(u)g(t?u)dudt=∫0∞?e?st(f?g)dt
性質
摘自wiki https://zh.wikipedia.org/wiki/%E5%8D%B7%E7%A7%AF
交換律 f?g=g?ff*g = g*ff?g=g?f
結合律 f?(g?h)=(f?g)?hf*(g*h) = (f*g)*hf?(g?h)=(f?g)?h
分配律 f?(g+h)=(f?g)+(f?h)f*(g+h) = (f*g) + (f*h)f?(g+h)=(f?g)+(f?h)
數乘結合律 a(f?g)=(af)?g=f?(ag)a(f*g) = (af)*g = f*(ag)a(f?g)=(af)?g=f?(ag) a為任意實數(或復數)
微分定理 D(f?g)=Df?g=f?Dg\mathcal{D}(f*g) = \mathcal{D}f*g = f*\mathcal{D}gD(f?g)=Df?g=f?Dg
理解卷積
可以理解為f(t)f(t)f(t)為系統的元素輸入率,g(t)g(t)g(t)理解為系統元素本身的變化比率
e.g. 核廢料填埋
填埋率為f(t)f(t)f(t),衰變率g(t)=e?ktg(t) = e^{-kt}g(t)=e?kt
某一時間uuu填進廢料f(u)duf(u)duf(u)du 這一堆放射性廢料在ttt時刻衰變為原來的e?k(t?u)e^{-k(t-u)}e?k(t?u)倍
所以ttt時刻的放射性物質為fff和ggg的卷積 即
∫0tf(u)e?k(t?u)du\int_0^tf(u)e^{-k(t-u)}du ∫0t?f(u)e?k(t?u)du
總結
以上是生活随笔為你收集整理的MIT_18.03_微分方程_Convolution_卷积_Notes的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: MIT_18.03_微分方程_Lapla
- 下一篇: C语言拾零(to be continue