范亞茹，亓宇欣

(1. 西南民族大學(xué)數(shù)學(xué)學(xué)院，四川成都 610065；2. 首都師范大學(xué)教師教育學(xué)院，北京100048)

數(shù)學(xué)源于生活，又運(yùn)用于生活. 數(shù)學(xué)建模運(yùn)用數(shù)學(xué)思想、理論和方法，以解決和解釋現(xiàn)實(shí)問(wèn)題為目的，被廣泛地應(yīng)用于各個(gè)領(lǐng)域[1-5]. 數(shù)學(xué)建模競(jìng)賽是高等院校培養(yǎng)大學(xué)生創(chuàng)新能力和實(shí)踐能力的一個(gè)重要手段.根據(jù)調(diào)查[6]，參加過(guò)數(shù)學(xué)建模競(jìng)賽的大學(xué)生的創(chuàng)新能力相對(duì)于沒(méi)有參加過(guò)數(shù)學(xué)建模競(jìng)賽的大學(xué)生較強(qiáng)，參加數(shù)學(xué)建模競(jìng)賽獲得獎(jiǎng)項(xiàng)的大學(xué)生的創(chuàng)新能力水平要明顯高于未參加數(shù)學(xué)建模競(jìng)賽大學(xué)生的創(chuàng)新能力水平，特別是獲得國(guó)家級(jí)獎(jiǎng)的大學(xué)生的創(chuàng)新能力.所以，數(shù)學(xué)建模活動(dòng)對(duì)于大學(xué)生的創(chuàng)新能力是有顯著影響的.

目前，數(shù)學(xué)應(yīng)用軟件(如 MATLAB，SPSS，Visual C ＋ ＋，Python，JAVA 等)在大學(xué)課程的教學(xué)中已廣泛開(kāi)展和實(shí)施[7-10].編程軟件具有強(qiáng)大的符號(hào)處理和繪圖功能，能輔助數(shù)學(xué)表達(dá)式的計(jì)算和推導(dǎo)，驗(yàn)證數(shù)學(xué)假設(shè)和猜想，求解數(shù)學(xué)模型的解析解，展示靜態(tài)和動(dòng)態(tài)的圖表分析結(jié)果.將編程軟件與大學(xué)課程的理論知識(shí)相融合，應(yīng)用于數(shù)學(xué)建模，創(chuàng)新創(chuàng)業(yè)等競(jìng)賽活動(dòng)，能讓學(xué)生真正地將書(shū)本知識(shí)用于求解實(shí)際問(wèn)題，切實(shí)感受到理論與實(shí)踐相結(jié)合的高效性.

本文以2021 年全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽中的B題為例，運(yùn)用統(tǒng)計(jì)軟件SPSS(Statistics Package for Social Sciences)和科學(xué)計(jì)算軟件MATLAB(MATrix LABoratory)高效求解部分問(wèn)題并進(jìn)行分析.SPSS 是世界流行的統(tǒng)計(jì)軟件之一，它采用類似Excel 表格的方式輸入與管理數(shù)據(jù)，有完整的編輯、統(tǒng)計(jì)分析、報(bào)表、圖像制作等功能，易學(xué)易用，操作簡(jiǎn)單[11-14].MATLAB 軟件將數(shù)值計(jì)算、符號(hào)計(jì)算、圖像處理、系統(tǒng)仿真和程序流程控制等功能集成在一個(gè)環(huán)境中，并能夠與其他程序設(shè)計(jì)語(yǔ)言通過(guò)應(yīng)用接口連接，具有許多面向特定應(yīng)用領(lǐng)域的工具箱，功能強(qiáng)大，開(kāi)放性強(qiáng)[15-18].如果學(xué)生在數(shù)學(xué)建模競(jìng)賽中，靈活運(yùn)用SPSS 和MATLAB 軟件統(tǒng)計(jì)和分析數(shù)據(jù)，能更高效地解決問(wèn)題，更深入地分析問(wèn)題，建立更優(yōu)的數(shù)學(xué)模型.

