王 靜， 劉 冰， 李 輝， 王 濤＊

（1.長(zhǎng)春工業(yè)大學(xué)計(jì)算機(jī)科學(xué)與工程學(xué)院，吉林長(zhǎng)春 130012；2.吉林省塑料研究院，吉林長(zhǎng)春 130022）

基于約束滿足的ERP測(cè)試計(jì)劃方法

王 靜1， 劉 冰1， 李 輝2， 王 濤1＊

（1.長(zhǎng)春工業(yè)大學(xué)計(jì)算機(jī)科學(xué)與工程學(xué)院，吉林長(zhǎng)春 130012；2.吉林省塑料研究院，吉林長(zhǎng)春 130022）

將約束編程與ERP產(chǎn)品自身特點(diǎn)相結(jié)合，提出了一種基于約束滿足的ERP測(cè)試計(jì)劃方法。通過(guò)對(duì)某飼料行業(yè)ERP中記賬憑證軟件進(jìn)行約束建模。該測(cè)試計(jì)劃方法已經(jīng)應(yīng)用于飼料行業(yè)ERP平臺(tái)中。

約束滿足；ERP；測(cè)試計(jì)劃；約束建模

0 引 言

企業(yè)資源計(jì)劃（Enterprise Resource Planning，ERP）以其信息化管理方式改進(jìn)企業(yè)的業(yè)務(wù)流程、完善供需鏈、提高企業(yè)經(jīng)濟(jì)效益，雖然目前ERP在我國(guó)的應(yīng)用并不成熟，但是這種信息化的管理方式具有充分的發(fā)展空間。對(duì)于ERP產(chǎn)品，其持續(xù)有效存活取決于產(chǎn)品的質(zhì)量，因此產(chǎn)品的測(cè)試工作尤為重要。據(jù)統(tǒng)計(jì)，在一個(gè)軟件產(chǎn)品開發(fā)過(guò)程中，軟件測(cè)試工作所占比重大致為軟件開發(fā)總工作量的40%［1］。ERP具有開發(fā)周期長(zhǎng)、流程化的特點(diǎn)，目前，ERP測(cè)試計(jì)劃的制定多以經(jīng)驗(yàn)為指導(dǎo)，只是根據(jù)產(chǎn)品模塊以及測(cè)試人員的時(shí)間等因素主觀進(jìn)行測(cè)試資源分配，文獻(xiàn)［1］給出一種軟件測(cè)試計(jì)劃方法，一定程度上提高了測(cè)試的效率，但對(duì)ERP產(chǎn)品而言并不完全適合。ERP具有模塊聯(lián)系緊密的特性，對(duì)業(yè)務(wù)操作流程不熟悉的測(cè)試人員而言，測(cè)試結(jié)果可能并不盡人意，遺留bug過(guò)多必然影響產(chǎn)品的正常實(shí)施［1］。因此對(duì)于ERP的測(cè)試，一般需要具有符合產(chǎn)品特點(diǎn)的實(shí)施測(cè)試的計(jì)劃和方法，對(duì)測(cè)試的模塊、內(nèi)容、分配的測(cè)試人員等方面進(jìn)行詳細(xì)描述，制定合理的測(cè)試資源分配方法是十分必要的。

約束程序（Constraint Programming，CP）是人工智能領(lǐng)域一個(gè)重要的研究分支，其在描述約束滿足問(wèn)題方面取得了很多成熟的研究結(jié)果，軟件測(cè)試的分配方法類似于約束滿足問(wèn)題中的一類調(diào)度問(wèn)題，因此，文中將ERP測(cè)試計(jì)劃問(wèn)題轉(zhuǎn)化為約束滿足問(wèn)題（Constraint Satisfaction Problem，CSP）進(jìn)行求解，具有一定的可行性與合理性?；诖?，文中提出了一種基于約束滿足的ERP測(cè)試計(jì)劃方法，對(duì)某飼料行業(yè)ERP中記賬憑證軟件約束建模，測(cè)試驗(yàn)證了文中方法的有效性，目前，該測(cè)試方法已經(jīng)應(yīng)用于飼料行業(yè)ERP平臺(tái)并取得較好效果。這對(duì)于指導(dǎo)ERP產(chǎn)品的測(cè)試工作具有很重要的現(xiàn)實(shí)意義［2－3］。

