D2D通信中一种资源分配与功率控制结合的方案*

王 蓉,江 帆,许腾驹,孙长印

(西安邮电大学通信与信息工程学院，西安　710121)

D2D通信中一种资源分配与功率控制结合的方案*

王 蓉**,江 帆,许腾驹,孙长印

(西安邮电大学通信与信息工程学院，西安710121)

以under1ay方式工作的D2D(Device-to-Device)通信通过资源共享复用蜂窝网络中的资源,在提高系统资源利用率的同时,对已有蜂窝链路带来了同频干扰。为了减小因无线资源复用带来的干扰,提出了一种资源分配与功率控制相结合的方法。以最小化系统干扰为目的,通过在D2D链路之间合理地分配资源,保证了蜂窝链路的通信质量;同时动态地调整D2D链路的发射功率,在保证相应蜂窝链路干扰可控的情况下,合理提高了D2D链路的通信质量。仿真结果表明,与现有方案相比,所提出的算法能够有效地提高D2D通信与蜂窝用户共存场景下的系统性能。

D2D通信;资源复用干扰;资源分配;功率控制

引用格式:王蓉，江帆，许腾驹，等.D2D通信中一种资源分配与功率控制结合的方案[J].电讯技术，2016，56(3)：295-301.[WANG Rong，JIANG Fan，XU Tengju，et a1.A joint method of resource a11ocation and Power contro1 for device-to-device(D2D)communication[J].Te1ecommunication Engineering，2016，56(3)：295-301.]

1　引 言

随着移动多媒体业务及应用的快速发展，蜂窝网络需要处理的数据流量急剧地上升，直连通信[1]通过合理利用现有可用频谱资源，成为提高系统性

能的一种重要技术，而D2D(Device-to-Device)通信技术[2]则是目前数据直通技术研究的重点。D2D通信是指在LTE系统控制下，物理距离较近的两终端之间，不经过基站(Evo1ved NodeB，eNB)中继而直接进行数据传输的技术，从而减轻eNB负载，并且减小终端与eNB之间的路径损耗，减少端到端的时延，尤其是使小区边缘用户通信质量得到很大改善，提高小区覆盖率。

目前D2D的工作方式主要分为over1ay和under1ay两种[3]，前者指D2D链路使用与原有蜂窝链路频谱资源正交的预留资源，链路间无干扰，但由于等待接入的D2D用户数目的不确定性，系统的频谱资源利用率较低；后者指D2D链路复用原有蜂窝链路的频谱资源，相对于over1ay方式，under1ay的优点是提高了频谱利用率，但信令开销较大[4]。为了提升有限频谱资源的利用率，目前的研究点主要集中在后者[5]。不同于目前已有的WLAN、Wi-Fi、Ad-Hoc等直连通信技术，以under-1ay方式工作的D2D通信技术是受eNB控制的，使用授权频带传输数据，因此可靠性、安全性更高；与传统的中继频谱资源复用相比，只需一些简单控制信令，不会对eNB产生大的负荷[6]，端到端时延低，频谱利用率更高[7]。本文将经过eNB进行中继通信的常规用户称为蜂窝用户，记作CUE(Ce11u1ar UE)，采用D2D模式通信的用户称为D2D用户对，记作DUE(D2D UE)。

为了减小由于D2D在under1ay方式下复用CUE的资源造成的同频复用干扰，目前已有很多相关研究，其中资源分配、功率控制是研究的重点[7-8]。文献[9]以空口为背景，提出了一种基于业务质量的D2D资源分配方案；文献[10]提出了一种只考虑个体干扰的资源分配方案，但忽视了系统整体受到的干扰；文献[11]提出在多层小区系统中通过高优先级用户发送命令进行功率控制的方案；文献[12]提出动态调整D2D发射功率的方案，但没有区分复用资源的上下行模式，影响了功率控制的效果。基于以上文献[10，12]存在的一些问题，并考虑到D2D通信以under1ay方式与蜂窝用户共存场景下的资源分配及干扰管理所面临的挑战，本文首先以最小化系统干扰为目的，通过为D2D链路合理地分配复用资源，保证蜂窝链路的通信质量；同时适当地调整D2D链路的发射功率，在保证蜂窝链路干扰可控的情况下，使D2D链路的通信质量得到合理提升。

2　系统模型及干扰分析

2.1系统模型