文章的結(jié)構(gòu)如下:接下來(lái)的第二部分是2021 年全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽B 題題目；第三部分是采用MATLAB 和SPSS 數(shù)學(xué)軟件對(duì)競(jìng)賽題目的求解與分析；第四部分是全文總結(jié).

1 數(shù)學(xué)建模競(jìng)賽B 題題目

B 題 乙醇偶合制備C4 烯烴

C4 烯烴廣泛應(yīng)用于化工產(chǎn)品及醫(yī)藥的生產(chǎn)，乙醇是生產(chǎn)制備C4 烯烴的原料.在制備過(guò)程中，催化劑組合(即:Co 負(fù)載量、Co/SiO2 和 HAP 裝料比、乙醇濃度的組合)與溫度對(duì)C4 烯烴的選擇性和C4 烯烴收率將產(chǎn)生影響.因此通過(guò)對(duì)催化劑組合設(shè)計(jì)，探索乙醇催化偶合制備C4 烯烴的工藝條件具有非常重要的意義和價(jià)值.

某化工實(shí)驗(yàn)室針對(duì)不同催化劑在不同溫度下做了一系列實(shí)驗(yàn)，結(jié)果如附件1 和附件2 所示. 請(qǐng)通過(guò)數(shù)學(xué)建模完成下列問(wèn)題:

(1)對(duì)附件1 中每種催化劑組合，分別研究乙醇轉(zhuǎn)化率、C4 烯烴的選擇性與溫度的關(guān)系，并對(duì)附件2中350 度時(shí)給定的催化劑組合在一次實(shí)驗(yàn)不同時(shí)間的測(cè)試結(jié)果進(jìn)行分析.

(2)探討不同催化劑組合及溫度對(duì)乙醇轉(zhuǎn)化率以及C4 烯烴選擇性大小的影響.附錄:名詞解釋

選擇性:某一個(gè)產(chǎn)物在所有產(chǎn)物中的占比.

乙醇轉(zhuǎn)化率:單位時(shí)間內(nèi)乙醇的單程轉(zhuǎn)化率，其值為100%*(乙醇進(jìn)氣量-乙醇剩余量)/乙醇進(jìn)氣量.

表1 附件1 部分?jǐn)?shù)據(jù)Table 1 Part of data of attachment 1

表2 附件2 數(shù)據(jù)Table 2 Data of attachment 2

2 求解與分析

我們先用MATLAB 軟件，將附件1 中的溫度、乙醇轉(zhuǎn)化率、C4 烯烴的選擇性數(shù)據(jù)轉(zhuǎn)化成圖形，便于觀測(cè).下圖分別是編號(hào) A1、A2、A5、A11 催化劑組合方式中溫度、乙醇轉(zhuǎn)化率、C4 烯烴的選擇性的關(guān)系圖.

圖2 相關(guān)性分析流程圖Fig.2 The flow-sheet of correlation analysis

圖1 表明隨著溫度的升高，乙醇轉(zhuǎn)化率、C4 烯烴的選擇性都在升高，呈正相關(guān)關(guān)系. 乙醇轉(zhuǎn)化率圖形更陡峭，而C4 烯烴的選擇性圖形更平緩，表明乙醇轉(zhuǎn)化率受溫度影響較大.

圖1 乙醇轉(zhuǎn)化率、C4 烯烴的選擇性與溫度間的關(guān)系圖Fig.1 The correlation between temperature and ethanol conversion，C4 alkene selectivity

其次，我們通過(guò)SPSS 軟件中的相關(guān)性分析進(jìn)一步研究問(wèn)題.由于附件1 中有21 種催化劑組合方式(A1-A14，B1-B7)，為了節(jié)省論文空間，我們以編號(hào)A1 催化劑組合方式為例. SPSS 軟件可以直接讀取Excel 數(shù)據(jù)文件，然后對(duì)三個(gè)變量進(jìn)行定義；在“分析”菜單欄中的“相關(guān)”子菜單下選擇“雙變量”，可進(jìn)行兩兩變量的相關(guān)性分析；具體操作步驟如下圖所示，其余催化劑組合方式的相關(guān)性分析步驟一致.