1 測(cè)試模塊的選取

軟件測(cè)試是軟件質(zhì)量保證的重要手段，研究表明：軟件產(chǎn)品開發(fā)費(fèi)用和軟件質(zhì)量與軟件產(chǎn)品中存在問(wèn)題發(fā)現(xiàn)時(shí)間早晚有直接關(guān)系。盡早發(fā)現(xiàn)產(chǎn)品中的問(wèn)題，后期維護(hù)費(fèi)用也相應(yīng)會(huì)減少。目前軟件測(cè)試的工作主要放在測(cè)試技術(shù)上，對(duì)軟件測(cè)試過(guò)程關(guān)注很少，軟件測(cè)試過(guò)程對(duì)軟件的有效性及軟件質(zhì)量具有很好的衡量作用，因此，軟件測(cè)試過(guò)程的改進(jìn)具有重要的現(xiàn)實(shí)意義。CMMI是針對(duì)軟件過(guò)程改進(jìn)提出的，軟件產(chǎn)品測(cè)試為軟件過(guò)程的一部分，因此，CMMI模型對(duì)于軟件產(chǎn)品測(cè)試具有一定的指導(dǎo)意義。度量由CMMI提供，主要是為了獲取過(guò)程的表征數(shù)據(jù)［45］。ERP對(duì)企業(yè)的采購(gòu)、生產(chǎn)、銷售進(jìn)行統(tǒng)一整合，每部分功能關(guān)聯(lián)緊密。軟件產(chǎn)品的測(cè)試人員并不一定是參與開發(fā)的人員，對(duì)產(chǎn)品的代碼結(jié)構(gòu)并不一定熟悉，因此，選取功能點(diǎn)數(shù)作為產(chǎn)品規(guī)模的度量元。

ERP產(chǎn)品中功能點(diǎn)數(shù)根據(jù)需求計(jì)劃書確定。最終功能點(diǎn)數(shù)Nnm即為測(cè)試模塊的個(gè)數(shù)。功能點(diǎn)之間存在相互制約和無(wú)約束兩種關(guān)系。存在制約關(guān)系的功能點(diǎn)記為：mi∧mj（i，j＜Nnm）；無(wú)約束關(guān)系記為：mx∨my（x，y＜Num）。子測(cè)試模塊集Mn包括多個(gè)測(cè)試模塊，至此得到測(cè)試模塊M：

由上式可得測(cè)試模塊及模塊間關(guān)系。

2 ERP測(cè)試方法CP建模

2.1 約束滿足問(wèn)題

約束滿足問(wèn)題是由一組預(yù)先定義的變量、變量的值域和變量間的約束關(guān)系組成，對(duì)約束滿足問(wèn)題的求解意味著在滿足所有約束條件下對(duì)每個(gè)變量指派一個(gè)合法的域值。其形式化為：?jiǎn)栴}P是一個(gè)三元組P＝（X，D，C）。其中，X＝｛x1，x2，…，xn｝是CSP中的變量集合；D＝｛D1，D2，…，Dn｝是論域的集合，其中，Di∈D對(duì)應(yīng)xi∈X，是變量xi的有限值域；C＝｛c1，c2，…，ck｝是一個(gè)有限約束集合，其中，任意cj∈C表示變量取值之間的制約關(guān)系［6］。其解為每個(gè)變量賦其論域上的一個(gè)值，使之同時(shí)滿足所有的約束［7－10］。

約束程序綜合了運(yùn)籌學(xué)和人工智能的知識(shí)，可以用來(lái)求解組合優(yōu)化問(wèn)題。約束指一個(gè)包含若干變量的關(guān)系表達(dá)式，用以表示這些變量必須滿足的條件［11－12］。求解CSP的基本思想是定義問(wèn)題領(lǐng)域的變量，每個(gè)變量預(yù)先設(shè)置各自的值域，通過(guò)描述各變量之間的約束，尋找滿足全局所有約束的一個(gè)或多個(gè)解［13－16］。ERP測(cè)試計(jì)劃方法問(wèn)題涉及到測(cè)試模塊的測(cè)試時(shí)間、測(cè)試人員對(duì)測(cè)試工作的能力值、測(cè)試人員對(duì)模塊的熟悉程度、測(cè)試模塊并發(fā)之間的平衡，具有NP-難的復(fù)雜性。

