杜 寧，張靜宇，張軍鋒

(1 海南大學材料與化工學院，海南 ?？?570228；2 海南政法職業(yè)學院，海南 ?？?571100)

?

現(xiàn)代計算機技術在物理化學實驗數(shù)據(jù)處理中的應用

杜 寧1,2，張靜宇1，張軍鋒1

(1 海南大學材料與化工學院，海南 海口 570228；2 海南政法職業(yè)學院，海南 ?？?571100)

隨著現(xiàn)代計算機技術的發(fā)展，將一些大型的軟件運用于物理化學實驗數(shù)據(jù)處理方面，會使得物理化學實驗數(shù)據(jù)處理更加方便、快捷、準確。本文運用Origin、Matlab、SAS、EViews以及Excel軟件對物理化學中幾個典型的實驗進行數(shù)據(jù)處理，并對部分軟件應用的優(yōu)缺點進行了分析總結，這四種數(shù)據(jù)處理法簡便、快捷、效率高，大大減小了計算和手工繪制相圖的誤差。旨在為物理化學實驗數(shù)據(jù)處理提供參考。

Origin；Matlab；SAS；EViews；Excel；物化數(shù)據(jù)處理

隨著現(xiàn)代計算機技術的發(fā)展，越來越多的軟件應用于科研工作等各個領域。在物理化學實驗數(shù)據(jù)處理方面，也有許多大型的軟件可以運用。

物理化學實驗往往獲得的實驗數(shù)據(jù)量大，計算比較繁瑣，有很多結果是要通過作圖得到的。因此，物理化學實驗的數(shù)據(jù)處理是很重要的，處理的準確性對實驗結果起重要作用。目前，學生數(shù)據(jù)處理方法主要是用結合計算器手工計算，然后用坐標紙手工繪圖求解，其弊端是易出現(xiàn)較大誤差，甚至錯誤，費時、費力。采用計算機軟件處理數(shù)據(jù)，可以避免繁瑣的處理過程，從而使得物理化學實驗數(shù)據(jù)處理更加方便、快捷、準確，還可以提高學生的學習效率和興趣。

1 SAS軟件

SAS(Statistical Analysis System)由美國SAS軟件研究所研制，是一套大型集成式應用軟件，具有完備的數(shù)據(jù)存款、管理和分析功能。在數(shù)據(jù)處理和統(tǒng)計分析領域，SAS被譽為國際權威軟件，廣泛應用于科研、教育、行政管理、生產(chǎn)和金融等不同領域[1]。下面以實驗《陽極極化曲線的測定》為例進行簡單運用介紹，原始實驗數(shù)據(jù)記錄見表1。

表1 陽極極化曲線原始數(shù)據(jù)記錄表

續(xù)表1

140 26040 84150 21040 8160 16040 86170 11040 86180 06040 68190 01040 6920-0 01040 3121-0 06040 3622-0 11040 3723-0 16040 424-0 21040 4425-0 26040 4926-0 31040 4827-0 36040 428-0 41040 3829-0 46040 3730-0 51040 3831-0 56040 3932-0 61040 3833-0 66040 434-0 71040 4135-0 76040 5436-0 81042 1637-0 860412 338-0 910429 239-0 960449 140-1 010471 4941-1 060492 742-1 1104114 57

運用SAS軟件以電流密度為縱坐標，電極電位(相對于參比電極)為橫坐標，繪制陽極極化曲線的步驟如下：

①雙擊圖標，啟動SAS。

②在編輯器窗口輸入如下代碼(其中電極電位數(shù)據(jù)u，電流數(shù)據(jù)i是根據(jù)實驗結果確定的)：

