朱 柯

一、研究背景

邢瓷以“類雪”“類銀”名震古今、蜚聲中外。上世紀(jì)八十年代初，邢窯遺址 一經(jīng)考古發(fā)現(xiàn)即引起各界關(guān)注。目前在河北省的臨城、內(nèi)丘和邢臺(tái)市區(qū)都有與制瓷有關(guān)的遺址被相繼發(fā)現(xiàn)。這些遺址大致沿著京廣線以西、作南北縱向分布。根據(jù)目前已發(fā)表的考古材料可知邢窯燒制始于北朝，盛于唐，衰于五代，延續(xù)至宋、元、明、清等朝代。從歷次考古發(fā)掘出土的遺物來看， 碗類是邢窯生產(chǎn)規(guī)模較大、數(shù)量較多以及延續(xù)時(shí)間較長(zhǎng)的一種器物，具有較強(qiáng)的代 表性。在對(duì)隋代邢窯碗類材料的整理時(shí)，筆者發(fā)現(xiàn)墓葬的發(fā)掘材料一般不能確定其是實(shí)用器抑或?yàn)槊髌?。而窯址發(fā)現(xiàn)的考古材料一般數(shù)量較大，不但包含實(shí)用器，也會(huì)包含明器。 因此以窯址為主體，能更好的反映社會(huì)的生產(chǎn)使用概況

3。

隋代邢窯遺址發(fā)掘材料比較豐富，1997年6月，邢臺(tái)市區(qū)順德路北段的一處施工工地中，意外地發(fā)現(xiàn)了一處邢窯窯址。邢臺(tái)市文物管理處隨即對(duì)該遺址進(jìn)行了搶救性發(fā)掘，出土了數(shù)量龐大的瓷器、窯具等遺物，其種類主要有粗白瓷和黑釉瓷等，器型有碗、杯、瓶、盂、缽、罐、盅等，與周邊地區(qū)已發(fā)掘的邢窯窯址中出土的遺物基本一致，如臨城縣陳劉莊和內(nèi)丘縣西關(guān)北、中豐洞窯址等。邢臺(tái)市區(qū)順德路北段的這處遺址，其發(fā)現(xiàn)的瓷碗數(shù)量較大，器型標(biāo)準(zhǔn)，本文對(duì)隋代邢窯碗類的研究主要采用這處遺址發(fā)掘的瓷碗。出土的碗類大致有白釉、黑釉、黃釉（青釉）三大類。用英文大寫字母W、B、Y 依次表示，編號(hào)為 W1-50、 B1-29、Y1-25。其它相關(guān)資料將會(huì)在下文中逐一介紹。

20世紀(jì)60年代以后，隨著新考古學(xué)的興起，美國類型學(xué)研究發(fā)生質(zhì)的變化，就是在類型學(xué)的研究中引入數(shù)學(xué)精確計(jì)算方法，如貝氏統(tǒng)計(jì)、聚類分析、主成分分析、要素分析等方法。相對(duì)于國外，我國考古學(xué)中使用定量分析較晚。90年代后，一些大學(xué)和地方考古教學(xué)研究部門逐漸同國外加強(qiáng)合作，通過多種抽樣調(diào)查方案對(duì)包括日照、赤峰地區(qū)、以及河南伊洛河流域、靈寶和洛陽盆地，山西垣曲盆地，陜西周原和周公廟遺址等地進(jìn)行考古調(diào)查。應(yīng)用多種定量統(tǒng)計(jì)方法建立了相應(yīng)地區(qū)的考古地理信息系統(tǒng)。這種應(yīng)用多種定量統(tǒng)計(jì)方法的研究模式揭示了相應(yīng)時(shí)空中人與自然環(huán)境的關(guān)系、考古學(xué)文化的演化以及不同文化間的交流和碰撞，進(jìn)一步推動(dòng)定量考古學(xué)的發(fā)展。使得考古工作者對(duì)定量研究和統(tǒng)計(jì)學(xué)觀念和方法的認(rèn)同日益增強(qiáng)。同時(shí)期的學(xué)術(shù)研究中出現(xiàn)一大批貫穿統(tǒng)計(jì)學(xué)觀念和定量分析方法的專著。

二、研究方法