所考虑的系统由7个正交正六边形小区组成，采用正交频分复用技术(Orthogona1 Frequency Division Mu1tiP1exing，OFDM)，以保证同小区内蜂窝链路之间不存在干扰。假定蜂窝网络是采用时分双工(Time Division DuP1exing，TDD)模式，整个系统是同步的，所有eNBs都是同时使用上行链路资源或者同时使用下行链路资源，以保证上下行传输之间不存在干扰。

每个小区所分配的频带是相同的，包含K个物理资源块(Physica1 Resource B1ock，PRB)。为提高系统频谱利用率，D2D通信采用基于under1ay的工作方式，即D2D用户复用CUE已用的PRBs。每个小区中有N个CUE组成集合C={CUE1，CUE2，…，CUEN}，CUE位置随机分布，且N≤K，等待接入的D2D通信双方以D2D用户对的形式组成集合D，D ={DUE1，DUE2，…，DUEM}，D2D对也是随机分布的，且M≤N。另外，为了控制D2D的资源复用对CUE带来的干扰，要求每条蜂窝链路最多只能由一对D2D进行复用，并且一个D2D对最多复用一条蜂窝链路的资源。为了清楚表示，下文将D2D对的发射端、接收端分别表示为D2D-Tx、D2D-Rx。

2.2干扰分析

首先分析DUE复用蜂窝上、下行链路不同场景下的干扰情况，如图1所示。其中eNB与用户CUEk间信道的信干噪比(Signa1 to Interference P1us Noise Ratio，SINR)估算公式为

式中：Pk为eNB与CUEk间蜂窝链路的发射功率；Pi为DUEi发射端D2D-Tx的发射功率；小区间干扰信号强度Iinter，由文献[13]进行干扰测量得到；N0为接收端接收到的热噪声；xi，k表示DUEi复用蜂窝链路k的资源复用因子，取值为0或1，即

图1　D2D复用蜂窝链路资源的干扰情况Fig.1 The interference between D-1ink and C-1ink

链路发射端到接收端的信道增益用G表示，计算公式为

式中：PL、shadowing分别代表收发端之间的路径损耗、阴影衰落[12]。因此，根据公式(3)可以计算得到公式(1)中eNB与CUEk间链路的信道增益Gk、D2D-Txi与蜂窝链路接收端间的干扰信道增益Gi，k。如图1所示，以DUE复用CUE的下行链路资源为例，如果DUE复用距离较近的CUE的链路资源，即使D2D-Tx的发射功率较小，CUE也会受到较大干扰[7]。同样地，复用上行链路时，若D2D对随机复用CUE的链路资源，D2D-Rx可能会受到较大的干扰。因此如何令D2D对合理地复用已有的蜂窝链路资源，从而协调由于资源复用所带来的干扰是D2D资源分配必须解决的问题。

文献[10]研究了上行链路复用的干扰情况，由于上、下行研究方法类似，只是受干扰方有所不同，因此本文仅以复用下行链路为例来进行研究。

3　资源分配

本文以复用下行蜂窝链路资源为例，首先分析目前已有的研究方案，并在此基础上提出一种基于整体干扰最小的资源分配策略，以保证D2D用户全部接入情况下，对原有蜂窝链路造成的干扰最小；最后比较不同的资源分配策略的性能。

3.1随机资源分配(Random Assignment,RA)策略

文献[14]介绍了一种随机资源复用算法，即完全不考虑资源复用带来的干扰对双方的影响，随机地为D2D对分配CUE的频谱资源进行复用。

3.2基于个体干扰最小的资源分配(Individual Optimization,IO)策略

文献[10]提出了一种基于个体干扰最小的资源分配策略，该策略的基本思想是先将本小区内D2D用户对的收发间信道增益进行估计，记为Gi(i =1，2，…，M)，优先为链路质量好的D2D对选择复用的资源，并且选择受到此D2D对干扰最小的CUE的资源进行复用，然后为链路质量次优的D2D对选择合适的复用资源，直到为所有等待接入的D2D对都分配到资源为止。因此，该方案可以表示为个体干扰最小的资源分配算法，即IO。

设Ii，k表示DUEi复用CUEk的PRB时，CUEk受到的D2D-Txi的干扰，xi，k表示D2Di对于CUEk链路资源的复用情况，IO算法如下：

(1)初始化等待接入的DUEs集合D'=D，处于通信状态的CUEs集合C'=C；

(2)for i=1：M(DUE用户对的序号1～M)

