諧振式巨磁阻微機(jī)械陀螺噪聲分析*

褚偉航，白曉曉，蔣孝勇，李孟委

(1. 中北大學(xué) 電子測(cè)試技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，山西 太原 030051； 2. 中北大學(xué) 儀器與電子學(xué)院，山西 太原030051；3. 中北大學(xué) 微系統(tǒng)集成研究中心，山西 太原 030051)

?

諧振式巨磁阻微機(jī)械陀螺噪聲分析*

褚偉航1,2,3，白曉曉1,2,3，蔣孝勇1,2,3，李孟委1,2,3

(1. 中北大學(xué) 電子測(cè)試技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，山西 太原 030051； 2. 中北大學(xué) 儀器與電子學(xué)院，山西 太原030051；3. 中北大學(xué) 微系統(tǒng)集成研究中心，山西 太原 030051)

摘要:巨磁阻效應(yīng)微陀螺具有高靈敏、 低噪聲等優(yōu)點(diǎn). 本文詳細(xì)闡述了巨磁阻效應(yīng)微陀螺的工作原理，設(shè)計(jì)了微陀螺結(jié)構(gòu)，其驅(qū)動(dòng)方向頻率為5 553.7 Hz，檢測(cè)方向頻率為5 556.7 Hz，頻率差僅為3 Hz，實(shí)現(xiàn)了模態(tài)的良好匹配. 設(shè)計(jì)了應(yīng)用于該微陀螺中的檢測(cè)磁場(chǎng)，測(cè)試了一款巨磁阻器件的靈敏度，基于上述結(jié)果完成了巨磁阻微陀螺的系統(tǒng)級(jí)仿真，并分析了巨磁阻微陀螺的噪聲組成以及噪聲的計(jì)算. 仿真結(jié)果表明，將巨磁阻效應(yīng)應(yīng)用于陀螺的角速率信號(hào)檢測(cè)中，可以實(shí)現(xiàn)靈敏度為46 mV·((°)/s)-1，噪聲為 的信號(hào)檢測(cè).

關(guān)鍵詞:巨磁阻； 微機(jī)械陀螺； 系統(tǒng)級(jí)仿真； 噪聲分析

0引言

1工作原理及理論分析

1.1巨磁阻微陀螺的工作原理

如圖 1 所示為巨磁阻效應(yīng)微陀螺的工作原理. 首先，陀螺放置在由永磁體產(chǎn)生的勻強(qiáng)磁場(chǎng)中，然后，在驅(qū)動(dòng)導(dǎo)線上加載交變驅(qū)動(dòng)電流以產(chǎn)生交變洛倫茲力，使質(zhì)量塊在該驅(qū)動(dòng)力的作用下沿驅(qū)動(dòng)方向(X軸)往復(fù)振動(dòng). 當(dāng)有Y軸方向的角速率輸入時(shí)，質(zhì)量塊在柯氏力的作用下沿檢測(cè)方向(Z軸)運(yùn)動(dòng)，并帶動(dòng)磁體遠(yuǎn)離或靠近磁敏電阻. 由于二者間距發(fā)生了變化，使磁敏電阻敏感到的磁場(chǎng)發(fā)生變化，從而引起磁敏電阻中自旋相關(guān)的隧穿電流發(fā)生變化，進(jìn)一步導(dǎo)致磁敏電阻的阻值發(fā)生劇烈變化，通過測(cè)量阻值變化能夠?qū)崿F(xiàn)對(duì)微弱柯氏力地檢測(cè).

圖 1　磁阻隨間距變化原理圖Fig.1　The principle diagram of the magnetic resistance changing with distance

1.2微機(jī)械陀螺的動(dòng)力學(xué)方程

根據(jù)微機(jī)械陀螺的工作原理，它的簡(jiǎn)化動(dòng)力學(xué)模型如圖 2 所示.

圖 2　微陀螺的動(dòng)力學(xué)模型Fig.2　Dynamic model of micro gyroscope

微機(jī)械陀螺的動(dòng)力學(xué)方程為

(1)

(2)

式中：F0為驅(qū)動(dòng)力幅值；mx,mz分別為微陀螺驅(qū)動(dòng)質(zhì)量與檢測(cè)質(zhì)量；kx，cx分別為驅(qū)動(dòng)模態(tài)的剛度與阻尼；kz，cz分別為檢測(cè)方向的剛度與阻尼. 求解動(dòng)力學(xué)方程，可以得到微陀螺驅(qū)動(dòng)方向的穩(wěn)態(tài)幅值為

