2025年3月12日 星期三

解一維熱傳導方程

眾所周知,Black-Sholes 方程可以化為一維的熱傳導方程,而這個方程可以用傅立葉變換來解。但是一般財經書籍不會去解 Black-Sholes 方程,而是用了 risk neutral 方法直接解出。不過這有種天外飛來一筆的感覺,因為當初是先解了 Black-Sholes 方程,解出來以後才知道可以用 risk neutral 方法來解。所以導致一個奇怪的現象,明明列了 Black-Sholes 方程,卻不去解它。傅立葉變換的代數操作複雜些,以下就來介紹。

首先我們需要一個傅立葉積分的結果,這裡就不推導怎麼算出來的:$$\int_0^\infty e^{-ax^2} \cos bx\ dx=\frac{1}{2} \sqrt{\frac{\pi}{a}}e^{\frac{-b^2}{4a}}\tag{1}$$另外,傅立葉變換不同的作者會採用不同的定義,通常是只是差一個係數,只要變換和反變換係數的乘積為$1/2\pi$就好。這裡採用的定義是:$$\mathcal{F}\{f\}=\hat f(\omega)=\int_{-\infty}^{\infty}f(x)e^{-i\omega x}\ dx\\ \mathcal{F}^{-1}\{\hat f\}=f(x)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\hat f(\omega)e^{i\omega x}\ d\omega$$這裡 $\mathcal{F}\{f\}$ 就是對 $f(x)$ 做傅立葉變換。
我們還需要一些傅立葉變換的性質。首先是 $f(x)$ 的微分性質,在傅立葉變換下就變成簡單的代數操作,也就是乘以 $i\omega$。若 $f$ 是多元函數,而傅立葉變換只對 $x$ 變換,則 $f$ 關於其他變數的偏微分可以直接提到積分號外。若只關心被提到積分號外的變數,可以寫成常微分。$$ (a)\ \ \ \  \mathcal{F}\left\{\frac{\partial f}{\partial x}\right\}=i\omega\mathcal{F}\{f\},\ \ \ \ (b)\ \ \ \  \mathcal{F}\left\{\frac{\partial^2 f}{\partial x^2}\right\}=-\omega^2\mathcal{F}\{f\},\ \ \ \ (c)\ \ \ \  \mathcal{F}\left\{\frac{\partial f}{\partial t}\right\}=\frac{\partial}{\partial t}\mathcal{F}\{f\}$$其次是卷積性質。首先兩個函數 $f(x)$ 和 $g(x)$ 的卷積定義為$$f\ast g = \int_{-\infty}^{\infty}f(u)g(x-u)\ du$$傅立葉變換的卷積性質為$$\mathcal{F}\{f\ast g\}=\mathcal{F}\{f\}\mathcal{F}\{g\}$$傅立葉變換傅立葉變換所需性質已經條列完畢,正式開始解一維熱傳導方程:$$\frac{\partial f}{\partial t}=\kappa\frac{\partial^2 f}{\partial x^2},\ \ \ \ f(x,0)=g(x),\ \ \ \ |f(x,t)|<M$$求解 $f$ 。對偏微分式兩端做傅立葉變換,並利用其微分性質可得$$\frac{d}{dt}\mathcal{F}\{ f\}=-\kappa\omega^2 \mathcal{F}\{ f\}$$這是一基本的對自變數 $t$ 的微分方程,其解為$$\mathcal{F}\{ f\}=Ce^{-\kappa\omega^2 t}$$根據邊界條件可得$$\mathcal{F}\{f(x,0)\}=C=\mathcal{F}\{g(x)\}$$因此$$\mathcal{F}\{ f\}=\mathcal{F}\{g(x)\}e^{-\kappa\omega^2 t}\tag{2}$$接下來我們利用 $(1)$ 找出 $e^{-\kappa\omega^2 t}$是哪一個函數的傅立葉變換。利用歐拉公式 $e^{i\theta}=\cos\theta+i\sin\theta$ 和奇偶函數的積分性質可得 $$\int_{-\infty}^{\infty}e^{-ax^2}e^{-ibx}dx=2\int_0^{\infty}e^{-ax^2}\cos bx\ dx=\sqrt{\frac{\pi}{a}}e^{\frac{-b^2}{4a}}$$令 $b=\omega$得 $$\int_{-\infty}^{\infty}e^{-ax^2}e^{-i\omega x}dx=\mathcal{F}\{e^{-ax^2}\}=\sqrt{\frac{\pi}{a}}e^{\frac{-\omega^2}{4a}}$$令 $\frac{-\omega^2}{4a}=-\kappa\omega^2 t$ 得 $a=\frac{1}{4\kappa t}$ 代入上式得$$\mathcal{F}\{e^{-x^2/4\kappa t}\}=\sqrt{4\kappa\pi t}e^{-\kappa\omega^2t}$$利用傅立葉變換的線性性質調整係數得$$\mathcal{F}\{\sqrt{\frac{1}{4\kappa\pi t}}e^{-x^2/4\kappa t}\}=e^{-\kappa\omega^2t}$$代入 $(2)$ 式得$$\mathcal{F}\{ f\}=\mathcal{F}\{g\}\mathcal{F}\{\sqrt{\frac{1}{4\kappa\pi t}}e^{-x^2/4\kappa t}\}$$因此$$f(x,t)=g\ast\sqrt{\frac{1}{4\kappa\pi t}}e^{-x^2/4\kappa t}=\int_{-\infty}^{\infty}g(u)\sqrt{\frac{1}{4\kappa\pi t}}e^{-(x-u)^2/4\kappa t}du$$