得到的所有xi，k(i=1，2…，M，k=1，2…，N)组成X矩阵对应相应的资源分配情况。

3.3启发式算法:基于整体干扰最小(Overall Op

timization,OO)的资源分配策略

文献[10]提出的方案在为D2D对分配资源时，只考虑了此D2D对造成的干扰最小化，而不是从整体的角度进行DUEs与CUEs之间的配对。参考文献[10]，本文提出了基于整体干扰最小的资源分配策略(OO)，该策略的基本思想是最小化整体蜂窝链路所受干扰，具体实现思路如下：

(1)首先遍历得到每个D2D对若复用不同CUE的链路资源时造成的所有干扰值Ii，k(i=1，2…，M；k=1，2…，N)，组成干扰矩阵IM×N作为效率矩阵；

(2)所有D2D对接入时，以蜂窝链路受到的总干扰最小为目标，将干扰最小化问题转化为对资源复用因子xi，k(i=1，2，…，M；k=1，2，…，N)进行0-1最优指派问题[15-16]，在所有配对组合中找到最合适的一种作为最佳分配策略；

(3)以蜂窝系统受到的总干扰最小为目标，根据得到的干扰矩阵IM×N，利用匈牙利算法[16]，对干扰矩阵变换得到各行、各列均有零元素的矩阵，然后判断并得到M个独立零元素(在干扰矩阵中，一组处在不同行、不同列的零元素称为独立零元素组，其中的每个元素称为独立零元素)，从而得到最优指派方案，即DUE的资源复用的决策矩阵XM×N。

具体算法描述如下：

(1)初始化等待接入的DUEs对数m=M，处于通信状态的CUEs数目n=N；

(2)遍历计算所有情况下CUEs受到的D2D用户对的同频干扰Ii，k：

(3)用Hungarian算法进行资源分配：

4　功率控制(Power Control,PC)

包含D2D通信在内的蜂窝系统中，由于各类用户及eNB的发射功率对于系统链路之间的干扰较大，需要考虑如何在DUE以及CUE和eNB之间合理地分配功率。因此，为D2D对分配频谱资源后，需要对D2D-Tx的发射功率Pd进行调整，以提高D2D对的通信质量，同时避免对相应链路造成太大干扰。

文献[12]提出了一种动态功率控制算法，通过周期性地对每对D2D-Tx的发射功率Pd进行调整，以保证D2D链路对相应CUEs的干扰可控。其中，eNB调整每个Pd的标准就是保证相关的eNB，以及本小区和邻小区内相关CUEs受到的该D2D-Tx的干扰值都在可忍受范围内，但是该方法没有区分D2D复用的资源是蜂窝上行链路资源还是下行的。本文采用的功率控制方法是基于文献[12]所提方法的改进，对D2D通信复用的链路资源加以区分，这样对Pd的限制条件也就减少了。

同样以复用下行链路资源为例，由图2可知，D2D-Tx对本小区CUE造成干扰，同时对邻小区也造成干扰，需要将频率复用所引起的干扰值限定在CUE可忍受的范围内[12]。

图2　D2D功率控制的干扰限制条件Fig.2 The restrictions of the transmit Power of D2D-Tx Pd

假设第i对D2D复用本小区CUEk的下行链路资源，则要保证本小区CUEk接收到D2D-Txi的干扰功率小于预先设定的干扰门限值；同时由于D2D -Txi对邻小区相应蜂窝链路也可能造成干扰，因此，选择邻小区有可能接收到该D2D-Tx干扰功率最大的位置(该边缘位置记作edge)作为参考(在实际应用场景中，并不需要得到edge的具体位置，而只需要得到D2D-Tx到邻小区边缘位置的最小距离即可)，需要满足该位置接收到的D2D-Tx的干扰功率值也小于设定的门限值。设D2D-Txi的发射功率记为Pdi，首先由公式(3)得到D2D-Tx到两个受干扰位置的信道增益分别为Gi，k、Gi，edge，则本小区CUEk及边界edge处接收到的该D2D对的干扰功率需要满足

式中：γth为设定的门限值，该值表示接收端可以忍受的最大干扰功率。为使相应蜂窝链路的通信质量得到保证，D2D-Txi的发射功率需要同时满足

因此，在此取D2D-Txi的发射功率为

5　性能仿真及分析

系统模型为7小区模型，仿真参数选用3GPP LTE所规定的OFDMA系统的参数[17]，如表1所示。为了作比较，对采用不同的频率资源复用和功率控制方案时，系统中用户的SINR和吞吐量分别进行仿真。