(3)

式中：ωx為驅(qū)動(dòng)模態(tài)角頻率；ξx為驅(qū)動(dòng)模態(tài)阻尼比.

檢測(cè)方向的穩(wěn)態(tài)幅值為

(4)

式中：ωz為驅(qū)動(dòng)模態(tài)角頻率；ξz為驅(qū)動(dòng)模態(tài)阻尼比.

微陀螺的結(jié)構(gòu)靈敏度為

(5)

1.3巨磁阻微陀螺的靈敏度分配

根據(jù)隧道磁阻微陀螺原理及其系統(tǒng)構(gòu)成，我們可以將隧道磁阻微陀螺的靈敏度分為三部分： 結(jié)構(gòu)靈敏度、 磁場(chǎng)靈敏度以及磁阻靈敏度. 信號(hào)轉(zhuǎn)換過程分別為角速率信號(hào)轉(zhuǎn)換為位移信號(hào)，位移信號(hào)轉(zhuǎn)換為磁信號(hào)，磁信號(hào)轉(zhuǎn)換為電信號(hào). 微陀螺靈敏度可按如下方式進(jìn)行分解，三者相乘可得微陀螺的靈敏度，即

(6)

通過對(duì)圖 3 中所示三部分靈敏度的規(guī)劃，可以實(shí)現(xiàn)對(duì)器件總靈敏度的規(guī)劃.

圖 3　靈敏度分配圖Fig.3　The diagram of distribution of the sensitivity

2巨磁阻微陀螺的靈敏度設(shè)計(jì)

2.1結(jié)構(gòu)靈敏度設(shè)計(jì)

利用有限元仿真軟件對(duì)微陀螺結(jié)構(gòu)進(jìn)行模態(tài)分析與動(dòng)力學(xué)分析，結(jié)果如圖 4 和表 1 所示.

圖 4　微陀螺模態(tài)仿真Fig.4　Mode simulation of micro gyro

模態(tài)頻率理論值/Hz頻率仿真值/Hz誤差/%一階模態(tài)5553.753064.46二階模態(tài)5556.753543.65

模態(tài)分析結(jié)果顯示，微陀螺模態(tài)頻率的理論值與仿真值誤差均在5%以內(nèi)，在可接受范圍內(nèi)，證明了所設(shè)計(jì)結(jié)構(gòu)的頻率匹配性較好. 仿真結(jié)果與理論計(jì)算結(jié)果之間的誤差主要產(chǎn)生原因?yàn)椋珹NSYS有限元分析軟件網(wǎng)格劃分的疏密會(huì)影響計(jì)算結(jié)果，若網(wǎng)格稀疏會(huì)致使有限元計(jì)算不精確； 再者，理論計(jì)算公式是對(duì)實(shí)際理論模型的抽象與簡(jiǎn)化，簡(jiǎn)化后的公式計(jì)算與實(shí)際模型相比會(huì)有一些誤差. 下面對(duì)微陀螺的振動(dòng)特性進(jìn)行分析，在驅(qū)動(dòng)方向加載與驅(qū)動(dòng)力等效的位移載荷，Y方向施加1 rad/s的角速率，利用ANSYS仿真得到檢測(cè)方向位移結(jié)果，如圖 5 所示.

由圖 5 可以得到，仿真得到微陀螺的檢測(cè)方向位移幅值為2.070 8×10-6m,利用式(5)可以求得

理論計(jì)算得到的微陀螺位移靈敏度為32.33 nm·((°)/s)-1，仿真值與理論值較為接近，說明結(jié)構(gòu)設(shè)計(jì)的合理性.

以上完成了微陀螺結(jié)構(gòu)靈敏度的設(shè)計(jì)，下面進(jìn)行磁場(chǎng)靈敏度以及磁阻靈敏度的設(shè)計(jì).

圖 5　檢測(cè)方向位移分析結(jié)果Fig.5　The displacement analysis results of Detection direction

2.2檢測(cè)磁場(chǎng)靈敏度的設(shè)計(jì)

