白金漢π定理

**白金汉π定理**是因次分析中的重要定理，在工程、應用數學及物理中都會用到。白金汉π定理可以視為是形式化的。簡單的說，白金汉π定理指出若有一個物理上有意義的方程，其中有*n*個物理量，而這些物理量共有*k*個獨立的因次，則原方程式可以寫成由*p* = *n* &minus; *k* 個無因次的參數1, 2, ..., p 組成的方程式（此處的*k*可以用特定矩陣的秩而得），而這些無因次的參數是由原方程式中的物理量所組成。

白金汉π定理可以視為是一種無因次化的框架，其中提供方法，從已知的物理量中找到一組無因次的參數，甚至此時方程式的具體形式還不清楚也沒有關係。

例如在流體中運動的物體，其阻力方程中包括以下五個物理量：速度 *u*、流體密度 *ρ*、
動黏滯係數 *ν*、 物體截面大小*A*以及阻力 *F<sub>D</sub>*

: f1(u, ρ, ν, A, FD)= 0，

利用白金漢π定理，可以將阻力方程簡化為由阻力係數 *C<sub>D</sub>*及雷諾數 *R*e組成的方程

: f2(CD, Re)= 0，

而這二個物理量是由上述物理量組合而成。

## 歷史
白金汉π定理得名自美國物理學家，不過最早是由法國數學家約瑟·伯特蘭證明，伯特蘭在1878年時只考慮到一些由電動力學及熱傳導而來的一些問題，不過他的文章中包括所有現在證明此定理會用到的基本概念，而且清楚的指出此定理可以用在物理現象的建模。使用此定理的相關技術（因次法）因為約翰·斯特拉特，第三代瑞利男爵的貢獻而廣為人知（他是第一位將π定理「應用在通用情形上」的人，他在研究流體通過管子的壓降的統御參數時用到此方法，日期可能可以追溯到1892年，而在1894年有一個啟發式的證明，再配合數列展開）。

可以針對任意多參數的π定理是由A. Vaschy在1892年提出，1911年分別由A. Federman及俄國物理學家獨立發現，1914年又由埃德加·白金汉發現。白金汉在他的文章中用符號*i*表示無因次量，這也成為這定理的名稱。

## 定理描述
白金漢π定理中，無因次量的個數*p*等於方程式中因次矩陣的零化度（nullity），而*k*是因次矩陣的秩。在實驗上，若不同的系統其無因次量的描述相同，則這二個系統是等效的。

在數學上，若有以下在物理上有意義的方程式

:f(q_1,q_2,\ldots,q_n)=0,

其中*q**i*為*n*個物理量，而這些物理量以*k*個獨立變數表示，則上式可以重組為以下的方程式

:F(\pi_1,\pi_2,\ldots,\pi_p)=0,

其中*i*為無因次的參數，由*q**i*組合而成，而會有 *p* = *n* − *k*個無因次的方程式（稱為Pi群）如下

:\pi_i=q_1^{a_1}\,q_2^{a_2}\cdots q_n^{a_n},

其中的指數*a**i*為有理數，多半情形會是整數。

## 重要性
白金漢π定理提供方式可以針對一些物理現象的物理量計算無因次量，就算還不知方程式的具體形式也沒關係。不過無因次量的選擇不是唯一的，白金漢π定理只提供一個方式找出一組無因次量，不表示這些無因次量是「物理上最有意義」的。

二個系統若有相同的參數，則稱為「相似」（就算相似三角形一樣，相似三角形之間只是大小不同而已）。在方程式上是等效的，因此想要確認方程式形式的實驗者可以選擇最合適，最方便的系統進行。而且白金漢π定理提供了變數及基礎因次個數之間的關係。

## 定理說明
### 簡述
會假設基本物理量及衍生物理量的空間會形成一個有理數的向量空間，基本物理量為單位向量，而物理量之間的相乘變成向量的加法，而物理量的次幂變成向量的純量乘法：