表1　仿真参数设定Tab.1 Simu1ation Parameters

在仿真中，before表示蜂窝系统中未引入D2D通信的情况，其中OONPC、IONPC、RANPC分别表示不考虑功率控制时整体干扰最小资源分配方案、个体干扰最小资源分配方案以及随机资源分配方案；OOPC、IOPC、RAPC则分别表示以上3种不同的资源分配方案与改进的功率控制结合的方案；而PC为Gu等人提出的功率控制方案[12]，该方案并没有考虑资源分配的情况。

5.1不同资源分配策略性能比较

考虑无功率控制情况下系统传输的性能。资源分配方式采用第3节中的3种不同的资源分配方法，D2D-Tx以恒定的发射功率进行通信，蜂窝链路的信干噪比分布情况如图3所示。从图3中可以看出，因D2D通信的引入对蜂窝链路造成同频干扰，使蜂窝用户的链路质量变差，并且OONPC的蜂窝用户SINR分布大概比IONPC高5 dB，比RANPC高15 dB，这是由于OONPC算法找到了使整个系统中CUEs受到的干扰最小的最佳配对策略，而IONPC算法只考虑个体，不考虑其他D2D对造成的干扰，CUE整体所受干扰相对较大，对于RANPC算法由于采用随机复用策略，CUEs通信性能更差。

图3　无功率控制时不同资源分配方案CUEs性能比较Fig.3 SINR of C-1inks using different resource a11ocation methods without Power contro1

5.2功率分配与不同资源分配策略结合

考虑功率控制之后，采用3种资源分配算法的蜂窝用户SINR分布情况仿真结果如图4所示。从图4中可以看出，与图3相比，加功率控制后OOPC和IOPC的分布都变差了一些，这是由于系统为DUE分配了较远的CUE，DUE对CUE产生的干扰较小，会适当地增大Pd，因此CUE的SINR有所降低。而RAPC因为给DUE分配的CUE可能距离比较近，DUE造成干扰较大，进行功率控制会适当降低Pd，因而CUE的SINR分布会有所提高。而RAPC与PC几乎重合，原因是虽然RAPC对Pd的限制少了，但由于经过路径损耗，Pd的提高对相应CUE的影响较小。总体来看，采用OOPC后，CUE的SINR分布更好，虽然相比OONPC差一些，但是影响较小。

图4　资源分配与功率控制结合后CUEs性能比较Fig.4 SINR of C-1inks with different strategies

图5给出了不考虑功率控制的条件下采用OONPC资源分配方法为D2D传输分配频谱资源时，D2D链路的信干噪比分布情况。由于在未考虑功率控制之前，D2D复用下行蜂窝链路资源，所以受到的干扰来自eNB，与所分配资源无关。因此当D2D-Tx端采用恒定发射功率时，3种资源分配算法D2D的SINR分布是相同的，在此画出了对应曲线OONPC。

图5　采用不同方案DUEs的性能比较Fig.5 SINR of D-1inks with different strategies

而结合所提出的功控算法之后的3种方法，DUE的SINR都变得更好，而且，OOPC要比IOPC的分布高3～4 dB，要比RAPC高3 dB，这是由于采用OOPC策略时DUE对CUE造成的干扰较小，调整Pd值后，高于其他两种算法，直接影响了DUE的SINR分布。Gu[12]等提出的功率控制方案PC，由于对D2D发射功率限制较多，Pd值较小，导致DUE的SINR分布较差。总体来看，采用OOPC之后，DUE的SINR分布最好。

5.3系统吞吐量

通过对系统的吞吐量进行仿真，结果如图6所示。可以看出采用3种不同资源分配方案在加入功率控制后，系统吞吐量与之前相比都有提高。这是因为加入功率控制后，虽然对蜂窝链路造成一些影响，但D2D链路质量较好，整体吞吐量得到提高。而且采用了所提出的OOPC方案后，系统吞吐量分布最优。因此，可以得出引入D2D通信之后，通过采用合理的资源分配与功率控制结合的策略，可以有效提高系统吞吐量。

图6　采用不同方案时系统吞吐量分布情况Fig.6 System throughPut with different methods

6　结束语

D2D通信的引入显著地减轻了eNB的负载压力，并且工作在under1ay方式的D2D通信提高了系统的频谱利用率，但由于复用原有蜂窝链路的资源又引入了新的干扰。为减小此复用干扰，本文提出了一种将资源分配与功率控制结合的算法。通过性能仿真，验证了该算法通过基于整体干扰最小的资源分配方法在一定程度上协调了D2D引入的干扰，并结合所提出的功率控制方法有效地提高了D2D用户的通信性能，同时不会对蜂窝链路带来较大的影响，而且提高了系统的吞吐量。所提出的算法可以有效解决D2D与蜂窝用户共存场景下的资源分配和干扰控制问题。与已有算法相比，所提出的算法能更有效地提高系统性能。但目前的研究仅考虑了D2D对分布较少的场景，在今后的工作中会考虑D2D用户大规模分布的场景，进一步降低算法的复杂度。

[1] LEE N，LIN X Q，ANDREWS J G，et a1.Power contro1 for D2D under1aid ce11u1ar networks：mode1ing，a1gorithms，and ana1ysis[J].IEEE Journa1 on Se1ected Areas in Communications，2015，33(1)：1-13.

[2] MARCO B，GABOR，DEMIAD P，et a1.Benchmarking Practica1 RRM a1gorithms for D2D communications in LTE advanced[J].Wire1ess Persona1 Communications，2015，82(2)：883-910.

[3] HAKOLA S，TAO C，LEHTOMAKI J，et a1.Device-todevice(D2D)communication in ce11u1ar network-Performance ana1ysis of oPtimum and Practica1 communication mode se1ection[C]//Proceedings of 2010 IEEE Wire1ess CommunicationsandNetworkingConference (WCNC).Sydney，Austra1ia：IEEE，2010：1-6.

[4] CAO Y，JIANG T，WANG C G.CooPerative device-todevice communications in ce11u1ar networks[J].IEEE Wire1ess Communications，2015，22(3)：124-129.

[5] TEHRANI M N，UYSAL M，YANIKOMEROGLU H.Device-to-device communication in 5G ce11u1ar networks：cha11enges，so1utions，and future directions[J].IEEE Communications Magazine，2014，52(5)：86-92.

[6] HUAN T，ZHI D，LEVY B C.Enab1ing D2D communications through neighbor discovery in LTE ce11u1ar networks [J].IEEE Transactions on Signa1 Processing，2014，62 (19)：5157-5170.

[7] MALANDRINO F，LIMANI Z，CASETTI C.Interferenceaware down1ink and uP1ink resource a11ocation in HetNets with D2D suPPort[J].IEEE Transactions on Wire1ess Communications，2015，14(5)：2729-2741.

[8] FENG D Q，LU L，YI Y W，et a1.Device-to-device communications in ce11u1ar networks[J].IEEE Communications Magazine，2014，52(4)：49-55.

[9] ZHU X Y，WEN S，CAO G，et a1.QoS-based resource a1-1ocation scheme for device-to-device(D2D)radio under-1aying ce11u1ar networks[C]//Proceedings of 19th Internationa1 Conference on Te1ecommunications.Jounieh：IEEE，2012：1-6.

[10] 韩涛.LTE终端直通通信系统的资源复用与干扰协调[D].杭州：浙江大学，2013. HAN Tao.Resource reusing and interference coordination of D2D under1aying LTE[D].Hangzhou：Zhejiang University，2013.(in Chinese)

[11] HOSSAIN E，RASTI M，TABASSUM H，et a1.Evo1ution towards 5G mu1ti-tier ce11u1ar wire1ess networks：an interference management PersPective[J].IEEE Wire1ess Communications Magazine，2014，21(3)：118-127.

[12] GU J，BAE S J，CHOI B G.Dynamic Power contro1 mechanism for interference coordination of device-to-device communication in ce11u1ar networks[C]//Proceedings of Third Internationa1 Conference on Ubiquitous and Future Networks(ICUFN).Da1ian：IEEE，2011：71-75.

[13] Nokia CorPoration.Method and aPParatus for Providing interference measurement for D2D communication：US，A1，2010/0093364[P].2010.

[14] SUN H G，SHENG M，WANG X J，et a1.Resource a11ocation for maximizing the device-to-device communications under1aying LTE-Advanced networks[C]//Proceedings of 2013 IEEE Internationa1 Conference on Communications.Xi'an：IEEE，2013：60-64.

[15] 宋巨龙，王香柯，冯晓慧.最优化方法[M].西安：西安电子科技大学出版社，2012：184-192.