datajihuaqux?ian；inputu1i；u=-u1；cards；0 91040 090 86041 920 81044 470 76048 380 710412 650 66046 250 61041 310 56040 310 51040 480 46040 560 41040 880 36040 980 31040 800 26040 840 21040 800 16040 860 11040 860 06040 680 01040 69-0 01040 31-0 06040 36-0 11040 37-0 16040 40-0 21040 44-0 26040 49-0 31040 48-0 36040 40-0 41040 38-0 46040 37-0 51040 38-0 56040 39-0 61040 38-0 66040 40-0 71040 41-0 76040 54-0 81042 16-0 860412 30-0 910429 20-0 960449 10-1 010471 49-1 060492 70-1 1104114 57；procgplotdata=jihuaquxian；ploti?u；symbolc=blackv=stari=spline；title金屬極化曲線；run；

③輸入代碼后按提交快捷鍵，運行結果如圖1所示，得到極化曲線。

圖1 用SAS得到的極化曲線

2 EViews 軟件

EViews(Econometric Views)是基于Windows平臺的統(tǒng)計分析、模型估計、預測、作圖、模擬及數(shù)據(jù)管理工具。EViews能夠提供統(tǒng)計分析(各種假設檢驗、方差分析、主成分分析、因子分析等)、回歸建模分析(包括單方程線性、非線性模型，聯(lián)立方程模型、離散選擇模型、時間序列模型、分布滯后模型、向量自回歸模型、誤差修正模型、ARCH模型、GARCH模型、面板數(shù)據(jù)模型等)、預測、作圖、模型的求解和模擬、數(shù)據(jù)庫管理等強大功能，在系統(tǒng)數(shù)據(jù)分析與評價、金融分析、經(jīng)濟模型的估計和仿真、宏觀經(jīng)濟預測、成本分析和銷售預測等眾多領域有著廣泛的應用。下面以實驗《旋光法測定蔗糖水解反應速率常數(shù)》為例介紹其使用方法。每隔5 min取一次讀數(shù)，將不同時刻讀取的旋光度值記入表2中；取Δ=30 min，整理有關數(shù)據(jù)也記入表2中。

操作步驟：下載軟件并打開，新建一個變量，選擇變量類型為其他，變量1～6，然后點擊OK。建立一個空的變量，再輸入數(shù)據(jù)，最后得到數(shù)據(jù)組。點擊quick菜單下的Equation Estimation，在彈出對話框中輸入y c t, Method選擇LS，然后點擊“確定”，(或者直接在命令窗口中輸入datayt回車，ls y c t)得到回歸直線結果如圖2所示。以group的形式打開數(shù)據(jù)，出現(xiàn)一個對話框，點擊左上方的 View，選擇 graph，出現(xiàn) Graph Options對話框，在Specific欄選擇Scatter，并在Fitlines欄中選擇Regression Line，然后點擊“確定”，得到回歸直線圖(圖3)。得到線性回歸方程為y=2.269957-0.045591x，擬合優(yōu)度R2=0.998192，說明擬合效果較好。從線性回歸方程可以得到蔗糖水解反應的速率常數(shù)值為k=0.0456 min-1。

表2 不同時刻溶液旋光度

圖2 回歸直線結果

圖3 回歸曲線圖

3 MATLAB軟件

Matlab是當今國際上公認的在科技領域方面最為優(yōu)秀的應用軟件和開發(fā)環(huán)境。在歐美及我國部分高等院校，Matlab已經(jīng)成為應用線性代數(shù)、自動控制理論、數(shù)據(jù)統(tǒng)計、數(shù)字信號處理、時間序列分析、動態(tài)系統(tǒng)仿真、圖形處理等高級課程的基本數(shù)學工具，是攻讀學位的大學生、碩士生、博士生必須掌握的基本技能。Matlab有可靠的數(shù)值和符號計算能力、強大的圖形和可視化功能、簡單易學的程序語言、為數(shù)眾多的應用工具包，大大降低了對使用者的數(shù)學基礎和計算機語言知識的要求，所以它的確為一個高效的科學助手。以實驗《二元液系平衡相圖》為例說明其使用方法，實驗數(shù)據(jù)記錄于表3和表4中。

表3 環(huán)己烷-乙醇標準溶液的折射率

表4 不同組成的環(huán)己烷-乙醇溶液的折射率

下面是詳細操作過程：

