張青娟,潘東海,張藝萍,楊 凡
(1.中水北方勘測(cè)設(shè)計(jì)研究有限責(zé)任公司,天津 300222;2.天津市新華職工大學(xué),天津 300041;3.天津城市職業(yè)學(xué)院,天津 300250)
張青娟1,潘東海1,張藝萍2,楊 凡3
(1.中水北方勘測(cè)設(shè)計(jì)研究有限責(zé)任公司,天津 300222;2.天津市新華職工大學(xué),天津 300041;3.天津城市職業(yè)學(xué)院,天津 300250)
本文論述了二項(xiàng)式系數(shù)計(jì)算的四種方法和它們?cè)诮Y(jié)構(gòu)力學(xué)中應(yīng)用的一個(gè)例子。這四種方法是楊輝(Pascal)三角法、二項(xiàng)式定理、系數(shù)遞推法和線性方程組法,其中線性方程組法已脫離了楊輝(Pascal)三角系統(tǒng)成為一個(gè)新的方法,簡(jiǎn)單明了。它們?cè)诮Y(jié)構(gòu)力學(xué)中的應(yīng)用也是初步嘗試。
楊輝(Pascal)三角法;二項(xiàng)式定理;系數(shù)遞推法;線性方程組法
圖1
基于楊輝(Pascal)三角法,而建立起來的二項(xiàng)式定理為:
①
其系數(shù):
②
這一定理說明二項(xiàng)式展開式的系數(shù)與a、b的取值無關(guān)。同時(shí)說明①式為齊次多項(xiàng)式,a為降冪,b為升冪。其它特性亦可由②得出。如:
系數(shù)的對(duì)稱性:
i與n-i為對(duì)稱項(xiàng),于是
因此,ci=cn-i(這是組合數(shù)學(xué)公式之一)
這就是二項(xiàng)式系數(shù)的對(duì)稱性。
二項(xiàng)式定理把楊輝(Pascal)三角的特性通過數(shù)學(xué)歸納法,擴(kuò)充到了任意n值,這樣在理論上很有價(jià)值,尤其在組合數(shù)學(xué)中,基于二項(xiàng)式定理導(dǎo)出許多很有用的公式。二項(xiàng)式定理不僅能計(jì)算二項(xiàng)式的系數(shù),它更是計(jì)算許多問題的一種重要方法和工具。
筆者之一在文獻(xiàn)③中給出了一個(gè)計(jì)算二項(xiàng)式系數(shù)的簡(jiǎn)便方法——系數(shù)遞推法,該法可用一句話來說明:“前一項(xiàng)的系數(shù)值,乘以前一項(xiàng)a的指數(shù),再除以前一項(xiàng)b的指數(shù)加1,即為下一項(xiàng)的系數(shù)”。寫成公式如下:
③
這一公式非常容易記憶。因?yàn)橛散偈娇芍猚k的k正好與b的指數(shù)一一對(duì)應(yīng),因此也就知道了a的指數(shù)為n-k,而前一項(xiàng)ck本身也為已知,這樣就可以從k=0,1,2,3,…,n逐項(xiàng)計(jì)算出來。
該式比原二項(xiàng)式定理要便捷、易記。其實(shí)ci+1與ci相比,即為此式。
倘若計(jì)算某n值的指定k項(xiàng)的系數(shù)時(shí),可按下式計(jì)算:
④
k0=1,k1=n
④式容易證明,根據(jù)系數(shù)遞推法,ck由ck-1組成,ck-1又由ck-2組成,等等,即可導(dǎo)出④式。但這算式有不理想的地方,通用性太差,k0和k1都計(jì)算不出來。實(shí)際上,④式分子分母都乘以(n-k)!后:
此即二項(xiàng)式定理。這也就是為什么二項(xiàng)式定理不把分子、分母的共同因子(n-k)!消去的原因。
至此,可以看到系數(shù)遞推法主要優(yōu)點(diǎn)是比原定理——二項(xiàng)式定理簡(jiǎn)單、易記,實(shí)用性更強(qiáng)。
從上面所述,已經(jīng)知道,二項(xiàng)式的系數(shù)與a、b的取值無關(guān)。但利用a、b的不同取值,卻可算出各系數(shù)。
令a=1,2,3,···,n+1;b≡1,則①式可寫成如下形式:
該行列式顯然不為0,故可唯一地解出c0,c1…cn諸系數(shù)。
該法利用了①式,而①式可由(a+b)n的定義直接導(dǎo)出。它就完全脫離了楊輝(Pascal)三角系統(tǒng),一目了然,成為一個(gè)新法。
N=7時(shí),求二項(xiàng)式的各系數(shù)。
當(dāng)n=7時(shí),有4個(gè)邊界條件已知,即1、7和7、1,又有對(duì)稱性,故僅由兩個(gè)方程式即可求得?,F(xiàn)選擇n=7的前兩個(gè)方程:
c0+c1+c2+c3+c4+c5+c6+c7=27
27c0+26c1+25c2+24c3+23c4+22c5+2c6+c7=37
去掉邊值,利用對(duì)稱性有:
c2+c3=56
3c2+2c3=133
∴c2=21
c3=35
故其解為:1,7,21,35,35,21,7,1。
這同楊輝(Pascal)三角的第七級(jí)是一樣的。
二項(xiàng)式定理在數(shù)學(xué)及其他領(lǐng)域中多有應(yīng)用。本文在結(jié)構(gòu)力學(xué)中進(jìn)行了初步探索,其方法如下。
一個(gè)由等腰對(duì)稱三角形組成的剛架,頂點(diǎn)受垂直單位荷載壓力,計(jì)算各層的壓力分布。頂點(diǎn)將單位垂直壓力一分為二的向下傳遞,類似楊輝(Pascal)三角法的傳播規(guī)律,如圖2所示。提出各層的分母后,留下的分子即為楊輝(Pascal)三角,見圖3。單位荷載在各層的分布可視為擴(kuò)散系數(shù),擴(kuò)散系數(shù)乘實(shí)際頂點(diǎn)壓力即為各層所受壓力。如實(shí)際頂點(diǎn)垂直壓力為500噸,在第五層的分布為:
它們之和為500噸。
圖2
圖3(提出分母,即為楊輝(Pascal)三角)
本文提供了多種形式的有關(guān)二項(xiàng)式系數(shù)的計(jì)算方法。從本文看,前三種方法大同小異,系數(shù)遞推比較簡(jiǎn)便、快捷、易記,是它的一大特點(diǎn)。當(dāng)n相當(dāng)大時(shí),可由線性方程第一分支進(jìn)行初步檢查,如不滿足,一定計(jì)算有誤。它只是必要條件。對(duì)于重要的計(jì)算,也可代入線性方程組,看是否都滿足,如都能滿足,則計(jì)算是成功的。
線性方程組法在現(xiàn)代電子計(jì)算機(jī)技術(shù)如此發(fā)達(dá)時(shí)代解大型方程組不成問題,況且邊界點(diǎn)有四個(gè)已知,利用對(duì)稱性再減去一半的工作量。所以也有實(shí)用價(jià)值。該法的意義是脫離了傳統(tǒng)方法,另成一種簡(jiǎn)單明了的方法。上述應(yīng)用的例子,只是結(jié)構(gòu)力學(xué)的一個(gè)單元,如何應(yīng)用于真實(shí)結(jié)構(gòu),將另文探討。
[1]華羅庚.從楊輝三角談起[M].北京:人民教育出版社,1964.
Brief Review on Calculation of Coefficient of Binomial (a+b)n
ZHANG Qing-juan1, PAN Dong-hai1, ZHANG Yi-ping2, YANG Fan3
(1.RecycledWaterNorthSurveyDesignResearchCo.,LTD,Tianjin, 300222;2.TianjinXinhuaStaffUniversity,Tianjin, 300052;3.TianjinCityVocationalCollege,Tianjin, 300250)
the papers discusses four methods to calculate the coefficient of binomial (a+b)n. An example applied in structural mechanics is also cited. The four methods are Pascal triangle method, binomial theorem, coefficient recursion method and system of linear equations. The last method - system of linear equations has broken away from the Pascal triangle method and developed to be a new method, simple and easy to use. It is also the first time to attempt to apply the four methods in structural mechanics.
Pascal triangle method; binomial theorem; coefficient recursion method; system of linear equations
2014-09-12
張青娟(1963-),女,陜西人,中水北方勘測(cè)設(shè)計(jì)研究有限責(zé)任公司高級(jí)工程師,主要從事計(jì)算機(jī)應(yīng)用方面研究。
O13
A
1673-582X(2015)02-0089-04
天津職業(yè)院校聯(lián)合學(xué)報(bào)2015年2期