2.2 CP建模

2.2.1.ERP程序的變量描述

根據(jù)飼料行業(yè)ERP記賬憑證模塊測(cè)試的實(shí)際情況，確定測(cè)試分配方法中的變量如下：

1）測(cè)試工作開始時(shí)間：Ts。

2）測(cè)試模塊交付時(shí)間：Te。

3）測(cè)試模塊集合：M。

4）每個(gè)測(cè)試模塊在理想狀態(tài)下的測(cè)試工時(shí)：W＝｛w1，w2，…，wi｝。

5）測(cè)試人員集合：P＝｛p1，p2，…，pz｝。

6）測(cè)試人員P對(duì)測(cè)試中關(guān)注項(xiàng)T＝｛t1，t2，…，tk｝的能力值Q＝｛q1，q2，…，qz｝。

7）測(cè)試人員P對(duì)測(cè)試模塊mi熟悉程度：Fzi。

下面對(duì)上述變量進(jìn)行化簡(jiǎn)。每個(gè)測(cè)試模塊mi的測(cè)試開始時(shí)間Tsi有兩種情況：

1）整個(gè)測(cè)試工作的開始時(shí)間Ts即為模塊測(cè)試的開始時(shí)間Tsi＝Ts；

2）根據(jù)前置測(cè)試模塊的測(cè)試結(jié)束時(shí)間Te

i－1確定Tsi。每個(gè)測(cè)試模塊的測(cè)試工時(shí)并不確定，工時(shí)的長(zhǎng)短取決于測(cè)試人員對(duì)該測(cè)試模塊的熟悉程度Fzi，以及測(cè)試人員的個(gè)人能力值T。

根據(jù)飼料行業(yè)ERP實(shí)際情況，測(cè)試人員對(duì)測(cè)試主要關(guān)注項(xiàng)T的能力值見表1（完全熟悉記作1.0；本測(cè)試團(tuán)隊(duì)人員為5）。

表1 測(cè)試人員對(duì)測(cè)試過(guò)程主要關(guān)注項(xiàng)能力值

測(cè)試人員P的測(cè)試能力值用對(duì)各關(guān)注項(xiàng)T的加權(quán)平均表示。

測(cè)試工時(shí)W是指測(cè)試人員能力值以及測(cè)試人員對(duì)測(cè)試模塊的熟悉程度均為1.0時(shí)的測(cè)試用時(shí)。根據(jù)表1以及測(cè)試人員pz對(duì)測(cè)試模塊的熟悉程度Fzi按照公式wi／（qz×Fzi）可以得到測(cè)試人員P對(duì)每個(gè)測(cè)試模塊M進(jìn)行測(cè)試的工時(shí)，記為Wzi。

至此，變量3），5），6）化簡(jiǎn)為測(cè)試人員對(duì)每個(gè)測(cè)試模塊的測(cè)試工時(shí)：Wzi。

2.2.2.ERP程序的約束描述

測(cè)試模塊時(shí)間順序約束。模塊測(cè)試開始時(shí)間Tsi應(yīng)該等于或晚于整個(gè)模塊測(cè)試開始時(shí)間Ts，記作：Tsi≤Ts；測(cè)試模塊的結(jié)束時(shí)間Tei應(yīng)該早于或等于模塊交付時(shí)間Te，記作：Te≤Tei；測(cè)試模塊mj應(yīng)該在測(cè)試模塊mi測(cè)試結(jié)束之后進(jìn)行測(cè)試，即滿足：Tsi＞Tsj∧Tei＞Tsj，其中i，j為模塊集個(gè)數(shù)。

測(cè)試模塊間并發(fā)約束。測(cè)試中的部分測(cè)試模塊的測(cè)試是可以同時(shí)進(jìn)行的。用M′＝｛m′1，m′2，…，m′i｝表示，其中M′模塊的個(gè)數(shù)大于等于2，并發(fā)模塊M′是依據(jù)測(cè)試過(guò)程不斷變化的。

個(gè)人能力約束。測(cè)試模塊人員分配應(yīng)根據(jù)測(cè)試人員測(cè)試關(guān)注項(xiàng)的能力值和測(cè)試人員對(duì)測(cè)試模塊的熟悉程度確定。

同一時(shí)刻、同一測(cè)試模塊只能由一位測(cè)試人員進(jìn)行測(cè)試。