①雙擊軟件圖標，啟動MATLAB。

②在命令窗口輸入如下代碼 (n為所測對應濃度的折射率數(shù)據(jù)，c為所測溶液濃度數(shù)據(jù)，ng為測定氣相冷凝的折射率所有數(shù)據(jù)，nl為測定液相的折射率所有數(shù)據(jù)，T1為各溶液沸點，dp為當時氣壓與標準大氣壓之差) ：

n=[1.42338 1.4149 1.4032 1.3908 1.3750 1.35935]; c=[100 80 60 40 20 0]; p=polyfit(n,c,2);%擬合，濃度(c)與折射率(n)的關系p； n1=1.365:0.001:1.414%為了畫曲線將自變量定義更精細； c1=polyval(p,n1);%按擬合的多項式p從n1計算c1； plot(n1,c1,n,c,′*′);%畫圖，從圖上比較擬合的直線與測量點，觀察擬合效果； ng=[1.35935 1.3714 1.3931 1.3938 1.3966 1.4000 1.4030 1.4072 1.4181 1.42338];%氣相冷凝的折射率數(shù)據(jù)； nl=[1.35935 1.3615 1.3708 1.3809 1.3937 1.4043 1.4187 1.4227 1.4230 1.42338];%液相的折射率數(shù)據(jù)； cg=polyval(p,ng);%折射率數(shù)據(jù)計算氣相濃度； cl=polyval(p,nl);%折射率數(shù)據(jù)計算液相濃度； T1=[80.7 75.4 67.7 65.5 64.4 64.5 65.4 67.6 73.5 78.5]; dp=0;%P-PO T=T1+T1/10.*dp;%校正壓力對沸點的影響； plot(cg,T,′*′,cg,T,′b′,cl,T,′o′,cl,T,′black′);%用圓圈標出實驗點。 title(′環(huán)己烷-乙醇氣液平衡溫度-組成圖′); xlabel(′環(huán)己烷的摩爾分數(shù)c/%′);ylabel(′溫度T/℃′);

③輸入代碼后按Enter鍵，得到圖像如圖4，在編輯菜單中選擇復制圖形，然后就可粘貼至所寫的實驗報告中了。④利用數(shù)據(jù)游標確定最低沸點的溫度和組成，單擊圖標，將光標移至曲線最低點，單擊可得最低點的溫度約為64.4 ℃，組成為49.28%，如圖5所示。

圖4 氣相與液相的平衡曲線

圖5 最低恒沸點的溫度和組成

4 Origin軟件

Origin是由美國MicroCall Inc.總公司(世界一流的高靈敏熱量計設計公司)于1991年3月首次推出的基于Windows平臺下用于數(shù)據(jù)分析和工程繪圖的軟件。該軟件不僅包括計算、統(tǒng)計、直線和曲線擬合等各種完善的數(shù)據(jù)分析功能，而且提供了幾十種二維和三維繪圖模板，并將高質(zhì)量科技圖形繪制、C語言編程和NAG數(shù)學統(tǒng)計功能庫集成為一體，其功能強大，是當今世界上最著名的繪圖和數(shù)據(jù)處理軟件之一。Origin軟件是一個多文檔界面應用程序，在使用上，采用直觀的、圖形化的、面向對象的窗口菜單和工具欄操作，容易上手，是公認的簡單易學、操作靈活快速的工程制圖軟件，可以滿足一般用戶及高級用戶的制圖、數(shù)據(jù)分析和函數(shù)擬合的需要。因此在世界各國科技工作者中使用較為普遍。下面以《最大泡壓法測定溶液的表面張力》實驗數(shù)據(jù)處理為例介紹其使用方法。表5中記錄了測定溶液表面張力實驗測得的最大壓力差。

表5 表面張力測定數(shù)據(jù)記錄表

圖6 σ-c擬合曲線圖

圖7 Γ-c擬合曲線圖

圖8 c /Γ-c擬合曲線圖

5 Excel軟件