用基本物理量及其次幂表示現象中有關的物理量（若某一基本物理量沒有出現在物理量的表示式中，基本物理量的次幂寫0即可）。例如標準重力 *g* 的單位為D/T^2=D^1T^{-2}（距離除以時間平方），因此若基本物理量為（距離, 時間），標準重力會用向量(1,-2)表示。

因此找到物理量可組成的無因次量，就變成在物理單位的向量空間中處理線性限制的問題。

### 因次矩陣及無因次量
假設一個系統，其中有*n*個有因次的參數，其基本因次有*k*個。即可定義其因次矩陣，其各-{zh-hans:行; zh-hant:列;}-為基本因次，而各-{zh-hans:列; zh-hant:行;}-是有因次的參數：(*i*,&nbsp;*j*)個元素是第*j*個參數在第*i*個基本因次中的幂次項。矩陣可以視為是將所有變數在各基本因次上的幂次項整合到一個矩陣中。

無因次變數是指在各基本因次上的幂次均為0（向量空間中的零向量），也等於矩陣的核。

根據秩－零化度定理，有*n*向量（矩陣的-{zh-hans:列; zh-hant:行;}-）且有*k*個線性獨立向量的系統，其零化度p會滿足(*p*&nbsp;=&nbsp;*n*&nbsp;&minus;&nbsp;*k*)，其中零化度是額外因次的個數，可以選為無因次量。

無因次量可以用因次量的整數次乘積及除法表示。在數學上沒有哪一種無因次量的選擇是比較自然的，不過有些無因次量會比較有物理意義，理想上也會以這類選擇為主。

## 例子

### 速度
以下的例子很簡單，但已足以說明白金漢π定理。

假設一部車以時速100&nbsp;km/h的速度行駛，行駛200&nbsp;km需要多少時間？

問題中有三個因次：距離*D*、時間*T*及速度*V*。這些因次可以用二個基本因次組成：時間*T*及距離*D*，因此會有3&nbsp;&minus;&nbsp;2&nbsp;=&nbsp;1個無因次量。

其因次矩陣為
:M=\begin{bmatrix}
1 & 0 &  1\\
0 & 1 & -1
\end{bmatrix}.

其各-{zh-hans:行; zh-hant:列;}-對應基本因次*D*及*T*，而-{zh-hans:列; zh-hant:行;}-對應三個因次*D*、*T*及*V*。例如第三-{zh-hans:列; zh-hant:行;}-(1,&nbsp;&minus;1)，所以因次*V*（速度）用基本因次表示時，為 D^1 T^{-1} = D/T 。

為了要求得無因次的常數\pi=D^{a_1}T^{a_2}V^{a_3}，需找到一個向量a=[a_1,a_2,a_3]使得矩陣*M*和向量*a*的乘積是零向量[0,0]。在線性代數中，向量*a*稱為矩陣*M*的核，可以生成矩陣的零空間，在此例中，零空間是一維的。以上寫的矩陣是為-{zh-hans:行; zh-hant:列;}-規範形矩陣，因此可以計算得到向量*a*如下，頂多再加一個和向量相乘的係數：

:a=\begin{bmatrix} -1\\ 1\\ 1\\ \end{bmatrix}.

若因次矩陣還不是-{zh-hans:行; zh-hant:列;}-規範形矩陣，也可以用高斯消去法得到因次矩陣的核，因此要求的無因次數可以寫成下式：

:\pi = D^{0}T^0V^0 = D^{-1}T^1(D^1T^{-1}) = D^{-1}T^1V^1 = TV/D .

根據因次分析，可以得到有關上述三個物理量之間的關係式：

:f(\pi)=0,

可以寫成：

:T=\frac{C D}{V},

其中*C*是使C=f^{-1}(0)成立的常數。若剛剛的向量*a*有再加上係數，所求得的*C*也會隨之不同，是原來*C*值的乘幂。

這三個變數之間的關係其實是D=VT，因此實際的無因次量方程（f(\pi)=0）為：