3 模型求解算法

3.1 約束傳播

約束傳播技術(shù)是通過(guò)在搜索過(guò)程中不斷縮小搜索空間來(lái)簡(jiǎn)化問(wèn)題的求解難度［17］。該案例ERP憑證部分的測(cè)試方法通過(guò)使用約束傳播技術(shù)，可以及時(shí)、有效地去除部分無(wú)用值，進(jìn)而簡(jiǎn)化求解過(guò)程。

3.2 回溯

當(dāng)某變量被賦值后，則進(jìn)行約束傳播，此時(shí)約束傳播的值域縮減，這就使得未賦值集合中某些變量的值域?yàn)榭?，這就說(shuō)明此時(shí)的賦值并沒(méi)有可行解［18－20］。這個(gè)時(shí)候需采用回溯方法取得可行解。

在變量的集合和約束已知的前提下，CSP技術(shù)可以找到滿足所有約束的合理的解集。文獻(xiàn)［1］根據(jù)需求樹確定測(cè)試模塊集，文中的測(cè)試模塊則是根據(jù)功能點(diǎn)確定，這更符合ERP產(chǎn)品的行業(yè)特色。對(duì)模型求解算法如下：

輸入部分根據(jù)測(cè)試人員對(duì)測(cè)試模塊的熟悉程度和測(cè)試人員對(duì)測(cè)試關(guān)注項(xiàng)的能力值，得到測(cè)試人員對(duì)各測(cè)試模塊的測(cè)試工時(shí)；輸入模塊中包含由式（1）得到的測(cè)試模塊及關(guān)系。

算法開始，若測(cè)試模塊中包含未測(cè)試的模塊mi，則對(duì)其分配測(cè)試人員。若mi包含在并發(fā)測(cè)試模塊中，則對(duì)其進(jìn)行并發(fā)處理。為當(dāng)前模塊mi分配測(cè)試人員后需要進(jìn)行約束傳播，若發(fā)生沖突，則進(jìn)行回溯操作，直至找到可行解并輸出。

3.3 并發(fā)測(cè)試模塊

并發(fā)問(wèn)題在測(cè)試工作中比較常見［20］，本ERP憑證部分的測(cè)試顯然存在并發(fā)問(wèn)題。測(cè)試開始時(shí)，對(duì)于手工填制憑證、采購(gòu)轉(zhuǎn)憑證、生產(chǎn)轉(zhuǎn)憑證、銷售轉(zhuǎn)憑證是可以同時(shí)發(fā)生的，對(duì)于并發(fā)測(cè)試模塊集的處理如下。

選取并發(fā)測(cè)試模塊中理想狀態(tài)下測(cè)試工時(shí)最大的子測(cè)試模塊Max（W），在測(cè)試人員中選取對(duì)當(dāng)前測(cè)試模塊測(cè)試工時(shí)最小Min（Wzi）的測(cè)試人員，若該測(cè)試人員資源被占用，則回溯。若沒(méi)有未分配任務(wù)的測(cè)試人員，則未分配測(cè)試人員的測(cè)試模塊處于等待狀態(tài)，等待最近一次結(jié)束的測(cè)試人員，其它模塊則開始測(cè)試，此時(shí)測(cè)試過(guò)程轉(zhuǎn)化為順序執(zhí)行測(cè)試。

4 模型求解方法示例與分析

