MathStudio使用方法入門(mén)7

1、2022年3月6日星期日1/23MathStudioMathStudio for iPadfor iPad 使用方法入門(mén)使用方法入門(mén) （7 7）楊輝三角楊輝三角 2014年4月18日2022年3月6日星期日2/23早年工作之余，喜歡涉獵一些科普讀物，物理、數(shù)學(xué)的都有；這些讀物即開(kāi)闊了視野，增長(zhǎng)了知識(shí)，又不失為排遣休閑的好方法?？吹饺A羅庚的從楊輝三角談起，用現(xiàn)代數(shù)學(xué)方法探究中國(guó)古算的奧秘，深入淺出，饒有興趣。 退休后，自由支配的時(shí)間充裕了，閑暇無(wú)事，找一些數(shù)學(xué)題琢磨，對(duì)于不抽煙、不喝酒、不打牌的“枯燥型”老人，既能延緩大腦退化，又能打發(fā)時(shí)間，經(jīng)濟(jì)實(shí)惠，何樂(lè)不為！ 近日，重溫從楊輝三角談起，又選讀

2、了一些關(guān)于楊輝三角的課件，頗有啟發(fā)，遂萌發(fā)了用MathStudio為工具，闡發(fā)中國(guó)古算題奧秘的念頭。原來(lái)覺(jué)得深?yuàn)W難算的高階等差級(jí)數(shù)題，用MathStudio的 Binomial、list、loopend等函數(shù)，竟而簡(jiǎn)捷明快、輕而易舉地得到答案。高興之余，整理成文，與同好者共享。 本文引用了一些百度文庫(kù)里課件的圖形資料，給我省了很多時(shí)間，未查到作者，在此敬向無(wú)名作者深致謝意。2022年3月6日星期日3/23 左圖 現(xiàn)稱(chēng)楊輝三角，本名是開(kāi)方作法本源或稱(chēng)乘方求廉圖，是古人對(duì)任何數(shù)求任意次方根操作過(guò)程的運(yùn)算工具。 設(shè)想如果在沒(méi)有電腦、計(jì)算器 、計(jì)算尺、對(duì)數(shù)表等現(xiàn)代計(jì)算工具的情況下，只用紙、筆（甚至紙筆