:f(\pi)=\pi-1 = VT/D - 1 = 0.

換句話說*C*只有一個值，就是1，不過這個資訊不是靠因次分析所能得到的。

### 單擺

第二個問題是要確認單擺在小幅度振動情形下的週期*T*，假設週期是長度*L*、質量*M*以及地表重力加速度*g*的函數。其模型如下

:f(T,M,L,g) = 0.\,

（此處寫成四個物理量之間的關係，不是直接將週期*T*寫成另*M*, *L*及*g*的函數）

在上式中有三個基本物理量：時間*t*、質量*m*以及長度*l*，以及四個有因次的物質量。因此只會有4&nbsp;−&nbsp;3&nbsp;=&nbsp;1個無因次量，稱為π，模型可以改寫為

:f(\pi) = 0,

其π為

:\pi =T^{a_1}M^{a_2}L^{a_3}g^{a_4}

而*a*1, ..., *a*4為待確定的數值

上述物理量的因次為：

:T = t, M = m, L = \ell, g = \ell/t^2.

因次矩陣為：

:M=\begin{bmatrix}
1 & 0 & 0 & -2\\
0 & 1 & 0 &  0\\
0 & 0 & 1 &  1
\end{bmatrix}.

（-{zh-hans:行; zh-hant:列;}-表示基本因次*t*、*m*及*l*，而-{zh-hans:列; zh-hant:行;}-表示物理量*T*, *M*, *L*和*g*。例如第四行為(&minus;2,&nbsp;0,&nbsp;1)表示變數*g*的因次為t^{-2}m^0 \ell^1）

計算的目的是要找核向量*a* = [*a*1,&nbsp;*a*2,&nbsp;*a*3,&nbsp;*a*4] 使得矩陣*M*和向量*a*的乘積為零向量[0,0,0]。上述的因次矩陣為列規範形矩陣，因此可以得到核向量*a*如下：

:a=\begin{bmatrix}2\\ 0 \\ -1 \\ 1\end{bmatrix}.

因此無因次量為：

:\begin{align}\pi &= T^2M^0L^{-1}g^1\\
                        &= gT^2/L\end{align}.

若只考慮基本因次：

:\pi=(t)^2(m)^0(\ell)^{-1}(\ell/t^2)^1 = 1,

也是無因次量。

這個例子比較簡單，因為四個物理量中有三個是基本單位，因此最後一個（*g*）的單位可以從前三個單位轉換而來。注意若*a*2不為零，沒有其他的方式可以將*M*的單位消去，因此*a*2需為零。因次分析允許單擺的週期不是質量的函數。（在質量、時間及距離三維組成的三維空間中，質量的向量和其他三個單位的因次向量是線性獨立的，因此若不考慮比例的係數，只有\vec g + 2 \vec T - \vec L是唯一可組成無因次參數的非顯然解。）

因此模型可以用下式表示：

:f(gT^2/L) = 0.\ 

假設*f*只有有限個零點，即可令*gT*2/*L*&nbsp;=&nbsp;*C**n*，其中*C<sub>n</sub>*為第*n*個零點。若只有一個零點，則*gT*2/*L*&nbsp;=&nbsp;*C*。需要更多物理的概念或是實驗才能確認此函數只有一個零點，因此*C*&nbsp;=&nbsp;4π2。

若單擺的振幅較大，分析會比較複雜，需要導入額外的因次變數（最大擺動角度）。上述的分析在擺動角度接近零時有很好的近似效果。

## 参考文献

## 延伸閱讀
* 
* 
*
*
* 
* Mike Sheppard, 2007 Systematic Search for Expressions of Dimensionless Constants using the NIST database of Physical Constants 
* 
* 
* 
*
*
*
*

## 外部連結
* Some reviews and original sources on the history of pi theorem and the theory of similarity (in Russian) 

## 参见

* 衝擊波
* 因次分析
* 無因次
* 自然單位
* 無因次化
* 
*