編號(hào)A1 的催化劑組合方式的乙醇轉(zhuǎn)化率、C4 烯烴的選擇性與溫度間的相關(guān)性結(jié)果分別如表3 和表4所示，其中皮爾遜相關(guān)系數(shù)(Pearson correlation coefficient)反映兩個(gè)隨機(jī)變量之間的線性相關(guān)程度；顯著性(雙尾)值反映兩個(gè)隨機(jī)變量之間的相關(guān)性的顯著性，其值介于0 到1 之間，越小表明相關(guān)性越顯著.

表3 乙醇轉(zhuǎn)化率與溫度間的相關(guān)性(A1)Table 3 The correlation between temperature and ethanol conversion (A1)相關(guān)性

表4 C4 烯烴選擇性與溫度間的相關(guān)性(A1)Table 4 The correlation between temperature and C4 alkeneselectivity(A1)相關(guān)性

皮爾遜相關(guān)系數(shù)r的定義如下:

其中X，Y是兩個(gè)隨機(jī)變量，n是樣本數(shù)量(X，Y的維度)，分別是X，Y的樣本平均數(shù).r的取值介于-1到1 之間，取值為1 時(shí)，表示兩個(gè)隨機(jī)變量之間呈完全正相關(guān)關(guān)系；取值為-1 時(shí)，表示兩個(gè)隨機(jī)變量之間呈完全負(fù)相關(guān)關(guān)系；取值為0 時(shí)，表示兩個(gè)隨機(jī)變量之間線性無(wú)關(guān).

顯著性(雙尾)值p在相關(guān)性分析中，與皮爾遜相關(guān)系數(shù)是相輔相成的. 當(dāng)顯著性水平一定時(shí)(通常為5%)，p≥0.05 表明兩個(gè)隨機(jī)變量的相關(guān)性不顯著，p＜0.05 表明兩個(gè)隨機(jī)變量的相關(guān)性顯著，且p越小表明相關(guān)性越顯著.

從輸出結(jié)果表3 可以看出乙醇轉(zhuǎn)化率與溫度間的皮爾遜相關(guān)系數(shù)為0.965，顯著性(雙尾)值是0.008( ＜0.05).表明當(dāng)溫度升高時(shí)，乙醇轉(zhuǎn)化率在增加，且乙醇轉(zhuǎn)化率與溫度間有較強(qiáng)的線性相關(guān)性. 同理，表4 顯示C4 烯烴選擇性與溫度間的皮爾遜相關(guān)系數(shù)為0.887，顯著性(雙尾)值是0.045( ＜0.05).當(dāng)溫度升高時(shí)，C4 烯烴選擇性也在增加，C4 烯烴選擇性與溫度間的線性相關(guān)性一般.通過(guò)比較表3 和表4中的皮爾遜相關(guān)系數(shù)和顯著性(雙尾)值，發(fā)現(xiàn)溫度與乙醇轉(zhuǎn)化率的線性相關(guān)性要比與C4 烯烴的選擇性的線性相關(guān)性較強(qiáng)，正如圖1(b)所示.

對(duì)附件1 中的每一種催化劑組合方式，我們都進(jìn)行了乙醇轉(zhuǎn)化率、C4 烯烴的選擇性與溫度間的相關(guān)性分析，結(jié)果如表5 所示.對(duì)于編號(hào) A5、A8、A9、A11、A12、A13、B1、B2、B3、B5、B6、B7 的催化劑組合方式，乙醇轉(zhuǎn)化率與溫度間的皮爾遜相關(guān)系數(shù)小于C4 烯烴的選擇性與溫度間的皮爾遜相關(guān)系數(shù)，表明C4 烯烴的選擇性與溫度間的線性相關(guān)性較強(qiáng).而其余的催化劑組合方式 A1、A2、A3、A4、A6、A7、A10、A14、B4，它們的乙醇轉(zhuǎn)化率與溫度間的皮爾遜相關(guān)系數(shù)大于C4烯烴的選擇性與溫度間的皮爾遜相關(guān)系數(shù)，表明乙醇轉(zhuǎn)化率與溫度間的線性相關(guān)性較強(qiáng). 正如圖1(a)和(b)所示，A1 和A2 催化劑組合方式中，乙醇轉(zhuǎn)化率與溫度存在較強(qiáng)的線性相關(guān)性；而A5 和A11 催化劑組合方式中，C4 烯烴的選擇性與溫度存在較強(qiáng)的線性相關(guān)性，如圖1(c)和(d)所示.