本文旨在通過使用數(shù)學(xué)定量方法對(duì)邢臺(tái)市區(qū)順德路北段隋代邢窯的部分碗類進(jìn)行研究，探究它們的輪廓造型特征。以往的陶瓷研究中對(duì)器物參數(shù)的選擇一般包括口徑、腹徑、通高等?！督y(tǒng)計(jì)分析在商代陶鬲研究中的應(yīng)用》一文中就是選擇商代陶鬲最大腹徑、通高、襠高、實(shí)足高作為參數(shù)，而吳雋在《古瓷類文物器型結(jié)構(gòu)的數(shù)字化特征鑒定》中研究五代、宋的撇口碗則選取了口徑、通高外還選擇了口沿處斜率、拐點(diǎn)。

在本文中，不僅采用了以往類型學(xué)研究中常用的器物參數(shù)度量，包括：口徑、底徑、高等，還有基于Matlab軟件對(duì)研究對(duì)象邊緣檢測(cè)和函數(shù)擬合獲得斜率方面相關(guān)參數(shù)，如：平均斜率參數(shù)以及斜率比值的關(guān)系等。利用這一參數(shù)，對(duì)不同種類的器物進(jìn)行細(xì)致的研究，以實(shí)現(xiàn)不同種類器物的數(shù)學(xué)定量描述?；谶@些描述，可望探討不同碗類瓷器制作工藝以及當(dāng)時(shí)手工業(yè)的行業(yè)標(biāo)準(zhǔn)。關(guān)于這一軟件，已有學(xué)者在對(duì)陶瓷器型進(jìn)行分析時(shí)使用過。

由于碗類的造型相對(duì)簡(jiǎn)單，可選擇的數(shù)據(jù)比較少，本文盡可能多的選擇度量數(shù)據(jù)。其中，基礎(chǔ)數(shù)據(jù)中以口徑、底徑、通高為主，然而這些參數(shù)只能孤立地代表碗的口沿、底、高度特征，不能反映出器物輪廓線的變化特征。因此基于上述分析筆者考慮將研究對(duì)象輪廓線的平均斜率也就是（口徑-底徑）/高度納入本文的研究范圍。此參數(shù)在一定程度上反映了考察碗類造型的平均外撇情況。為了更好地反映出器物輪廓線的具體特征，我們擬合輪廓線的函數(shù)解析式，考察它們輪廓線的斜率變化。研究中先求出它們的拐點(diǎn)，觀察它們的輪廓線總體的凹凸變化情況，最后通過對(duì)輪廓線上下半段斜率的對(duì)比情況來描述輪廓線形狀特征。（圖1）

圖1 參數(shù)示意圖

本文研究對(duì)象的參數(shù)來源于兩方面：一種是通過考古報(bào)告中的基礎(chǔ)統(tǒng)計(jì)獲得，包括口徑、底徑、通高、足高等，這些數(shù)據(jù)一般在考古報(bào)告或其它資料中都會(huì)出現(xiàn)。本文對(duì)這些數(shù)據(jù)進(jìn)行了收集整理。另一種是通過使用軟件處理獲得。我們使用excel軟件計(jì)算平均斜率；采用Matlab軟件對(duì)研究對(duì)象的電子線圖進(jìn)行了邊緣檢測(cè)和函數(shù)擬合，通過對(duì)擬合函數(shù)進(jìn)行求解獲得研究對(duì)象輪廓線不同點(diǎn)的斜率值，基于平均斜率以及輪廓線的不同斜率值的統(tǒng)計(jì)，進(jìn)一步進(jìn)行數(shù)學(xué)定量分析。簡(jiǎn)要示意可見圖2。

對(duì)于處理數(shù)據(jù)的獲得，這里要介紹一下Matlab邊緣檢測(cè)的算法。我們使用sobel算子算法進(jìn)行邊緣檢測(cè)。該方法檢測(cè)出的圖像的邊緣的圖像，對(duì)噪聲具有平滑作用，可以提供比較精確的邊緣方向的信息。至于更為具體的數(shù)學(xué)分析過程，因篇幅原因，在此略去。

三、數(shù)據(jù)分析

（一）類型學(xué)初步劃分

邢臺(tái)市區(qū)順德路北段窯址出土的碗類，按釉色主要分為三類：白釉、黑釉、黃釉。原報(bào)告將三類瓷碗做了以下類型劃分。白釉瓷碗基本分為兩個(gè)類型：A型為深腹碗，高8厘米左右，口大，腹淺，腹壁較直，體形顯瘦高；B型為淺腹碗，高在6厘米左右，體形較A型矮，口大，腹淺，弧腹，略成變長(zhǎng)方體。深腹碗的編號(hào)為1-40，淺腹碗為41-50。黑釉瓷碗也被劃為兩個(gè)大的類型：編號(hào)1到23，是A型深腹碗，高在8厘米左右，口大，腹深，體形瘦高。編號(hào)24到29是B型淺腹碗，高在6厘米左右，體形比A矮，口大，腹淺，略成扁長(zhǎng)方體。其中A型深腹碗根據(jù)底足的不同分為兩個(gè)亞型Aa和Ab，主體底足以上輪廓差別不大。黃釉（青釉）瓷碗，同樣劃分為兩個(gè)類型：1到19是A型深腹碗，高在8厘米左右，體形顯得瘦高，尖圓唇，斜弧腹，餅狀足，組緣斜削一周。20到25是B型淺腹碗，高在6厘米左右，器型大體與A相似，但比A矮，口大，腹淺，略成扁長(zhǎng)方體。

（二）輪廓線研究

輪廓線是指器物的輪廓曲線，由于隋代邢窯碗、碟、罐等對(duì)稱性器物基本采用輪盤拉坯成型，屬于所謂的“旋轉(zhuǎn)體”，所以，輪廓線包含了豐富的器物特征信息的。輪廓曲線的平均斜率表示該曲線在在所有點(diǎn)上的平均彎曲程度，可以在大體上表征其口、腹的外撇程度；而輪廓線的擬合則可以在細(xì)節(jié)上描述輪廓線的具體走勢(shì)。輪廓線的具體研究主要考慮的是基于類型學(xué)同類型間的關(guān)系，以及在輪廓線研究的基礎(chǔ)上類型的進(jìn)一步劃分。利用Matlab對(duì)瓷碗的進(jìn)行定量研究，并非意圖否定傳統(tǒng)類型學(xué)研究成果，而是希望通過定量分析的嘗試，以更客觀的方式為瓷器的斷代、演變及制作工藝提供一些思考。

本文所選碗類的分析標(biāo)本，其白、黑、黃釉三類瓷碗的基礎(chǔ)數(shù)據(jù)見附表1。附表1中的平均外撇程度為（口徑-底徑/（通高-足高）的值。

對(duì)于以上三種釉色的瓷碗首先按照單一釉色進(jìn)行研究，之后相互比較綜合討論。

由圖3可以看出白釉瓷碗中編號(hào)前40的平均外撇值基本上在1上下浮動(dòng)，除了個(gè)別異常點(diǎn)最大值都不超過1.2；40-50的平均外撇值則全部在1.2以上。很明顯白釉瓷碗在外撇值上面可以分為兩類。這種分類與原報(bào)告中類型的劃分相一致。白釉瓷碗數(shù)量巨大而且表現(xiàn)出較好的規(guī)律性，這一點(diǎn)可以在一定程度表征隋朝白釉瓷碗的規(guī)范化生產(chǎn)情況。規(guī)范化生產(chǎn)代表了精細(xì)化的一方面，精細(xì)化的方向預(yù)示著白瓷的興盛。

圖3 白釉瓷碗平均斜率分布圖

根據(jù)圖4可以看出，原報(bào)告中劃分的兩種類型在分布圖中差別不夠明顯。所有的碗的平均斜率大致都分布在1.15線上下。在以口徑和高為參考的劃分中，淺腹碗高度在6厘米左右，相比深腹碗的8厘米左右這個(gè)差值是非常明顯的。當(dāng)碗的高度降低時(shí)平均斜率理想狀況下是應(yīng)該變大的。

圖4 黑釉瓷碗平均斜率分布圖

從圖5可以看出黃釉瓷碗的平均斜率大致分布在“1.1線”上下，沒有呈現(xiàn)出規(guī)律的分類情況，也與原報(bào)告在高和口徑中劃分的類別是不一致的。它和黑釉瓷碗一樣沒有在深腹和淺腹碗的差別上表現(xiàn)出理想的規(guī)律性。

圖5 黃釉瓷碗平均斜率分布圖

以上以平均斜率考察時(shí)都沒能很好的描述輪廓線的具體情況，因此進(jìn)一步擬合輪廓線的函數(shù)考察輪廓線的變化。對(duì)輪廓線采用輪廓線邊緣檢測(cè)后的四次多項(xiàng)式擬合。在獲得相應(yīng)函數(shù)的情況下計(jì)算參數(shù)進(jìn)行統(tǒng)計(jì)分析。

對(duì)三類瓷碗的擬合函數(shù)系數(shù)列表見附表2。由附表2中數(shù)據(jù)可知，三類碗的第一種類型都是數(shù)量比較多且多少不同。為了保證嚴(yán)謹(jǐn)性，統(tǒng)一取以上三類碗中A型碗的編號(hào)前十五的標(biāo)本進(jìn)行研究。B型碗相對(duì)較少，全部都進(jìn)行統(tǒng)計(jì)研究。統(tǒng)一對(duì)這些碗類輪廓線擬合的多項(xiàng)式函數(shù)進(jìn)行一次、二次、三次求導(dǎo)，得到口沿處斜率、中點(diǎn)斜率、起始點(diǎn)斜率、拐點(diǎn)。對(duì)相關(guān)參數(shù)統(tǒng)計(jì)見附表3。

對(duì)附表3中的數(shù)據(jù)進(jìn)行篩選，按照上下半段比值從小到大的順序排列見附表4。

根據(jù)重新的排序發(fā)現(xiàn)其輪廓線的類別按數(shù)值劃分與之前的傳統(tǒng)類型劃分大不相同。

圖6中四只碗按行排列分別是W12、W46、W49、B1，它們的上段對(duì)下段比值介于0-1之間，說明其上半段的斜率值變化不明顯，也就是擬合曲線各點(diǎn)斜率之間比較接近，所以比較趨于直線化。如果單以輪廓線論，而不考慮碗的高矮口徑問題，這四只碗更符合劃分為同一類的標(biāo)準(zhǔn)。之所以存在高矮胖瘦的不同在于它們的口徑高不同，以上四只碗考察口徑與高的比值依次為1.707317073、2.4、2.259259259、1.575。顯然W12與B1的值比較接近，W46與W49的值比較接近。通過圖中可以明顯看出W12和B1的輪廓基本相同，W46與W49的輪廓基本相同。因此可以說在輪廓線上下段比值在0-1時(shí)而口徑與高的比值比較接近時(shí)，它們的輪廓造型基本相同，可以歸為一類。而且這種劃分也是與考古報(bào)告的的劃分相一致的，只是顯得比較細(xì)致。

圖6 斜率比值0-1階段碗類圖

圖7中四只碗按行排列分別是W1、B3、B10、W13，它們的上下段斜率比值介于1-2之間，可以明顯的看出這幾只碗的腹部不明顯，輪廓曲線顯得上下較為一致、舒緩。由于口沿處斜率不是很大，大多數(shù)可以稱為敞口形狀，或者說是直口微敞。同樣考慮它們的口徑與高的比值依次為1.794871795、1.581395349、2、1.707317073。顯 然 除 了B10的值外，其余三者的差別不是很大，因此在輪廓線相似，高矮胖瘦差別不大的情況下，其余三者可以很好地歸為一類，同樣與報(bào)告中的型式劃分一致又顯得比較細(xì)致。

圖7 斜率比值1-2階段碗類圖

在上下段斜率差值比處于2以上并不斷增大時(shí)，它們的輪廓線造型也表現(xiàn)出明顯的變化趨勢(shì)。這里選擇B4、B7、B14、B2、Y20、B28、Y24、B6這幾只碗，其上下段比值依次為 2.3750、 3.6412、 4.8013、 5.5501、6.6250、 7.4881、8.4000、9.4918。

明顯可見，圖8中的瓷碗的口沿逐漸向直口靠攏，腹部逐漸明顯，由于比值增大，上部的斜率在同樣的橫坐標(biāo)長(zhǎng)度展開的話需要一定高度保證，所以對(duì)比之后造型會(huì)向深腹表現(xiàn)。可以大膽地預(yù)測(cè)隨著這一比值的增大，這些碗類輪廓線將基出現(xiàn)直口、甚至直口微斂，球形腹、鼓腹的情況。例如Y15、W47基本上就成這種趨勢(shì)。（圖8）

圖8 斜率比值在2以上的碗類圖

圖9 斜率階段分布驗(yàn)證示意圖

同樣考慮它們的口徑和高的比值依次為：1.910447761、 1.934210526、 無 、1.642857143、 1.538461538、 2、2.028985507、1.666666667。從中可以看出在輪廓線向著直口甚至微斂的情況變化時(shí)，當(dāng)標(biāo)本的口徑與高度比值接近時(shí)它們的造型是非常相似的，可以很好地歸為一類。我們也從中發(fā)現(xiàn)了一個(gè)異常的歸類對(duì)象Y20，在輪廓線相似的情況下它的口徑與高比值1.538461538與Y24的2.028985507相 差 比 較大，而在考古報(bào)告中它們是歸為一類的。對(duì)于其它在輪廓線相似的變化趨勢(shì)，口徑與高的比值接近時(shí)可以很好地歸為一類。

四、小結(jié)

僅以口徑和高簡(jiǎn)單劃分以上碗類造型的類別時(shí)，可以比較明顯地分為深腹和淺腹類型。用平均斜率來考察，發(fā)現(xiàn)白釉瓷碗表現(xiàn)出很好的規(guī)律性，深腹碗的斜率基本上在1上下浮動(dòng)，除了個(gè)別異常點(diǎn)最大值都不超過1.2；淺腹碗的平均外撇值則全部在1.2以上。黑釉和黃釉瓷碗在平均斜率上都沒有表現(xiàn)出原報(bào)告中型式劃分的分類趨勢(shì)。黑釉瓷碗所有的碗的平均斜率大致都分布在1.15線上下波動(dòng)，黃釉瓷碗所有碗的平均斜率大致分布在1.1線上下波動(dòng)。使用邊緣檢測(cè)和函數(shù)擬合考察他們輪廓線的性質(zhì)。發(fā)現(xiàn)所有碗輪廓線的上下半段斜率差的比值呈現(xiàn)出分段的特點(diǎn)。大致可以分為0-1，1-2以及2以上，它們的輪廓線隨著這個(gè)分段呈現(xiàn)規(guī)律性的變化。當(dāng)這個(gè)比值在0-1時(shí)，處于這一區(qū)間段的碗類的輪廓線上半段的斜率比較接近，因此表現(xiàn)出趨直線化，下半段的腹部表現(xiàn)不明顯，顯示出敞口的特點(diǎn)。而處于這一區(qū)間的同色瓷碗在分類上當(dāng)口徑與高的比值（高矮胖瘦）確定時(shí)，這種階段性的劃分與原報(bào)告中的型式劃分也是一致的。當(dāng)這個(gè)比值處于1-2這一區(qū)間時(shí)，輪廓線表現(xiàn)出上下半段的輪廓線變化較為一致得特點(diǎn)，沒有明顯腹部，口部敞；但相對(duì)于0-1的區(qū)間輪廓線略向直口靠攏。當(dāng)確定它們的口徑與高的比值比較接近時(shí)，這一分類段內(nèi)標(biāo)本的輪廓造型（忽略底足釉色等，下同）基本是相同的，歸類也與原報(bào)告中的劃分一致。當(dāng)這個(gè)比值處于2以上時(shí)，所有碗類在口沿和腹部呈現(xiàn)一致的變化趨勢(shì)：口沿處趨于直口，甚至斂口；腹部明顯，趨于鼓腹，球形腹。如果考慮它們的口徑與高的比值比較接近時(shí)，這一階段的標(biāo)本的外輪廓造型基本相同，可以歸為一類。此時(shí)的劃分發(fā)現(xiàn)原報(bào)告中出現(xiàn)異常的劃分Y20，在輪廓線相似的情況下它的口徑與高比值1.538461538與Y24的2.028985507 相差比較大，視覺上它們?cè)谠煨蜕鲜遣煌模谠瓐?bào)告中它們是歸為一類的。

對(duì)于這一偏差，筆者認(rèn)為，黑釉和黃釉瓷碗的平均斜率存在波動(dòng)情況，推測(cè)是由于生產(chǎn)過程是人為操作的，存在不可預(yù)知的誤差，一定程度上說明了這些碗類生產(chǎn)的精度和規(guī)律性比較差，產(chǎn)生了在肉眼視覺上歸屬于同一類，而在數(shù)學(xué)精確度量時(shí)又存在與預(yù)期不完全相同的情況。

