[傅里葉變換及其應用學習筆記] 五. 傅里葉級數連續性討論,熱方程


這份是本人的學習筆記,課程為網易公開課上的斯坦福大學公開課:傅里葉變換及其應用。

 

熱方程后續

上節課推導出熱方程的傅里葉系數:

$C_k(t) = C_k(0)e^{-2\pi ^2 k^2t}$

那么$C_k(0)$是什么?

 

上節課有提到溫度有如下關系式:

$U(x,t) = \displaystyle{\sum_{k=-\infty}^{\infty}C_k(t)e^{2\pi ikx} }$

當$t=0$,代表初始時刻圓環上的溫度分布

$f(x) = U(x,0) = \displaystyle{\sum_{k=-\infty}^{\infty}C_k(0)e^{2\pi ikx} }$

則,$C_k(0)$為$f(x)$的傅里葉系數

$C_k(0) = \hat{f}(k)$

因此,溫度分布公式(熱方程)如下:

$U(x,t) = \displaystyle{\sum_{k=-\infty}^{\infty}\hat{f}(k)e^{-2\pi^2k^2t}e^{2\pi ikx} }$

溫度$U$與時間$t$的關系為:當$t \to \infty$,$-2\pi^2k^2t \to –\infty$,$e^{-2\pi^2k^2t} \to 0$,$U \to 0$。因此,圓環的溫度最終會變為0。

 

熱方程進一步推導,引入卷積

我們可以對熱方程中的$\hat{f}(k)$進行進一步分解

$\hat{f}(k) = \displaystyle{\int_0^1 e^{-2\pi iky}f(y)dy}$

考慮到初始時刻的溫度分布$f(x)$與熱方程$U(x,t)$中的位置變量$x$可能會取不同的值,我們在此把$f(x)$寫成$f(y)$。

把$\hat{f}(k)$代入熱方程后,得

$\begin{align*}
U(x,t) &=\displaystyle{\sum_{k=-\infty}^{\infty}(\int_0^1 e^{-2\pi iky}f(y)dy) e^{-2\pi^2k^2t}e^{2\pi ikx}} \\
&=\displaystyle{\int_0^1(\sum_{k=-\infty}^{\infty}e^{-2\pi iky}e^{2\pi ikx}e^{-2\pi^2k^2t})f(y)dy } \\
&=\displaystyle{\int_0^1(\sum_{k=-\infty}^{\infty}e^{2\pi ik(x-y)}e^{-2\pi^2k^2t})f(y)dy }
\end{align*}$

$g(x,t) = \displaystyle{\sum_{k=-\infty}^{\infty}e^{2\pi ikx}e^{-2\pi^2k^2t} }$

上面的等式被稱為熱核方程(heat kernel),則

$U(x,t) = \displaystyle{\int_0^1g(x-y,t)f(y)dy }$

如上面的等式,熱方程被轉換成了卷積的表現形式

 

 

從傅里葉級數到傅里葉變換

傅里葉級數到傅里葉變換是從周期現象到非周期現象的轉變,我們可以將非周期函數看做是周期函數的一種特殊情況:周期趨於無窮。

 

對於周期為1的函數

$C_k = \displaystyle{\hat{f}(k) = \int_0^1e^{-2\pi ikt}f(t)dt }$

$f(t) = \displaystyle{\sum_{k=-\infty}^{\infty}\hat{f}(k)e^{2\pi ikt} }$

頻譜圖如下

image

由於$C_k$為復數形式,因此我們無法在圖上畫出,因此只能畫出$\left | C_k \right | = \left | a + bi \right | = \sqrt{a^2 + b^2}$。另外我們在第二節課的時候也學過,$C_k$是y軸對稱的。

 

對於周期為T的函數

$\begin{align*}
C_k = \displaystyle{\hat{f}(k) } &= \displaystyle{\frac{1}{T}\int_0^1e^{-\frac{2\pi }{T}ikt}f(t)dt } \\
    &= \displaystyle{\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}e^{-2\pi i\frac{k}{T}t}f(t)dt }
\end{align*}.$

$f(t) = \displaystyle{\sum_{k=-\infty}^{\infty}\hat{f}(k)e^{2\pi i\frac{k}{T}t} }$

頻譜圖如下

image

由於周期為$T$,因此頻率為$\frac{1}{T}$。當$T \to \infty$,$\frac{1}{T} \to 0$,此時頻譜會變得連續了。

 

 

$T \to \infty$

但是是否僅僅讓$T \to \infty$就能得到傅里葉變換?答案是否定的,下面來看一個例子

有一個函數f(t)如下圖

image

該函數的傅里葉系數求解過程如下

$\begin{align*}
C_k = \displaystyle{\hat{f}(k) }
&= \displaystyle{\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}e^{-2\pi i\frac{k}{T}t}f(t)dt } \\
&= \displaystyle{\frac{1}{T}\int_a^b e^{-2\pi i\frac{k}{T}t}f(t)dt } \\
&\leqslant \displaystyle{\frac{1}{T}\int_a^b \left | e^{-2\pi i\frac{k}{T}t}\right | \left |f(t) \right |dt } \\
&= \displaystyle{\frac{1}{T}\int_a^b Mod(e^{-2\pi i\frac{k}{T}t}) \left |f(t) \right |dt } \\
&= \displaystyle{\frac{1}{T}\int_a^b Mod(cos(-2\pi \frac{k}{T}t) + isin(-2\pi \frac{k}{T}t)) \left |f(t) \right |dt } \quad spread \ with \ Eular \ Formula \\
&= \displaystyle{\frac{1}{T}\int_a^b \sqrt{cos^2(-2\pi \frac{k}{T}t) + sin^2(-2\pi \frac{k}{T}t)} \left |f(t) \right |dt } \\
&= \displaystyle{\frac{1}{T}\int_a^b 1\left |f(t) \right |dt } \\
&= \displaystyle{\frac{1}{T}\int_a^b \left |f(t) \right |dt } \\
&= \frac{M}{T}
\end{align*}.$

即對於所有$C_k$都有$C_k \leqslant \frac{M}{T}$。

$M$是該函數絕對值的積分,是有限值,如果$T \to \infty$,則所有的$C_k \to 0$。所有傅里葉系數為0則該傅里葉變換毫無意義。


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM