• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    HetNets Under Decoupled Uplink and Downlink Access with UE Random Discontinuous Transmission: Local Delay and Energy Efficiency

    2024-04-01 02:08:40ZhangYueZhengFuchunLuoJingjing
    China Communications 2024年3期

    Zhang Yue,Zheng Fuchun,Luo Jingjing

    School of Electronic and Information Engineering,Harbin Institute of Technology(Shenzhen),Shenzhen 518055,China

    Abstract: Due to the limited uplink capability in heterogeneous networks(HetNets),the decoupled uplink and downlink access(DUDA)mode has recently been proposed to improve the uplink performance.In this paper,the random discontinuous transmission (DTX)at user equipment(UE)is adopted to reduce the interference correlation across different time slots.By utilizing stochastic geometry,we analytically derive the mean local delay and energy efficiency(EE)of an uplink HetNet with UE random DTX scheme under the DUDA mode.These expressions are further approximated as closed forms under reasonable assumptions.Our results reveal that under the DUDA mode,there is an optimal EE with respect to mute probability under the finite local delay constraint.In addition,with the same finite mean local delay as under the coupled uplink and downlink access(CUDA)mode,the HetNets under the DUDA mode can achieve a higher EE with a lower mute probability.

    Keywords: decoupled uplink and downlink access(DUDA);heterogeneous networks(HetNets);random discontinuous transmission(DTX);stochastic geometry

    I.INTRODUCTION

    With the rapid growth of global mobile data and increased demand for a better quality of service (QoS),the heterogeneous networks (HetNets),where multiple low-cost and low-power small base stations(SBSs)are deployed over the existing macro base stations tier(MBSs),have become a key architecture for the future communications networks [1,2].Recently,the importance of uplink (UL) performance in HetNets has arisen along with the development of social networking,UAV surveillance and cloud solutions [3].However,under the conventional coupled UL and downlink (DL) access (CUDA) mode,where both the UL and DL of a user are associated with the same BS,the UL performance may be limited by the DL parameters such as different transmission powers of the BSs in different tiers.Although several techniques under CUDA,e.g.,cell range extension(CRE)[4],have been proposed to mitigate the impact of inter-tier differences,the association criterion based on the best DL performance generally does not result in the best UL performance as well.Therefore,the concept of decoupled uplink and downlink access (DUDA) has been proposed in [5],which enables users or UEs to be connected with different BSs for their UL and DL.Under the DUDA mode,without affecting the DL association,a user is allowed to be associated with the BS which receives the strongest signal in the UL,thus ensuring a high-quality UL transmission.In [6],a UL performance comparison was carried out,showing that the DUDA mode leads to better load balance for HetNets and has a better or equal average spectral efficiency(SE)and energy efficiency(EE)than the CUDA mode.Based on DUDA,the impact of fractional power control(FPC)on the delay and EE of UL HetNets was investigated in [7],but the influence of the interference correlation was not considered.

    Apart from SE and EE,local delay,defined as the average time taken until a data packet has been transmitted successfully,is a key component of transmission delay and is considered an important performance metric reflecting the QoS and network reliability [8].In [9],an approximation of the interfering user point process was used to derive the mean local delay in a single-tier uplink network.For a static low-mobility HetNet,the time correlation in the spatial distribution of the nodes will cause interference correlation between different time slots [10],resulting in an increase of local delay (i.e.,the same set of interfering users may be faced by a BS across different time slots).Particularly,the local delay may become an infinity for certain poor interferer distribution,referred to as network-contention phase transition [11].As such,[12] employed random discontinuous transmission (DTX) at BSs to reduce DL mean local delay.[13] investigated the local delay of DL HetNets with random DTX where user and BS locations are correlated.Inspired by this,we in this paper also adopt the UE random DTX scheme,by allowing UEs to be muted independently and randomly in each time slot,to weaken the interference correlation and reduce the local delay.

    On the other hand,the uplink DTX scheme,in the context of the explosive growth of mobile smart devices,is also beneficial for saving energy consumption on the UE side [14].Yet the decrease in energy consumption may not imply a rise in energy efficiency(another key performance metric in this paper),since throughput is also affected by device sleeping.In contrast with the legacy cycle-based DTX or discontinuous reception (DRX) [15],the impact of UE random DTX on uplink energy efficiency remains to be studied.

    In fact,the trade-off between the delay and EE is one of the important trade-off issues in HetNets and has been extensively investigated [16-18].The Het-Nets,especially DL HetNets with DTX at BSs,have been well analyzed.In [19],considering the random BS-DTX,the mean local delay and EE of DL HetNets based on the Poisson cluster process were analyzed.

    However,the relationship between the mean local delay and EE in the UL HetNets with UE random DTX has not been adequately studied.On the one hand,due to the usual orthogonal access mode within each BS cell,only one UE per BS transmits on a particular resource block and forms interference only to the UE in other cells.The uplink interfering point process therefore is more difficult to model than downlink interference due to the correlation between the Voronoi cells of the Poisson-Voronoi tessellation and the channel access scheme in each cell[9].On the other hand,the analysis of downlink performance does not need to consider the direct impact of uplink parameters(as in most literature),whereas the uplink performance is directly subject to the downlink parameters such as downlink power.This is another motivation for the DUDA analysis in this paper.Furthermore,the respective impact of UE DTX on DUDA-based HetNets and the performance differences between DUDA and CUDA modes,have yet to be explored.In this paper,we therefore focus on a DUDA-based UL HetNet and investigate the impact of UE random DTX on the mean local delay and EE by comparing the performance of the DUDA and CUDA modes.The main contributions can be summarized as follows:

    ? Utilizing stochastic geometry,we derive thebth moments of the conditional successful transmission probability(STP)with UE random DTX,leading to general results for the mean local delay and EE in a UL DUDA-basedK-tier HetNet with UE random DTX.Furthermore,simplified expressions of the mean local delay and EE under the same cross-tier path loss exponent are obtained for DUDA and CUDA mode,respectively.

    ? Based on a reasonable assumption on the distribution of interfering users,we further obtain some approximation results with closed forms for local delay and energy efficiency in DUDA-based and CUDA-based HetNets,respectively.Then,the range of the signal-to-interference ratio(SIR)threshold that makes the local delay a finite value and the bounds of optimal SIR threshold that maximizes the energy efficiency are obtained.In addition,the effect of the mute probability is also analyzed.

    ? With the simulation results,we reveal the effect of the key parameters (e.g.,SIR and mute probability) on local delay and EE.Under the finite mean local delay constraint,the relationship between the maximum EE and the mute probability is investigated.Finally,the performance difference between the DUDA and CRE-based CUDA is also analyzed.

    II.SYSTEM MODEL

    2.1 Heterogeneous Network Model

    We consider an uplink HetNet with DUDA comprisingKindependent tiers of BSs.DenotingK={1,...,K},for anyk ∈K,the positions of BSs in thekth tier follow a homogeneous Poisson point process (PPP) Φkwith intensityλk.Moreover,the BSs across tiers differ in terms of transmit powerPk.We assume that each BS schedules only one user on each resource block(RB)at a given time,i.e.,orthogonal access under the same BS is assumed.The user point process is also modeled as a homogeneous PPP Φuwith densityλuand independent across{Φk}k∈K.Note that the user density is assumed to be sufficiently high such that there is always at least one interfering user served by another BS.All users have the same transmit powerPu.Furthermore,the standard path loss model with the path loss exponents{αk}k∈K >2 for tierkis used.Letht,xdenote the power fading coefficient between the user atxand the tpical user’s serving BS(termed tagged BS)at the origino ∈R2in time slott.We assume that{ht,x}are spatially independent over the different links,temporally independent across time slots,and follow the exponential distribution with unit mean (i.e.,the Rayleigh fading model).In addition,since the interfering user density is typically high in an ultra-dense HetNet,we ignore the noise and focus on the interference-limited regime.

    In this paper,we focus on the uplink transmission performance of the DUDA mode versus the CUDA mode.A simple case of transmissions in the DUDA and CUDA modes is depicted in Figure 1.In the conventional CUDA mode based on cell range extension(CRE),each user is connected to the same BS for both UL and DL transmissions,i.e.,the BS which offers the maximum long-term averaged biased received-power(BRP)in DL transmission.In the DUDA mode,however,the UL association can be different from the DL association.Each user accesses the BS which receives the strongest long-term averaged signal (transmitted by this user)in UL transmission,and therefore obtains the strongest uplink signal strength.DenoteDkas the association distance between a UE and its serving BS.According to[20,Lemma 3],when the biased transmit power ratio is equal to 1,the probability density function(pdf)ofDkin the DUDA mode can be expressed as

    Figure 1.Illustration of CUDA and DUDA in a two-tier HetNet.

    whereAkis the probability that a UE is associated with thekth tier BS in the DUDA mode [20],which is given by

    2.2 UE Random Discontinuous Transmission Scheme

    We assume that the time is divided into equal slots and each transmission attempt requires one time slot.In the random DTX scheme,each UE is in either of the two transmission states for each time slot,i.e.,the mute state and the active state.We model the transmission state in each time slot as a Bernoulli trial with a tunable parameterζ(0≤ζ <1)called mute probability.Namely,a UE in each time slot stays mute temporally independently with probabilityζand stays active with probability 1-ζ.

    Denoting ΦI,j ?ΦuandIt,j,respectively,as the set of interfering users associated with thejth tier and their aggregate interference at the tagged BS in the time slott,then the overall interference at the tagged BS in the time slottis

    where1(·)is the indicator function,Φu,tdenotes the set of active UEs in the time slott,andxj,idenotes the location of theith interfering user connected within tierj.The received SIR at the tagged BS in thekth tier and time slottis therefore

    wherexk,0is the location of the typical user andRkdenotes the distance between the typical user under consideration and its serving BS at the tierk(i.e.,the tagged BS).

    2.3 Performance Metrics

    In this paper,we consider two performance metrics:the mean local delay and energy efficiency.

    Mean Local Delay:We assume that retransmission will occur if a transmission fails in a certain time slot.The mean local delay is the mean number of transmission attempts(or time slots)needed before the first success if the transmitter is allowed to keep transmitting.

    We consider the event that the received uplink SIR at the tagged BS in thekth tier is higher than a specific valueθ(i.e.the SINR threshold)in a time slot as a successful transmission event,which is denoted byCk|ΦIand conditioned on the distanceRkand the set of interfering users ΦI=∪j∈KΦI,j.The corresponding probability averaged over fading is termed conditional successful transmission probability(STP)and is given by

    Given the point process,the needed number of time slotsLk|ΦIuntil a successful transmission is geometrically distributed with Pr(Ck|ΦI),thus the conditional mean ofLk|ΦIis Pr-1(Ck|ΦI).Therefore,the mean local delay in thekth tier can be obtained by the expectation with respect toRkand ΦI,i.e.,

    From the law of total probability,the mean local delay of a HetNet can be expressed as

    Energy Efficiency:We define the uplink energy efficiency(with the unit of nats/J/Hz)as the ratio of uplink area spectrum efficiency to average area user power consumption.Assuming that a fixed-rate transmission scheme is applied,i.e.,data packets are transmitted at the rate of log(1+θ)nat/sec/Hz in each time slot,the uplink spectrum efficiency can then be expressed by

    Furthermore,the total power consumption of a UE in this paper is modeled as

    wherePcandPsdenote the circuit power expenditure in the active state and mute state,ρis the UE power amplifier (PA) efficiency.Therefore,the energy efficiency can be calculated by

    From (6) and (9),it can be observed that the mean local delay and the STP of a tier respectively are the-1-st moment and the 1-st moment of the conditional STP.In the next section,we will derive theb-th moment of the conditional STP,thus obtaining the expressions for the mean local delay and EE of aK-tier Het-Net with UE random DTX.

    III.MEAN LOCAL DELAY AND ENERGY EFFICIENCY

    Here,under the UE random DTX,we first analyze theb-th moment of the conditional STP in thekth tier(i.e.,then derive the general result of the mean local delay and EE in a HetNet under DUDA and further give some simpler expressions for a special case.

    3.1 General Case and Main Results

    Theorem 1.In a K-tier HetNet in DUDA mode with a UL random DTX scheme,the b-th moment of the conditional STP in the k-th tier is given by

    Proof.The conditional STP in thekth tier can be evaluated as

    where(a)and(b)follow the fact thatandare exponentially distributed with unit mean.Due to the uplink orthogonal access within each cell,only one UE per cell transmits on a resource block and there is no intra-BS or intra-cell interference.Therefore,ΦI,jis a non-stationary process with the intensity depending on the interfering distance[9],and the exact statistics of the interference are hard to derive.In this paper,we adopt the following approximation for the intensity function of ΦI,jaccording to[9],which is

    where(a)is obtained by utilizing the probability generating functional(PGFL)of the general PPP[21]andfRk(r) follows (1).Then by employing the substitutiont=x2/r2we can obtain(12).

    3.2 Special Case and Results

    The expression of the mean local delay and EE can be simplified when the links in each tier are assumed to have the same path loss exponent,i.e.,{αk}k∈K=α.Note that under this assumption a UE will access the geographically closest BS in UL transmission when the DUDA mode is adopted.Then we have the results in the following corollary.

    Corollary 1.When {αk}k∈K=α,the uplink mean local delay and EE of a K-tier HetNet in DUDA modewith UE random DTX are given by

    In addition,based on a similar analytical procedure,we can obtain the same analysis to obtain the mean local delay and energy efficiency of the CRE-based Het-Net with CUDA according to [20,Lemma 3],which are given by(18)and(19).

    Bjdenotes the bias factor of thejth tier,and=Bj/Bk,=Pj/Pk.

    It can be seen from (18) and (19) that the performance of CUDA mode under=1 will be equivalent to that of DUDA mode.This is because,under such a condition,the BS selection in the UL will not be decided based on the difference of transmission power of BSs in different tiers-the BS is chosen according to the closest distance in the UL transmission,which ensures the strongest uplink signal strength.

    3.3 Approximate Analytical Expressions

    Since the equations(16)-(19)still contain complex integration terms,it is hard to directly obtain useful conclusions about the effects of network parameters on local delay and energy efficiency.Therefore,in this part,based on a reasonable assumption,we are to obtain some approximate closed-form equations of the analytical results and give some insightful propositions based on them.

    For further simplification,we approximate the interfering users as following a homogeneous PPP distribution.Since the users further away from the tagged base station receive little influence from the tagged BS when choosing cells,their locations can be regarded as independent from the tagged BS.Therefore,it is reasonable to consider the interfering users as a homogeneous PPP distribution outside a certain range centered around the tagged BS.Moreover,in the area close to the tagged BS,the assumption may exaggerate the density of interfering users due to orthogonal intra-cell access,thus the results obtained based on this assumption will be slightly more pessimistic than the real situation (hence they can serve as a bound).As such,the interference density of users in thej-th tier is assumed to beλj.We then have the following conclusions.

    Corollary 2.The mean local delay of DUDA-based and CRE-based CUDA HetNets can be approximated,respectively,as

    Proof.For the DUDA mode,based on our assumption,the local delay in(16)can be rewritten as

    It is worth noting that since(BjPj)2/α >1,we can conclude thatDC≥DD,where the equality holds taken when=1 forj ∈K.

    The analytical results of (16) and (20) are given in Figure 2,and it can be seen that they have similar trends,with(20)performing slightly worse than(16),which confirms our previous analysis.It can also be seen that they converge to infinity at a certain SIR threshold,which we term the critical thresholdθc:

    Figure 2.The analytical results of local delay for α=4 dB,ζ=0.7.

    Determining the value of the critical thresholdθcis essential to allow the network to achieve a finite mean local delay,thus next we give the upper bound(i.e.,θc)on the value ofθthat enables local delay to maintain a finite value in DUDA and CUDA modes,respectively.

    Proposition 1.In DUDA mode and CRE-based CUDA mode,the mean local delay is finite when the SIR thresholds θ satisfy

    for CRE-based CUDA,respectively.

    Sinceα >2,θcis an increasing function with respect toζ,illustrating thatθccan be increased by setting a larger mute probability,which leads to a lower mean local delay in the severe interference scenarios.

    Utilizing a similar analytical approach,we can also obtain the approximations for EE.

    Corollary 3.The energy efficiency of DUDA-based and CRE-based CUDA HetNets can be approximated,respectively,as

    Similarly,it can be observed that≥since (BjPj)2/α >1.It also can be seen that both log(1+θ) andψ(θ,α) are increasing functions with respect toθ,and log(1+θ) has a larger slope whenθis small and gradually decreases and eventually falls below the slope ofψ.Therefore,there should be an optimal valueθ*for energy efficiency,i.e.,the solution for

    For the DUDA-based HetNets,we give the approximate range of the value ofθ*below.

    Proposition 2.In DUDA-based HetNets,the bounds of the optimal SIR threshold θ*that maximizes energy efficiency are given as

    where γ=α2sin(2π/α)(2π(1-ζ))-1.

    Proof.The derivative of(26)with respect toθyields

    Sinceθ/(1 +θ) ≤log(1+θ) ≤θ,we have the following result

    Then we can obtain the bounds ofθ*in(29).

    It can be seen that the upper and lower bounds ofθ*are both decreasing functions ofζ,i.e.,asζincreases,the optimal valueθ*taking interval will move positively along the coordinate axis,which is not conducive to taking a finite local delay,so there is a tradeoff between local delay and energy efficiency.

    IV.SIMULATION RESULTS

    In this section,we present some numerical and simulation results to investigate the effect of the critical network parameters (i.e.,the SIR thresholdθand mute probabilityζ) on the mean local delay and energy efficiency.Moreover,some comprehensive comparisons of the performance for the DUDA mode and CRE-based CUDA mode are also provided.Here,we consider an interference-limited two-tier HetNet composed of MBSs (Tier 1) and SBSs (Tier 2) withα1=α2=α,and the system parameters are chosen according to[6],i.e.,P1=43 dBm,P2=21 dBm,Pu=10 dBm,Pc=Ps=13 dBm,ρ=1,λ1=1/(π5002),λ2=5λ1andλu=200λ1.In addition,we perform the Monte Carlo simulations for a circular area with a radius of 1 km,and the results are obtained from an average of 10,000 independent realizations.

    4.1 Effect of SIR Threshold

    Figure 3 shows the mean local delay and EE as a function of the SIR thresholdsθwith different values of the UE mute probabilityζ.First,in Figure 3(a),the mean local delayDincreases withθ,and when the SIR threshold exceeds a critical valueθc(ζ),a phase transition occurs(i.e.,the mean local delay becomes infinite from being finite).Namely,due to the interference correlation in different time slots,consecutive unsuccessful transmission events are more likely to appear at a highθ(i.e.,high-rate regime).Moreover,it is observed thatθc(ζ)increases asζincreases.The reason is that the interference correlation is further reduced by the higher mute probability,and the finite mean local delay therefore can be achieved at a high SIR threshold.However,the mean local delay increases withζwhenθis low (i.e.,a low-rate regime whereθ →0),which means that the mean local delay is mainly influenced by the active state of the typical user rather than interference at a low thresholdθ.In addition,whenζ=0.7,Dhas a better performance in the DUDA mode compared to the CUDA mode with a 15 dB factor.On the one hand,the advantage of DUDA mode over CUDA is more obvious as the SIR threshold increases,e.g.,atθ=-1 dB,the former has a 2 dB advantage over the latter.On the other hand,it can be seen that DUDA mode has a higherθc(ζ),which means that users can transmit at a higher rate with a guaranteed finite mean local delay in DUDA mode.

    Figure 3.The mean local delay and EE as a function of θ,where α=4.

    In Figure 3(b),the analytical (ana) and simulation(sim)results of EE are shown and their gaps lie within a tolerable range,and the reason for the existence of gaps is mainly due to the approximation of the interference intensity in (14).It is observed that for differentζand association rules,there exist optimal SIR thresholds making the energy efficiency optimal as a result of the combined effect of the transmission rate log(1+θ)and the STP.We also observe that DUDA mode has a higher EE compared to CUDA mode for the sameζ,because it can be ensured in DUDA mode that the connected BS can receive the strongest signal,while it cannot be achieved under CUDA whenP2B21according to(19).In addition,the energy efficiency decreases due to the decrease in throughput caused by the increase in the mute probabilityζ.So there exists a trade-off between reducing the mean local delay (or increasing the critical SIR thresholdθc(ζ)) and improving the EE.Here,however,we focus on the impact ofζon the maximum EE under the restriction of finite mean local delay since it is meaningless to discuss the performance of EE at the infinite local delay (i.e.,the message cannot be transmitted successfully).As can be seen from Figure 3(a),in the DUDA mode,whenζchanges from 0.3 to 0.5 and 0.7,the correspondingθc(ζ)increases,but the corresponding maximum EEηEE(θc) in Figure 3(b) slightly increases and decreases respectively.Clearly,there is an optimalζbetween 0.3 and 0.7 that optimizesηEE(θc)under the constraint of the finite mean local delay.(i.e.,increasingζdoes not always reduce the maximum EE at the finiteD).Moreover,the mild changes inηEE(θc) also show that the maximum EE may be achieved with an nearly stable level by increasing theθ(or transmission rate)when a higherζis adopted to increaseθc(ζ).

    4.2 Effect of Mute Probability

    Figure 4 shows the mean local delay and EE as a function of the mute probabilityζ.In Figure 4(a),first,for different values of SIR thresholdθ,there exists an optimal mute probabilityζoptthat minimizes the mean local delay in DUDA mode or CUDA mode,and it can be observed thatζoptand the corresponding optimal mean local delay increase with the increase inθ.Whenζ <ζopt,random DTX can decrease the mean local delay by reducing the interference correlation.Whenζ >ζopt,however,increasing the mute probability reduces the number of time slots that a typical user is active,i.e.,fewer transmissions are performed,which in turn increases the mean local delay.As theθincreases and tends to infinity,the probability of unsuccessful transmission also rises,then a higherζwill be needed to reduce the interference correlation (i.e.,ζopt →1).In addition,by comparing with CRE-based CUDA mode at the same SIR threshold,DUDA only needs a lowerζto have a lower bound ofD,which is certainly beneficial to achieving higher EE,as is also demonstrated in the following energy efficiency curves.

    Figure 4.The mean local delay and EE as a function of ζ,where α=4.

    In Figure 4(b),we also give analytical and simulation curves of energy efficiency with respect to the mute probability atθ=0 dB,and they can be matched well.Like our previous analysis in Figure 3,the EE degrades due to the mute probability increase that reduces the throughput of the system.We show the performance of EE in DUDA mode compared to the performance in CUDA mode based on different bias factors.It can be seen that the bias factor does have the ability to improve EE due to its ability to attenuate the effect of transmit power differences between different tiers of BSs and to balance the network load.As the bias factor increases (i.e.,P2B2/P1→1),it gets closer to the performance of DUDA.Furthermore,we compare the maximum EE that can be achieved with the same finite local delay in the DUDA and CUDA modes.WhenD=10 time slots in Figure 4(a),we first obtain the lower corresponding mute probabilityζ1in DUDA mode andζ2in 15 dB bias-based CUDA mode,respectively.Next,however,as can be seen in Figure 4(b),ζ1corresponds to an EE with a significant advantage of 3 dB relative toζ2.This shows that for the sameθ(or transmission rate),the DUDA mode,compared to CUDA mode,can achieve the same mean local delay but significantly higher energy efficiency with a lower mute probability.

    V.CONCLUSION

    First,we have analytically derived the mean local delay and energy efficiency of an uplinkK-tier HetNet with UE random DTX scheme in DUDA mode via the derivation of theb-th moments of the conditional STP.Furthermore,based on a reasonable assumption for interfering users,the simplified expressions with closed forms for the mean local delay and EE in the DUDA HetNet and the CRE-based CUDA HetNet were obtained,respectively.The analysis,along with the simulation results,provides some insight into the impact of the SIR threshold and mute probability on the mean local delay and energy efficiency.The simulation results revealed that in the DUDA mode,a finite mean local delay can be achieved via the random DTX scheme at the high SIR threshold,and the maximum EE is not always sacrificed by mute probability under the constraint of the finite local delay.In addition,comparing the performance of DUDA Het-Net and CRE-based CUDA HetNet,it has been found that the DUDA HetNet has a lower lower-bound of the mean local delay and a higher upper-bound of EE under a certain mute probability.With the same finite mean local delay as in CUDA mode,the HetNet in DUDA mode can achieve a significantly higher EE with a smaller mute probability.Further works will include clustered user distribution and overall HetNet performance with time division duplexing.

    ACKNOWLEDGEMENT

    This work was supported in part by the National Key R&D Program of China under Grant 2021YFB 2900304 and the Shenzhen Science and Technology Program under Grants KQTD20190929172545139 and ZDSYS20210623091808025.

    如何舔出高潮| 国产一区二区激情短视频 | 看十八女毛片水多多多| 午夜福利网站1000一区二区三区| 国产精品人妻久久久影院| 91精品国产国语对白视频| 97精品久久久久久久久久精品| 国产精品麻豆人妻色哟哟久久| www.av在线官网国产| 18在线观看网站| 成人无遮挡网站| 999精品在线视频| 久久精品熟女亚洲av麻豆精品| 熟女电影av网| 精品午夜福利在线看| 观看美女的网站| 侵犯人妻中文字幕一二三四区| 男女无遮挡免费网站观看| 精品国产露脸久久av麻豆| 少妇人妻 视频| 亚洲精品久久午夜乱码| 一二三四在线观看免费中文在 | 国产一区二区在线观看日韩| 国产综合精华液| 亚洲av免费高清在线观看| 成年av动漫网址| 99国产精品免费福利视频| 丝袜人妻中文字幕| 丰满乱子伦码专区| 麻豆精品久久久久久蜜桃| 久久精品aⅴ一区二区三区四区 | 中文字幕人妻熟女乱码| 人人妻人人添人人爽欧美一区卜| 香蕉精品网在线| 欧美人与善性xxx| 1024视频免费在线观看| 国产亚洲欧美精品永久| 在线 av 中文字幕| 欧美变态另类bdsm刘玥| 亚洲精品,欧美精品| 97在线视频观看| a级毛片黄视频| 国产色爽女视频免费观看| 90打野战视频偷拍视频| 亚洲av.av天堂| 王馨瑶露胸无遮挡在线观看| av在线老鸭窝| 日韩av不卡免费在线播放| 午夜激情久久久久久久| 少妇熟女欧美另类| 尾随美女入室| 大香蕉久久成人网| av在线播放精品| 纯流量卡能插随身wifi吗| 国产一区二区三区综合在线观看 | 老女人水多毛片| 最近中文字幕2019免费版| 日本vs欧美在线观看视频| 99久久综合免费| 久久精品国产a三级三级三级| 女性生殖器流出的白浆| 我的女老师完整版在线观看| 成年人午夜在线观看视频| 在线观看一区二区三区激情| 一级片免费观看大全| 精品久久国产蜜桃| 日韩一区二区三区影片| 黄色一级大片看看| 日日撸夜夜添| 国产精品久久久久久精品电影小说| 夜夜爽夜夜爽视频| 亚洲一区二区三区欧美精品| 观看美女的网站| 欧美成人精品欧美一级黄| 亚洲精品美女久久av网站| 欧美人与性动交α欧美精品济南到 | 一区在线观看完整版| av国产精品久久久久影院| 黄色 视频免费看| 亚洲精品美女久久久久99蜜臀 | 涩涩av久久男人的天堂| 免费观看av网站的网址| 亚洲精品国产色婷婷电影| 少妇熟女欧美另类| videosex国产| 亚洲国产av新网站| 五月天丁香电影| 精品国产一区二区三区四区第35| 最新的欧美精品一区二区| 久久久久人妻精品一区果冻| 18禁动态无遮挡网站| 女人久久www免费人成看片| 免费播放大片免费观看视频在线观看| 春色校园在线视频观看| av视频免费观看在线观看| 精品人妻偷拍中文字幕| 91在线精品国自产拍蜜月| 女的被弄到高潮叫床怎么办| 天堂中文最新版在线下载| 五月玫瑰六月丁香| 亚洲精品成人av观看孕妇| 亚洲,欧美精品.| 香蕉国产在线看| 国精品久久久久久国模美| 久久人人97超碰香蕉20202| 亚洲av电影在线进入| av片东京热男人的天堂| 热99国产精品久久久久久7| 亚洲国产看品久久| 91国产中文字幕| 午夜免费鲁丝| 丁香六月天网| 亚洲一区二区三区欧美精品| 黄色怎么调成土黄色| 纯流量卡能插随身wifi吗| 夫妻性生交免费视频一级片| 国产激情久久老熟女| 亚洲欧美一区二区三区国产| 久久av网站| 哪个播放器可以免费观看大片| 亚洲欧美日韩卡通动漫| tube8黄色片| 中文天堂在线官网| 久久精品久久精品一区二区三区| 日本91视频免费播放| 最近最新中文字幕大全免费视频 | 国产又色又爽无遮挡免| 男人操女人黄网站| 亚洲国产最新在线播放| h视频一区二区三区| 丝袜人妻中文字幕| 热re99久久国产66热| 久久久久久久亚洲中文字幕| 日韩av不卡免费在线播放| 人人澡人人妻人| 久久久久久人人人人人| 国产福利在线免费观看视频| 各种免费的搞黄视频| 王馨瑶露胸无遮挡在线观看| 日韩欧美一区视频在线观看| 美女主播在线视频| 春色校园在线视频观看| 午夜福利,免费看| 日韩精品有码人妻一区| 精品一区二区三区视频在线| 日韩,欧美,国产一区二区三区| 婷婷色综合大香蕉| 18禁在线无遮挡免费观看视频| 少妇精品久久久久久久| 制服诱惑二区| 又黄又粗又硬又大视频| 日韩成人伦理影院| 中国三级夫妇交换| 欧美丝袜亚洲另类| 老司机影院毛片| 中文精品一卡2卡3卡4更新| 亚洲精品成人av观看孕妇| 国产欧美亚洲国产| 日韩三级伦理在线观看| 国产成人欧美| 久久99一区二区三区| 国产高清国产精品国产三级| 久久久国产欧美日韩av| 制服诱惑二区| 日韩免费高清中文字幕av| 精品视频人人做人人爽| 精品国产一区二区久久| 人妻系列 视频| www.av在线官网国产| 黑人高潮一二区| 99热6这里只有精品| 久久久久国产网址| 国产精品三级大全| 老熟女久久久| av在线观看视频网站免费| 男男h啪啪无遮挡| 十八禁高潮呻吟视频| 最近最新中文字幕大全免费视频 | av黄色大香蕉| 免费人妻精品一区二区三区视频| 一本大道久久a久久精品| 欧美老熟妇乱子伦牲交| 少妇人妻久久综合中文| 天堂俺去俺来也www色官网| videosex国产| 亚洲av电影在线观看一区二区三区| 内地一区二区视频在线| 亚洲av国产av综合av卡| av电影中文网址| 少妇的逼水好多| 免费人妻精品一区二区三区视频| 亚洲精品日本国产第一区| 亚洲精品一二三| 视频区图区小说| 丝袜美足系列| av在线app专区| 丝袜人妻中文字幕| www日本在线高清视频| 久久av网站| 极品人妻少妇av视频| 亚洲色图综合在线观看| 久久精品国产鲁丝片午夜精品| 亚洲成色77777| 亚洲天堂av无毛| 亚洲精品国产av成人精品| 久久亚洲国产成人精品v| 久久这里只有精品19| 日韩 亚洲 欧美在线| 亚洲精品美女久久久久99蜜臀 | 少妇的逼水好多| 久久这里只有精品19| 国产黄色视频一区二区在线观看| 国产精品一区二区在线观看99| 大片电影免费在线观看免费| 高清在线视频一区二区三区| 国产精品一国产av| 最近的中文字幕免费完整| 看免费成人av毛片| 男女国产视频网站| 亚洲欧美色中文字幕在线| 久久99蜜桃精品久久| 夜夜骑夜夜射夜夜干| 亚洲欧美色中文字幕在线| 亚洲色图 男人天堂 中文字幕 | 亚洲经典国产精华液单| 国产激情久久老熟女| 成年av动漫网址| 男女啪啪激烈高潮av片| 国产精品熟女久久久久浪| 日韩,欧美,国产一区二区三区| 午夜激情久久久久久久| 丁香六月天网| 一区在线观看完整版| a级毛片黄视频| 性高湖久久久久久久久免费观看| 亚洲国产精品一区三区| h视频一区二区三区| 国产精品秋霞免费鲁丝片| 免费黄色在线免费观看| 国产成人精品无人区| 精品少妇久久久久久888优播| 国产日韩欧美在线精品| 欧美日韩视频精品一区| 少妇 在线观看| 久久久久久久久久人人人人人人| 啦啦啦中文免费视频观看日本| 精品一区二区三卡| 18禁裸乳无遮挡动漫免费视频| 丰满饥渴人妻一区二区三| 国产精品一区二区在线观看99| 成人影院久久| 久久精品熟女亚洲av麻豆精品| 亚洲美女视频黄频| 国产精品国产三级国产专区5o| 日韩熟女老妇一区二区性免费视频| 男女边吃奶边做爰视频| 久久精品国产a三级三级三级| 免费人妻精品一区二区三区视频| 最近2019中文字幕mv第一页| 最近中文字幕2019免费版| 青青草视频在线视频观看| 天天影视国产精品| 色婷婷av一区二区三区视频| 久久99热6这里只有精品| 国产免费一级a男人的天堂| 国产成人a∨麻豆精品| 97超碰精品成人国产| 最近最新中文字幕免费大全7| 热re99久久精品国产66热6| 寂寞人妻少妇视频99o| 亚洲色图综合在线观看| 99国产精品免费福利视频| 国产激情久久老熟女| 国产欧美日韩一区二区三区在线| 成人毛片60女人毛片免费| 视频中文字幕在线观看| 亚洲精华国产精华液的使用体验| 各种免费的搞黄视频| 男人舔女人的私密视频| 成人黄色视频免费在线看| 亚洲国产色片| 午夜福利在线观看免费完整高清在| 久久国产亚洲av麻豆专区| 久久亚洲国产成人精品v| 少妇的逼好多水| 97在线人人人人妻| 日韩制服骚丝袜av| 国产精品蜜桃在线观看| 欧美日韩视频精品一区| 热re99久久精品国产66热6| 永久免费av网站大全| 嫩草影院入口| 亚洲,欧美,日韩| 国产片内射在线| 丝袜喷水一区| 极品人妻少妇av视频| 韩国高清视频一区二区三区| 91午夜精品亚洲一区二区三区| 国内精品宾馆在线| 国产成人精品在线电影| av在线播放精品| 久久亚洲国产成人精品v| 有码 亚洲区| 国产精品麻豆人妻色哟哟久久| 26uuu在线亚洲综合色| 日韩一区二区三区影片| 亚洲精品久久久久久婷婷小说| av女优亚洲男人天堂| 99国产综合亚洲精品| 中文字幕av电影在线播放| 亚洲精华国产精华液的使用体验| 中文字幕最新亚洲高清| 巨乳人妻的诱惑在线观看| 久久这里只有精品19| 建设人人有责人人尽责人人享有的| 亚洲精品aⅴ在线观看| 久久 成人 亚洲| 亚洲在久久综合| 国产日韩一区二区三区精品不卡| 97超碰精品成人国产| 国产亚洲欧美精品永久| 亚洲第一av免费看| 国产精品人妻久久久影院| 免费看av在线观看网站| 一级,二级,三级黄色视频| 成年人午夜在线观看视频| 中文字幕另类日韩欧美亚洲嫩草| 夜夜骑夜夜射夜夜干| 老司机影院毛片| 少妇猛男粗大的猛烈进出视频| 欧美激情国产日韩精品一区| 日本av手机在线免费观看| 18+在线观看网站| √禁漫天堂资源中文www| 国产一区有黄有色的免费视频| 亚洲国产av影院在线观看| 欧美日韩亚洲高清精品| 又大又黄又爽视频免费| 欧美日韩亚洲高清精品| 久久久久精品久久久久真实原创| 自线自在国产av| 久久精品国产自在天天线| 啦啦啦啦在线视频资源| 一区在线观看完整版| 热re99久久精品国产66热6| 日产精品乱码卡一卡2卡三| 国产精品熟女久久久久浪| 日韩成人伦理影院| 精品第一国产精品| 国产探花极品一区二区| 久久97久久精品| 国产成人一区二区在线| 国产综合精华液| 国产成人一区二区在线| 久久青草综合色| 人妻一区二区av| 中文字幕最新亚洲高清| 久久av网站| 中国国产av一级| 久久久欧美国产精品| 91久久精品国产一区二区三区| h视频一区二区三区| 国产男女内射视频| 精品人妻在线不人妻| 午夜久久久在线观看| 成人综合一区亚洲| 精品一区二区三区四区五区乱码 | 久久99精品国语久久久| 激情视频va一区二区三区| 丝袜人妻中文字幕| 亚洲成人一二三区av| 亚洲精品,欧美精品| 亚洲五月色婷婷综合| 久久久久久久久久久久大奶| 99热全是精品| 大香蕉久久成人网| 欧美日韩亚洲高清精品| 亚洲国产av影院在线观看| 成年av动漫网址| 中文字幕另类日韩欧美亚洲嫩草| av一本久久久久| 黑人巨大精品欧美一区二区蜜桃 | 自拍欧美九色日韩亚洲蝌蚪91| 久久久久国产精品人妻一区二区| 插逼视频在线观看| 久久热在线av| 又黄又爽又刺激的免费视频.| 国产 一区精品| 美国免费a级毛片| 国产成人午夜福利电影在线观看| 久久久久精品久久久久真实原创| 亚洲欧美一区二区三区国产| 国产av精品麻豆| 在线观看美女被高潮喷水网站| 成年女人在线观看亚洲视频| 免费播放大片免费观看视频在线观看| 亚洲人与动物交配视频| 免费av不卡在线播放| 欧美精品av麻豆av| 久久久a久久爽久久v久久| 午夜av观看不卡| 丝袜脚勾引网站| 人体艺术视频欧美日本| 啦啦啦中文免费视频观看日本| 日韩一区二区视频免费看| 草草在线视频免费看| 欧美3d第一页| 男女边摸边吃奶| 亚洲av免费高清在线观看| 免费大片黄手机在线观看| 久久综合国产亚洲精品| 国产精品.久久久| 王馨瑶露胸无遮挡在线观看| 男女免费视频国产| 国产精品99久久99久久久不卡 | 日本免费在线观看一区| 精品一区二区三区视频在线| 蜜臀久久99精品久久宅男| 麻豆乱淫一区二区| 欧美人与性动交α欧美精品济南到 | 我要看黄色一级片免费的| 99视频精品全部免费 在线| av不卡在线播放| 国产女主播在线喷水免费视频网站| 最近2019中文字幕mv第一页| 国产 一区精品| 久久人人97超碰香蕉20202| 美女国产视频在线观看| 国产又爽黄色视频| av一本久久久久| 99国产精品免费福利视频| 国产免费视频播放在线视频| 你懂的网址亚洲精品在线观看| 国产片特级美女逼逼视频| 亚洲精品自拍成人| 18禁动态无遮挡网站| www.色视频.com| 国产精品久久久av美女十八| 伦理电影免费视频| 黄色一级大片看看| 欧美 亚洲 国产 日韩一| 纵有疾风起免费观看全集完整版| 日韩 亚洲 欧美在线| 久久99蜜桃精品久久| 飞空精品影院首页| 99久久中文字幕三级久久日本| 日本欧美视频一区| 成人18禁高潮啪啪吃奶动态图| 捣出白浆h1v1| 国产在视频线精品| a级毛片在线看网站| 妹子高潮喷水视频| 久久久久精品久久久久真实原创| 久久久久久久亚洲中文字幕| 久久午夜综合久久蜜桃| 久久精品人人爽人人爽视色| 大香蕉97超碰在线| 丝瓜视频免费看黄片| av福利片在线| 国产精品麻豆人妻色哟哟久久| 男女边吃奶边做爰视频| 国产在线视频一区二区| 中文字幕精品免费在线观看视频 | 免费看av在线观看网站| 极品人妻少妇av视频| 夜夜爽夜夜爽视频| 国产精品女同一区二区软件| 亚洲欧美精品自产自拍| 少妇熟女欧美另类| 少妇人妻 视频| av在线app专区| 国产激情久久老熟女| 国产午夜精品一二区理论片| 在线观看美女被高潮喷水网站| 日本黄色日本黄色录像| 99久久精品国产国产毛片| 亚洲精品aⅴ在线观看| 欧美xxxx性猛交bbbb| 欧美精品一区二区免费开放| 熟妇人妻不卡中文字幕| 少妇人妻久久综合中文| 久久婷婷青草| 亚洲 欧美一区二区三区| 99国产综合亚洲精品| 国产精品99久久99久久久不卡 | a级毛片在线看网站| 国产一级毛片在线| 久久久久视频综合| 18禁国产床啪视频网站| 卡戴珊不雅视频在线播放| 亚洲欧美精品自产自拍| 热re99久久精品国产66热6| 久久这里只有精品19| 久久精品熟女亚洲av麻豆精品| 综合色丁香网| av.在线天堂| 韩国av在线不卡| 搡女人真爽免费视频火全软件| 精品国产一区二区三区四区第35| 午夜免费鲁丝| 日本-黄色视频高清免费观看| 人妻系列 视频| 精品一区二区三区四区五区乱码 | 99热网站在线观看| 免费少妇av软件| 免费大片18禁| 日本免费在线观看一区| 天堂俺去俺来也www色官网| 极品人妻少妇av视频| 久久久精品免费免费高清| 久久久久久久精品精品| 国产免费福利视频在线观看| 黄网站色视频无遮挡免费观看| 99热国产这里只有精品6| 国产精品久久久av美女十八| 欧美日韩视频高清一区二区三区二| 99国产综合亚洲精品| 亚洲欧美日韩另类电影网站| 搡女人真爽免费视频火全软件| 夜夜爽夜夜爽视频| 国产视频首页在线观看| 成人亚洲欧美一区二区av| 狂野欧美激情性bbbbbb| 久久这里只有精品19| 欧美变态另类bdsm刘玥| 丝袜美足系列| 婷婷色麻豆天堂久久| 精品99又大又爽又粗少妇毛片| 精品午夜福利在线看| 精品少妇黑人巨大在线播放| 亚洲四区av| 国产 一区精品| 国产成人精品一,二区| 观看美女的网站| 欧美xxⅹ黑人| 欧美成人精品欧美一级黄| 亚洲av福利一区| 亚洲国产日韩一区二区| 中文乱码字字幕精品一区二区三区| 哪个播放器可以免费观看大片| 国产成人精品一,二区| 校园人妻丝袜中文字幕| 天堂中文最新版在线下载| 成人二区视频| 五月玫瑰六月丁香| 看免费成人av毛片| 人妻少妇偷人精品九色| 99热国产这里只有精品6| 日韩一区二区视频免费看| 在线观看免费高清a一片| 日韩人妻精品一区2区三区| 三上悠亚av全集在线观看| 国产精品三级大全| 人人妻人人澡人人爽人人夜夜| 亚洲五月色婷婷综合| 美女国产视频在线观看| 深夜精品福利| 国产老妇伦熟女老妇高清| 日本色播在线视频| 在线观看免费视频网站a站| 国产精品久久久久久av不卡| 男女边摸边吃奶| 十八禁高潮呻吟视频| 精品第一国产精品| 国产成人精品久久久久久| 久久久久久久大尺度免费视频| 欧美精品国产亚洲| 成人免费观看视频高清| 51国产日韩欧美| 国产欧美日韩综合在线一区二区| 精品视频人人做人人爽| 视频中文字幕在线观看| 制服丝袜香蕉在线| 晚上一个人看的免费电影| 又大又黄又爽视频免费| 久久精品熟女亚洲av麻豆精品| 国产成人a∨麻豆精品| 婷婷色麻豆天堂久久| 日本91视频免费播放| 天美传媒精品一区二区| av网站免费在线观看视频| 亚洲av在线观看美女高潮| 丰满迷人的少妇在线观看| 国产免费福利视频在线观看| 精品99又大又爽又粗少妇毛片| 精品一区在线观看国产| 99精国产麻豆久久婷婷| 亚洲精品中文字幕在线视频| 欧美精品国产亚洲| videosex国产| 建设人人有责人人尽责人人享有的| 黑人高潮一二区| 大香蕉久久网| 最近的中文字幕免费完整| 中国三级夫妇交换| www.av在线官网国产| 亚洲五月色婷婷综合| 爱豆传媒免费全集在线观看| 大香蕉久久网| 国产精品国产三级国产专区5o| 亚洲国产欧美在线一区| 黄色一级大片看看| 韩国av在线不卡| 国产男人的电影天堂91| 亚洲经典国产精华液单| 欧美丝袜亚洲另类| 日本猛色少妇xxxxx猛交久久| 国产av精品麻豆| 欧美亚洲日本最大视频资源| 免费av中文字幕在线| 99国产综合亚洲精品| 亚洲国产欧美日韩在线播放| 色哟哟·www| 国产av精品麻豆| 18禁动态无遮挡网站|