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

    Exploiting the Direct Link in IRS Assisted NOMA Networks with Hardware Impairments

    2023-02-17 03:13:48ZiweiLiuXinweiYueShuoChenXuliangLiuYafeiWangandWanweiTang

    Ziwei Liu,Xinwei Yue,*,Shuo Chen,Xuliang Liu,Yafei Wang and Wanwei Tang

    1Key Laboratory of Modern Measurement&Control Technology,Ministry of Education,School of Information and Communication Engineering,Beijing Information Science and Technology University,Beijing,100101,China

    2The Party Committee Office,Yellow River Conservancy Technical Institute,Kaifeng,475001,China

    3Tangshan University,Tangshan,063000,China

    ABSTRACT Hardware impairments(HI)are always present in low-cost wireless devices.This paper investigates the outage behaviors of intelligent reflecting surface (IRS) assisted non-orthogonal multiple access (NOMA) networks by taking into account the impact of HI. Specifically, we derive the approximate and asymptotic expressions of the outage probability for theIRS-NOMA-HI networks. Based on the asymptotic results, the diversity orders under perfect self-interferencecancellation and imperfect self-interference cancellation scenarios are obtained to evaluate the performance of the considered network. In addition, the system throughput of IRS-NOMA-HI is discussed in delay-limited mode. The obtained results are provided to verify the accuracy of the theoretical analyses and reveal that:1)The outage performance and system throughput for IRS-NOMA-HI outperforms that of the IRS-assisted orthogonalmultiple access-HI(IRS-OMA-HI)networks;2)The number of IRS elements,the pass loss factors,the Rician factors,and the value of HI are pivotal to enhancing the performance of IRS-NOMAHI networks;and 3)It is recommended that effective methods of reducing HI should be used to ensure system performance,in addition to self-interference cancellation techniques.

    KEYWORDS Hardware impairments; imperfect SIC; intelligent reflecting surface; non-orthogonal multiple access; outage probability

    1 Introduction

    As a promising technology to achieve significant improvement in system throughput and energy efficiency, non-orthogonal multiple access (NOMA) has recently attracted growing attention for the sixth-generation (6G) communication networks [1,2]. Different from the traditional orthogonal multiple access (OMA), the main advantage of NOMA technology is that it allows multiple users to utilize the same resource and realize the correct demodulation by introducing the successive interference cancellation (SIC) technology, which can improve the utilization of spectrum resources [3]. It has been demonstrated that power domain NOMA can ensure much higher throughput in underlay cognitive radio networks than conventional OMA [4]. NOMA also can contribute significantly to improving the quality of service, which provides a balanced service comparatively for serving multiple devices [5,6]. In order to enhance the system’s performance, the integration of NOMA and cooperative communication has been proposed in [7,8]. In particular, the authors of [9] concluded that the outage behaviors and ergodic rate of cooperative NOMA systems are better compared to traditional cooperative communication. To achieve enhanced spectrum efficiency, the authors of [10,11] analyzed the outage performance of users for amplify-and-forward and decode-and-forward relays in cooperative NOMA systems.

    Meanwhile, intelligent reflecting surface (IRS) is deemed as another new transmission technology for future wireless networks [12-14], which consists of numerous precisely designed arrangements of electromagnetic units, and the signal is controlled by changing the phase and amplitude of the incident wave [15]. There are some advantages of IRS:i) IRS can enrich the channel scattering conditions, enhance the system’s multiplexing gain, and improve the received signal strength [16];ii) The lightweight and low power requirements of the IRS make it easy to deploy on various scattering surfaces in wireless propagation environment, such as buildings,roadside billboards, windows, and interior walls [17]; andiii) Unlike the traditional active relay,IRS reflects signals in full-duplex mode through passive beamforming and does not impose additional energy consumption and noise [18]. These features have attracted researchers to apply it to various new wireless communication scenarios. Specifically, the authors of [19] showed that IRS deployment could effectively improve wireless networks’throughput, energy efficiency and coverage performance. In [20], the IRS-assisted wireless communications were introduced in an end-toend transmission scenario, which demonstrated that deploying an IRS structure can substantially improve the pairwise error probability metric. The aforementioned works manifest that IRS yields significant performance in communication scenarios.

    Inspired by the appealing advantages of IRS and NOMA, some recent works have investigated the interplay between IRS and NOMA, which can achieve enhanced massive connectivity with higher network spectral efficiency and user fairness [21-23]. The authors of [24] demonstrated that combining IRS and NOMA is an effective strategy to enhance communication coverage and energy efficiency. Furthermore, the authors of [25] proved that IRS-NOMA could achieve superior performance compared to the counterpart IRS-OMA schemes. With the help of a 1-bit coding scheme, the outage probability and ergodic rate for IRS-NOMA networks under Rayleigh fading channels have been analyzed considering the impact of perfect SIC (pSIC) and imperfect SIC(ipSIC) scenarios [26]. Moreover, in [27], the authors assessed the impact of IRS phase-shifting design on the outage probability through two-phase configurations, i.e., a random phase and coherent phase shifting. The authors of [28] minimized the transmit power of the IRS-NOMA system by optimizing the beamforming vector at base station (BS) and the phase shifters at IRS. As a further development, the interplay between IRS and NOMA is widely investigated for various technologies like unmanned aerial vehicle [29], deep reinforcement learning [30], mobile edge computing [31] and physical layer security [32]. The previous works were applied to the ideal hardware conditions.

    However, it is challenging to achieve the perfect performance in practice with the existence of hardware impairments (HI), and residual interference [33,34]. Specifically, HI causes errors in the transmitted signal at the transmitter and reduces the detection quality of the received signal at the receiver [35]. Under the realistic assumptions of residual HI, channel estimation errors and ipSIC,the authors of [36] investigated the NOMA systems with an emphasis on reliability and security. In[37], the performance of full-duplex NOMA networks over Rician fading channels was described and shown the influence of ipSIC on the system is more pronounced than that of HI. The authors employed multiple relays with finite energy storage capability to assist the BS to two NOMA device users’transmission transmission of the NOMA network in [38], and the outage probability expressions of users were derived to evaluate the system performance. Considering the IRS-assisted wireless communications, the existence of HI will affect both the source transmit beamforming and the IRS reflecting beamforming designs [39]. To compensate for the losses caused by HI,various analog and digital signal processing schemes are discussed in [40]. Despite these efforts,the HI still degrades the performance of wireless communication networks. As a further advance,the authors of [41,42] analyzed the system performance under the impact of HI at both BS and users on various NOMA and multiple input multiple output (MIMO) networks.

    1.1 Motivations and Contributions

    Motivated by the aforementioned observations, we study the performance of an IRS-NOMA network by incorporating signal distortions caused by HI. Until now, most of the existing research contributions for IRS-NOMA with ipSIC/pSIC over Rician fading channels have considered the optimization questions. Motivated by the observations mentioned earlier [43,44], we considered an IRS-NOMA communication system over Rician channels by incorporating signal distortions caused by HI in the transmitter and receiver. In contrast to the literature [26], we extend our analytical results from the Rayleigh channels to the more general Rician fading channels. Compared to [24], the direct communication link from the BS to the users is taken into account in this paper. This is another motivation to develop this treatise. More specifically, we investigate the performance of paired users, i.e., the usernand userffor IRS-NOMA-HI networks, in terms of outage probability and system throughput in delay-limited transmission mode. Moreover, the collective impact of HI and residual interference is also considered in the considered networks.IRS-OMA-HI is shown as a comparison benchmark for better evaluation of system performance.Against the previous works, the contributions of the current paper can be summarized as follows:

    1. We derive the expressions of outage probability for the usernwith ipSIC/pSIC in IRSNOMA-HI networks over Rician fading channels. To gain the corresponding diversity order,we also derive the asymptotic outage probability expressions in the high SNR region and obtain the diversity orders of the usernwith ipSIC/pSIC schemes. We observe an error floor for the outage probability of the usernwith ipSIC in the high SNR regime.

    2. We derive the expressions of outage probability for the userfin IRS-NOMA-HI networks.Based on the asymptotic results, the diversity orders are obtained. We further confirm that the diversity order is connected with the number of IRS elements and the Rician factorK.We also derive the exact expressions of outage probability for IRS-OMA-HI.

    3. We evaluate the system throughput of IRS-NOMA-HI and IRS-OMA-HI networks in delaylimited transmission modes. We further compare the outage behaviors and system throughput of users for IRS-NOMA-HI with IRS-OMA-HI. Through extensive simulation results, it is shown that the outage performance and system throughput of IRS-NOMA-HI is superior to that of IRS-OMA-HI.

    1.2 Organization

    The remainder of this paper is structured as follows. In Section 2, the system model for IRSNOMA-HI is introduced in detail. The expressions of outage probability are derived and the system throughput is studied in Section 3. Furthermore, Section 4 gives simulation results and corresponding performance analysis, followed by conclusions and future directions in Section 5.

    2 System Model

    We consider an IRS-NOMA-HI network, as depicted in Fig. 1, where one BS communicates with two users via the assistance of an IRS. The two users are classified as the nearby user n and the distant user f, which can receive both the superimposed signal from the BS and the signal reflected by IRS. To simplify the design, the BS and users are equipped with a single antenna while an IRS consists of M configurable elements. Anddenote the phase shifting matrixes, where θm∈[0,2π) is the phase shift of the m-th reflection element.1Note that multiple antennas equipped by the BS and users will further suppress the self-interference and enhance the performance of the IRS-NOMA networks,which are set aside for our future work.All wireless communication channels for the IRS-NOMA-HI networks are assumed to follow the baseband equivalent Rician fading channel models. Note that the channel coefficients from the BS to users, from BS to IRS, and from IRS to users are denoted by hsi, hsr∈CM×1and hri∈CM×1with i ∈(n,f), respectively.2It is worth noting that the nearby user and distant user are distinguished based on the distance from the users to BS.For example,user n are near to BS,while user f are far away from BS.More specifically, letand ?hsi~C N (0,1), where αsidenotes the large-scale fading coefficients. dsi, ?and η denote the distance from BS to the user i, path loss exponent and frequency-dependent factor, respectively.The vectors of channel coefficients from the BS to IRS and from the IRS to users are given aswhereand~C N (0,1) with ι ∈(sr,ri). The Rician factor K is the ratio of the power of the direct path signal to the variance of the other multipath power, which reflects the influence of multipath scattering on signal distribution. In particular, the power allocation for the user n and the user f are anand af, satisfying the relationship an+af=1 and an<af. The coherent phase shifting has the ability to improve the performance of IRS-NOMA-HI network, where the phase shift of each reflecting element is matched with the phases of its incoming and outgoing fading channels. This paper chooses coherent phase shifting to simplify computational complexity and provide distinct analytical results. It is worth pointing out that the random phase-shifting of IRS can also affect the outage behaviors, which will be set aside in our future work. To characterize the optimal performance of the IRS-NOMA-HI networks, we assume that the IRS can obtain perfect channel state information of the users.

    2.1 IRS-NOMA-HI

    The BS simultaneously transmits the signals of two users by adopting superposition coding,and the user n receives both the superimposed signal from the BS and the signal reflected by IRS.Accordingly, the received signal is expressed by

    where Psdenotes the normalized transmission power, nndenotes the complex additive white Gaussian noise (AWGN) sample with the mean power N0at user n, and w is HI term and can be modeled by a complex Gaussian distributed random variable with zero mean and finite variance,It is worth noting that kTand kRrepresent the levels of HI at the transmitter and the receiver, respectively. Especially, they are modeled as a complex Gaussian process with zero mean and variancesrespectively.

    Figure 1: System model for IRS-NOMA-HI networks

    According to the NOMA principle, SIC is adopted at usernto detect the information of userf. Therefore, the received signal-to-interference-plus-noise ratio (SINR) at the usernfor detectingxfis given by

    whereis the transmit SNR. After informationxfis detected, it is eliminated from the received signal by performing the SIC process. Thus, the SINR at usernfor detectingxnis given by

    where?=0 and 0<?≤1 denotes the situations of pSIC and ipSIC. Without loss of generality, the complex channel coefficient of residual interference from the ipSIC is denoted ashI~C N (0,ΩI).

    With the help of IRS, the received signal at usernis given by

    wherenfis the AWGN at userfwith the mean powerN0.

    The received SINR to decodexffrom userfcan be given by

    2.2 IRS-OMA-HI

    In this subsection, the IRS-OMA-HI scheme is used as a benchmark for comparison with the NOMA scheme. The BS transmit a superimposed signal that propagates from the BS for the users through the IRS, the SINR at OMA users can be expressed as

    3 Outage Probability

    Outage probability metrics are critical to the reliability of transmissions in 6G systems. Hence,this section evaluates the outage performance of the IRS-NOMA-HI network. Specifically, we derive the expressions of outage probability for each user in the pSIC and ipSIC scenarios. Then,we obtain the diversity orders to provide further insights into the communication performance.

    3.1 Statistical Channel Characterization

    In order to facilitate the derivation of outage probabilities, we characterize the statistical properties of the Rician channel in this subsection.is an random variables (RVs), where the cumulative distribution function (CDF) of Z should be obtained. Letbe an RV, with the probability density function (PDF) can be given by

    whereK0(·)is the modified Bessel function of the second kind with order zero.

    His caretakers said he knew when it was Thursday, even though he knew little else and could not communicate how he knew. He could see, but not speak, could not even sit up unassisted. Yet, he knew when it was his day to go ride. He was only ten, and he didn t live to his teens.

    The PDF ofX=is more difficult to calculate. Based on the characteristics of the Rician fading channels, the PDF is obtained using the Laguerre series approximation, which can be approximated as

    3.2 Outage Probability for IRS-NOMA-HI

    In wireless communication networks, the theoretical analyses of outage probability are the crucial works, which can guide the design and performance optimization of wireless communication systems’practical. In the next part, the outage behaviors of usernand userffor IRS-NOMA-HI networks are investigated in details. Furthermore, it is worth pointing out that implementing MRC in IRS-NOMA-HI transmission will further enhance the performance of IRS-NOMA-HI, which are set aside for our future work.

    3.2.1 User n of Outage Probability

    For the usern, system communication is interrupted when one of the following two events occurs. Casei)Usernis unable to detect the signalxfof userf; Caseii)Userncannot decode its own signalxnafter detectingxf. Thus, the outage probability of usernfor RIS-NOMA-HI networks is expressed as

    whereγthf=2Rf-1andγthn=2Rn-1denote the SNR threshold of usernand userf. Specially,RnandRfis the target rate at the users to detectxnandxf. The outage probability expression of usernfor the IRS-NOMA-HI isis approximated by the following theorem.

    Theorem 1.The expression of outage probability for user f with ipSIC in IRS-NOMA-HI networks is approximated by

    Proof:SeeAppendix A.

    Corollary 1.For the particular case ε=0,the expression of outage probability for user f with pSIC in IRS-NOMA-HI networks is approximated by

    3.2.2 User f of Outage Probability

    Proof:SeeAppendix B.

    3.3 Outage Probability for IRS-OMA-HI

    3.4 Diversity Analysis

    To gain more insights, we select the diversity order to evaluate the outage behaviors in the high SINR region [45] based on the above analytical results. In the IRS-NOMA-HI networks, the diversity order is mathematically defined as

    whereρsis the transmit SNR, andP∞(ρs)denotes the asymptotic outage probability of the users in the high SNR regime for the considered IRS-NOMA-HI networks, which are summarized in the following corollaries.

    Remark 2.Upon substituting (19) into (17) and after some manipulations, the diversity order for the user f with pSIC in IRS-NOMA-HI is obtained as+1.It can be observed that the diversity order of user f with pSIC is in connection with the number of IRS elements and Rician factor K.

    Corollary 4.When ρs tends to infinity, the asymptotic outage probability of the user n in IRSNOMA-HI networks is written by

    Remark 3.Upon substituting(20)into(17),the diversity order of user n is equal towhich is also related to the number of IRS elements and Rician factor K.

    Remark 4.Similar to the proof of the diversity order for IRS-NOMA-HI,the diversity order of user f and user n for IRS-OMA-HI are equal torespectively,which are also related to the number of IRS elements and Rician factor K.

    We plotted the diversity order of IRS-NOMA-HI and IRS-OMA-HI networks in a table, as shown in Table 1.

    Table 1: Diversity order for IRS-NOMA-HI and IRS-OMA-HI networks

    3.5 Throughput Analysis

    The system throughput metric is paramount to evaluating the reliability of the transmission in 6G systems. Next, by using the derived results of outage probability above, the delay-limited system throughput of IRS-NOMA-HI networks can be formulated by

    whereδ∈(pSIC,ipSIC),andPfcan be obtained from (11), (12) and (14),respectively.

    4 Numerical Results

    In this section, we aim to confirm through Monte-Carlo simulations the accuracy of our previous mathematical analysis and illustrate the outage performance of the IRS-NOMA-HI networks. To verify the feasibility of the IRS-NOMA-HI networks, IRS-OMA-HI is shown as the comparison benchmarks, and the performance of the transmission schemes is evaluated through computer simulation. The simulation results are averaged over 106realizations. More specifically,the parameters adopted are presented in Table 2 to highlight the performance of the proposed scheme, which is similar to many NOMA research contributions in [10,26,27]. BPCU denotes the short for a bit per channel use and the unit of normalized distancedfor dimensionless physical quantity is generally taken as one. Furthermore, we show the impact of residual interference,target rate, Rician factor and IRS elements on the performance of IRS-NOMA-HI networks. It is worth pointing out that these dots ·are the Monte Carlo simulation values, and the lines are the theoretically derived values.

    Table 2: The parameters for numerical results

    Fig. 2 plots the outage probability of two users for IRS-NOMA-HI networksvs.the transmit SNR, withM=5,K=5,=-30 dB,Rn=0.03 andRf=0.02 BPCU. For comparison,the outage performance of the traditional IRS-OMA-HI schemes is presented as well. The exact black and red solid curves for the outage probability of usernwith pSIC/ipSIC are plotted according to (11) and (12), respectively. The blue dotted curve for asymptotic outage probability of usernaccording to (18) and (19), respectively. Furthermore, this figure depicts the outage probability of userfin (14)vs.the SNR obtained from the approximate expression in (20) and Monte-Carlo simulations from (13). The exact outage probability cures for the IRS-OMA-HI scheme are plotted according to (16). The simulation is identical to analytical curves across the entire SNR range, which verifies the outage performance. As a benchmark, the outage behavior for IRS-OMA-HI shows the worst performance. The reason is that NOMA differs from conventional OMA by using non-orthogonal communication at the transmitter side and achieving correct demodulation at the receiver side utilizing SIC techniques. The IRS-NOMA-HI has been able to provide higher user fairness relative to IRS-OMA-HI. One can observe that the outage probability with ipSIC converges to an error floor in the high SNR region and obtain a zero diversity order.This is due to the fact that the residual interference from ipSIC for IRS-NOMA-HI, which is also confirmed in Remark 1. Another observation is that the outage behavior of usernwith pSIC outperforms that of usernwith ipSIC and userffor IRS-NOMA-HI networks. This is due to the fact that the usernwith pSIC has access to a larger diversity order than userf, in line with the insights of Remark 2. Due to residual interference, the outage probability of usernwith ipSIC converges to an error floor and the diversity gain is zero. It can also be observed that the outage performance of usernwith ipSIC gets worse as the residual interference value increases.In addition, the outage performance of userncan be improved in the high SNR region, although this user has HI.

    Fig. 3 illustrates the outage probability of two users for IRS-NOMA-HI networksvs.the transmit SNR, for a simulation setting withK=5,Rn=0.03 andRf=0.02 BPCU. For high values of the number of IRS elements, the outage probability for IRS-NOMA-HI gains a steeper slope and achieves enhanced outage performance. It is noted that the diversity order found in Remark 2 and Remark 3 depends on the number of IRS elements and the values of HI.This phenomenon indicates the promising integration of IRS technologies in cooperative NOMA communications. Fig. 4 shows the outage probability of two users for IRS-NOMA-HI networks vs. the transmit SNR with different values of target rate. More specifically, the values of Rnare increased from 0.02 BPCU to 0.16 BPCU, and the values of Rfare increased from 0.01 BPCU to 0.12 BPCU. It can be observed that the outage probability decreases with the increase of target rate values.

    Figure 2: Outage probability vs. the transmit SNR, with Rf =0.03, Rn=0.02 BPCU, M=5, K=5 and =-5 dB

    Figure 3: Outage probability vs. the transmit SNR, with the different M

    Figure 4: Outage probability vs. the transmit SNR, with the different target rates

    Fig. 5 plots the outage probability of two users for IRS-NOMA-HI networksvs.the transmit SNR, withM= 5,K= 5 and= -30 dB. It can be seen that as the Rician factor decreases, the IRS-NOMA-HI network is capable of achieving enhanced outage performance. This phenomenon can be explained by the larger the Rician factorK, the more significant the impact on channel capacity and the increasing damage to the IRS-NOMA-HI network, especially as the SNR increases. As a further advance, Fig. 6 depicts the outage probabilityvs.the transmit SNR, with the different values of pass loss factors. We can observe the influence of the values of pass loss factors on the performance of IRS-NOMA-HI; the higher the values are, the better the performance of the outage probability is. This is due to the fact that the path loss is determined by the radiated spread of the transmit power and the propagation characteristics of the channel,and the channels are mainly influenced by line-of-sight component when the Rician factor is large.

    Figure 5: Outage probability vs. the transmit SNR, with the different Rician factors

    Figure 6: Outage probability vs. the transmit SNR, with the different pass loss factors

    Fig. 7 plots the outage probability of two users for IRS-NOMA-HI networksvs.the transmit SNR, withM=5,K=5 and=-30 dB. It can be observed from the figure that the bad impact of IRS-NOMA-HI networks becomes more severe than the value of HI. The perfect value of HIkT=kR=0 is shown as a benchmark, and the best outage performance can be observed for the ideal case in which both the transmitter and receiver do not experience the impact of HI.It is worth noting the importance of accurately modeling HI at the transmitter and receivers when evaluating the performance of IRS-NOMA-HI networks.

    Figure 7: Outage probability vs. the transmit SNR, with the different value of HI

    Fig. 8 plots the curve of system throughputvs.the transmit SNR in delay-limited transmission mode, for a simulation setting withRn= 1,Rf= 1 BPCU,M= 5 andK= 5. In this figure, the curves of system throughput with the transmit SNR in IRS-OMA-HI and IRSNOMA-HI schemes are drawn according to (21). It can be observed from the figure that the system throughput for IRS-NOMA-HI in delay-limited transmission mode is better than that of IRS-OMA-HI. It can be explained that the system throughput depends on achieved the outage probability of users. This finding facilitates the design of IRS in wireless systems to cater to mass connections.

    Figure 8: System throughput vs. the transmit SNR in delay-limited transmission mode

    5 Conclusion

    Although NOMA communication can raise the outage performance and capacity of the traditional OMA wireless system, the residual interference generated by ipSIC and HI imposes a detrimental effect on system performance. This paper examined the joint impact of HI and residual interference on the outage probability and system throughput for the IRS-NOMA-HI and IRS-OMA-HI networks over Rician channels. In particular, we derived the exact and asymptotic expressions of outage probability in both pSIC and ipSIC cases and analyzed the diversity orders to evaluate their performance under various scenarios. Based on the numerical results, it illustrated that the outage performance and system throughput of IRS-NOMA-HI outperformed that of IRS-OMA-HI networks, which can be a promising candidate for future systems. It is recommended that more efforts should be made to suppress the effects of HI and residual interference when deploying IRS-NOMA-HI networks in realistic scenarios.

    Funding Statement:This work was supported by the National Natural Science Foundation of China under Grants 62071052, 61901043 and theR&DProgram of Beijing Municipal Education Commission under Grant KM202011232003. The work of W. Tang was supported by Scientific Research Projects supported by Talent Engineering Training Funds of Hebei Province under Grant A202101106 and Science and Technology Project of Hebei Education Department under Grant QN2020508.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    Appendix A.Proof of Theorem 1

    According to the definition of outage probability, the outage probability of usernfor IRSNOMA-HI networks can be expressed by

    By using the Gauss-Chebyshev quadrature [46], the complex definite integral in the above equation can be approximated as

    Appendix B.Proof of Theorem 2

    Base on (13), the outage probability of userffor IRS-NOMA-HI networks can be rewritten as

    By using the Gauss-Chebyshev quadrature, the outage probability of usernfor IRS-NOMAHI networks can be finally expressed as

    国产伦在线观看视频一区| 欧美极品一区二区三区四区| 夜夜爽夜夜爽视频| 色网站视频免费| 日产精品乱码卡一卡2卡三| 国产成人精品婷婷| 亚洲四区av| 国产乱人偷精品视频| 伊人久久精品亚洲午夜| 午夜福利在线在线| 99在线人妻在线中文字幕| 一边摸一边抽搐一进一小说| 蜜桃久久精品国产亚洲av| 秋霞伦理黄片| 成人毛片a级毛片在线播放| 久久久久久伊人网av| 最近中文字幕高清免费大全6| 久久精品夜色国产| 97超碰精品成人国产| 国产极品精品免费视频能看的| 精品久久久久久电影网 | 国产精品99久久久久久久久| 国产极品天堂在线| 成人漫画全彩无遮挡| 国产在视频线精品| 久99久视频精品免费| 久久久久久久亚洲中文字幕| 在线观看66精品国产| 亚洲欧洲国产日韩| 嫩草影院新地址| 啦啦啦啦在线视频资源| 亚洲最大成人中文| 热99在线观看视频| 日韩高清综合在线| 国产乱人偷精品视频| 免费av不卡在线播放| 国内揄拍国产精品人妻在线| 丰满少妇做爰视频| 夜夜爽夜夜爽视频| 午夜福利在线在线| 波多野结衣巨乳人妻| .国产精品久久| 精品人妻视频免费看| 99九九线精品视频在线观看视频| 久热久热在线精品观看| 国产伦在线观看视频一区| 国内少妇人妻偷人精品xxx网站| 亚洲av福利一区| 美女cb高潮喷水在线观看| 中文字幕制服av| 天天一区二区日本电影三级| 国产精品av视频在线免费观看| 天美传媒精品一区二区| 内地一区二区视频在线| 欧美日韩综合久久久久久| 亚洲av成人精品一二三区| 毛片一级片免费看久久久久| 日本三级黄在线观看| 国产精品国产高清国产av| 黄片无遮挡物在线观看| 激情 狠狠 欧美| 韩国av在线不卡| 91精品一卡2卡3卡4卡| 国产女主播在线喷水免费视频网站 | 国产精品一区二区性色av| 欧美激情久久久久久爽电影| 白带黄色成豆腐渣| 麻豆久久精品国产亚洲av| 欧美日韩在线观看h| 久久精品影院6| av在线播放精品| 久久久久网色| 美女高潮的动态| 十八禁国产超污无遮挡网站| 亚洲无线观看免费| 国产极品精品免费视频能看的| av在线观看视频网站免费| 亚洲综合色惰| 亚洲精品成人久久久久久| 亚洲精品亚洲一区二区| 国产精品一区二区性色av| av天堂中文字幕网| 日本与韩国留学比较| 看非洲黑人一级黄片| av免费在线看不卡| 日本爱情动作片www.在线观看| 精品久久久久久久末码| 日韩成人av中文字幕在线观看| 亚洲精品久久久久久婷婷小说 | 国产三级在线视频| 日韩,欧美,国产一区二区三区 | 久久精品夜夜夜夜夜久久蜜豆| 我的老师免费观看完整版| 成年版毛片免费区| 日韩在线高清观看一区二区三区| 天堂影院成人在线观看| 久久草成人影院| 久久精品夜色国产| 日本欧美国产在线视频| 国产精品国产三级国产专区5o | 国产精品女同一区二区软件| 蜜桃久久精品国产亚洲av| 波野结衣二区三区在线| 亚洲内射少妇av| 51国产日韩欧美| 国产成人午夜福利电影在线观看| 国产不卡一卡二| 久久午夜福利片| 亚洲国产精品久久男人天堂| ponron亚洲| 国产视频首页在线观看| 婷婷色av中文字幕| 99视频精品全部免费 在线| 亚洲一级一片aⅴ在线观看| 国内精品宾馆在线| 国产精品人妻久久久久久| 国产三级中文精品| av卡一久久| av在线老鸭窝| 国产亚洲5aaaaa淫片| 精品免费久久久久久久清纯| 真实男女啪啪啪动态图| 婷婷色麻豆天堂久久 | 老司机影院毛片| 18禁在线无遮挡免费观看视频| 亚洲欧美清纯卡通| 久久精品熟女亚洲av麻豆精品 | 亚洲精品成人久久久久久| 男女视频在线观看网站免费| 亚洲在线观看片| 亚洲av日韩在线播放| 一级毛片电影观看 | 永久免费av网站大全| 男女边吃奶边做爰视频| 美女高潮的动态| 免费无遮挡裸体视频| 亚洲av男天堂| 亚洲成人久久爱视频| 在线播放无遮挡| 国产一区二区在线观看日韩| 一本一本综合久久| 精品久久久久久久人妻蜜臀av| 国产精品福利在线免费观看| 成人无遮挡网站| 久久久精品大字幕| av在线观看视频网站免费| av在线播放精品| videos熟女内射| 边亲边吃奶的免费视频| 欧美性感艳星| 最近最新中文字幕大全电影3| 成人毛片a级毛片在线播放| 听说在线观看完整版免费高清| 中文字幕亚洲精品专区| 边亲边吃奶的免费视频| 国产高清视频在线观看网站| 99久久精品国产国产毛片| 国产在线男女| 久久婷婷人人爽人人干人人爱| 午夜精品在线福利| 日韩,欧美,国产一区二区三区 | 小蜜桃在线观看免费完整版高清| 久久久成人免费电影| 久久久色成人| 一个人看视频在线观看www免费| 免费播放大片免费观看视频在线观看 | 亚洲欧美精品综合久久99| 黄色日韩在线| a级毛片免费高清观看在线播放| 欧美丝袜亚洲另类| 国语自产精品视频在线第100页| 精品久久久久久电影网 | 国产三级在线视频| 最近手机中文字幕大全| 中文资源天堂在线| 七月丁香在线播放| 午夜激情福利司机影院| h日本视频在线播放| 久久久久久久久久黄片| 三级国产精品片| 色5月婷婷丁香| 国产精品综合久久久久久久免费| 永久免费av网站大全| 久久久国产成人精品二区| 国产亚洲一区二区精品| 久久久国产成人精品二区| 成人午夜精彩视频在线观看| 国产成人午夜福利电影在线观看| 欧美日韩综合久久久久久| 亚洲欧洲国产日韩| 亚洲aⅴ乱码一区二区在线播放| 干丝袜人妻中文字幕| kizo精华| 国产极品天堂在线| 国产精品一区二区三区四区久久| 能在线免费看毛片的网站| 久久久久久久久久久丰满| av天堂中文字幕网| 高清日韩中文字幕在线| 一级毛片我不卡| 老师上课跳d突然被开到最大视频| av福利片在线观看| 七月丁香在线播放| 老司机影院成人| 一级毛片aaaaaa免费看小| 国产午夜福利久久久久久| 国产欧美另类精品又又久久亚洲欧美| 午夜福利在线观看吧| 天堂√8在线中文| 我的女老师完整版在线观看| 国模一区二区三区四区视频| 乱系列少妇在线播放| 亚洲在线自拍视频| 久久精品91蜜桃| 丰满人妻一区二区三区视频av| 久久这里有精品视频免费| 国产成人freesex在线| 一区二区三区高清视频在线| 国产真实伦视频高清在线观看| 国产日韩欧美在线精品| 国产精品麻豆人妻色哟哟久久 | 女人被狂操c到高潮| 尾随美女入室| 亚洲自拍偷在线| 嫩草影院新地址| 晚上一个人看的免费电影| 国产成人freesex在线| 精品久久久久久久久亚洲| 欧美激情在线99| 精品久久久久久成人av| 黄色配什么色好看| 热99在线观看视频| 午夜老司机福利剧场| 国产色爽女视频免费观看| 国内精品美女久久久久久| 久久人人爽人人爽人人片va| 国产老妇伦熟女老妇高清| 成人综合一区亚洲| 中文天堂在线官网| 亚洲欧洲日产国产| 国产精品美女特级片免费视频播放器| 小蜜桃在线观看免费完整版高清| 欧美日本亚洲视频在线播放| 日日摸夜夜添夜夜爱| 久久久精品大字幕| 只有这里有精品99| 女人十人毛片免费观看3o分钟| 中文字幕免费在线视频6| 久久久欧美国产精品| 亚洲精品一区蜜桃| 18禁动态无遮挡网站| 又爽又黄a免费视频| 免费av不卡在线播放| 精品国产三级普通话版| 欧美xxxx黑人xx丫x性爽| 99久久人妻综合| 久久久欧美国产精品| 可以在线观看毛片的网站| 国产高清不卡午夜福利| 久久久久久久久大av| 国产免费视频播放在线视频 | 精品欧美国产一区二区三| 亚洲丝袜综合中文字幕| 久久欧美精品欧美久久欧美| 69人妻影院| 久久久久久久亚洲中文字幕| 在线观看av片永久免费下载| 淫秽高清视频在线观看| 三级毛片av免费| 久久精品国产亚洲av涩爱| 亚洲成人中文字幕在线播放| 性插视频无遮挡在线免费观看| 内地一区二区视频在线| 日本免费在线观看一区| 少妇的逼好多水| 青春草国产在线视频| 国产成人freesex在线| 精品久久久久久电影网 | av在线蜜桃| 91久久精品国产一区二区成人| 国产欧美日韩精品一区二区| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美成人综合另类久久久 | 菩萨蛮人人尽说江南好唐韦庄 | 91午夜精品亚洲一区二区三区| 亚洲国产最新在线播放| 亚洲成色77777| 99九九线精品视频在线观看视频| 日韩欧美在线乱码| 国产精品三级大全| 久久久精品大字幕| 97超视频在线观看视频| 高清毛片免费看| 国产真实伦视频高清在线观看| 2021少妇久久久久久久久久久| 三级经典国产精品| 免费搜索国产男女视频| 国产精品一区二区三区四区免费观看| 狠狠狠狠99中文字幕| 91精品一卡2卡3卡4卡| 精品久久久久久久久av| 日韩视频在线欧美| 亚洲18禁久久av| 99久久人妻综合| 亚洲成人中文字幕在线播放| 免费看日本二区| 日韩精品青青久久久久久| 好男人在线观看高清免费视频| av福利片在线观看| 好男人视频免费观看在线| 看免费成人av毛片| 国产精品精品国产色婷婷| 中文字幕人妻熟人妻熟丝袜美| 午夜福利高清视频| 91在线精品国自产拍蜜月| 亚洲乱码一区二区免费版| 亚洲欧美日韩无卡精品| 久久综合国产亚洲精品| 国产视频内射| 精品酒店卫生间| 久久久久久久久久黄片| 精品久久国产蜜桃| 热99在线观看视频| 五月玫瑰六月丁香| 国产精品一区www在线观看| 欧美成人午夜免费资源| 免费无遮挡裸体视频| 亚洲伊人久久精品综合 | 久热久热在线精品观看| 九九在线视频观看精品| 国国产精品蜜臀av免费| 极品教师在线视频| 亚洲国产高清在线一区二区三| 欧美性猛交黑人性爽| 99热这里只有精品一区| 亚洲国产高清在线一区二区三| 最近视频中文字幕2019在线8| 99热这里只有精品一区| 一区二区三区高清视频在线| 嫩草影院新地址| 乱码一卡2卡4卡精品| 亚洲真实伦在线观看| 国产人妻一区二区三区在| 欧美三级亚洲精品| 亚洲最大成人av| 九九爱精品视频在线观看| 大香蕉97超碰在线| 成人鲁丝片一二三区免费| 欧美激情久久久久久爽电影| 久久99热这里只有精品18| 最近2019中文字幕mv第一页| 一个人看视频在线观看www免费| 亚州av有码| 国产乱来视频区| 成人亚洲欧美一区二区av| 九色成人免费人妻av| 国产亚洲av嫩草精品影院| 日本一二三区视频观看| 亚洲丝袜综合中文字幕| 村上凉子中文字幕在线| 成人无遮挡网站| 国产激情偷乱视频一区二区| 日本色播在线视频| 国产在视频线精品| 亚洲成av人片在线播放无| 久久99热这里只有精品18| 久久人人爽人人片av| 精品国产一区二区三区久久久樱花 | 一级爰片在线观看| 91av网一区二区| 国产真实乱freesex| 伦理电影大哥的女人| 久久精品综合一区二区三区| 亚洲不卡免费看| 日韩av不卡免费在线播放| 亚洲精品影视一区二区三区av| av专区在线播放| 亚洲内射少妇av| 免费一级毛片在线播放高清视频| 久久久久久久午夜电影| 在线播放无遮挡| 99九九线精品视频在线观看视频| 免费看光身美女| av视频在线观看入口| 国产精品一区二区三区四区免费观看| 午夜福利在线观看吧| 成人亚洲欧美一区二区av| 午夜老司机福利剧场| 熟女电影av网| 亚洲欧美中文字幕日韩二区| 男人舔奶头视频| 免费黄色在线免费观看| 人人妻人人澡欧美一区二区| 极品教师在线视频| 日韩制服骚丝袜av| 国产成人a区在线观看| 久久精品影院6| 免费一级毛片在线播放高清视频| 国产精品精品国产色婷婷| 麻豆成人av视频| 搡老妇女老女人老熟妇| 欧美不卡视频在线免费观看| videossex国产| 国产老妇女一区| 男插女下体视频免费在线播放| 国产极品精品免费视频能看的| 亚洲欧美日韩高清专用| 久久欧美精品欧美久久欧美| 免费看日本二区| 在线观看一区二区三区| 又爽又黄a免费视频| 亚洲精品色激情综合| 亚洲在线观看片| 国产一区二区在线av高清观看| 小蜜桃在线观看免费完整版高清| 日本免费a在线| 亚洲欧美日韩东京热| 一二三四中文在线观看免费高清| 精品国产露脸久久av麻豆 | 麻豆乱淫一区二区| av国产免费在线观看| 亚洲国产精品久久男人天堂| 亚洲自偷自拍三级| 欧美高清成人免费视频www| 青青草视频在线视频观看| 少妇熟女aⅴ在线视频| 嫩草影院入口| 寂寞人妻少妇视频99o| 国产精华一区二区三区| 国产白丝娇喘喷水9色精品| av国产免费在线观看| 亚洲国产欧洲综合997久久,| 女人被狂操c到高潮| 日本黄大片高清| 嫩草影院新地址| 免费人成在线观看视频色| 男人舔女人下体高潮全视频| 免费在线观看成人毛片| 国产免费一级a男人的天堂| 亚洲高清免费不卡视频| 亚洲,欧美,日韩| 九九热线精品视视频播放| 国产精品麻豆人妻色哟哟久久 | 麻豆乱淫一区二区| 午夜福利网站1000一区二区三区| 美女内射精品一级片tv| 亚洲欧美成人精品一区二区| 日本猛色少妇xxxxx猛交久久| 久久人妻av系列| 亚洲av成人精品一二三区| 国产精品久久久久久久久免| 日本与韩国留学比较| 国产乱来视频区| 亚洲精品亚洲一区二区| 老司机影院成人| 久久精品91蜜桃| 18禁裸乳无遮挡免费网站照片| 在线观看66精品国产| 中文精品一卡2卡3卡4更新| 蜜桃久久精品国产亚洲av| 久久久午夜欧美精品| 一边摸一边抽搐一进一小说| 日本熟妇午夜| 深爱激情五月婷婷| 天堂网av新在线| 日本wwww免费看| 免费看av在线观看网站| 亚洲中文字幕一区二区三区有码在线看| 久久综合国产亚洲精品| 看片在线看免费视频| 五月玫瑰六月丁香| 亚洲国产精品国产精品| 91av网一区二区| 简卡轻食公司| 亚洲成人av在线免费| 校园人妻丝袜中文字幕| 精品国内亚洲2022精品成人| 精品午夜福利在线看| 精品99又大又爽又粗少妇毛片| 久久99蜜桃精品久久| 国产中年淑女户外野战色| 国产老妇女一区| 美女内射精品一级片tv| 国产成人freesex在线| 精品午夜福利在线看| 国产一区有黄有色的免费视频 | 成人毛片60女人毛片免费| 女人十人毛片免费观看3o分钟| 内射极品少妇av片p| 嫩草影院精品99| 秋霞伦理黄片| 国产精品三级大全| 联通29元200g的流量卡| 亚洲伊人久久精品综合 | 精品国产一区二区三区久久久樱花 | 超碰av人人做人人爽久久| av视频在线观看入口| 成人综合一区亚洲| 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产亚洲网站| 天堂中文最新版在线下载 | 内地一区二区视频在线| 国产三级在线视频| 亚洲人与动物交配视频| 亚洲国产高清在线一区二区三| 激情 狠狠 欧美| 久久亚洲国产成人精品v| 中文精品一卡2卡3卡4更新| 亚洲精品一区蜜桃| 婷婷色麻豆天堂久久 | 国产极品天堂在线| 亚洲最大成人手机在线| 国产av码专区亚洲av| 岛国在线免费视频观看| 欧美日本视频| 久久精品综合一区二区三区| 熟妇人妻久久中文字幕3abv| 国内精品宾馆在线| 精品人妻视频免费看| 亚洲久久久久久中文字幕| 国产亚洲91精品色在线| 中文资源天堂在线| 成人无遮挡网站| 神马国产精品三级电影在线观看| 久久久久久久国产电影| 久久人人爽人人爽人人片va| 成人高潮视频无遮挡免费网站| 亚洲av免费高清在线观看| 欧美日韩在线观看h| 日本一本二区三区精品| 偷拍熟女少妇极品色| 一级黄色大片毛片| 51国产日韩欧美| 草草在线视频免费看| 夜夜爽夜夜爽视频| 两个人视频免费观看高清| 亚洲中文字幕一区二区三区有码在线看| 麻豆精品久久久久久蜜桃| 高清日韩中文字幕在线| 日韩在线高清观看一区二区三区| 性色avwww在线观看| av在线观看视频网站免费| 国产av一区在线观看免费| 国产高清国产精品国产三级 | 男人和女人高潮做爰伦理| 久久99蜜桃精品久久| 亚洲天堂国产精品一区在线| 欧美性感艳星| 日本免费在线观看一区| 亚洲精品456在线播放app| 国产成年人精品一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 午夜免费激情av| 国产一区二区三区av在线| 真实男女啪啪啪动态图| 嫩草影院新地址| 日本黄大片高清| 午夜福利成人在线免费观看| 国产成人福利小说| 亚洲av二区三区四区| 乱系列少妇在线播放| 草草在线视频免费看| kizo精华| 成人无遮挡网站| 中文欧美无线码| 国产一区有黄有色的免费视频 | 97热精品久久久久久| 国产乱人视频| 日韩欧美国产在线观看| 在线观看美女被高潮喷水网站| 日本黄色视频三级网站网址| 99久国产av精品| 国产亚洲午夜精品一区二区久久 | 欧美日韩国产亚洲二区| 一区二区三区乱码不卡18| 欧美日本亚洲视频在线播放| 久久久a久久爽久久v久久| 日本爱情动作片www.在线观看| 精品一区二区三区人妻视频| 亚洲国产成人一精品久久久| 久久精品综合一区二区三区| .国产精品久久| 色视频www国产| 国产精品人妻久久久影院| 少妇熟女aⅴ在线视频| 国产黄色视频一区二区在线观看 | 久久这里有精品视频免费| 久久韩国三级中文字幕| 日日啪夜夜撸| videos熟女内射| 91在线精品国自产拍蜜月| 国产精品综合久久久久久久免费| 国产在线一区二区三区精 | 99热精品在线国产| 一二三四中文在线观看免费高清| 日韩精品青青久久久久久| videos熟女内射| 免费观看性生交大片5| 午夜福利高清视频| 国产91av在线免费观看| 舔av片在线| 麻豆久久精品国产亚洲av| av线在线观看网站| 我要搜黄色片| 午夜福利网站1000一区二区三区| 色尼玛亚洲综合影院| 午夜亚洲福利在线播放| 国产精品综合久久久久久久免费| 国产av一区在线观看免费| 天美传媒精品一区二区| 欧美精品一区二区大全| 日本三级黄在线观看| 我要搜黄色片| 中文天堂在线官网| 色哟哟·www|