對(duì)于ERP系統(tǒng)而言，財(cái)務(wù)是整個(gè)系統(tǒng)的核心，企業(yè)的采購(gòu)、生產(chǎn)、庫(kù)存狀況最終都將體現(xiàn)在財(cái)務(wù)系統(tǒng)，財(cái)務(wù)數(shù)據(jù)的準(zhǔn)確性可直觀地體現(xiàn)企業(yè)的經(jīng)濟(jì)效益，因此，財(cái)務(wù)系統(tǒng)成為企業(yè)人員最為關(guān)心的部分，對(duì)于測(cè)試人員而言，在業(yè)務(wù)流程穩(wěn)定的前提下，測(cè)試的重點(diǎn)理應(yīng)放在財(cái)務(wù)系統(tǒng)上。因此，這里選取飼料行業(yè)ERP財(cái)務(wù)系統(tǒng)的核心記賬憑證模塊進(jìn)行測(cè)試。該財(cái)務(wù)模塊分為兩部分：手工填制的憑證和業(yè)務(wù)數(shù)據(jù)的轉(zhuǎn)憑證。手工填制憑證主要對(duì)憑證中各部分內(nèi)容進(jìn)行填制、對(duì)憑證進(jìn)行審核簽字、憑證的打印設(shè)置、打印及刪除功能；業(yè)務(wù)數(shù)據(jù)的轉(zhuǎn)憑證主要是將業(yè)務(wù)部分的采購(gòu)單、生產(chǎn)單、銷售單轉(zhuǎn)為記賬憑證，根據(jù)需求，這里的轉(zhuǎn)憑證是一單一轉(zhuǎn)，因此，對(duì)于這三大部分的轉(zhuǎn)憑證應(yīng)分為3個(gè)模塊，而每部分所轉(zhuǎn)的記賬憑證也具有審核簽字、打印設(shè)置、打印以及刪除功能。至此得到測(cè)試模塊集M：

飼料行業(yè)ERP有測(cè)試人員5人，模塊測(cè)試標(biāo)準(zhǔn)工時(shí)為48h，測(cè)試人員對(duì)測(cè)試模塊的熟悉程度Fzi見表2。

表2 測(cè)試人員對(duì)測(cè)試模塊熟悉度

每個(gè)測(cè)試模塊的標(biāo)準(zhǔn)測(cè)試工時(shí)W見表3。

表3 每個(gè)測(cè)試模塊的標(biāo)準(zhǔn)測(cè)試工時(shí) W／h

由表1～表3可得測(cè)試人員對(duì)每個(gè)測(cè)試模塊的工時(shí)，見表4。

表4 測(cè)試人員對(duì)測(cè)試模塊的工時(shí) W／h

如此得到一個(gè)可行解，如圖1所示。

圖1 基于約束滿足的測(cè)試分配方法圖

圖中，橫坐標(biāo)表示測(cè)試時(shí)間，縱坐標(biāo)表示測(cè)試模塊。由圖1可知，記賬憑證部分測(cè)試完成為22h。

較之基于經(jīng)驗(yàn)的測(cè)試分配方法，文中提出的方法不僅縮短了完成工期，而且測(cè)試更全面。測(cè)試人員基于經(jīng)驗(yàn)的測(cè)試方法和文中測(cè)試方法在時(shí)間和測(cè)試全面性（bug數(shù)目為依據(jù)）的對(duì)比情況分別如圖2和圖3所示。

圖2 測(cè)試人員測(cè)試用時(shí)對(duì)比圖

圖3 測(cè)試人員測(cè)試全面性對(duì)比圖

每個(gè)測(cè)試人員的測(cè)試工時(shí)縮短的同時(shí)，測(cè)試更具全面性、科學(xué)性，這就降低了模塊中出現(xiàn)異常的可能性，為整個(gè)產(chǎn)品的交付提供了保障。

5 結(jié) 語(yǔ)

根據(jù)ERP本身測(cè)試的特點(diǎn)，將測(cè)試方法問(wèn)題模為約束求解問(wèn)題，并提出了一種全新的基于束滿足的測(cè)試方法。將該方法應(yīng)用于飼料行業(yè)P記賬憑證模塊的測(cè)試方法取得了較好的效。目前，有很多不同類型、不同性質(zhì)的軟件產(chǎn)，由于不同軟件產(chǎn)品對(duì)軟件測(cè)試所關(guān)注點(diǎn)以及試人員能力值差異等問(wèn)題，約束集合的確定還要根據(jù)具體情況而定，這就需要進(jìn)一步對(duì)約束足問(wèn)題進(jìn)行研究，根據(jù)自身特點(diǎn)構(gòu)造出合適的件測(cè)試方法。

［1］李悅，李娟，陳偉.基于約束滿足的測(cè)試計(jì)劃方法［J］.計(jì)算機(jī)工程與設(shè)計(jì)，2008（3）：535-539.

［2］柏淑琴.高校排課問(wèn)題的約束滿足優(yōu)化模型與算法［J］.科技視界，2012（18）：39-42.

［3］Freuder E C，Mackworth A K.Constraint Satisfaction：An Emerging Paradigm.RossiF，van Beek P，Walsh T，eds.Handbook of Constraint Programming［C］／／Amsterdam：Elservier，2006：13-28.