表5 乙醇轉(zhuǎn)化率、C4 烯烴的選擇性與溫度的皮爾遜相關(guān)系數(shù)Table 5 Pearson correlation coefficients between temperatureand ethanol conversion，C4 alkene selectivity

我們繼續(xù)用MATLAB 和SPSS 軟件研究附件2 中350 度時(shí)給定的催化劑組合在一次實(shí)驗(yàn)中乙醇轉(zhuǎn)化率、C4 烯烴的選擇性與時(shí)間之間的關(guān)系.圖3 表明隨著時(shí)間的增長(zhǎng)，乙醇轉(zhuǎn)化率圖形迅速下降，而C4 烯烴的選擇性圖形平緩波動(dòng)，表明時(shí)間對(duì)乙醇轉(zhuǎn)化率影響較大.表6 顯示乙醇轉(zhuǎn)化率與時(shí)間的皮爾遜相關(guān)系數(shù)為-0.966，負(fù)相關(guān)，即乙醇轉(zhuǎn)化率與時(shí)間存在較強(qiáng)的線性相關(guān)性，當(dāng)時(shí)間增長(zhǎng)，乙醇轉(zhuǎn)化率減少. C4 烯烴的選擇性與時(shí)間的皮爾遜相關(guān)系數(shù)為0.215，正相關(guān)，即C4 烯烴的選擇性與時(shí)間不存在顯著的線性相關(guān)性.乙醇轉(zhuǎn)化率與C4 烯烴的選擇性的皮爾遜相關(guān)系數(shù)為-0.150，負(fù)相關(guān)，表明乙醇轉(zhuǎn)化率與C4 烯烴的選擇性線性相關(guān)性不明顯.

圖3 乙醇轉(zhuǎn)化率、C4 烯烴的選擇性與時(shí)間的關(guān)系圖Fig.3 The correlation between time and ethanol conversion，C4 alkene selectivity

表6 乙醇轉(zhuǎn)化率、C4 烯烴的選擇性與時(shí)間的相關(guān)性Table 6 The correlation between time and ethanol conversion，C4 alkene selectivity相關(guān)性

綜上所述，乙醇轉(zhuǎn)化率與催化劑組合方式、溫度和時(shí)間都有較強(qiáng)的相關(guān)性，而C4 烯烴的選擇性與催化劑組合方式和溫度有較強(qiáng)的相關(guān)性，時(shí)間和乙醇轉(zhuǎn)化率對(duì)C4 烯烴的選擇性影響較小.

3 結(jié)論

上述例子表明大學(xué)生數(shù)學(xué)建模競(jìng)賽題目通常求解難度大，求解方法不唯一，且有復(fù)雜的應(yīng)用背景.解決這樣的難題不僅需要實(shí)際問(wèn)題的背景知識(shí)和足夠的數(shù)學(xué)理論知識(shí)，計(jì)算機(jī)軟件應(yīng)用能力及編程能力也非常重要.在數(shù)學(xué)建模過(guò)程中，在思路和方法確定的情況下，復(fù)雜的運(yùn)算可由數(shù)學(xué)軟件完成(如MATLAB，C ＋＋和Python)；在方法思路不明的情況下，數(shù)學(xué)軟件可以測(cè)試構(gòu)思、指引方向并驗(yàn)證猜想.總之，數(shù)學(xué)軟件是數(shù)學(xué)建模必不可少的重要工具，運(yùn)用數(shù)學(xué)軟件高效地解決社會(huì)問(wèn)題應(yīng)是當(dāng)代大學(xué)生必備的技能.