Excel軟件集數(shù)據(jù)的編輯、整理、統(tǒng)計分析、圖表繪制于一身，利用它處理物理實驗數(shù)據(jù)可以減少枯燥的數(shù)據(jù)運算，防止運算中的錯誤；它的圖表功能使作圖非常容易,并且可以對實驗結果進行分析。物理實驗中實驗數(shù)據(jù)的處理、不確定度的計算、繪制表格利實驗數(shù)據(jù)的圖示，這些工作可以利用高版本Excel中的內(nèi)置工作表函數(shù)得到很方便地解決。以實驗《乙酸乙酯皂化反應速率常數(shù)及活化能的測定》為例介紹Excel 2010處理數(shù)據(jù)過程。

表6 不同時間t下的kt值以及值隨時間t變化的記錄表

注：25 ℃時，測得k0=11.22μS·cm-1，40 ℃時，測得k0=13.98μS·cm-1。

圖9 25 ℃時擬合曲線

6 結 語

Origin功能強大，可以擬合曲線，剔除粗差，尋求經(jīng)驗公式；但較少用于圖像的仿真和圖像的再現(xiàn)，且對計算機的要求相對較高。Matlab具有功能全面，數(shù)據(jù)處理精確，繪制的圖形可從不同角度觀看，曲線擬合不受限制等優(yōu)點；缺點是部分功能的使用需要編程，對使用者的要求較高。Excel功能強大，易學易用，無需編程；但它僅限于直線,對曲線的擬合，誤差較大。EViews和SAS是兩款較多在特定領域應用的軟件。

計算機技術的運用使物理化學實驗數(shù)據(jù)的處理變得快速、準確，大大節(jié)省了處理時間。在教學中運用計算機技術來處理實驗數(shù)據(jù)已經(jīng)得到了大多數(shù)教師的肯定和推廣，能夠提高學生的學習效率和興趣。

[1] 張曉冉.統(tǒng)計分析及其SAS實現(xiàn)[M].北京:清華大學出版社,2011:1-85.

[2] 唐典勇.計算機輔助物理化學實驗[M].北京:化學工業(yè)出版社, 2013:19-30,34-36.

[3] 謝祖芳.物理化學實驗及其數(shù)據(jù)處理[M].成都:西南交通大學出版社, 2014:106-107.

[4] 原安娟,王吉有.幾種軟件在物理實驗數(shù)據(jù)處理中的應用比較[J].大學物理實驗,2007,20(3):82-85.

[5] 肖厚貞,庾名槐.SAS在物理化學實驗數(shù)據(jù)處理中的應用[J].實驗室研究與探索,2009,28(5):75-76.

[6] 盧森杰.基于EViews6.0軟件下的普朗克常量測定[J].物理通報,2012(10):76-80.

[7] 攸頻,張曉峒.EViews6實用教程[M].北京:中國財政經(jīng)濟出版社,2008.

Application of Modern Computer Technology in Data Processing of Physical Chemistry Experiment

DUNing1,2,ZHANGJing-yu1,ZHANGJun-feng1

(1 Materials and Chemical Engineering, Hainan University, Hainan Haikou 570228; 2 Hainan Vocational College of Political Science and Law, Hainan Haikou 571100, China)

With the development of modern computer technology, some large software used in data processing of physical chemistry experiment would make the physical chemistry experimental data processing more convenient, fast and accurate. Origin, Matlab, SAS, EViews and Excel software were applied to experimental chemistry of several typical physical data processing and optimization. The disadvantages of some software applications were analyzed and summarized. These methods were simple, rapid and efficient, which greatly reduced the error of calculation and manual drawing of phase diagram. The aim was to provide reference for physical and chemical experiment data processing.

Origin; Matlab; SAS; EViews; Excel; physical and chemical data processing

杜寧(1981-)，女，助教，主要從事基礎教學及研究。

張軍鋒(1976-)，男，副教授，主要從事天然藥物開發(fā)與基礎化學教學及物理化學實驗教學的研究。

O647 2

A

1001-9677(2016)024-0117-05