［4］李宏博，李占山，王濤.改進(jìn)求解約束滿足問(wèn)題粗粒度弧相容算法［J］.軟件學(xué)報(bào)，2012，23（7）：1816-1823.

［5］Minton S，Johnston M D，Laird P.Minimal conflicts：A heuristic repair method for constraint satisfaction and scheduling problems［J］.Artificial Intelligence，1992，58（1）：161-206.

［6］Bessière C，Cardon S，Debruyne R，et al.Efficient algorithms for singleton arc consistency［J］.Constraints，2011，16（1）：25-50.

［7］Bessière C.Constraint Propagation.Rossi F，van Beek P，Walsh T，eds.Handbook of Constraint Programming［C］／／Amsterdam：Elservier，2006：29-84.

［8］Freuder E C，Mackworth A K.Constraint Satisfaction：An Emerging Paradigm［C］／／Rossi F，van Beek P，Walsh T，eds.Handbook of Constraint Programming.Amsterdam：Elservier，2006：13-28.

［9］Bessière C，Regin J C，Yap RHC，et al.An optimal coarse-grained arc consistency algorithm［J］.Artificial Intelligence，2005，165（2）：165-185.

［10］Jones C.Software assessments，benchmarks，and best practices［M］.Beijing：China Machine Press，2003.

［11］Ivar Jacobsen，Grady Booch，James Rumbaugh.The unified software development process［M］.Beijing：China Machine Press，2002.

［12］James Rumbaugh，Ivar Jacobsen，Grady Booch.The unified modeling language reference manual［M］.USA：Addison Wesley Longman Inc，1999.

［13］Lecoutre C，Boussemart F，Hemery F.Exploiting multidirectionality in coarse-grained arc consistency algorithms［C］／／In：Proc.of the CP 2003.Kinsale，2003.

［14］Mehta D，Van Dongen MRC.Reducing checks and revisions in coarse-grained MAC algorithms［C］／／In：Proc.of the IJCAI，2005.

［15］Wallace R J，F(xiàn)reuder E C.Ordering heuristics for arc consistency algorithms［C］／／In：Proc.of the 9th Canadian Conf.on Artificial Intelligence.Vancouver，1992.

［16］Van Dongen M R C.Saving support checks does not always save time［J］.Artificial Intelligence Review，2004，21（3／4）：125-132.

［17］Li Z S，Li H B，Zhang Y G，et al.An approach of solving constraint satisfaction problem based on cycle-cut［J］.Chinese Journal of Computers，2011，34（8）：235-242.

［18］Boussemart F，Hemery F，Lecoutre C，et al.Boosting systematic search by weighting constraints［C］／／In：Proc.of the 16th European Conf.on Artificial Intelligence.Valencia，2004.

［19］Smith B M，Dyer M E.Locating the phase transition in binary constraint satisfaction problems［J］.Artificial Intelligence，1996，81（1）：565-582.

［20］SEI Software.Engineering measurement and analysis［EB／OL］.http：／／www.sei.cmu.edu／sema／welcome.html，2001.

Constraint-satisfaction based method for ERP test plan

WANG Jing1， LIU Bing1， LI Hui2， WANG Tao1＊
（1.School of Computer Science ＆Engineering，Changchun University of Technology，Changchun 130012，China；2.Jili Province Plastics Research Institute，Changchun 130022，China）

Combing the constraint programming with the characteristics of the ERP products，we put forward a ERP test plan based on constraint-satisfaction.The mathematical model of the vouchers in ERP of a feed industry is established.The test plan has be applied in the ERP platform of the feed industry.

constraint satisfaction；Enterprise Resource Planning（ERP）；test plan；constraint modeling.

N 39

A

1674-1374（2014）01-0071-06

2013-05-21

吉林省科技廳重點(diǎn)項(xiàng)目（20110361）

王 靜（1987－），女，漢族，山東萊蕪人，長(zhǎng)春工業(yè)大學(xué)碩士研究生，主要從事約束求解方向研究，E-mail：wj1grace＠163.com.＊通訊作者：王 濤（1969－），女，漢族，吉林長(zhǎng)春人，長(zhǎng)春工業(yè)大學(xué)副教授，碩士，主要從事約束求解方向研究，E-mail：wangtao＠m(xù)ail.ccut.edu.cn.