2025年3月11日 星期二

解偏微分方程的特徵曲線法(characteristic curve)

 在大學生數學競賽習題精講(第3版)的 4.2 節有如下題目:$$解偏微分方程\  y\frac{\partial z}{\partial x}-x\frac{\partial z}{\partial y}=0\ \tag{1}$$解法是設 $u=x,\ v=x^2+y^2$。這個解法挺無厘頭,因為等於是事先知道了答案再代回去。實際上解這類微分方程有標準方式,就是特徵曲線法,以下介紹,用的符號也改用一般講偏微分方程的書的符號。

如下的一階線性偏微分方程:$$\frac{\partial u}{\partial x}+p(x,y)\frac{\partial u}{\partial y}=0$$,可以用特徵曲線法來解。首先來研究這個方程的幾何意義。這個方程可以改寫成$$ (1,\ p(x,y)) \cdot \nabla u=0$$也就是 $u$ 的梯度和沿著 $(1,\ p(x,y))$ 這個方向的內積等於零。多元微積分告訴我們這就是
$u$ 在 $(1, \ p(x,y))$ 這個方向的變化為零,或說是 $u$ 在這方向為定值。而 $(1, \ p(x,y))$ 就是 $dy/dx = p(x,y)$ 的解曲線的切線方向。我們把 $dy/dx = p(x,y)$ 的解曲線寫成這個形式:$\varphi(x,y)=C$ 。所以 $u$ 在 $\varphi(x,y)$ 這條曲線 (稱之為特徵曲線) 上為定值,也就是 $u$ 是 $\varphi$ 的函數,可以寫成 $u=f(\varphi(x,y))$ ,$f$ 為任意函數(當然要符合微分條件),而這就是解了。

總結一下解法。首先解出:$$\frac{dy}{dx}=p(x,y)$$這個微分方程,解為 $\varphi(x,y)=C$ 。$u$ 的解就是:$$ u = f(\varphi(x,y))$$

作為範例,我們試著來解 $(1)$ 式。將 $(1)$ 式改寫成$$  \frac{\partial z}{\partial x}-\frac{x}{y}\frac{\partial z}{\partial y}=0$$所以特徵曲線的微分方程就是$$\frac{dy}{dx}=-\frac{x}{y}$$此式可以用變數分離法解出$$\varphi(x,y)=x^2+y^2=C$$,答案就是$$ z=f(x^2+y^2)$$。也因此我們可以發現變數代換 $v=x^2+y^2$ 本身就已經是答案了,所以用這種變數代換去解沒有意義。