磁場(chǎng)設(shè)計(jì)是巨磁阻微陀螺的關(guān)鍵內(nèi)容之一，本文設(shè)計(jì)被應(yīng)用于巨磁阻效應(yīng)微陀螺中的磁場(chǎng)，分驅(qū)動(dòng)磁場(chǎng)與檢測(cè)磁場(chǎng)兩部分. 驅(qū)動(dòng)磁場(chǎng)為勻強(qiáng)磁場(chǎng)，磁場(chǎng)強(qiáng)度可根據(jù)所需驅(qū)動(dòng)力大小來進(jìn)行選擇. 檢測(cè)磁場(chǎng)的設(shè)計(jì)為本文分析的重點(diǎn)，它應(yīng)具有兩個(gè)特點(diǎn)，首先在離面方向上應(yīng)具有較高的磁場(chǎng)變化率，可以實(shí)現(xiàn)微陀螺的高靈敏檢測(cè)； 其次驅(qū)動(dòng)時(shí)檢測(cè)磁場(chǎng)應(yīng)盡量接近為勻強(qiáng)磁場(chǎng)，避免與檢測(cè)時(shí)的磁場(chǎng)信號(hào)發(fā)生耦合干擾. 結(jié)合工藝條件與磁場(chǎng)分布特性，采用矩形磁體應(yīng)用于巨磁阻微陀螺中，圖 6 所示為矩形永磁體磁場(chǎng)分布的建模[7].

圖 6　矩形永磁空間外磁場(chǎng)分布模型Fig.6　External magnetic field distribution model of rectangular permanent magnet

定義Z方向?yàn)槲⑼勇莸臋z測(cè)方向，X方向?yàn)轵?qū)動(dòng)方向. 設(shè)矩形磁體的長(zhǎng)、 寬、 高分別為a、b、h，永磁體內(nèi)部任一點(diǎn)(x0,y0,z0)的z0與dz0構(gòu)成的磁體薄膜在其空間外任一點(diǎn)P(x,y,z)處，由電流環(huán)A′B′C′D′產(chǎn)生的磁場(chǎng)為dB，則空間中的磁場(chǎng)由3個(gè)方向的分量表示為[8]

(7)

(8)

(9)

其中

式中：Γ是以γ1，γ2，γ3為自變量的一個(gè)函數(shù)記號(hào)；Js為與平面xOy平行的任一平面上的面電流密度；μ0為真空磁導(dǎo)率；Ψi是以ψ1，ψ2，ψ3為自變量的一個(gè)函數(shù)記號(hào). 利用以上這三個(gè)式子，可以將一個(gè)磁體周圍空間中任意一點(diǎn)的磁場(chǎng)強(qiáng)度分解到X、Y、Z三個(gè)方向上，從而計(jì)算出磁場(chǎng)在各個(gè)方向上的分量.

針對(duì)所得磁場(chǎng)在三個(gè)方向分量公式，設(shè)計(jì)磁體尺寸為長(zhǎng)340 μm，寬200 μm，厚5 μm. 利用MATLAB軟件對(duì)矩形永磁體磁場(chǎng)分布特性進(jìn)行仿真分析，檢測(cè)磁場(chǎng)在驅(qū)動(dòng)方向分量的仿真結(jié)果如圖 7，圖 8 所示.

圖 7　驅(qū)動(dòng)時(shí)檢測(cè)磁場(chǎng)分布圖Fig.7　Magnetic field distribution of detection direction when driving

圖 8　檢測(cè)磁場(chǎng)在驅(qū)動(dòng)方向的磁場(chǎng)變化率Fig.8　The magnetic field change rate of detection magnetic field in the driving direction

從圖 7 中可以看出，在x方向除邊緣區(qū)域，檢測(cè)磁場(chǎng)在驅(qū)動(dòng)方向分量的磁場(chǎng)分布近乎為勻強(qiáng)，圖8所示其勻強(qiáng)范圍為100～240 μm. 所設(shè)計(jì)微陀螺驅(qū)動(dòng)梁的運(yùn)動(dòng)間隙為80 μm，所以巨磁阻器件在驅(qū)動(dòng)方向可布置在磁體的120～180 μm這個(gè)范圍內(nèi)，如圖 7 所示. 圖 8 中取接近勻強(qiáng)磁場(chǎng)的區(qū)域進(jìn)行線性擬合，得到在圖 7 所示勻強(qiáng)磁場(chǎng)范圍內(nèi)，磁場(chǎng)的變化率為0.005 4 Oe/nm，磁場(chǎng)變化率很小，下面分析檢測(cè)磁場(chǎng)在檢測(cè)方向的磁場(chǎng)變化率情況.