SONG Ju1ong，WANG Xiangke，FENG Xiaohui.The oPtimization method[M].Xi'an：Xidian University Press，2012：184-192.(in Chinese)

[16] 常庭懋，韩中庚.用“匈牙利算法”求解一类最优化问题[J].信息工程大学学报，2004，5(1)：60-62.

CHANG Tingmao，HAN Zhonggeng.Using the Hungarian a1gorithm to so1ve a c1ass of oPtimization Prob1ems[J]. Journa1 of Information Engineering University，2004，5 (1)：60-62.(in Chinese)

[17] 3GPP TR 25.814V7.1.0，Physica1 1ayer asPects for evo1ved UTRA[S].

王 蓉(1990—)，女，河北衡水人，2013年于石家庄学院获工学学士学位，现为硕士研究生，主要研究方向为宽带无线通信、D2D通信的干扰协调技术；

WANG Rong was born in Hengshui，Hebei Province，in 1990.She received the B.S.degree from Shijiazhuang University in 2013.She is now a graduate student.Her research concerns wideband wire-1ess communication，interference coordination techno1ogy of device-to-device communications under1aying ce11u1ar network.

Emai1：18729224395@163.com

江 帆(1982—)，女，2010年于北京邮电大学获工学博士学位，现为副教授、硕士生导师，主要研究方向为超密异构网络、D2D通信、协作中继网络等；

JIANG Fan was born in 1982.She received the Ph.D.degree from Beijing University of Posts and Te1ecommunications in 2010.She is now an associate Professor and a1so the instructor of graduate students.Her research concerns u1tra dense heterogeneous network architecture，D2D communications，cooPerative re-1ay network.

许腾驹(1990—)，男，2013年于西安电子科技大学获工学学士学位，现为硕士研究生，主要研究方向为未来5G系统中的D2D通信的中继接入技术；

XU Tengju was born in 1990.He received the B.S.degree from Xidian University in 2013.He is now a graduate student. His research concerns the re1ay access techno1ogy of device-todevice communications in 5G te1ecommunication network.

孙长印(1963—)，男，2000年于西安电子科技大学雷达信号处理重点实验室获博士学位，现为副教授，主要研究方向为下一代无线通信技术。

SUN Changyin was born in 1963.He received the Ph.D. degree from Xidian University in 2000.He is now an associate Professor.His research concerns interference management and radio resource management in the evo1ved mobi1e communication system.

A Joint Method of Resource Allocation and Power Control for Device-to-Device(D2D)Communication

WANG Rong，JIANG Fan，XU Tengju，SUN Changyin
(Schoo1 of Communication and Information Engineering，Xi'an University of Posts and Te1ecommunications，Xi'an 710121，China)

By sharing the resource with ce11u1ar users，device-to-device(D2D)communication under1aying ce11u1ar network can imProve the system resource uti1ization.But unfortunate1y，frequency reuse interference wi11 be introduced between D2D users and ce11u1ar users.In order to coordinate the co-channe1 interference，a combined strategy of resource a11ocation with Power contro1 is ProPosed in this PaPer.The Pro-Posed strategy a11ocates resource to D2D 1inks based on the overa11 system oPtimization for minimum interference.Furthermore，each D2D transmitter can dynamica11y adjust the transmit Power so as to imProve the qua1ities of D2D Pairs whi1e ensuring the interference to ce11u1ar users to remain in a to1erab1e range.Simu1ation ana1ysis demonstrate that comPared with the existing methods，the ProPosed method can effective1y imProve overa11 system Performance.

D2D communication；frequency reuse interference；resource a11ocation；Power contro1

The Nationa1 Natura1 Science Foundation of China(No.61501371);The Nationa1 High-tech R&D Program of China(863 Program)(2014AA01A703);The Natura1 Science Foundation of Shaanxi Province(2011JM8027);The Science and Techno1ogy Research and Deve1oPment Program of Shaanxi Province(2014K09-14)

TN929.5

A

1001-893X(2016)03-0295-07

10.3969/j.issn.1001-893x.2016.03.011

2015-07-29;

2015-11-20 Received date:2015-07-29;Revised date:2015-11-20

国家自然科学基金资助项目(61501371);国家高技术研究发展计划(863计划)项目(2014AA01A703);陕西省自然科学基金资助项目(2011JM8027);陕西省科技攻关计划项目(2014K09-14)

**通信作者:18729224395@163.com Corresponding author:18729224395@163.com