3、也不用）進(jìn)行大數(shù)值的開(kāi)高次方運(yùn)算，我們?cè)撛趺醋?？?huì)嗎？ 在這里，讓我們用MathStudio為工具，探索、演示楊輝三角的一些有趣的特性；這就好比騎著摩托車(chē)在古人開(kāi)辟的古道上徜徉，感受一下前賢先哲們的精微睿智，體驗(yàn)他們的艱辛苦樂(lè)。2022年3月6日星期日4/23上圖是擴(kuò)展的楊輝三角，增加上圖是擴(kuò)展的楊輝三角，增加了通項(xiàng)表達(dá)式。了通項(xiàng)表達(dá)式。每行數(shù)字對(duì)應(yīng)于（每行數(shù)字對(duì)應(yīng)于（a+ba+b）n n 的的 展開(kāi)式的分離系數(shù)。展開(kāi)式的分離系數(shù)。楊輝三角有很多有趣的特性1.三角形的兩側(cè)斜邊上都是數(shù)字1，而其余的數(shù)都等于它肩上的兩個(gè)數(shù)字其余的數(shù)都等于它肩上的兩個(gè)數(shù)字相加相加 (C(n,0)=1,C(n,n)=

4、0;C(n,r)=C(n-1,r-1)+C(n-1,r)2.具有對(duì)稱(chēng)性，每行與首末兩端“等距離 ”的兩個(gè)數(shù)相等(C(n,r)=C(n,n-r) 3.每一行的第二個(gè)數(shù)就是這行的行數(shù)（C(n,1)=1）4.第n行包含n+1個(gè)數(shù)，每行數(shù)字和為2n。5.所有行的第二個(gè)數(shù)構(gòu)成等差數(shù)列； 所有行的第三個(gè)數(shù)構(gòu)成二階等差數(shù)列； 所有行的第r個(gè)數(shù)構(gòu)成r-1階等差數(shù)列； 6.左下向右上的斜線(xiàn)上數(shù)字之和構(gòu)成斐波那契數(shù)列能用能用 MathStudio MathStudio 演繹這些有趣的特演繹這些有趣的特性嗎？請(qǐng)接著看性嗎？請(qǐng)接著看2022年3月6日星期日5/23nCr(n,r) 從n件物品中取出r件物品的組合數(shù)，只

5、作數(shù)值計(jì)算。Binomial(n,r)= n! / r! * (n-r)! , 二項(xiàng)式系數(shù)，可有限度地進(jìn)行符號(hào)運(yùn)算，r可以是非整數(shù)。Binomial 比 nCr 有更大的計(jì)算范圍因?yàn)橐驗(yàn)锽inomial(n,0)=1Binomial(n,0)=1 Binomial(n,n)=1 Binomial(n,n)=1所以楊輝三角每行兩側(cè)斜線(xiàn)數(shù)字都是所以楊輝三角每行兩側(cè)斜線(xiàn)數(shù)字都是1 1 Binomial(n,0),Binomial(n,n+1)以及nCr(n,0), nCr(n,n+1) 本無(wú)意義，為便于表達(dá)，作為記號(hào)分別定義為1，0 2022年3月6日星期日6/23 楊輝三角的基本性質(zhì)楊輝三角形的每

6、行數(shù)列，除兩側(cè)斜邊上的數(shù)字1外，其余任何一個(gè)數(shù)都其余任何一個(gè)數(shù)都等于它肩上的兩個(gè)數(shù)字相加之和等于它肩上的兩個(gè)數(shù)字相加之和 上圖上圖 上式，稍加整理、化簡(jiǎn)，就得下式上式，稍加整理、化簡(jiǎn)，就得下式Binomial(n-1,r-1)+Binomial(n-1,r)=Binomial(n,r)2022年3月6日星期日7/23對(duì) 稱(chēng) 性C(n，r)-C(n，n-r)=0Binomial(n，r)-Binomial(n，n-r)0即即C(n，r)=C(n，n-r)Binomial(n，r)=Binomial(n，n-r)對(duì)稱(chēng)中心（最大值項(xiàng)）對(duì)稱(chēng)中心（最大值項(xiàng)）n=n=偶數(shù)時(shí)偶數(shù)時(shí) r=n/2r=n/2n=

7、n=奇數(shù)時(shí)奇數(shù)時(shí) r1=(n-1)/2r1=(n-1)/2 r2=(n+1)/2 r2=(n+1)/2 注：b=-1.42109E-14 可能是Binomial計(jì)算過(guò)程的截尾誤差，極小，視為02022年3月6日星期日8/23MathStudio可以列出 指定行的整行數(shù)據(jù)或任意項(xiàng)的數(shù)據(jù)根據(jù)楊輝三角的通項(xiàng)表達(dá)式第n+1行C(n,0),C(n,1),C(n,2)C(n,r)C(n,n)C(n,0),C(n,1),C(n,2)C(n,r)C(n,n)左圖 n=16 的整行數(shù)據(jù)，a數(shù)列有17個(gè)數(shù)據(jù)，向左拖曳a數(shù)列表可看到全部數(shù)據(jù)。這行17個(gè)數(shù)據(jù)之和為 216 = 65536 2022年3月6日星期日9/

8、23 對(duì) 稱(chēng) 性N=17 整行數(shù)據(jù)當(dāng)n為奇數(shù)時(shí)，最大值有2個(gè)， r1=(n-1)/2 r2=(n+1)/2 ncr(17,8)=ncr(17,9)=24310 n ncr(n,r)=2n r=02022年3月6日星期日10/23第第n n行，逐項(xiàng)平方和行，逐項(xiàng)平方和 等于等于第第2n2n行的第行的第n n項(xiàng)項(xiàng) n c (n,r)2= c (2n,n) r=02022年3月6日星期日11/23 楊輝三角 與高階等差級(jí)數(shù)C(r,r)+C(r+1,r)+C(r+2,r)+C(r+3,r)+C(r+m-1,r)=C(r+m,r+1) nrMathStudioMathStudio可以列出首項(xiàng)為可以列出首

9、項(xiàng)為1 1的的r r階等差級(jí)數(shù)數(shù)列及其前階等差級(jí)數(shù)數(shù)列及其前m m項(xiàng)之和項(xiàng)之和r r階等差級(jí)數(shù)階等差級(jí)數(shù) 首項(xiàng)首項(xiàng)1 1，第，第m m項(xiàng)是項(xiàng)是 C C(r+m-1,r(r+m-1,r）, ,前前m m項(xiàng)之和是項(xiàng)之和是 C C(r+m,r+1(r+m,r+1）2022年3月6日星期日12/23楊輝三角的斜線(xiàn)部分楊輝三角的斜線(xiàn)部分 構(gòu)成等差級(jí)數(shù)構(gòu)成等差級(jí)數(shù)r=2 1，2，3，4 等差級(jí)數(shù)r=3 1，3，6，10 二階等差級(jí)數(shù)r=4 1，4，10，20 三階等差級(jí)數(shù)r=5 1，5，15，35 四階等差級(jí)數(shù)左上圖 四階等差級(jí)數(shù)前6項(xiàng)，及6項(xiàng)之和左下圖 五階等差級(jí)數(shù)前10項(xiàng)，及10項(xiàng)之和 n=r+m-1

10、2022年3月6日星期日13/23MathStudio怎么計(jì)算任意高階等差級(jí)數(shù)從楊輝三角圖的斜線(xiàn)上觀(guān)察到的高階等差級(jí)數(shù)數(shù)列都是首項(xiàng)a1=1，dr=1的“原始型”高階等差級(jí)數(shù)，如果a11，dr1，怎么辦？下面列出運(yùn)用MathStudio計(jì)算任意高階等差級(jí)數(shù)的通式，式子怎么來(lái)的？不贅述了，有興趣的朋友自己找資料探究吧an=an=a a1+C(n-1,1)1+C(n-1,1)d d1+C(n-1,2)1+C(n-1,2)d d2+C(n-1,r)2+C(n-1,r)d dr rSn=C(n,1)Sn=C(n,1)a a1+C(n,2)1+C(n,2)d d1+C(n,3)1+C(n,3)d d2+C

11、(n,r+1)2+C(n,r+1)d dr rnrnr如果要用上式計(jì)算r階高階等差級(jí)數(shù)的第n項(xiàng)，我們必須知道a1,d1,d2,dr的數(shù)值,共r+1個(gè)數(shù)據(jù)，用到式子的r+1項(xiàng)常見(jiàn)的二階等差級(jí)數(shù)，只需知道3個(gè)數(shù)字，用到式子的3項(xiàng)，很容易。至此，我們可以探討堆垛術(shù)了2022年3月6日星期日14/23 堆垛術(shù) 長(zhǎng)方垛按左上圖題意，酒壇堆垛類(lèi)似左下圖，上層49，下面每層長(zhǎng)寬各遞加1，共8層。各層酒壇數(shù)如以下數(shù)列：3636，5050，6666，8484，104104，126126，150150，176176各層酒壇數(shù)之差為：1414，1616，1818，2020，2222，2424，2626各層酒壇數(shù)之差

12、之差為：2 2，2 2，2 2，2 2，2 2，2 2顯然，這是個(gè)二階等差級(jí)數(shù)a1=36,d1=14,d2=2左圖，算出各層酒壇數(shù)，累加數(shù)為792.按前面的通式計(jì)算：n=8,a1=ab=36 d1=(a+1)(b+1)-a1=50-36=14,d2=(a+2)(b+2)-(a+1)(b+1)-d1=66-50-14=2Sn=C(n,1)a1+C(n,2)d1+C(n,3)d2Sn=C(n,1)a1+C(n,2)d1+C(n,3)d2 =8 =836+2836+2814+5614+562=7922=792無(wú)誤無(wú)誤2022年3月6日星期日15/23長(zhǎng)方垛在華羅庚科普著作從楊輝三角談起里，給出了計(jì)算

13、長(zhǎng)方垛的通式（底層計(jì)算長(zhǎng)方垛的通式（底層14141010，由，由底層向上遞減）底層向上遞減）：左圖 底層a=14 b=10起，每層長(zhǎng)寬各遞減1，10層長(zhǎng)方垛的計(jì)算結(jié)果。a數(shù)列 由頂?shù)降祝鲗觽€(gè)數(shù)c 各層累計(jì)總和m 按長(zhǎng)方垛計(jì)算通式的計(jì)算結(jié)果按前述通式計(jì)算：n=10,a1=5,d1=12-5=7,d2=21-12-7=2Sn=C(n,1)a1+C(n,2)d1+C(n,3)d2Sn=C(n,1)a1+C(n,2)d1+C(n,3)d2 =10 =105+455+457+1207+1202=6052=605不同方法計(jì)算結(jié)果相同不同方法計(jì)算結(jié)果相同 S2022年3月6日星期日16/23長(zhǎng)方垛 Mat

14、hStudioMathStudio用通式的計(jì)算程序用通式的計(jì)算程序an=an=a a1+C(n-1,1)1+C(n-1,1)d d1+C(n-1,2)1+C(n-1,2)d d2 Sn=C(n,1)2 Sn=C(n,1)a a1+C(n,2)1+C(n,2)d d1+C(n,3)1+C(n,3)d d2 22022年3月6日星期日17/23正方垛如果每層長(zhǎng)寬酒壇數(shù)相等，就是正方垛各層酒壇數(shù)為1, 4, 9, 16, 25, 36, 49, 64各層酒壇數(shù)之差為3，5，7，9，11，13，15各層酒壇數(shù)之差之差為2，2，2，2，2，2Sn=C(n,1)a1+C(n,2)d1+C(n,3)d2Sn

15、=C(n,1)a1+C(n,2)d1+C(n,3)d2 =8 =81+281+283+563+562=2042=204在華羅庚科普著作從楊輝三角談起里，給出了計(jì)算正方垛的通式計(jì)算正方垛的通式： n=8 代入上式 S=8917/6 = 204計(jì)算結(jié)果相同S=2022年3月6日星期日 18/23C(0,0)=1, F（1）C(1,0)=1, F（2）C(2,0)+C(1,1)=2 F（3）C(3,0)+C(2,1)=3 F（4）C(4,0)+C(3,1)+C(2,2)=5 F（5）C(5,0)+C(4,1)+C(3,2)=8 F（6） C(n,0)+C(n-1,1)+C(n-2,2)+C(n-ce

16、il(n/2),floor(n/2) = Fibonacci(n+1) n3楊輝三角 與斐波那契數(shù)列從左圖可看出從左圖可看出Fibonacci(n)Fibonacci(n)是是4545斜線(xiàn)上數(shù)字之斜線(xiàn)上數(shù)字之和，數(shù)字的個(gè)數(shù)是和，數(shù)字的個(gè)數(shù)是ceil(n/2)ceil(n/2)因?yàn)镃(5,0)=C(4,0)=C(3.0)=C(2,0)=1 C(4,1)=C(3,0)+C(3,1) C(3,2)=C(2,1)+C(2,2)所以F(6)=1+C(4,1)+C(3,2) =1+C(3,0)+C(3,1)+C(2,1)+C(2,2) =1+C(3,1)+C(2,2)+1+C(2,1) =F(5)+F(4

17、) 符合斐波那契數(shù)列定義，余類(lèi)推2022年3月6日星期日19/23MathStudio用二項(xiàng)式系數(shù)計(jì)算斐波那契數(shù)列左圖 分別是Fibonacci(9)和Fibonacci(20)的計(jì)算結(jié)果a數(shù)列 楊輝三角圖上相應(yīng)斜線(xiàn)上的數(shù)字b a數(shù)列各數(shù)字之和2022年3月6日星期日20/23關(guān)于數(shù)值范圍的限制前面說(shuō)到運(yùn)用MathStudio計(jì)算任意高階等差級(jí)數(shù)的通式，似乎說(shuō)MathStudio法力無(wú)邊、神奇無(wú)比，如果不說(shuō)明它的數(shù)域范圍，就未免有夸大之嫌了。事實(shí)上，我們的計(jì)算能力是受到軟硬件限制的；MathStudio這個(gè)只有1MB多個(gè)頭的精干軟件，它的計(jì)算能力究竟怎樣呢？請(qǐng)看左圖：階乘階乘 171171！溢

18、出！溢出 (7.2(7.2E E306)306)nCr (360,125)nCr (360,125)溢出溢出 (2.1(2.1E E99)99)BinomialBinomial（10301030，500500）溢出）溢出 (7.1 (7.1E E307)307)2n 10242n 1024溢出溢出 (8.9(8.9E E307)307)FibonacciFibonacci（605605）溢出）溢出 (7.5(7.5E E125)125)（猜想：觀(guān)察上述情況，在一般條件下，（猜想：觀(guān)察上述情況，在一般條件下，計(jì)算過(guò)程允許計(jì)算到計(jì)算過(guò)程允許計(jì)算到 E E308,308,在高精度條在高精度條件下，允

19、許計(jì)算到件下，允許計(jì)算到 E E100100）Binomial比nCn有大得多的數(shù)值范圍MathStudioMathStudio內(nèi)置函數(shù)內(nèi)置函數(shù)BinomialBinomial只能計(jì)算到第只能計(jì)算到第604604項(xiàng)，項(xiàng)，能不能突破這個(gè)能不能突破這個(gè)“瓶頸瓶頸”呢？能！呢？能！上面左圖是按斐波那契數(shù)列的定義，用上面左圖是按斐波那契數(shù)列的定義，用loopendloopend逐項(xiàng)計(jì)算到第逐項(xiàng)計(jì)算到第14761476項(xiàng)項(xiàng)上面右圖是用斐波那契數(shù)列的二項(xiàng)式系數(shù)表達(dá)式（楊輝三角上面右圖是用斐波那契數(shù)列的二項(xiàng)式系數(shù)表達(dá)式（楊輝三角4545斜線(xiàn)），逐項(xiàng)計(jì)算到斜線(xiàn)），逐項(xiàng)計(jì)算到 第第14761476項(xiàng)，注意：用項(xiàng)

20、，注意：用BinomialBinomial，不要用，不要用nCrnCr兩種方法都能把斐波那契數(shù)列的可計(jì)算項(xiàng)數(shù)擴(kuò)大兩種方法都能把斐波那契數(shù)列的可計(jì)算項(xiàng)數(shù)擴(kuò)大2.442.44倍，答案數(shù)值擴(kuò)大了倍，答案數(shù)值擴(kuò)大了1.71.7E E182182倍倍2022年3月6日星期日22/231 2 3 4 5 6 7 8 910 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 3637 38 39 40 41 42 43 44 45 46 47 48 49數(shù)字金字塔左圖 網(wǎng)上流傳的數(shù)字金字塔求解 第n行的項(xiàng)數(shù)，及其首項(xiàng)、 中項(xiàng)、末項(xiàng)、每行數(shù)字之和應(yīng)用前面已有知識(shí)很容易得到答案各行的項(xiàng)數(shù)：各行的項(xiàng)數(shù)： 1，3，5，7，9（2n-12n-1）各行的首項(xiàng)： 1，2，5，10，17，26，37各行首項(xiàng)之差： 1，3，5，7， 9， 11各行首項(xiàng)之差之差： 2，2，2，2， 2這是2階等差級(jí)數(shù)，