查看“︁白金漢π定理”︁的源代码

{{NoteTA
|G1 = Math
|1 = zh:量綱; zh-hans:量纲; zh-hant:因次
|2 = zh-hans:无量纲量; zh-hant:無因次量;
}}

'''白金汉π定理'''是[[因次分析]]中的重要[[定理]]，在工程、應用數學及物理中都會用到。白金汉π定理可以視為是形式化的{{le|雷諾因次分析法|Rayleigh's method of dimensional analysis}}。簡單的說，白金汉π定理指出若有一個物理上有意義的方程，其中有''n''個物理量，而這些物理量共有''k''個獨立的因次，則原方程式可以寫成由''p'' = ''n'' &minus; ''k'' 個[[無因次]]的參數{{pi}}<sub>1</sub>, {{pi}}<sub>2</sub>, ..., {{pi}}<sub>p</sub> 組成的方程式（此處的''k''可以用特定[[矩陣]]的[[秩_(线性代数)|秩]]而得），而這些無因次的參數是由原方程式中的物理量所組成。

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

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

: f<sub>1</sub>(u, ρ, ν, A, F<sub>D</sub>)= 0，

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

: f<sub>2</sub>(C<sub>D</sub>, R<sub>e</sub>)= 0，

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

== 歷史 ==
白金汉π定理得名自美國物理學家{{le|埃德加·白金汉|Edgar Buckingham}}，不過最早是由法國數學家[[約瑟·伯特蘭]]證明<ref>{{cite journal|last=Bertrand|first=J.|year=1878|title=Sur l'homogénéité dans les formules de physique|url=http://gidropraktikum.narod.ru/Bertrand-1878.djvu|journal=Comptes rendus|volume=86|number=15|pages=916–920|access-date=2016-08-17|archive-date=2016-04-09|archive-url=https://web.archive.org/web/20160409033832/http://www.gidropraktikum.narod.ru/Bertrand-1878.djvu|dead-url=no}}</ref>，伯特蘭在1878年時只考慮到一些由電動力學及熱傳導而來的一些問題，不過他的文章中包括所有現在證明此定理會用到的基本概念，而且清楚的指出此定理可以用在物理現象的建模。使用此定理的相關技術（因次法）因為[[約翰·斯特拉特，第三代瑞利男爵]]的貢獻而廣為人知（他是第一位將π定理「應用在通用情形上」的人<ref>When in applying the pi–theorem there arises an ''arbitrary function'' of dimensionless numbers.</ref>，他在研究流體通過管子的壓降的統御參數時用到此方法，日期可能可以追溯到1892年<ref>{{cite journal|last=Rayleigh|year=1892|title=On the question of the stability of the flow of liquids|url=http://gidropraktikum.narod.ru/Rayleigh-1892.djvu|journal=Philosophical magazine|volume=34|pages=59–70|doi=10.1080/14786449208620167|access-date=2016-08-17|archive-date=2021-04-13|archive-url=https://web.archive.org/web/20210413193327/http://gidropraktikum.narod.ru/Rayleigh-1892.djvu|dead-url=no}}</ref>，而在1894年有一個啟發式的證明，再配合數列展開<ref>Second edition of ``The Theory of Sound’’({{cite book|last=Strutt|first=John William|year=1896|title=The Theory of Sound|url=http://archive.org/details/theorysound05raylgoog|publisher=Macmillan|volume=2}}).</ref>）。

可以針對任意多參數的π定理是由A. Vaschy在1892年提出<ref>Quotes from Vaschy’s article with his statement of the pi–theorem can be found in: {{cite journal|last=Macagno|first=E. O.|year=1971|title=Historico-critical review of dimensional analysis|url=http://gidropraktikum.narod.ru/Macagno-1971.djvu|journal=Journal of the Franklin Institute|issue=6|volume=292|pages=391–402|doi=10.1016/0016-0032(71)90160-8|access-date=2016-08-17|archive-date=2021-04-13|archive-url=https://web.archive.org/web/20210413193345/http://gidropraktikum.narod.ru/Macagno-1971.djvu|dead-url=no}}</ref>，1911年分別由A. Federman<ref>{{cite journal|last=Федерман|first=А.|year=1911|title=О некоторых общих методах интегрирования уравнений с частными производными первого порядка|url=http://gidropraktikum.narod.ru/Federman.djvu|journal=Известия Санкт-Петербургского политехнического института императора Петра Великого. Отдел техники, естествознания и математики|issue=1|volume=16|pages=97–155|access-date=2016-08-17|archive-date=2021-04-13|archive-url=https://web.archive.org/web/20210413193256/http://gidropraktikum.narod.ru/Federman.djvu|dead-url=no}} (Federman A., On some general methods of integration of first-order partial differential equations, Proceedings of the Saint-Petersburg polytechnic institute. Section of technics, natural science, and mathematics)</ref>及俄國物理學家{{le|<!--迪米特里·-->Dimitri Riabouchinsky|Dimitri Riabouchinsky}}獨立發現<ref>{{cite journal|last=Riabouchinsky|first=D.|year=1911|title=Мéthode des variables de dimension zéro et son application en aérodynamique|url=http://gidropraktikum.narod.ru/Riabouchinsky-Aerophile-1911.djvu|journal=L'Aérophile|pages=407–408|access-date=2016-08-17|archive-date=2021-04-13|archive-url=https://web.archive.org/web/20210413193053/http://gidropraktikum.narod.ru/Riabouchinsky-Aerophile-1911.djvu|dead-url=no}}</ref>，1914年又由埃德加·白金汉發現<ref>{{Cite web |url=http://gidropraktikum.narod.ru/Buckingham.djvu |title=Original text of Buckingham’s article in ''Physical Review'' |access-date=2016-08-17 |archive-date=2014-04-26 |archive-url=https://web.archive.org/web/20140426233716/http://gidropraktikum.narod.ru/Buckingham.djvu |dead-url=no }}</ref>。白金汉在他的文章中用符號{{pi}}<sub>''i''</sub>表示無因次量，這也成為這定理的名稱。

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

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

:<math>f(q_1,q_2,\ldots,q_n)=0,</math>

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

:<math>F(\pi_1,\pi_2,\ldots,\pi_p)=0,</math>

其中{{pi}}<sub>''i''</sub>為無因次的參數，由''q''<sub>''i''</sub>組合而成，而會有 ''p'' = ''n'' − ''k''個無因次的方程式（稱為Pi群）如下

:<math>\pi_i=q_1^{a_1}\,q_2^{a_2}\cdots q_n^{a_n},</math>

其中的指數''a''<sub>''i''</sub>為有理數，多半情形會是整數。

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

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

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

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

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

=== 因次矩陣及無因次量 ===
假設一個系統，其中有''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個無因次量。

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

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

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

:<math>a=\begin{bmatrix} -1\\ 1\\ 1\\ \end{bmatrix}.</math>

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

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

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

:<math>f(\pi)=0,</math>

可以寫成：

:<math>T=\frac{C D}{V},</math>

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

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

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

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

===單擺===
[[File:Pendel PT.svg|thumb]]
第二個問題是要確認單擺在小幅度振動情形下的週期''T''，假設週期是長度''L''、質量''M''以及地表[[重力加速度]]''g''的函數。其模型如下

:<math>f(T,M,L,g) = 0.\,</math>

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

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

:<math>f(\pi) = 0,</math>

其π為

:<math>\pi =T^{a_1}M^{a_2}L^{a_3}g^{a_4}</math>

而''a''<sub>1</sub>, ..., ''a''<sub>4</sub>為待確定的數值

上述物理量的因次為：

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

因次矩陣為：

:<math>M=\begin{bmatrix}
1 & 0 & 0 & -2\\
0 & 1 & 0 &  0\\
0 & 0 & 1 &  1
\end{bmatrix}.</math>

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

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

:<math>a=\begin{bmatrix}2\\ 0 \\ -1 \\ 1\end{bmatrix}.</math>

因此無因次量為：

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

若只考慮基本因次：

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

也是無因次量。

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

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

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

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

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

== 参考文献 ==
{{Reflist|30em}}

== 延伸閱讀 ==
* {{cite web | author=Hanche-Olsen, Harald | year=2004 | url=http://www.math.ntnu.no/~hanche/notes/buckingham/buckingham-a4.pdf | title=Buckingham's pi-theorem | format=PDF | publisher=NTNU | accessdate=April 9, 2007 | archive-date=2008-09-20 | archive-url=https://web.archive.org/web/20080920234609/http://www.math.ntnu.no/~hanche/notes/buckingham/buckingham-a4.pdf | dead-url=no }}
* {{cite book
| last = Hart
| first = George W.
| date = March 1, 1995
| title = Multidimensional Analysis: Algebras and Systems for Science and Engineering
| publisher = Springer-Verlag
| isbn = 0-387-94417-6
| url = http://www.georgehart.com/research/multanal.html
| access-date = 2016-08-17
| archive-date = 2021-01-26
| archive-url = https://web.archive.org/web/20210126041621/http://www.georgehart.com/research/multanal.html
| dead-url = no
}}
*{{cite book
| last =  Kline
| first =  Stephen J.
| year = 1986
| title = Similitude and Approximation Theory
| publisher = Springer-Verlag, New York
| isbn = 0-387-16518-5
}}
*{{cite book
| last =  Wan
| first =  Frederic Y.M.
| year = 1989
| title = Mathematical Models and their Analysis
| url =  https://archive.org/details/mathematicalmode0000wanf
| publisher = Harper & Row Publishers, New York
| isbn = 0-06-046902-1
}}
* {{cite web
 | author=Vignaux, G.A.
 | year=1991
 | url=http://www.mcs.vuw.ac.nz/~vignaux/docs/maxent.pdf
 | title=Dimensional analysis in data modelling
 | format=PDF
 | publisher=Victoria University of Wellington
 | accessdate=December 15, 2005
 | archive-date=2007-07-20
 | archive-url=https://web.archive.org/web/20070720112213/https://www.mcs.vuw.ac.nz/~vignaux/docs/maxent.pdf
 | dead-url=no
 }}
* Mike Sheppard, 2007 [http://www.msu.edu/~sheppa28/constants/constants.html Systematic Search for Expressions of Dimensionless Constants using the NIST database of Physical Constants] {{Wayback|url=http://www.msu.edu/~sheppa28/constants/constants.html |date=20160308094715 }}
* {{cite book
| last =  Gibbings
| first =  J.C.
| year = 2011
| title = Dimensional Analysis
| publisher = Springer
| isbn = 1-84996-316-9
}}
* {{cite journal
| last=Vaschy | first=A.
| year=1892
| title=Sur les lois de similitude en physique
| journal=Annales Télégraphiques
| volume= 19 | pages=25–28
}}
* {{cite journal
| last=Buckingham | first=E.
| year=1914
| title=On physically similar systems; illustrations of the use of dimensional equations
| url=https://archive.org/details/sim_physical-review_1914-10_4_4/page/345 | journal=Physical Review
| volume=4 | issue=4 | pages=345–376
|doi = 10.1103/PhysRev.4.345
|bibcode = 1914PhRv....4..345B
}}
*{{cite journal
|last=Buckingham |first=E.
|year=1915
|title=The principle of similitude
|url=https://archive.org/details/sim_nature-uk_1915-12-09_96_2406/page/n4 |journal=Nature
|volume=96 |pages=396–397
|doi=10.1038/096396d0
|bibcode = 1915Natur..96..396B
 |issue=2406
}}
*{{cite journal
|last=Buckingham |first=E.
|year=1915
|title=Model experiments and the forms of empirical equations
|journal=Transactions of the American Society of Mechanical Engineers
|volume=37 |pages=263–296
}}
*{{cite journal
 | first = Sir G. | last = Taylor
 | year = 1950
 | title = The Formation of a Blast Wave by a Very Intense Explosion. I. Theoretical Discussion
 | journal = Proceedings of the Royal Society A
 | volume =  201 | pages = 159–174
 | doi = 10.1098/rspa.1950.0049
|bibcode = 1950RSPSA.201..159T
 | issue = 1065
}}
*{{cite journal
 |first        = Sir G.
 |last         = Taylor
 |year         = 1950
 |title        = The Formation of a Blast Wave by a Very Intense Explosion. II. The Atomic Explosion of 1945
 |journal      = Proceedings of the Royal Society A
 |volume       = 201
 |pages        = 175–186
 |url          = http://rspa.royalsocietypublishing.org/content/201/1065/175
 |doi          = 10.1098/rspa.1950.0050
 |bibcode      = 1950RSPSA.201..175T
 |issue        = 1065
 |access-date  = 2016-08-17
 |archive-date = 2018-06-02
 |archive-url  = https://web.archive.org/web/20180602163328/http://rspa.royalsocietypublishing.org/content/201/1065/175
 |dead-url     = no
}}

== 外部連結 ==
* [http://www.gidropraktikum.narod.ru/pi-theorem-history.htm Some reviews and original sources on the history of pi theorem and the theory of similarity (in Russian)] {{Wayback|url=http://www.gidropraktikum.narod.ru/pi-theorem-history.htm |date=20210413193235 }}

== 参见 ==
{{Div col}}
* [[衝擊波]]
* [[因次分析]]
* [[無因次]]
* [[自然單位]]
* [[無因次化]]
* {{le|相似 (模型)|Similitude (model)}}
* {{le|雷諾因次分析法|Rayleigh's method of dimensional analysis}}
{{Div col end}}

[[Category:因次分析]]
[[Category:物理定理]]