日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

MIT_18.03_微分方程_Convolution_卷积_Notes

發布時間:2025/4/5 编程问答 29 豆豆
生活随笔 收集整理的這篇文章主要介紹了 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=0n?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
稱為fffggg卷積,寫作f?gf*gf?g


文章目錄

  • Convolution
    • @[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?st0t?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時刻的放射性物質為fffggg的卷積 即
∫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的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。