圖 9 所示為檢測(cè)磁場(chǎng)在檢測(cè)方向的磁場(chǎng)分布情況，從圖中可以看出，在磁體表面，磁場(chǎng)變化率為零，磁場(chǎng)強(qiáng)度最大. 接下來對(duì)檢測(cè)磁場(chǎng)在檢測(cè)方向分量進(jìn)行分析，磁場(chǎng)檢測(cè)方向分量Hz，在檢測(cè)方向分布曲線如圖 10 所示.

取磁場(chǎng)強(qiáng)度分布線性度較好的范圍(-130～-70 μm)對(duì)其進(jìn)行線性擬合，得到結(jié)果如圖 10 所示，擬合曲線的斜率即磁場(chǎng)靈敏度約為0.1 Oe/nm.

圖 9　檢測(cè)時(shí)檢測(cè)磁場(chǎng)分布圖Fig.9　Magnetic field distribution of detection direction when detecting

圖 10　檢測(cè)方向磁場(chǎng)變化率Fig.10　The magnetic field change rate of detection direction

至此，完成了對(duì)磁場(chǎng)靈敏度的設(shè)計(jì). 關(guān)于磁阻靈敏度的設(shè)計(jì)，本文選擇了一款巨磁阻器件對(duì)其靈敏度進(jìn)行測(cè)試，得到其磁阻靈敏度. 以實(shí)測(cè)磁阻靈敏度水平作為本文微陀螺的磁阻靈敏度.

2.3磁阻靈敏度的設(shè)計(jì)及測(cè)試

通過搭建實(shí)驗(yàn)，在磁屏蔽環(huán)境下分別測(cè)試了巨磁阻器件與各向異性磁阻器件的靈敏度并進(jìn)行了對(duì)比. 測(cè)試環(huán)境及使用儀器如圖 11 所示.

圖 11　磁阻測(cè)試環(huán)境及使用儀器Fig.11　Magnetic resistance test environment and the use of instruments

無磁屏蔽桶提供實(shí)驗(yàn)所需的磁屏蔽環(huán)境，它提供的磁場(chǎng)范圍為(5.57×10-7～0.890 68)Oe，可以滿足測(cè)試需求. 磁屏蔽桶的電源由信號(hào)發(fā)生器輸出電壓提供，通過精密萬用表的反饋電壓值即可計(jì)算桶內(nèi)磁場(chǎng)強(qiáng)度. 傳感器的供電電壓由5 V 的穩(wěn)壓源提供，傳感器的輸出電壓由另外一個(gè)精密萬用表檢測(cè). 其測(cè)試系統(tǒng)示意圖如圖 12 所示.

圖 12　測(cè)試系統(tǒng)示意圖Fig.12　Schematic diagram of test system

將傳感器測(cè)試電路板固定在安裝臺(tái)上，保證傳感器的敏感軸方向與桶軸方向平行. 給傳感器供電，由屏蔽桶提供磁場(chǎng)環(huán)境，記錄傳感器的輸出電壓值. 最后得到實(shí)驗(yàn)測(cè)試結(jié)果如圖 13 與表 2 所示.

圖 13　實(shí)驗(yàn)測(cè)試結(jié)果Fig.13　Experimental test results

(mV/Oe)

通過測(cè)試結(jié)果可以得到，巨磁阻傳感器的靈敏度比各向異性磁阻傳感器的靈敏度高出7倍，優(yōu)勢(shì)比較明顯，在工作電壓為5 V的情況下，巨磁阻的靈敏度為14.151 mV/Oe.

3微陀螺系統(tǒng)級(jí)仿真

3.1巨磁阻微陀螺的靈敏度仿真

上述部分完成了微陀螺結(jié)構(gòu)靈敏度、 磁場(chǎng)靈敏度、 磁阻靈敏度的設(shè)計(jì)，根據(jù)圖 3 所示，可以由式(6)算得微陀螺的靈敏度為

36nm/((°)/s)×0.09Oe/nm×14.151mV/Oe=

