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

    Radio Optical Network Simulation Tool(RONST)

    2022-08-24 03:30:44YasmineAbdelhakFadyKamelMoustafaHafezHusseinKotbHaithamOmranTawfikIsmailandHassanMostafa
    Computers Materials&Continua 2022年5期

    Yasmine I.Abdelhak,F(xiàn)ady Kamel,Moustafa Hafez,Hussein E.Kotb,Haitham A.Omran,Tawfik Ismailand Hassan Mostafa,

    1Department of Electronics and Communications,F(xiàn)aculty of Engineering,Mansoura University,Dakahlia,35516,Egypt

    2Nanotechnology and Nanoelectronics Program,Zewail City of Science and Technology,Giza,12578,Egypt

    3Department of Electronics and Communications,F(xiàn)aculty of Engineering,Cairo University,Giza,12613,Egypt

    4Department of Electronics and Communications,F(xiàn)aculty of Engineering,Ain Shams University,Cairo,11566,Egypt

    5Laboratory of Micro Optics,F(xiàn)aculty of Information Engineering and Technology,German University in Cairo,11835,Cairo,Egypt

    6National Institute of Laser Enhanced Sciences,Cairo University,Giza,12613,Egypt

    7Wireless Intelligent Networks Center(WINC),Nile University,Giza,12677,Egypt

    Abstract: This paper presents a radio optical network simulation tool(RONST) for modeling optical-wireless systems.For a typical optical and electrical chain environment,performance should be optimized concurrently before system implementation.As a result,simulating such systems turns out to be a multidisciplinary problem.The governing equations are incompatible with co-simulation in the traditional environments of existing software(SW)packages.The ultra-wideband (UWB) technology is an ideal candidate for providing high-speed short-range access for wireless services.The limited wireless reach of this technology is a significant limitation.A feasible solution to the problem of extending UWB signals is to transmit these signals to endusers via optical fibers.This concept implies the need for the establishment of a dependable environment for studying such systems.Therefore,the essential novelty of the proposed SW is that it provides designers, engineers, and researchers with a dependable simulation framework that can accurately and efficiently predict and/or optimize the behavior of such systems in a single optical-electronic simulation package.Furthermore, it is supported by a strong mathematical foundation with integrated algorithms to achieve broad flexibility and low computational cost.To validate the proposed tool,RONST was deployed on an ultra-wideband over fiber(UWBoF)system.The bit error rate(BER)has been calculated over a UWBoF system,and there is good agreement between the experimental and simulated results.

    Keywords: Optical-wireless systems; mathematical modelling techniques;opto-electronic software tools;ultra-wide band over fiber systems(UWBoF)

    1 Introduction

    Among the computer-aided design(CAD)tools,software developers,and end-users,it is postulated that no simulation tool can bring a wide range of physically diverse systems onto a single fully comprehensive simulation platform.The optical-wireless system is a real example of such a principle consisting of different carriers’phenomena.Simulink, developed by MathWorks, and COMSOL approach this principle;however,both simulation platforms failed to provide a completely satisfactory solution.The same is true for optical CAD tools, which face a similar issue due to the diversity of optical sciences themselves.They are extremely diverse according to the most straightforward classification of engineering optics that lies under one of the four categories:geometrical optics,wave optics,electromagnetic optics,and quantum optics[1].For example,CAD tools developed by LINOS Products are devoted to simulating geometrical optical systems,in which the primary concern of photo design is the development of bulk-optics systems.

    On the other hand, various SW packages have been developed to simulate radio frequency(RF) and optical systems individually.However,only a few of these design tools are simultaneously concerned with the simulation of both RF and optical systems.It is readily apparent that the developers of RF CAD design tools such as the High-Frequency Structure Simulator (HFSS) and the Applied Computational Sciences tool (ACS) are unconcerned with integrating optical components, devices,and/or systems into their RF libraries.Later, we will explore some critical differentiation in the underlying principles of several of these hybrid design CAD tools.Such a research problem is addressed in the proposed radio optical network simulation tool by introducing a fully SW package for modeling optical-wireless systems over a single and effective platform without any integration process with the other tools.The physical and mathematical background has been processed using a multi-step procedure.

    At first, we have started with a semi-analytical approach to identify the key design parameters and establish the mathematical framework from which the theoretical guidelines are naturally derived.A comprehensive simulation study based on existing SW packages was carried out to verify the mathematical modeling techniques.As a result, for the first time, we were able to develop a reliable CAD tool with an adaptive mathematical solver based on different numerical modeling techniques that could be used to serve a wide range of different systems in a more applicable and efficient manner.Finally,a series of verification experiments,depending on the available equipment,have been performed to validate and check the precision of the proposed simulation tool.These experiments have been implemented to be close to realistic scenarios usually encountered in access and in-building networks.Additionally, the developed SW package will be used in conjunction with the proposed module, UWBoF, to evaluate its standalone efficiency from the perspective of an end-user.The RONST layout is developed after accumulating all of the information necessary about the components and including it in the constituting libraries.The modeling and simulation of the components’libraries are conducted via MATLAB,Python,and C++programming languages.

    The rest of this paper is organized as follows.A brief description of an UWBoF module is presented in Section II.The characteristics of the presented CAD tool are explained in Section III.The simulation and experimental results are highlighted in Section IV to validate the workability of the presented CAD tool.Finally,a conclusion has been drawn.

    2 Ultra-Wide Band Over Fiber System

    2.1 System Description

    UWB is a fast-growing successful technology and a suitable candidate for many communications applications due to the numerous advantages offered by this technology, which perfectly suits the needs of high data rate-short-range wireless communication networks[2,3].UWB signals show high robustness in dense multipath fading environments [4,5].Moreover, low power consumption and very low levels of spectral emission are some of the advantages of this technology [6].However,these properties may cause mutual interference problems due to coexistence with other conventional narrowband radio systems.Therefore, the US Federal Communications Commission (FCC) has licensed the use of the spectrum for UWB communications, with a maximum effective isotropic radiated power(EIRP)spectral density(PSD)of-41.3 dBm/MHz to avoid coexistence problems with currently existing narrowband systems[7].

    This upper spectral bound can be reached throughout a wide band of 7.5 GHz of the spectrum ranging from 3.1 to 10.6 GHz.According to the FCC report, a UWB signal should occupy at least 500 MHz of the spectrum or equivalently has a minimal fractional bandwidth of 20%[7].In the design of UWB systems,efficient integration is targeted with the other wireless technologies[8,9].Moreover,considerable utilization of the acceptable power level and the available bandwidth is one of the key challenges to generate UWB pulses with high spectral and power efficiencies while respecting the FCC spectral constraints.

    The band-limited signals that exactly follow the FCC mask within the mentioned useful UWB band possess an EIRP that does not exceed 0.56 mW or equivalently -2.5 dBm.Such extremely low radiation transmission power is further decreased due to the path loss inherent in wireless channels even if multipath propagation is prohibited and free space propagation is hypothesized.This extremely low received power effectively limits the reach of a UWB signal to about 4~10 m preventing geographically remote users from accessing such systems.A feasible solution to this problem is to transmit UWB signals to those users via optical fibers leading to the concept of the UWBoF system[10–14].Fig.1 depicts a general block diagram describing the photonic generation, transmission,and detection processes by the optical generation techniques utilized in UWBoF systems.A set of communication processing techniques such as pulse shaping and modulation are applied on a basis function,usually having a Gaussian profile,that drives the system in the optical domain by employing microwave photonic(MWP)techniques[15–19].Then,the signal can be transmitted over a standard single-mode fiber (SMF) to a remote area and finally detected in the electrical domain through a photodetector(PD).

    The detected waveform is then transmitted to the target location,which may be several kilometers away from the original access point,via a UWB transmitting antenna(Tx),where it is picked up by the remote user’s UWB receiving antenna(Rx).In the presented module,the UWB signal is derived from an impulse radio(IR)signal,a sequence of modulated impulses like waveforms with a short duration,typically less than 2 ns,that is transmitted over a radio channel.Another important consideration that has a significant impact on the overall performance of a UWBoF system is the careful design of the IR-UWB waveforms used to encode UWB information signals[20–22].

    Figure 1:Typical block diagram of an UWBoF system with optical and electrical chains

    2.2 System Background Analysis

    In order to identify the design parameters, a mathematical analysis for an UWBoF system has been performed by our research group which is documented in a series of publications[22–26].It is started with the development of closed-form expressions for the IR-UWB waveforms to maximize the received electrical signal-to-noise ratio(SNR)and ended with the appropriate analysis and design of radiative-efficient UWB antennas.Such analysis has been considered in the source codes of RONST.The source code of each component, in Fig.1, has been encapsulated in its own transfer function.The multiplication of the cascaded transfer functions of the system components leads to the overall input-output relationship.

    3 Radio Optical Network Simulation Tool

    3.1 Software Architecture

    The suggested CAD tool is a block diagram SW package with a multi-layered architecture for research or industrial development, such as testing, design, analysis, and simulation.It also works as an optimization platform for long-haul links or the last-mile connection before installing at the end-user.It consists of three layers, as shown in Fig.2.The first layer is where all the programming tasks take place.Due to the powerful capabilities offered by MATLAB and Python as a programming environment and computing tools, the core programming of RONST relies on a MATLAB-based Python source code,which is responsible for the modeling and the characterization phase.However,in our case,a collaborative simulation using C++is employed to assess the interfacing for the hardwired modules with the source codes.This high-level programming language is masked by a further attractive and easy-to-use graphical user interface(GUI),as shown in Fig.3.

    Figure 2:Multi-layer software architectural design of the proposed CAD tool.Lower layer:high level programming languages.Higher layer:presentation layer viewed to the user via a graphical user interface

    Figure 3:RONST GUI;1)Workspace:includes the system design,2)Components:optical,electrical and wireless,3)Configuration settings:for each used component in the system,and 4)Messages and logs:include errors,warnings,reports and the run status

    In order to develop a standalone SW package,which usually works under a Windows environment,the operating system interface and the installation shield shall be developed by establishing the necessary dynamic link library(.dll)files and Java scripts.The main role of this layer is to present easyto-access facilities for the user while masking most of the potential of the programming environment.The GUI workbench layer is where all the tools and component libraries exist, and the user can establish the project workspace.Typically,this layer consists of the toolbars,menus,and component libraries.Further, taskbars can be created according to the individual design of this space.The component libraries are proposed to develop a wealth of components classified under the following categories:1) Optical components, 2) RF components, 3) Microwave Photonics components, 4)Wireless Transmission, 5) Communications Functions, 6) Digital Signal Processing components, 7)General Purpose components,8)General Purpose Mathematics,9)Data Analysis and Visualization,and 10)Optimization Tools.

    3.2 Numerical Modeling Approach

    In general,the analytical solution of the full-vectorial wave equation for a nonlinear optical system with multiple discontinuities does not exist.The wave equation is a second-order partial differential equation(PDE)with four dimensions,three spatial and one temporal.It is represented as follows:

    whereEis the electric field,μ0is the permeability of free space,εis the medium permittivity andPis the induced electric polarization[27].Furthermore,the numerical solutions are extremely difficult to be implemented due to the dimensionality problem.Therefore,the integration between numerical modeling techniques is a demanding need to achieve a reliable and flexible environment and to increase the applicability of the proposed SW package.Furthermore,in order to solve a more accurate scalar form of the wave equation, assumptions and approximations based on propagation conditions and experimental results are required.

    The finite-difference time-domain(FDTD)is a numerical technique that has proved its efficiency in solving PDEs and characterizing many problems [28–31].However, in the case of increasing the dimensions or the complexity of the problem domain, the solution will need more computational resources in terms of time and memory.Further, the FDTD is not accurate enough, especially with the existence of the stair-casing problem[32,33].The finite element method(FEM)is a very accurate numerical technique.It consists of various mesh elements interpolated with different order of shape functions to cover the solution domain in the most efficient and less computational method[34,35].It is widely used even in other scientific fields due to its reliability and flexibility as a modeling method[36,37].The FEM is involved in RONST for the modal solution of the wave equation and to get the spatial distribution of the optical fields.

    For the standard telecommunication silica fibers where the nonlinear effects are relatively weak,the wave equation has been derived in the frequency domain after ignoring the nonlinear polarization as follows,

    wherek0=ω/cis the free space wave number.The separation of variables method is applied over(3)to get the radial componentthat represents the solution of the next equation,

    whereβis the propagation constant andmis the azimuthal number (=0, 1, 2,...).A scalar FEM analysis is applied by considering the Galerkin’s approach and discretizing the cross section by the proper order of mesh elements.The resulted system of equations is represented as:

    in which the matrices for each mesh element are assembled inK,the characteristic matrix,andM,the global mass matrix.The FEM solution has been developed based on the weighted-residual technique assumption[34],therefore,the matrices[K]and[M]are represented as follows:

    where {N} is the vector of the FEM weight functions and the vector transpose {N}Trepresents the FEM shape functions.

    For the propagation analysis,a non-iterative bidirectional beam propagation method(NI-BiBPM)is adopted.The NI-BiBPM is a numerical technique that proves an excellent efficiency in solving systems with multiple discontinuities.The non-iterative algorithm has achieved accurate results in only one sweep analysis reducing the dimensionality of the problem and the time of calculations[38–40].The NI-BiBPM represents such problems by transition and propagation operators.The transition operators link between the transmitted and reflected field components around each discontinuityiin the system.The propagation operators map between the peers of these components over the different discontinuities along the propagation direction,taking into account the effect of the medium in-between[38].The formulas are summarized and represented in(8),(9),and(10).

    In order to avoid the lengthy and complicated process of the modal analysis while running RONST source codes,some techniques have been considered.At first,the lumping rule was adopted to convert the global mass matrix,M, into a diagonal matrix [41].Then, the inverse of the resulted matrix is multiplied intoK.The lumped mass matrix and the modified characteristic matrix are represented in(11)and(12),respectively.

    Secondly,the square root operator in the transition and propagation formulas has been calculated using blocked Schur(BS)algorithm[42].BS has proved a reliable stabilityvs.the iterative techniques such as Taylor series expansion and Padé approximation[38].The characteristic matrix is decomposed into an upper triangular matrix representing all the eigenmodes in its diagonal entries as shown in(13),

    whereQis a unitary matrix andUis an upper triangular matrix.As a result, thedepends onand the required square root is represented by another triangular matrixLwhere:

    The elements of[L]2are equated with[U]leading to the following system of equations:

    The formulas in (15) and (16) can be solved either a column or a super-diagonal at a time.A standard blocking scheme is applied over[L]and[U]that allows more efficient use of cache memory and a fast calculation process[42].The algorithm has an ability to calculate the square root of nonsymmetrical matrices directly and to suppress the small negative and semi-simple zero eigenvalues.Therefore, all the unwanted modes such as evanescent, radiated and even the non-physical modes are treated properly.Moreover,the condition of positive definite matrices is achieved by applying the branch-cut technique in which the negative eigenvalues,representing the evanescent modes,are rotated by an angleαextends from[0°,180°]as represented in(17)[43].

    The proposed analysis guarantees more stability especially in cases with very high index contrast discontinuities such as in plasmonics structures[43].

    Most of assumptions and approximations usually depend on the situation at hands.For example,in the nonlinear photonic crystal fibers, the linear propagation assumption resulting from the weak variations in the envelope approximation may be violated.Therefore,in order to build a generalized and powerful SW package, such phenomena have been considered in RONST.The inclusion of nonlinear effects provides a study of phenomena such as the four-wave mixing and self-phase modulation which exist in practice and are ignored by most current SW packages.Starting with the wave equation in(1),the nonlinear polarizationPNLis included.Then,we have followed the derivation of the nonlinear Schr?dinger equation(NLSE),a PDE that describes the propagation in dispersive and nonlinear mediums[27].It has been solved numerically before by using the split-step Fourier method(SSFM)[27,44]and this solution is developed and considered in RONST too[45–50].The SSFM is a numerical technique that behaves much faster when compared to the conventional finite-difference methods[44].

    3.3 Capabilities and Limitations

    The purpose of this section is devoted to exploring and comparing the newly developed features of RONST with the capabilities and limitations of the three top similar and familiar SW packages that have gained a wide trust among the research community.These SW packages are VPI Transmission Maker,developed by VPI Photonics,INTERCONNECT,developed by Lumerical,and Opti-System,developed by Optiwave.In Tab.1, the properties of these tools are summarized and compared with the proposed RONST SW package.

    Table 1:Comparison between RONST and commercial software packages

    Moreover, the above simulators do not usually take into account the transmission of realistic signals.Even those which deal with transmission standards like WiFi,WiMax,LTE and LTE-A,do not assume a realistic signal.

    3.4 Software Reliability

    It is important in the development process of any simulation and modeling SW to ensure the reliability of the developed package.Several metrics can be used as a reliable measurement depending on the inspector’s performance evaluation viewpoint.In what follows are some of these reliability metrics that are considered in the proposed CAD tool,

    1.The accuracy of the results obtained in comparison with those obtained from experimental set-ups representing real world applications

    2.The simulation speed of a typical system model of moderate complexity compared to other similar SW simulation tools

    3.The diversity of the library components,analysis and visualization tools

    4.The integrity with other software packages for the purpose of co-simulation

    5.Inter-operability of the developed SW under different operating system environments

    4 Simulation Results and Discussion

    The experiments have been designed and performed with the aid of the necessary hardware tools and apparatus.These experiments simulated real-life scenarios and were compared with simulation results obtained from RONST to verify the results and adjust any fitting parameters.

    4.1 Transfer Function

    In this experiment,we have measured the s-parameters as a representation for the transfer function of the whole system implemented in Fig.4.The instruments’details of this setup are provided in Tab.2.As shown in Fig.5,a simulation setup was built by RONST based on the specifications of the listed components.The experimental setup aims to assess the capability of the RONST simulator to give results close to the experimental outputs.

    Figure 4:Schematic diagram used for measuring the transfer function of the optoelectronic link.PC:polarization controller.It is used to align the laser light polarization with the electro-optical modulator.SMF:single mode fiber

    The previous setup was performed using the following scenario.The optical wavelength,1550 nm,was applied to a single mode optical fiber of with 10 km long.The frequency of the RF signal has been changed from 10 MHz to 20 GHz.High speed photodetector was used to extract the RF signal.The transfer function of the opto-electronic link was measured by using a vector network analyzer(VNA) [51].Fig.6 shows quite matched results for the transfer function between the experimental output and RONST simulations that predict the same system behavior.Carrier suppression effect causes a dip in the transfer function of the optoelectronic link at around 17.4 GHz.

    Table 2:Components and instruments of experiment set-up in Fig.4

    Figure 5:An UWBoF module implemented by RONST,inset:the transfer function variation of the photodetector output over the frequency of the input RF signal

    Figure 6:The magnitude of the link transfer function at the output of photodetector vs.the frequency of RF signal

    4.2 Bit Error Rate vs.Input Power

    The variation in the input power of a laser source at 1550 nm, has been controlled by adding attenuators with different values at the input of an optical fiber,1 km long.The frequency of an RF signal was settled at 10 Gbps.The bit error rate (BER) was measured for the output signal using a SONET/SDH CMA5000a device(XTA module).A RONST simulation has been performed using the same parameters to calculate the BER that was validated with the experimental results as shown in Fig.7.The expected inverse proportionality between the BER and laser power has been noted in the system output results, experimental and simulation.In addition, there is a good agreement between the two sets of results,which indicates the reliability of the presented tool for modeling such systems.

    4.3 Bit Error Rate vs.Fiber Length

    The effect of the optical fiber length over the BER measurements has been studied.The singlemode fiber length was changed for four distances,3 m,1 km,35 km and 36.850 km and the BER has been calculated for each fiber length.The increase in fiber length leads to a degradation in the received power at the photodetector end,thus an increase in the BER will be recognized.The same fiber length values have been considered in the RONST module and the resulting BER was compared with the experiment output in Fig.8.The proposed CAD tool continues to demonstrate accurate results along with the efficient computations benefit.

    We have designed and implemented an efficient Tx and Rx antennas over three separation distances, 10 cm, 20 cm and 30 cm as we reported in [25].The representative transfer functions of these antennas were involved in RONST and we choose the first two separations to represent them in Fig.9.The transfer function has been recalculated again for the whole system after considering the Tx and Rx antennas.The results are represented in Fig.10 over two separation distances,10 cm and 20 cm.

    Figure 7:The BER at the photodetector output signal vs.the input power of a laser source

    Figure 8:The bit error rate in the photodetector output signal vs.the optical fiber length

    In addition, Fig.11 shows the updated results for the BERvs.input power of the laser after considering the previous setup of antennas.It can be seen that RONST is able to depict the increase in the bit error rate arising from the increase in the channel length between the Tx and Rx antennas.

    Figure 9:RONST interface shows the response of Tx and Rx antennas at separation distances of 10 cm and 20 cm

    Figure 10:The simulation results for the system transfer function over the input frequency of RF signal after considering Tx-Rx antennas with a separation distance.(a)10 cm and(b)20 cm

    Figure 11:The bit error rate vs.the input power of a laser source after considering Tx-Rx antennas at distances,10 cm and 20 cm

    5 Conclusion

    RONST is a proposed CAD tool that strongly influences various aspects of the public research community and the telecommunications sector.On the research community level,RONST will provide a low-cost alternative to the various and incompatible simulators in the fields of optical and RF engineering.Therefore,the unique features offered by RONST will encourage the research community to engage in the newly emerging field of microwave photonics.On the telecommunications service provider’s level, RONST offers better services to millions of end-users through the optimal design of their last hop connections using the developed software package and its associated modules.This leads to consumers’satisfaction and thus;increases revenues and reduces users’complaints.This paper introduces a brief description of RONST numerical techniques,besides the background mathematical analysis.Several approaches have been integrated to cover the physical phenomena of such opticalwireless systems.For the first time,the FEM-based NI-BiBPM approach is developed to reduce the dimensionality of the solution and hence,the computational resources.Mathematical techniques,such as BSA and the lumping method have been adopted to speed up the computation process.The RONST simulations have been validated by the experimental results of an UWBoF system achieving a very good agreement and reliable expectations of the behavior of such systems.

    Acknowledgement:The authors would like to thank the Opto-Nano-Electronics laboratory (ONE Lab),Cairo University,and National Telecom Regulatory Authority(NTRA),Egypt.

    Funding Statement:The authors received a funding for this study by National Telecom Regulatory Authority(NTRA),Egypt“https://www.tra.gov.e.g.,/en/”.

    Conflicts of Interest:The authors whose names are listed immediately below certify that they have not affiliations with or involvement in any organization or entity with any financial interest(such as honoraria;educational grants;participation in speakers’bureaus;membership,employment,consultancies,stock ownership,or other equity interest;and expert testimony or patent-licensing arrangements),or non-financial interest(such as personal or professional relationships,affiliations,knowledge or beliefs)in the subject matter or materials discussed in this manuscript.

    丝袜人妻中文字幕| 女性生殖器流出的白浆| 黑人巨大精品欧美一区二区mp4| 欧美日韩亚洲国产一区二区在线观看 | 91老司机精品| 国产人伦9x9x在线观看| 夫妻午夜视频| 午夜成年电影在线免费观看| 国产男女内射视频| 2018国产大陆天天弄谢| 国产精品免费大片| 成年女人毛片免费观看观看9 | 精品人妻一区二区三区麻豆| 乱人伦中国视频| 欧美日韩亚洲高清精品| av不卡在线播放| 我的亚洲天堂| 免费av中文字幕在线| 亚洲国产中文字幕在线视频| 国产成人av教育| 波多野结衣一区麻豆| 麻豆国产av国片精品| 亚洲五月婷婷丁香| 视频区图区小说| 国产国语露脸激情在线看| e午夜精品久久久久久久| 淫妇啪啪啪对白视频 | 精品久久久久久久毛片微露脸 | 亚洲五月色婷婷综合| 午夜福利视频在线观看免费| 亚洲精品中文字幕一二三四区 | 日韩免费高清中文字幕av| 嫩草影视91久久| 精品人妻熟女毛片av久久网站| 亚洲精品久久久久久婷婷小说| 亚洲精品国产一区二区精华液| 无限看片的www在线观看| 午夜免费成人在线视频| 大陆偷拍与自拍| 亚洲人成电影免费在线| 国产成人a∨麻豆精品| 婷婷成人精品国产| www.av在线官网国产| 99精品久久久久人妻精品| 中国国产av一级| 国产精品一区二区精品视频观看| 1024香蕉在线观看| 成人手机av| 伦理电影免费视频| 一边摸一边做爽爽视频免费| 国产精品二区激情视频| 成人手机av| 国产精品一区二区免费欧美 | 精品国产一区二区久久| 亚洲精品av麻豆狂野| 婷婷成人精品国产| 国产91精品成人一区二区三区 | 国产淫语在线视频| 三上悠亚av全集在线观看| 另类精品久久| 亚洲一区中文字幕在线| 9热在线视频观看99| 亚洲熟女毛片儿| 国产精品久久久久久精品古装| 成人av一区二区三区在线看 | 亚洲国产av影院在线观看| 国产97色在线日韩免费| 大码成人一级视频| 中国国产av一级| 又黄又粗又硬又大视频| 国产精品久久久久久精品电影小说| a级毛片黄视频| 美国免费a级毛片| 日本av手机在线免费观看| 亚洲成国产人片在线观看| 亚洲国产欧美在线一区| 人妻久久中文字幕网| 欧美另类亚洲清纯唯美| 亚洲欧洲日产国产| 精品第一国产精品| 精品卡一卡二卡四卡免费| 日韩免费高清中文字幕av| 一级a爱视频在线免费观看| 国产一级毛片在线| 亚洲国产看品久久| 高清视频免费观看一区二区| av福利片在线| 亚洲欧美清纯卡通| 国产精品香港三级国产av潘金莲| 欧美日韩成人在线一区二区| 成人18禁高潮啪啪吃奶动态图| 黄片播放在线免费| 欧美另类亚洲清纯唯美| av网站在线播放免费| 欧美日韩一级在线毛片| 亚洲人成电影免费在线| 操出白浆在线播放| 99久久精品国产亚洲精品| 一区二区三区精品91| 日韩有码中文字幕| 少妇的丰满在线观看| 久久久精品94久久精品| 亚洲精品第二区| 久久中文看片网| 国产黄频视频在线观看| 国产免费视频播放在线视频| www.自偷自拍.com| 成在线人永久免费视频| 多毛熟女@视频| 久久久久国产精品人妻一区二区| 女人久久www免费人成看片| bbb黄色大片| 黑人操中国人逼视频| 999精品在线视频| 精品人妻一区二区三区麻豆| 中国美女看黄片| 操出白浆在线播放| 日日夜夜操网爽| 日韩制服丝袜自拍偷拍| 大陆偷拍与自拍| 婷婷丁香在线五月| 999精品在线视频| 久久女婷五月综合色啪小说| 人妻久久中文字幕网| 中文字幕精品免费在线观看视频| 一二三四在线观看免费中文在| 高清欧美精品videossex| 亚洲精品av麻豆狂野| 国产精品熟女久久久久浪| 女人爽到高潮嗷嗷叫在线视频| 久久精品国产综合久久久| 丝袜美足系列| 精品卡一卡二卡四卡免费| 亚洲精品第二区| www.自偷自拍.com| 免费女性裸体啪啪无遮挡网站| 2018国产大陆天天弄谢| 国产视频一区二区在线看| 深夜精品福利| 各种免费的搞黄视频| 777久久人妻少妇嫩草av网站| 国产一卡二卡三卡精品| 日本av免费视频播放| 日日爽夜夜爽网站| 国产又爽黄色视频| 老鸭窝网址在线观看| 精品第一国产精品| 国产精品影院久久| 人成视频在线观看免费观看| 精品一区二区三区av网在线观看 | 超碰97精品在线观看| 日韩三级视频一区二区三区| av网站在线播放免费| 国产亚洲一区二区精品| 国产免费现黄频在线看| 91大片在线观看| 悠悠久久av| 大陆偷拍与自拍| 2018国产大陆天天弄谢| 精品国产超薄肉色丝袜足j| 亚洲精品国产精品久久久不卡| 菩萨蛮人人尽说江南好唐韦庄| 久久久久网色| 国产免费现黄频在线看| 亚洲精品国产色婷婷电影| 亚洲精品中文字幕在线视频| 久久久国产精品麻豆| 亚洲精品乱久久久久久| 12—13女人毛片做爰片一| 精品久久久久久久毛片微露脸 | 久久久国产欧美日韩av| 汤姆久久久久久久影院中文字幕| 99九九在线精品视频| 国产淫语在线视频| 人人妻人人澡人人看| 99国产精品一区二区三区| 欧美老熟妇乱子伦牲交| e午夜精品久久久久久久| 777久久人妻少妇嫩草av网站| 国产欧美日韩精品亚洲av| 国产区一区二久久| av免费在线观看网站| 日韩三级视频一区二区三区| 国产精品欧美亚洲77777| 国产成人欧美| 母亲3免费完整高清在线观看| 久久国产精品人妻蜜桃| 精品人妻1区二区| 大陆偷拍与自拍| 桃花免费在线播放| 熟女少妇亚洲综合色aaa.| 首页视频小说图片口味搜索| 亚洲人成电影观看| 狂野欧美激情性xxxx| 国产av国产精品国产| 亚洲自偷自拍图片 自拍| 少妇粗大呻吟视频| 久久久欧美国产精品| 男女之事视频高清在线观看| 久9热在线精品视频| 中文字幕最新亚洲高清| 丰满少妇做爰视频| 国产无遮挡羞羞视频在线观看| 国产真人三级小视频在线观看| 99精国产麻豆久久婷婷| 亚洲熟女精品中文字幕| 黄色视频在线播放观看不卡| 国产免费现黄频在线看| 午夜免费鲁丝| 一区福利在线观看| 性少妇av在线| 另类精品久久| 中国美女看黄片| 中文字幕高清在线视频| 1024香蕉在线观看| 欧美一级毛片孕妇| 成人av一区二区三区在线看 | 麻豆乱淫一区二区| 日韩视频一区二区在线观看| 亚洲av电影在线观看一区二区三区| 亚洲精品一区蜜桃| 两性夫妻黄色片| 美国免费a级毛片| 欧美激情 高清一区二区三区| 欧美人与性动交α欧美精品济南到| 亚洲va日本ⅴa欧美va伊人久久 | 国产精品一区二区免费欧美 | 伦理电影免费视频| 亚洲精品国产av蜜桃| 伊人亚洲综合成人网| 啦啦啦中文免费视频观看日本| 欧美97在线视频| 别揉我奶头~嗯~啊~动态视频 | 69精品国产乱码久久久| 国产精品99久久99久久久不卡| 亚洲久久久国产精品| 亚洲国产av新网站| 午夜激情av网站| 国产精品一区二区在线观看99| 国产片内射在线| 亚洲专区中文字幕在线| 久久人妻熟女aⅴ| 久久久久久久久免费视频了| 久久久精品免费免费高清| 女警被强在线播放| 一区二区三区激情视频| 亚洲av日韩在线播放| 午夜免费鲁丝| 91麻豆精品激情在线观看国产 | 欧美午夜高清在线| 超色免费av| 久久久久久亚洲精品国产蜜桃av| 亚洲第一青青草原| 91成年电影在线观看| 欧美精品人与动牲交sv欧美| 精品国产乱码久久久久久小说| 久久久精品94久久精品| 久久久久国产精品人妻一区二区| 80岁老熟妇乱子伦牲交| 啪啪无遮挡十八禁网站| 夜夜夜夜夜久久久久| 美国免费a级毛片| av天堂在线播放| 久久久水蜜桃国产精品网| 久久久久久久久久久久大奶| av又黄又爽大尺度在线免费看| 国产日韩一区二区三区精品不卡| 麻豆乱淫一区二区| 人妻 亚洲 视频| 欧美精品人与动牲交sv欧美| 搡老乐熟女国产| 丰满饥渴人妻一区二区三| 欧美黑人欧美精品刺激| 日本一区二区免费在线视频| 国产成+人综合+亚洲专区| 男人添女人高潮全过程视频| 久久九九热精品免费| 99国产精品一区二区三区| www.av在线官网国产| 正在播放国产对白刺激| 亚洲精品久久成人aⅴ小说| 亚洲欧美成人综合另类久久久| 水蜜桃什么品种好| 国产福利在线免费观看视频| 色播在线永久视频| 日韩,欧美,国产一区二区三区| 多毛熟女@视频| 欧美精品高潮呻吟av久久| 国产极品粉嫩免费观看在线| 国产无遮挡羞羞视频在线观看| 高清黄色对白视频在线免费看| 亚洲精品久久午夜乱码| 国产91精品成人一区二区三区 | 精品少妇久久久久久888优播| 日韩欧美一区视频在线观看| 日韩制服骚丝袜av| 午夜视频精品福利| 国产精品成人在线| 亚洲精品第二区| 在线看a的网站| 精品少妇内射三级| 国产欧美亚洲国产| 看免费av毛片| 99国产精品免费福利视频| 国产男女内射视频| 中文字幕制服av| 俄罗斯特黄特色一大片| 在线观看免费高清a一片| 丁香六月天网| 免费人妻精品一区二区三区视频| 又紧又爽又黄一区二区| 亚洲av欧美aⅴ国产| 亚洲欧美清纯卡通| 国产免费现黄频在线看| 国产99久久九九免费精品| av视频免费观看在线观看| 两人在一起打扑克的视频| 午夜精品国产一区二区电影| 最黄视频免费看| 国产亚洲av高清不卡| 欧美精品一区二区大全| 在线十欧美十亚洲十日本专区| 久久国产精品大桥未久av| 一级毛片女人18水好多| 深夜精品福利| 黄色怎么调成土黄色| 亚洲精品一二三| 中文精品一卡2卡3卡4更新| www.自偷自拍.com| 热99re8久久精品国产| 另类精品久久| tocl精华| 亚洲国产欧美日韩在线播放| 国产色视频综合| 欧美精品亚洲一区二区| 天天操日日干夜夜撸| 欧美乱码精品一区二区三区| 精品国产一区二区三区久久久樱花| 日本av手机在线免费观看| 欧美变态另类bdsm刘玥| 国产黄色免费在线视频| 亚洲精品第二区| 美女大奶头黄色视频| 两性午夜刺激爽爽歪歪视频在线观看 | 多毛熟女@视频| 欧美在线黄色| 久久久国产一区二区| 午夜福利视频在线观看免费| 桃花免费在线播放| 丝袜喷水一区| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久久国产电影| 别揉我奶头~嗯~啊~动态视频 | 一本色道久久久久久精品综合| 亚洲精品国产av蜜桃| 麻豆国产av国片精品| 在线观看www视频免费| 18在线观看网站| 一级黄色大片毛片| 1024视频免费在线观看| 97精品久久久久久久久久精品| 啪啪无遮挡十八禁网站| 一级a爱视频在线免费观看| 国产精品久久久久久人妻精品电影 | 美国免费a级毛片| 一个人免费在线观看的高清视频 | kizo精华| 欧美大码av| 国产精品久久久人人做人人爽| 老司机午夜福利在线观看视频 | av视频免费观看在线观看| 国产男女内射视频| 久久毛片免费看一区二区三区| svipshipincom国产片| 国产xxxxx性猛交| 亚洲免费av在线视频| 亚洲全国av大片| 日本一区二区免费在线视频| 亚洲伊人色综图| 国产国语露脸激情在线看| 亚洲 国产 在线| 国产欧美日韩一区二区精品| 亚洲精品中文字幕在线视频| 麻豆国产av国片精品| 黄色 视频免费看| 久久香蕉激情| 黑人巨大精品欧美一区二区蜜桃| 热99re8久久精品国产| 日韩制服丝袜自拍偷拍| 女人精品久久久久毛片| 巨乳人妻的诱惑在线观看| 少妇的丰满在线观看| 国产在线观看jvid| 日韩中文字幕视频在线看片| 成年动漫av网址| 老司机午夜福利在线观看视频 | 国产av一区二区精品久久| 自拍欧美九色日韩亚洲蝌蚪91| 精品人妻一区二区三区麻豆| 中文字幕色久视频| 亚洲欧美色中文字幕在线| 亚洲成国产人片在线观看| 国产成人系列免费观看| 国产精品免费大片| 91麻豆精品激情在线观看国产 | 丝袜美腿诱惑在线| 亚洲自偷自拍图片 自拍| 精品免费久久久久久久清纯 | 亚洲av成人一区二区三| 少妇裸体淫交视频免费看高清 | 中文字幕人妻丝袜制服| 又黄又粗又硬又大视频| 满18在线观看网站| 国产精品99久久99久久久不卡| 久久免费观看电影| 亚洲欧美精品综合一区二区三区| a级毛片在线看网站| 一级黄色大片毛片| 亚洲熟女精品中文字幕| 久久久久久久国产电影| 女人精品久久久久毛片| 天天影视国产精品| 久久久久国产一级毛片高清牌| 在线观看一区二区三区激情| 黑人猛操日本美女一级片| av视频免费观看在线观看| 欧美大码av| 99国产精品免费福利视频| av国产精品久久久久影院| 啦啦啦 在线观看视频| 精品高清国产在线一区| 黄频高清免费视频| 下体分泌物呈黄色| 国产精品偷伦视频观看了| 黄色 视频免费看| 日本一区二区免费在线视频| 亚洲成人免费av在线播放| 18禁裸乳无遮挡动漫免费视频| 欧美在线黄色| 欧美大码av| 在线天堂中文资源库| 看免费av毛片| 午夜福利在线免费观看网站| 欧美人与性动交α欧美软件| 国产精品久久久久久人妻精品电影 | 蜜桃国产av成人99| 国产麻豆69| 国产成人系列免费观看| 午夜成年电影在线免费观看| 欧美激情 高清一区二区三区| 欧美大码av| 捣出白浆h1v1| 丝袜美足系列| 人成视频在线观看免费观看| 亚洲精品中文字幕在线视频| 国产区一区二久久| 丁香六月天网| 如日韩欧美国产精品一区二区三区| 精品卡一卡二卡四卡免费| 丰满少妇做爰视频| 日韩,欧美,国产一区二区三区| 国产激情久久老熟女| 精品亚洲乱码少妇综合久久| 久久九九热精品免费| 侵犯人妻中文字幕一二三四区| 免费女性裸体啪啪无遮挡网站| 丁香六月欧美| 80岁老熟妇乱子伦牲交| 亚洲免费av在线视频| 啦啦啦啦在线视频资源| 一级,二级,三级黄色视频| 老司机影院成人| 黄频高清免费视频| 国产亚洲精品一区二区www | 亚洲色图综合在线观看| 精品少妇黑人巨大在线播放| 欧美久久黑人一区二区| 99国产精品免费福利视频| 精品卡一卡二卡四卡免费| 永久免费av网站大全| 精品国产超薄肉色丝袜足j| 国产欧美日韩综合在线一区二区| 欧美在线黄色| 黄色怎么调成土黄色| 99国产综合亚洲精品| 亚洲av片天天在线观看| 国产精品.久久久| 香蕉丝袜av| 国产日韩一区二区三区精品不卡| 久久狼人影院| 亚洲精品乱久久久久久| 黑人巨大精品欧美一区二区蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利影视在线免费观看| 99re6热这里在线精品视频| av网站在线播放免费| 国产日韩欧美在线精品| 国产精品欧美亚洲77777| 欧美97在线视频| 丝袜喷水一区| 亚洲免费av在线视频| 国产亚洲精品第一综合不卡| 欧美日韩一级在线毛片| 女人被躁到高潮嗷嗷叫费观| 国产一区二区三区av在线| 国产精品免费视频内射| 亚洲精品自拍成人| 欧美精品亚洲一区二区| 国产精品av久久久久免费| 午夜91福利影院| 91av网站免费观看| 色婷婷久久久亚洲欧美| 国产亚洲精品久久久久5区| 不卡av一区二区三区| 高清黄色对白视频在线免费看| 女性被躁到高潮视频| 免费观看av网站的网址| 亚洲精品国产av成人精品| 天堂中文最新版在线下载| a级毛片黄视频| 最近最新免费中文字幕在线| 九色亚洲精品在线播放| av超薄肉色丝袜交足视频| 在线十欧美十亚洲十日本专区| 人成视频在线观看免费观看| 秋霞在线观看毛片| 1024视频免费在线观看| 999精品在线视频| 国产国语露脸激情在线看| 国产一区有黄有色的免费视频| 国产精品欧美亚洲77777| 国产欧美日韩一区二区三 | 精品高清国产在线一区| 国产精品国产av在线观看| 亚洲第一欧美日韩一区二区三区 | 另类精品久久| 如日韩欧美国产精品一区二区三区| 成人国产一区最新在线观看| 十分钟在线观看高清视频www| 国产主播在线观看一区二区| 黄片播放在线免费| 午夜福利免费观看在线| 欧美精品啪啪一区二区三区 | 波多野结衣一区麻豆| 十八禁高潮呻吟视频| 性少妇av在线| av在线app专区| 首页视频小说图片口味搜索| 亚洲国产毛片av蜜桃av| 精品久久蜜臀av无| 男人添女人高潮全过程视频| 亚洲情色 制服丝袜| 精品少妇黑人巨大在线播放| 亚洲精品久久午夜乱码| 五月开心婷婷网| 国产在视频线精品| 伊人亚洲综合成人网| 老汉色∧v一级毛片| 色婷婷av一区二区三区视频| av欧美777| 久久热在线av| 人人妻人人爽人人添夜夜欢视频| 精品国产乱子伦一区二区三区 | 我的亚洲天堂| 俄罗斯特黄特色一大片| 久久久久国产精品人妻一区二区| 国产av又大| 欧美亚洲日本最大视频资源| 视频区图区小说| 老熟女久久久| 法律面前人人平等表现在哪些方面 | 亚洲av日韩精品久久久久久密| 国产免费一区二区三区四区乱码| 国产精品自产拍在线观看55亚洲 | 免费久久久久久久精品成人欧美视频| 黄片小视频在线播放| 9191精品国产免费久久| 亚洲精华国产精华精| 美女视频免费永久观看网站| 国产一区二区在线观看av| 亚洲第一av免费看| 亚洲国产中文字幕在线视频| 老汉色av国产亚洲站长工具| 欧美 亚洲 国产 日韩一| 国产精品 国内视频| 97精品久久久久久久久久精品| 18在线观看网站| 12—13女人毛片做爰片一| 99精品久久久久人妻精品| 精品亚洲成国产av| 国产精品 欧美亚洲| 免费在线观看完整版高清| 午夜日韩欧美国产| 电影成人av| 欧美日韩黄片免| 人人妻人人澡人人看| 日韩制服丝袜自拍偷拍| 精品国产一区二区久久| 亚洲第一av免费看| 精品第一国产精品| 色94色欧美一区二区| 欧美精品啪啪一区二区三区 | 国产99久久九九免费精品| 丝袜美腿诱惑在线| 欧美日韩精品网址| 免费观看a级毛片全部| videos熟女内射| www.熟女人妻精品国产| 99国产精品免费福利视频| 99九九在线精品视频| 亚洲国产欧美一区二区综合| 美女视频免费永久观看网站| 欧美黄色淫秽网站| 欧美午夜高清在线| 精品高清国产在线一区|