五、探討

綜上所述，可以清楚的看出白釉瓷碗的規(guī)律性比較強(qiáng)，黃釉瓷碗與黑釉瓷碗規(guī)律性不是特別強(qiáng)。從輪廓線考慮的話，所有瓷碗呈現(xiàn)出隨著上下段斜率比值分段，同時(shí)呈現(xiàn)出相應(yīng)的分段變化趨勢(shì)。而且這一分段趨勢(shì)大體上也是與原報(bào)告中類別的劃分是一致的。但是表現(xiàn)出數(shù)值劃分越精細(xì)，相對(duì)原報(bào)告中的分類也更精細(xì)的特點(diǎn)。利用Matlab等軟件進(jìn)行數(shù)學(xué)定量分析，可以作為傳統(tǒng)類型學(xué)研究的得力助手，更加精細(xì)化與系統(tǒng)地分辨與分析典型器物及其演變特征，以更好地完成考古類型學(xué)的任務(wù)。

對(duì)于以上變化背后的社會(huì)因素，有學(xué)者已經(jīng)就邢窯各時(shí)期生產(chǎn)青瓷、白瓷、黑釉瓷器型特征、燒造工藝以及變化原因進(jìn)行了諸多探討與研究。結(jié)合上文中的定量分析數(shù)據(jù)，筆者提出一點(diǎn)淺見，以求教大方。筆者認(rèn)為北朝時(shí)期北方瓷器的燒造來自于南方的傳入，此時(shí)北方瓷器燒造技術(shù)還未經(jīng)過長(zhǎng)時(shí)間的積淀，因此可能會(huì)在碗類上表現(xiàn)出較差的規(guī)范性；隋代的燒造工藝相對(duì)北朝已經(jīng)有了很大進(jìn)步，出現(xiàn)桶裝匣缽，而且可以控制窯爐火焰，生產(chǎn)上表現(xiàn)出了進(jìn)步，碗類的規(guī)范化就更加明顯。尤其是研究中，白瓷標(biāo)本的規(guī)范化顯得更加突出。而南青北白的原因是因?yàn)槟戏降那啻晒に囈呀?jīng)很成熟，北方要在這方面超越南方是很困難的，出于競(jìng)爭(zhēng)原因，促使北方白瓷進(jìn)一步發(fā)展，因此在白瓷碗上體現(xiàn)出很好的規(guī)律性。而隋代邢窯碗類之間在輪廓線上表現(xiàn)出很好的分段歸類趨勢(shì)，也在一定程度上說明了工藝的規(guī)范和進(jìn)步。比較遺憾的是隋代這批材料沒有具體年代劃分，不能從細(xì)致的類型劃分來同時(shí)表征相應(yīng)時(shí)代特征，但還是可以從整體上看出整個(gè)隋代的生產(chǎn)特點(diǎn)的。

注釋：

1.馮先銘.談邢窯有關(guān)諸問題[J].故宮博物院院刊,1981(4):49:50.

2.李恩瑋,張志忠,李軍.邢窯遺址研究[M].北京:科學(xué)出版社,2007.

3.楊文山.隋代邢窯遺址的發(fā)現(xiàn)和初步分析[J].文物,1984(12).

4.河北省文物管理所.邢臺(tái)隋代邢窯[M].北京:科學(xué)出版社,2006.

5.原報(bào)告中認(rèn)為是黃釉，但比對(duì)南方岳州窯、洪州窯等窯址的出土實(shí)物，筆者認(rèn)為該釉色可以歸類于青釉，為尊重原報(bào)告，這里仍沿用“黃釉”的說法。

6.湯惠生.考古類型學(xué)的誕生與發(fā)展——考古類型學(xué)系列研究之一[J].文博,2006(04):27-28.

7.陳鐵梅.定量考古學(xué)[M].北京大學(xué)出版社,2005.

8.同上.

9.吳雋等.古陶瓷類文物器型結(jié)構(gòu)的數(shù)字化特征鑒定[J].中國科學(xué):技術(shù)科學(xué),2012.09.

10.同上.

11.吳曦.基于MATLAB的圖像邊緣檢測(cè)算法的研究和實(shí)現(xiàn)[D].碩士論文.吉林大學(xué)，2014.05:14-15.

12.同4.

?

?

?

?

?

?

?

?