下面用MATLAB軟件中的SIMULINK模塊對(duì)巨磁阻微陀螺進(jìn)行系統(tǒng)級(jí)仿真，驗(yàn)證角速率信號(hào)向電學(xué)信號(hào)的轉(zhuǎn)化過程. 微陀螺系統(tǒng)級(jí)仿真模型如圖 14 所示.

圖 14　巨磁阻微陀螺系統(tǒng)級(jí)仿真模型Fig.14　The system level simulation model of GMR gyro

圖 14 所示模型為巨磁阻微陀螺的整體框架，在驅(qū)動(dòng)方向施加正弦驅(qū)動(dòng)力時(shí)，微陀螺在驅(qū)動(dòng)方向做簡(jiǎn)諧振動(dòng)，當(dāng)有角速率信號(hào)輸入時(shí)，在柯氏力作用下微陀螺在檢測(cè)方向做簡(jiǎn)諧振動(dòng)，振動(dòng)的幅值與輸入角速率信號(hào)成正比. 系統(tǒng)級(jí)仿真得到的結(jié)果如圖 15 所示.

如圖 15(a)所示為輸入的正弦角速率信號(hào)，以驅(qū)動(dòng)方向位移信號(hào)為調(diào)制信號(hào)即高頻載波信號(hào)，得到載波后的信號(hào)如圖15(b)所示； 同樣以驅(qū)動(dòng)方向位移信號(hào)作為解調(diào)信號(hào)，解調(diào)后的曲線如圖15(c)，此時(shí)信號(hào)中仍?shī)A雜有高頻信號(hào)，采用低通濾波器濾除高頻信號(hào)，還原原始信號(hào)，得到微陀螺的電壓靈敏度輸出，如圖15(d)所示.

通過系統(tǒng)級(jí)仿真對(duì)微陀螺的設(shè)計(jì)參數(shù)進(jìn)行分析，驗(yàn)證了角速率信號(hào)到電學(xué)信號(hào)的轉(zhuǎn)換過程. 仿真驗(yàn)證微陀螺的靈敏度約為46 mV·((°)/s)-1，與理論值相對(duì)應(yīng).

3.2巨磁阻微陀螺的輸出噪聲分析

根據(jù)文獻(xiàn)[9-10]中對(duì)微陀螺輸出噪聲的求解方法，得到計(jì)算微機(jī)械陀螺輸出噪聲的計(jì)算公式為

(10)

4結(jié)論

參考文獻(xiàn):

[1]曹慧亮, 李宏生, 王壽榮, 等. MEMS陀螺儀結(jié)構(gòu)模型及系統(tǒng)仿真[J]. 中國(guó)慣性技術(shù)學(xué)報(bào)，2013，21(4)： 524-529.

Cao Huiliang, Li Hongsheng, Wang Shourong, et al. Structure model and system simulation of MEMS gyroscope[J]. Journal of Chinese Inertial Technology, 2013, 21(4)： 524-529. (in Chinese)

[2]吳學(xué)忠，劉宗林，李圣怡. 微硅陀螺性能影響因素及其對(duì)策研究[J]. 微納電子技術(shù)，2003，40(7)： 274-276.

Wu Xuezhong, Liu Zonglin, Li Shengyi. Influence factors and enhancement on property of micro-silicon gyroscope[J]. Micronanoelectronic technology, 2003, 40(7)： 274-276. (in Chinese)

[3]李小卿, 常洪龍, 焦文龍, 等. 基于杠桿的微機(jī)械陀螺結(jié)構(gòu)設(shè)計(jì)與仿真[J]. 微納電子技術(shù)，2012，49(5)： 313-317.

Li Xiaoqing, Chong Honglong, Jiao Wenlong, et al. Structure design and simulation of the MEMS gyroscope based on leverage mechanism[J]. Micronanoelectronic technology, 2012, 49(5)： 313-317. (in Chinese)

[4]Jianbing Xie, Qiang Shen Yongcun Hao, et al. Design, fabrication and characterization of a low?noise Z?axis micromachined gyroscope[J]. Microsyst Technol, 2015, 21(2)： 625-630.

[5]陳志勇, 張嶸, 周斌. 一種低機(jī)械噪聲弱耦合的微機(jī)電陀螺敏感結(jié)構(gòu)[J]. 中國(guó)慣性技術(shù)學(xué)報(bào)，2015，23(3)： 373-378.