2025年1月23日 星期四

量綱函數的形式

$f(\alpha x)=\varphi(\alpha)f(x)$,其中 $f$ 和 $\varphi$ 是可微函數。現在欲證明:$\varphi(\alpha)={\alpha}^{d}$。這裏的 $\varphi$ 就是大名鼎鼎的量綱函數,現在要證明它的形式只能是冪函數。

這裏我們還需要一個 Euler 齊次函數定理:
$$xf_x+yf_y+zf_z=nf(x,y,z)$$滿足此偏微分方程的 $f$ 一定是n次齊次函數:
$$f(tx,ty,tz)=t^nf(x,y,z)$$
現在證明開頭的量綱函數形式。$$\begin{align*} f(\alpha x)=\varphi(\alpha)f(x) \tag{1}            \end{align*}$$令 $\alpha = 1$ 代入 $(1)$ 式可得 $\varphi(1)=1$。
令 $u=\alpha x $,$(1)$式對 $\alpha$ 微分: $$\begin{align*} f_u x &=\varphi '(\alpha)f(x) \\f_u \frac{u}{\alpha}&=\varphi '(\alpha)f(\frac{u}{\alpha}) \\ &=\varphi '(\alpha) \varphi(\frac{1}{\alpha})f(u)  \end{align*}$$令 $\alpha = 1$ 代入上式可得:$$\begin{align*} f_u u &= \varphi '(1)\varphi(1)f(u) \\ &=\varphi '(1)f(u)\end{align*}$$ 由 Euler 齊次函數定理可得$$f(\alpha x)=\alpha^{\varphi '(1)}f(x)$$ 因此$\varphi(\alpha)=\alpha^d$。核查正確性,$\varphi ' (1)$ 確實等於 $d$。

2023年6月13日 星期二

拉格朗日乘數法一個簡易的證明

先簡述一下拉格朗日乘數法:
$f$ 是定義在 $R^n$ 中一開集 $U$ 的一階連續可微的實函數,它被限制在由$$
g_1(x_1,\cdots,x_n)=c_1 \\
\vdots \\
g_m(x_1,\cdots,x_n)=c_m \\
$$所定義的緊致曲面 $S$ 上,並且 $S\subset U$ 。其中各 $g_i$ 都連續可微,並且它們的梯度 $\nabla g_i$ 都線性獨立。如果 $f$ 在 $S$ 中的 $p$ 點達到極值,那麼 $f$ 在 $p$ 點的梯度 $\nabla_p f$ 是 $\nabla_p g_1,\cdots,\nabla _p g_m$ 的線性組合。

這個定理在Charles Chapman Pugh所著的Real Mathematical Analysis一書中有很簡短的證明,以下採用該書的證明。

2023年6月6日 星期二

用微分形式為什麼可以導出散度,梯度,旋度等等的物理量

在這兩篇:
用微分形式導出正交曲線座標系的散度、梯度、旋度和拉普拉斯算符
續上篇,用微分形式導出拉普拉斯算符
當中,對微分形式作幾次d運算,並且把例如$\sigma_1$和$dx$對位起來,就能得到所需要的各種物理量,這中間完全沒提及物理量的定義,用機械式的方式就算出來了,這是為什麼?

2023年6月3日 星期六

續上篇,例題

做一些例題,看上兩篇(12)的公式怎麼使用。

2023年6月2日 星期五

續上篇,用微分形式導出拉普拉斯算符

 符號都沿用上一篇的設定。

首先來看歐氏座標的拉普拉斯算符如何跟微分形式對應。這裡我們需要另一個微分形式的運算,其實它就是說$\ *\ $運算具有分配律而已。$$\ast(\omega_1+\omega_2)={\omega_1}^\ast+{\omega_2}^\ast$$