Chen Zhiyong, Zhang Rong, Zhoubin. MEMS gyroscope structure with low mechanical-thermal noise and weak mechanical coupling[J]. Journal of Chinese Inertial Technology, 2015, 23(3)： 373-378. (in Chinese)

[6]石云波, 張瓊, 羅源源, 等. 基于隧道效應(yīng)的MEMS陀螺儀的研究[J]. 壓電與聲光，2010，32(5)： 750-753.

Shi Yunbo, Zhang Qiong, Luo Yuanyuan, et al. Research of MEMS gyroscope based on tunneling effect[J]. Piezoelectrics and acoustooptics. 2010, 32(5)： 750-753. (in Chinese)

[7]茍曉凡, 楊勇, 鄭曉靜. 矩形永磁體磁場(chǎng)分布的解析表達(dá)式[J]. 應(yīng)用數(shù)學(xué)和力學(xué)，2004，25(3)： 271-278.

Gou Xiaofan, Yang Yong, Zheng Xiaojing. Analytic expression of magnetic field distribution of rectangular permanent magnets[J]. Applied mathematics and mechanics, 2004, 25(3)： 271-278. (in Chinese)

[8]劉宏娟. 矩形永磁體三維磁場(chǎng)空間分布研究[D]. 北京： 北京工業(yè)大學(xué)，2006.

[9]Seltzer S J. Developments in Alkali-metalAtomic Magnetometry. Ph. D. Thesis, Princeton University, Princeton, NJ, USA, 2008.

[10]蔣慶華，苑偉政，李曉瑩. 電容式z軸微機(jī)械陀螺的噪聲抑制[J]. 微納電子技術(shù)，2009，46(10)： 616-620.

Jiang Qinghua, Yuan Weizheng, Li Xiaoying. Noise reduction of capacitivez-Axis microgyroscopes[J]. Micronanoelectronic technology, 2009, 46(10)： 616-620. (in Chinese)

[11]Pannetier-Lecoeur M, Fermon C, de Vismes A, et al.Low noise magnetoresistive sensors for current measurement and compasses[J]. Journal of Magnetism and Magnetic Materials, 2007, 31(6)： 246-248.

Design and Noise Analysis of a MEMS Gyroscope Based on Giant Magnetoresistance Effect Detection

CHU Wei-hang1,2,3, BAI Xiao-xiao1,2,3, JIANG Xiao-yong1,2,3, LI Meng-wei1,2,3

(1. Science and Technology on Electronic Test & Measurement Laboratory,North University of China, Taiyuan 030051, China;2. School of Instrument and Electronics, North University of China, Taiyuan 030051, China;3. Center for Microsystem Intergration, North University of China, Taiyuan 030051, China)

Abstract:MEMS gyroscope based on giant magnetoresistance (GMR) effect has the advantages of high sensitivity and low noise. In this paper, the working principle of the micro gyroscope with GMR effect is described in detail, and the micro gyroscope structure is designed. The frequency of driving direction is 5 553.7 Hz, and the frequency of detection direction is 5 556.7 Hz, so the frequency difference is only 3 Hz, this design having good matching mode. The detection magnetic field applied in the micro gyro is designed and the sensitivity of GMR sensor has been tested. Based on the above results the system level simulation of the GMR micro gyroscope has been completed. Finally, the noise composition and noise calculation of the GMR gyroscope are analyzed. The simulation results show that the sensitivity of the GMR gyroscope can achieve 46 mV·((°)/s)-1, the noise level of the GMR gyroscope can achieve

Key words:giant magnetoresistance (GMR); MEMS gyroscope; system level simulation; noise analysis

文章編號(hào):1673-3193(2016)02-0142-08

*收稿日期:2015-09-05

作者簡(jiǎn)介：褚偉航(1989-)，女，碩士生，主要從事微納器件及系統(tǒng)方面的研究.通信作者： 李孟委(1975-)，男，副教授，博士，主要從事開展新原理MEMSN慣性傳感器及導(dǎo)航研究.

中圖分類號(hào):TP212

文獻(xiàn)標(biāo)識(shí)碼:A

doi:10.3969/j.issn.1673-3193.2016.02.009