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

    Variational Bayesian Based IMM Robust GPS Navigation Filter

    2022-08-24 12:58:06DahJingJwoandWeiYehChang
    Computers Materials&Continua 2022年7期

    Dah-Jing Jwoand Wei-Yeh Chang

    1Department of Communications, Navigation and Control Engineering, National Taiwan Ocean University Keelung,202301, Taiwan

    2TDK Taiwan Corp., Yangmei, Taoyuan, 326021, Taiwan

    Abstract: This paper investigates the navigational performance of Global Positioning System (GPS) using the variational Bayesian (VB) based robust filter with interacting multiple model (IMM) adaptation as the navigation processor.The performance of the state estimation for GPS navigation processing using the family of Kalman filter (KF) may be degraded due to the fact that in practical situations the statistics of measurement noise might change.In the proposed algorithm, the adaptivity is achieved by estimating the timevarying noise covariance matrices based onVB learning using the probabilistic approach, where in each update step, both the system state and time-varying measurement noise were recognized as random variables to be estimated.The estimation is iterated recursively at each time to approximate the real joint posterior distribution of state using the VB learning.One of the two major classical adaptive Kalman filter (AKF) approaches that have been proposed for tuning the noise covariance matrices is the multiple model adaptive estimate (MMAE).The IMM algorithm uses two or more filters to process in parallel, where each filter corresponds to a different dynamic or measurement model.The robust Huber’s M-estimation-based extended Kalman filter (HEKF) algorithm integrates both merits of the Huber M-estimation methodology and EKF.The robustness is enhanced by modifying the filter update based on Huber’s M-estimation method in the filtering framework.The proposed algorithm, referred to as the interactive multi-model based variational Bayesian HEKF (IMM-VBHEKF), provides an effective way for effectively handling the errors with time-varying and outlying property of non-Gaussian interference errors, such as the multipath effect.Illustrative examples are given to demonstrate the navigation performance enhancement in terms of adaptivity and robustness at the expense of acceptable additional execution time.

    Keywords: GPS; variational bayesian; Huber’s M-estimation; interacting multiple model; adaptive; outlier; multipath

    1 Introduction

    Non-Gaussian noise is often encountered in many practical environments and the estimation performance deteriorates dramatically.Known to be one of the dominant error sources in high accuracy global navigation satellite systems (GNSS) [1] positioning systems, such as the Global Positioning System (GPS) [2,3], multipath effects [2,4] occur when GPS signals arrive at a receiver site via multiple paths due to reflections from nearby objects, such as the ground and water surfaces,buildings, vehicles, hills, trees, etc.Many estimation algorithms have been studied to eliminate the positioning error caused by multipath.Since multipath errors are among uncorrelated errors that are not cancelled out during observation differencing, the performance of high precision GPS receivers are mostly limited by the multipath induced errors.One of the most important issues in GPS system performance improvement is the interference suppression techniques.

    The well-known Kalman filter (KF) [5] provides optimal estimate, namely minimum mean square error (MSE), of the system state vector and has been recognized as one of the most powerful state estimation techniques.As a recursive minimum norm estimation technique, the KF exhibits sensitivity to deviations from the assumed underlying probability distributions.Due to its simple structure, stable performance and low computational complexity, the conventional adaptive filtering algorithm where the least MSE is involved has been widely used in a variety of applications in the fields of adaptive signal processing and machine learning.However, the MSE criterion is limited to the assumption of linearity and Gaussianity while most of the noise in real word is non-Gaussian.The performance deteriorates significantly in the non-Gaussian noise environment.The traditional Kalman-type filter provides the best filter estimate when the noise is Gaussian, but most noise in real life is unknown,uncertain and non-Gaussian.The extended Kalman filter (EKF) is a nonlinear version of the KF and has been widely employed as the GPS navigation processor.The fact that EKF highly depends on a predefined dynamics model forms a major drawback.Furthermore, it is not optimal when the system is disturbed by heavy-tailed (or impulsive) non-Gaussian noises.To solve the performance degradation problem with non-Gaussian errors or heavy-tailed non-Gaussian noises, some robust Kalman filters have been developed by using non-minimum MSE criterion as the optimality criterion.

    Essentially based on modification of the quadratic cost function in the filter framework, the Huber M-estimation methodology [6-9] provides between smooth norm properties for small residuals and robust norm properties for large residuals.The technique that relies on Huber’s generalized maximum likelihood (ML) methodology exhibits robustness against deviations from the commonly assumed Gaussian probability density functions and can solve the non-Gaussian distribution problem efficiently and has been successfully employed for robust state estimation, inertial navigation system and visual tracking applications.The robust Huber’s M-estimation-based EKF (HEKF) algorithm integrates both merits of the Huber M-estimation methodology and EKF.For the signals contaminated with non-Gaussian noises or outliers, a robust scheme combining the Huber M-estimation methodology and the EKF framework is beneficial where the Huber M-estimation methodology is used to reformulate the measurement information of the EKF to handle outliers and provide robustness against deviations from Gaussianity.The measurement update can be modified to enhance robustness using the Huber M-estimation methodology.

    The variational Bayesian (VB) approach [10-15] is an inference method that utilizes a simple distribution to approximate the true posterior distribution of hidden variables, usually assuming that the hidden variables are independent of each other.The VB has been developed for a wide range of models to perform approximate posterior inference at low computational cost in comparison with the sampling methods.These methods assume a simpler, analytically tractable form for the posterior.The purpose of the variation is to find a variational distribution of the posterior probability density function (pdf) to minimize the Kullback-Leibler divergence (KLD) between the true posterior pdf.Two main approaches are either to derive a factored free form distribution or to assume a fixed form posterior distribution (e.g., a multivariate Gaussian, with possibly a suitable parametrization of the model).

    The classical way of solving the problem of uncertain parameters is to use adaptive filters [16]where the model parameters or the noise statistics are estimated together with the dynamic states.The classical noise adaptive filtering approaches can be divided into Bayesian, maximum likelihood,correlation and covariance matching methods.The Bayesian approach is the most general of these and the other approaches can often be interpreted as approximations to the Bayesian approach.Examples of Bayesian approaches to noise adaptive filtering include state augmentation based methods,multiple model methods such as the interacting multiple model (IMM) algorithm [17-21].The IMM algorithm has the configuration that runs in parallel several model-matched state estimation filters, which exchange information (namely interact) at each sampling time.The IMM algorithm has been originally applied to target tracking, and recently extended to navigation application.A model probability evaluator calculates the current probability of the vehicle being in each of the possible modes.A global estimate of the vehicle’s state is computed using the latest mode probabilities.The algorithm carries out a soft-switching between the various modes by adjusting the probabilities of each mode, which are used as weightings in the combined global state estimate.The covariance matrix associated with this combined estimate takes into account the covariances of the mode-conditioned estimates as well as the differences between these estimates.

    The algorithm proposed in this paper intends to provide an effective way for effectively handling the errors with time-varying and outlying property of non-Gaussian interference errors.The method utilizes the VB learning to approximate the noise strength for time-varying noise covariances, the Huber M-estimation methodology to enhance robustness especially for overcoming the problem of contamination distribution or outliers, and the IMM algorithm to furtherly tune the noise covariance matrices.The remaining of this paper is organized as follows.In Section 2, preliminary background on the Huber’s M-estimation-based extended Kalman filter is reviewed.The variational Bayesian algorithm is discussed in Section3.In Section4,the interacting multiple model algorithm is presented.In Section 5, simulation experiments are carried out to evaluate the performance and effectiveness.Conclusions are given in Section 6.

    2 Huber’s M-Estimation-Based Extended Kalman Filter

    For the nonlinear single model equations in discrete-time form, we have

    where the state vector xk∈Rn, process noise vector wk∈Rm, measurement vector zk∈Rm,measurement noise vector vk∈Rm,Qkis the process noise covariance matrix and Rkis the measurement noise covariance matrix.The vectors wkand vkin Eqs.(1) and (2) are zero mean Gaussian white sequences has zero cross-correlation with each other:

    whereE[·]represents expectation, and superscript “ T”denotes matrix transpose.The symbol δikstands for the Dirac function:

    Based on robust strategy, the Huber M-estimation methodology possesses the ability to solve pollution distribution or outliers by improving filter update.The essence of Huber’s robust estimation method is to solve the filter estimation value when the minimum cost function is obtained.For further applications,based on the generalized maximum likelihood estimation, amore general object function is given by Huber

    where ρ(·) is an appropriate positive symmetric function,υk,iis the residual.The cost function given by Huber is

    where γ is the default parameter.Amounts of experiences show that the performance is best when γ= 1.345.For Huber’s filtering robustness framework, reweighted average is used by

    where υk,iis the residual,φ(υk,i) =ρ′(υk,i) is the influence function.

    The Huber M-estimation methodology makes use of the weighting matrix to recast the measurement information.The weighting matrix can be interpreted in terms of reweighting the residual error covariance and constructing the measurement process.The modified covariance has the form

    The HEKF algorithm is summarized as follow:

    Stage 1: Prediction steps/time update equations

    (1) Predict the state vector

    (2) Predict the state error covariance matrix

    Stage 2: Correction steps/measurement update equations

    (3) Compute the residual

    (4) Compute the reweighting covariance matrix for measurement noise

    (5) Compute Kalman gain matrix

    (6) Update state vector

    (7) Update state error covariance matrix

    3 The Variational Bayesian Approach

    The variational Bayesian approach is an inference method that utilizes a simple distribution to approximate the true posterior distribution of hidden variables to find a variational distribution of the posterior probability density function (pdf) to minimize the KLD between the true posterior pdf.The variable decibel parameter is estimated by iterative calculations in which the gradient descent algorithm is involved.The difference between the lower bounds of two adjacent repeated calculations is used as the basis for judging convergence.

    3.1 Variational Bayesian Extended Kalman Filter

    Based on the Bayesian rule, the posterior distribution of the system state and the measurement noise covariance can be represented as

    which can be simplified top(zk|xk,Rk).For the first-order Markov process

    the Chapman-Kolmogorov equation can be obtained through the integration.

    Given the next measurement zk, the predictive distribution is updated to a posterior distribution using the Bayes rule

    The conditional probability distribution can be approximated as a product of Gaussian and Inverse-Gamma distributions

    whereN(?k|k-1,Pk|k-1) is the Gaussian distribution for prediction xkwith parameters?k|k-1and Pk|k-1.αk|k-1,i,βk|k-1,iare two parameters of the Inverse-Gamma distribution for, thei-th diagonal element of Rk.

    The variational Bayesian extended Kalman filter (VBEKF) employes a heuristic approach in the calculation process and predicts the distribution parameter by the first-order approximation

    The parameters of the new distribution is updated through

    and the measurement noise variance can be estimated

    The VBEKF algorithm is summarized as follows:

    Stage 1: Prediction steps/time update equations

    (1) Predict the state vector:

    (2) Predict the state error covariance matrix

    (3) Predict the distribution parameter

    Stage 2: Correction steps/measurement update equations

    (4) Compute measurement noise variance

    (5) Compute Kalman gain matrix

    (6) Update state vector

    (7) Update state error covariance matrix

    (8) Update the measurement distribution parameter

    3.2 Variational Bayesian Huber’s M-Estimation Based Extended Kalman Filter

    The positioning accuracy is degraded due to the presence of outliers and time-varying noise strength.Incorporation of the VB and Huber M-estimation methodologies into the EKF yielding the variational Bayesian Huber’s M-estimation based extended Kalman filter (VBHEKF) algorithm can effectively overcome the outliers and time-varying variance in the measurement noise and improve the positioning accuracy.The VBHEKF algorithm is summarized as follows:

    Stage 1: Prediction steps/time update equations

    (1) Predict the state vector

    (2) Predict the state error covariance matrix

    (3) Predict the distribution parameters

    Stage 2: Correction steps/measurement update equations

    (4) Compute measurement noise variance:

    (5) Compute the residual

    (6) Compute the reweighting covariance matrix for measurement noise

    (7) Compute Kalman gain matrix

    (8) Update state vector

    (9) Update state error covariance matrix

    (10) Update the measurement distribution parameter

    4 Interacting Multiple Model Extended Kalman Filter

    Based on pseudo-Bayesian theory, the IMM algorithm employs two or more filters to process in parallel where each filter corresponds to a different state space or measurement model.It is an adaptive mechanism which dynamically identifies uncertainties or modeling errors can be adopted and is mainly an algorithm involving Markov chain switching coefficients between different models.In each model an EKF is running, and the IMM algorithm makes uses of model probabilities to weight the inputs and output of a bank of parallel filters at each time instant.

    Assuming that a target with r states, corresponding to r models, set the system state equation and system observation equation represented by thej-th model are as follows:

    where the state vector xk∈Rn, process noise vector wk∈Rm, measurement vector zk∈Rm,measurement noise vector vk∈Rm,Qkis the process noise covariance matrix and Rkis the measurement noise covariance matrix.The vectors wkand vkin Eqs.(35) and (36) are zero mean Gaussian white sequences has zero cross-correlation with each other:

    whereE[·]represents expectation, and superscript “ T”denotes matrix transpose.The symbol δikstands for the Dirac function:

    The IMM-EKF algorithm is summarized as follows:

    Stage 1: Model interaction/mixing

    Compute mixed state and corresponding covariances of modelj

    Stage 2: Model filtering

    Stage 3: Model probabilities update

    Compute the likelihood function

    wherevj= zk-and.The model probability is updated according to the model likelihood and model transition probability

    Stage 4: Combination of model estimates and covariances

    The algorithm performs real-time detection by setting model filters corresponding to the number of possible models, setting weight coefficients and model update probabilities for each filter, and finally performing weighted fusion to calculate the current optimal system state estimation to achieve the purpose of model adaptation.

    5 Illustrative Examples

    To validate the effectiveness of the proposed approaches, simulation experiments have been carried out to evaluate the performance of the proposed approach in comparison with the other conventional methods for GPS navigation processing.The computer codes were developed by the authors using the Matlab?software.The commercial software Satellite Navigation (SatNav) Toolbox by GPSoft LLC [22] was employed for generation of the GPS satellite orbits/positions and thereafter, the satellite pseudoranges, carrier phase measurement, and constellation, required for simulation.

    The simulated pseudorange error sources corrupting GPS measurements include ionospheric delay, tropospheric delay, receiver noise and multipath.It is assumed that the differential GPS(DGPS)mode is available and therefore most of the receiver-independent common errors can be corrected,while the multipath and receiver thermal noise cannot be eliminated.The multipath interferences are added into the GPS pseudorange observation data during the vehicle moving.Two scenarios dealing with two types of interferences in pseudorange observables are carried out, which covers pseudorange errors involving (1) outlier type of multipath interferences during the vehicle moving, and (2) timevarying variance in the measurement noise with additional outliers.

    Fig.1 shows the test trajectory for the simulated vehicle and the skyplot during the simulation.In the simulation, there are 8 GPS satellites available.Since the research focus on the mitigation of multipath errors, the influence of measurement noise is relatively critical.A vehicle is designed to perform the uniform accelerated motion to reduce the impact caused by unmodeling system dynamic errors.First of all, performance comparison for the four estimator/filter including EKF,VBEKF, HEKF and VBHEKF is shown.Secondly, performance enhancement based on the IMM configuration is presented.

    Figure 1: Test trajectory for the simulated vehicle (left) and the skyplot (right) during the simulation

    5.1 Scenario 1: Pseudorange Observable Involving Outlier Type of Multipath Errors

    In Scenario 1, mitigation of the errors for pseudorange observable involving outlier type of multipath interferences is investigated.There are totally seven time durations where additional randomly generated errors are intentionally injected into the GPS pseudorange observation data during the vehicle moving.Tab.1 provides the information of the standard deviation for the 5 outliers.

    Table 1: Information of the standard deviation for the five outliers

    5.1.1 Performance Comparison for EKF and its Variants

    Comparison of GPS navigation accuracy for the four schemes: EKF, VBEKF, HEKF and VBHEKF is shown in Fig.2.The results show that either the VB or Huber’s algorithms can assist EKF to effectively deal with the outliers in the pesudorange observables individually and combination of the two algorithms can furtherly enhance the performance.

    Figure 2: Comparison of positioning accuracy for extended Kalman filter (EKF) and its variants-Scenario 1

    5.1.2 Performance Enhancement Based on the Interacting Multiple Model Configuration

    Incorporation of the IMM into the various filters is utilized for further performance enhancement.Fig.3 illustrates the positioning accuracy when the IMM configuration is involved for various algorithms: IMM-EKF, IMM-VBEKF, IMM-HEKF and IMM-VBHEKF.For the case that the VB does not possess sufficient capability to resolve the outlier type of interference, the IMM-VBHEKF demonstrates substantial performance improvement in navigation accuracy with acceptable extra computational expense.Tab.2 summarizes the estimation performance and execution time for various algorithms.

    Figure 3: Comparison of positioning accuracy when the interacting multiple model (IMM) configuration is involved for Scenario 1

    Table 2: Performance comparison for various algorithms-Scenario 1

    5.2 Scenario 2: Errors Involving Time-Varying Variance in Measurement Noise with Additional Outliers

    Scenario 2 is designed for investigating the performance when dealing with the time-varying Gaussian measurement noise with additional outliers.Description of time varying measurement variances in the five time intervals is shown in Tab.3.The variances of time-varying measurement noiserρiin the five time intervals for this scenario are assumed to follows the variation:rρi= 1→10→1→17→1 (in unit of m2), where the‘a(chǎn)rrows (→)’is employed for indicating the time-varying trajectory of variances.The set of unknown time-varying statistical parameters of noise needs to be estimated with the system states and the error covariance.

    Table 3: Description of the variation of noise strength and number of additional outliers in the five time intervals for Scenario 2

    5.2.1 Performance Comparison for EKF and its Variants

    Comparison of GPS navigation accuracy for the EKF, VBEKF, HEKF and VBHEKF is shown in Fig.4 where comparison of positioning accuracies without outlier and with outlier, respectively is provided.The results show that both VB and Huber based EKF can effectively improve the positioning performance.As can be seen, both the VB and Huber algorithms can be adopted into the EKF to improve GPS navigation accuracy in the environment of time-varying Gaussian noise.When additional outliers are involved in addition to the time-varying interferences, neither VB nor Huber only algorithm possesses sufficient capability to resolve the problem for such kind of challenge.The VBHEKF demonstrates substantial performance improvement in navigation accuracy at the expense of acceptable extra computational cost.

    Figure 4: Comparison of positioning accuracy for EKF and its variants-Scenario 2 (a) without outlier(b) with outlier

    5.2.2 Performance Enhancement Based on the Interacting Multiple Model Configuration

    Comparison of positioning accuracies for the four algorithms: IMM-EKF, IMM-VBEKF, IMMHEKF and IMM-VBHEKF is shown in Fig.5.Positioning accuracies for the cases without and with outliers, respectively, are presented.The Huber’s approach does not catch up the variation of noise strength very well.With the assistance of VB, the VBHEKF can further improve the performance.From the other view point, the adaptation capability of noise variance for the VBHEKF has been improved with the assistance of the IMM mechanism.Tabs.4 and 5 summarize the performance for various algorithms without and with outliers, respectively.Among the various approaches, the IMMVBHEKF provides the best positioning accuracy with only limited amount of extra execution time.

    Table 4: Performance comparison for Scenario 2-without outliers

    Table 5: Performance comparison for Scenario 2-with outliers

    Figure 5: Comparison of positioning accuracy for Scenario 2 when the IMM configuration is involved when IMM algorithm is incorporated (a) without outlier (b) with outlier

    6 Conclusions

    In this paper, an interacting multiple model based variational Bayesian robust extended Kalman filter is adopted for processing GPS navigation states.The method utilizes the variational Bayesian learning to approximate the noise strength and to provide a strong tracking capability for time-varying noise covariances; the Huber M-estimation methodology based on the robust strategy is employed to enhance robustness especially for overcoming the problem of contamination distribution or outliers;the IMM algorithm is introduced for further tuning of the noise covariance matrices.Combination the three merits leads to the estimator referred to as the IMM-VBHEKF.Illustrative examples have been provided for validation of the method.Two scenarios are presented representing two possible peudorange error patterns, including (1) observable involving outlier type of multipath errors;(2) observable involving time-varying variance in measurement noise with additional outliers.The proposed method has demonstrated performance enhancement in terms of adaptivity and robustness at the expense of acceptable additional execution time.

    Funding Statement:This work has been partially supported by the Ministry of Science and Technology,Taiwan [Grant Numbers MOST 108-2221-E-019-013 and MOST 109-2221-E-019-010].

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

    九九在线视频观看精品| 人人妻人人澡欧美一区二区| 欧美三级亚洲精品| 国产午夜福利久久久久久| 91久久精品电影网| 99久久九九国产精品国产免费| 99久久九九国产精品国产免费| 99视频精品全部免费 在线| 麻豆久久精品国产亚洲av| 亚洲精品视频女| 亚洲成人中文字幕在线播放| eeuss影院久久| 日韩大片免费观看网站| 国产精品日韩av在线免费观看| 嫩草影院新地址| 大香蕉97超碰在线| 久久99热6这里只有精品| 亚洲精品国产av蜜桃| 尾随美女入室| 日韩av不卡免费在线播放| 国产精品一及| 亚洲欧美清纯卡通| 精品久久久噜噜| 插逼视频在线观看| 亚洲av福利一区| 国产精品久久视频播放| 国产亚洲av片在线观看秒播厂 | 男女啪啪激烈高潮av片| 国产精品三级大全| 日本一二三区视频观看| 久久精品夜色国产| 国内精品宾馆在线| 午夜精品国产一区二区电影 | 日韩欧美国产在线观看| 亚洲一区高清亚洲精品| 午夜老司机福利剧场| 观看美女的网站| 日韩成人av中文字幕在线观看| a级毛片免费高清观看在线播放| 国产成人91sexporn| 免费高清在线观看视频在线观看| 国产真实伦视频高清在线观看| 久久99热这里只有精品18| 2018国产大陆天天弄谢| 天美传媒精品一区二区| 国语对白做爰xxxⅹ性视频网站| 成人av在线播放网站| 亚洲精品456在线播放app| 亚洲精品一区蜜桃| 免费看不卡的av| 成人亚洲精品一区在线观看 | 亚洲国产精品成人综合色| 一级黄片播放器| 天堂av国产一区二区熟女人妻| 人妻夜夜爽99麻豆av| 亚洲自偷自拍三级| 日韩 亚洲 欧美在线| 啦啦啦韩国在线观看视频| 人妻夜夜爽99麻豆av| 亚洲av福利一区| 好男人在线观看高清免费视频| 亚洲av成人精品一区久久| 天堂中文最新版在线下载 | 黑人高潮一二区| 欧美日本视频| 亚洲欧美清纯卡通| 欧美性感艳星| 亚洲高清免费不卡视频| 老司机影院成人| 免费黄网站久久成人精品| a级一级毛片免费在线观看| 在线天堂最新版资源| 黄色一级大片看看| 免费看光身美女| h日本视频在线播放| 男的添女的下面高潮视频| 亚洲成人中文字幕在线播放| 国产av不卡久久| av在线播放精品| 精品欧美国产一区二区三| 国产亚洲一区二区精品| 国产伦在线观看视频一区| 肉色欧美久久久久久久蜜桃 | 一级毛片黄色毛片免费观看视频| 亚洲精品aⅴ在线观看| 久久99精品国语久久久| av在线蜜桃| 欧美极品一区二区三区四区| 丝袜美腿在线中文| 国产精品伦人一区二区| 欧美日韩视频高清一区二区三区二| 六月丁香七月| 亚洲国产成人一精品久久久| 精品人妻一区二区三区麻豆| 秋霞伦理黄片| 亚洲成色77777| 亚洲成人精品中文字幕电影| 一级二级三级毛片免费看| 嫩草影院入口| 少妇熟女欧美另类| 日日啪夜夜爽| 中文资源天堂在线| 永久网站在线| 免费播放大片免费观看视频在线观看| 欧美另类一区| 亚洲精品国产av成人精品| 韩国高清视频一区二区三区| 亚洲在线观看片| 一个人看视频在线观看www免费| 身体一侧抽搐| eeuss影院久久| or卡值多少钱| 水蜜桃什么品种好| 免费少妇av软件| 国产熟女欧美一区二区| 亚洲综合色惰| 一本久久精品| 99久久人妻综合| 午夜免费男女啪啪视频观看| 少妇人妻一区二区三区视频| videossex国产| 夜夜爽夜夜爽视频| 日本色播在线视频| 男人和女人高潮做爰伦理| 80岁老熟妇乱子伦牲交| av国产免费在线观看| 舔av片在线| 国精品久久久久久国模美| 97超碰精品成人国产| 久久久久久久午夜电影| 欧美激情在线99| 国产成人aa在线观看| 高清在线视频一区二区三区| 特级一级黄色大片| 国产爱豆传媒在线观看| 国产精品久久久久久精品电影| 午夜亚洲福利在线播放| 亚洲人与动物交配视频| 午夜久久久久精精品| 欧美激情国产日韩精品一区| ponron亚洲| 亚洲av成人av| 插阴视频在线观看视频| 18禁动态无遮挡网站| 老师上课跳d突然被开到最大视频| 亚洲欧洲国产日韩| 老司机影院毛片| 美女主播在线视频| 超碰av人人做人人爽久久| a级毛片免费高清观看在线播放| av在线天堂中文字幕| 97热精品久久久久久| 色尼玛亚洲综合影院| 九九在线视频观看精品| 日韩制服骚丝袜av| 久久久精品欧美日韩精品| 熟女人妻精品中文字幕| 日韩电影二区| 精品国产一区二区三区久久久樱花 | 亚洲精品日本国产第一区| 精品一区二区三区视频在线| 如何舔出高潮| 少妇的逼水好多| 中文天堂在线官网| 免费电影在线观看免费观看| 精品久久久久久久末码| 亚洲精品456在线播放app| 欧美人与善性xxx| 国产永久视频网站| 五月玫瑰六月丁香| 国产亚洲精品久久久com| 精品久久久久久久久久久久久| 成人高潮视频无遮挡免费网站| 国产 一区精品| 国产黄片视频在线免费观看| 乱系列少妇在线播放| 不卡视频在线观看欧美| 国产男人的电影天堂91| 亚洲精品中文字幕在线视频 | 成人漫画全彩无遮挡| 直男gayav资源| 一区二区三区高清视频在线| 欧美成人a在线观看| 免费观看精品视频网站| 国产精品av视频在线免费观看| 亚洲av免费高清在线观看| 高清视频免费观看一区二区 | 免费看不卡的av| 亚洲欧洲国产日韩| 国产视频内射| 国产免费一级a男人的天堂| 一个人免费在线观看电影| av在线老鸭窝| 又大又黄又爽视频免费| 国产美女午夜福利| 日韩中字成人| 国产免费视频播放在线视频 | 久久国产乱子免费精品| 亚州av有码| 97热精品久久久久久| 亚洲久久久久久中文字幕| 亚洲丝袜综合中文字幕| 亚洲国产高清在线一区二区三| 亚洲av成人av| 亚洲人与动物交配视频| 一级a做视频免费观看| 亚洲国产精品成人久久小说| 国产亚洲精品av在线| 国产精品.久久久| 插阴视频在线观看视频| 久久99热这里只频精品6学生| 人妻系列 视频| 波多野结衣巨乳人妻| 男女啪啪激烈高潮av片| 国产亚洲最大av| 久久精品夜色国产| 亚洲av中文字字幕乱码综合| 久久热精品热| 国产成人福利小说| 一本一本综合久久| 成人欧美大片| 22中文网久久字幕| 一级毛片电影观看| 国产69精品久久久久777片| 伊人久久精品亚洲午夜| 哪个播放器可以免费观看大片| 国产精品久久久久久av不卡| 欧美激情国产日韩精品一区| 91久久精品国产一区二区三区| 欧美高清性xxxxhd video| 亚洲精品乱码久久久久久按摩| 91精品伊人久久大香线蕉| 99久久精品一区二区三区| 婷婷色综合大香蕉| 亚洲人与动物交配视频| 国内精品美女久久久久久| 欧美zozozo另类| 欧美潮喷喷水| 国产精品人妻久久久影院| 女的被弄到高潮叫床怎么办| 久久久午夜欧美精品| 春色校园在线视频观看| 久久久精品免费免费高清| 边亲边吃奶的免费视频| 精品久久久久久电影网| 日韩亚洲欧美综合| 国产在视频线在精品| 三级男女做爰猛烈吃奶摸视频| 哪个播放器可以免费观看大片| 最后的刺客免费高清国语| 精品一区二区三卡| 亚洲性久久影院| 免费播放大片免费观看视频在线观看| 国产成人a∨麻豆精品| 亚洲不卡免费看| 色网站视频免费| 亚洲一级一片aⅴ在线观看| 成年人午夜在线观看视频 | 欧美三级亚洲精品| 精品久久国产蜜桃| a级毛色黄片| 日本熟妇午夜| 日韩制服骚丝袜av| 最后的刺客免费高清国语| 午夜精品一区二区三区免费看| 久久精品人妻少妇| 搡老乐熟女国产| 久热久热在线精品观看| 免费无遮挡裸体视频| 国产精品国产三级国产av玫瑰| 最近的中文字幕免费完整| av福利片在线观看| 免费观看无遮挡的男女| av在线亚洲专区| 国产国拍精品亚洲av在线观看| 久久精品国产鲁丝片午夜精品| 亚洲欧美精品专区久久| 国产精品国产三级国产专区5o| 国产免费福利视频在线观看| 美女cb高潮喷水在线观看| 少妇熟女aⅴ在线视频| 国国产精品蜜臀av免费| 精品久久久噜噜| 99久久精品热视频| 欧美成人午夜免费资源| 美女主播在线视频| 麻豆av噜噜一区二区三区| 91久久精品电影网| 免费观看的影片在线观看| 看免费成人av毛片| 在线观看美女被高潮喷水网站| 成人午夜精彩视频在线观看| 国产成人精品婷婷| 如何舔出高潮| 九九在线视频观看精品| 日本wwww免费看| 国产精品久久视频播放| 一级爰片在线观看| 国产精品一区二区性色av| 亚洲婷婷狠狠爱综合网| 黄色日韩在线| 精品久久久久久久末码| 黄色日韩在线| 午夜福利视频精品| 免费无遮挡裸体视频| 搡老妇女老女人老熟妇| 不卡视频在线观看欧美| 嫩草影院新地址| 少妇的逼水好多| 2018国产大陆天天弄谢| 精品一区二区三卡| 久久久国产一区二区| 丝袜喷水一区| 国产欧美另类精品又又久久亚洲欧美| 日本一本二区三区精品| 日韩三级伦理在线观看| 国产日韩欧美在线精品| 亚洲熟妇中文字幕五十中出| 99热这里只有是精品50| 国产成年人精品一区二区| av在线亚洲专区| 国产午夜精品久久久久久一区二区三区| 成人av在线播放网站| 精品熟女少妇av免费看| 2021少妇久久久久久久久久久| 国产高清不卡午夜福利| 国模一区二区三区四区视频| 三级男女做爰猛烈吃奶摸视频| 91久久精品国产一区二区成人| 国产探花极品一区二区| 啦啦啦中文免费视频观看日本| 日韩成人伦理影院| 男女国产视频网站| av专区在线播放| 国产成人91sexporn| 少妇熟女欧美另类| 全区人妻精品视频| 性色avwww在线观看| 国产精品一区二区在线观看99 | 联通29元200g的流量卡| 一级a做视频免费观看| 日韩视频在线欧美| 日本猛色少妇xxxxx猛交久久| 亚洲人成网站高清观看| 777米奇影视久久| 边亲边吃奶的免费视频| 日韩欧美三级三区| 在线免费观看不下载黄p国产| 啦啦啦韩国在线观看视频| 亚洲av中文av极速乱| 欧美zozozo另类| 日韩电影二区| 亚洲成人中文字幕在线播放| 国语对白做爰xxxⅹ性视频网站| 国产成人午夜福利电影在线观看| 亚洲欧美成人综合另类久久久| 国产精品美女特级片免费视频播放器| 欧美人与善性xxx| 777米奇影视久久| 久久久久久久大尺度免费视频| 91久久精品电影网| 久久韩国三级中文字幕| 欧美潮喷喷水| 久久久成人免费电影| 非洲黑人性xxxx精品又粗又长| 美女脱内裤让男人舔精品视频| 激情五月婷婷亚洲| 国产探花在线观看一区二区| 久久久久久久午夜电影| 中国美白少妇内射xxxbb| 欧美性猛交╳xxx乱大交人| av专区在线播放| 欧美另类一区| 色5月婷婷丁香| 国产精品女同一区二区软件| 中文资源天堂在线| 久久精品久久精品一区二区三区| 精品酒店卫生间| 一二三四中文在线观看免费高清| 欧美不卡视频在线免费观看| 好男人视频免费观看在线| 欧美不卡视频在线免费观看| 国模一区二区三区四区视频| 国产免费福利视频在线观看| 亚洲欧美一区二区三区国产| 国产精品福利在线免费观看| 亚洲国产成人一精品久久久| 色综合亚洲欧美另类图片| 国产在线一区二区三区精| 久久99热这里只有精品18| av黄色大香蕉| 亚洲熟妇中文字幕五十中出| 小蜜桃在线观看免费完整版高清| 国产综合精华液| 天天一区二区日本电影三级| 69av精品久久久久久| 色综合色国产| 一本久久精品| 国产精品综合久久久久久久免费| 三级国产精品片| 欧美日韩在线观看h| 欧美3d第一页| 免费播放大片免费观看视频在线观看| 中文在线观看免费www的网站| 亚州av有码| 日韩电影二区| 成人漫画全彩无遮挡| 内射极品少妇av片p| 欧美zozozo另类| 国产男女超爽视频在线观看| 婷婷色综合www| 两个人的视频大全免费| 三级国产精品片| 美女大奶头视频| 熟妇人妻不卡中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 日韩人妻高清精品专区| 欧美3d第一页| 99久久精品热视频| 日韩视频在线欧美| 中国国产av一级| 天天一区二区日本电影三级| 能在线免费看毛片的网站| 亚洲第一区二区三区不卡| 成人av在线播放网站| 日日啪夜夜撸| 亚洲av不卡在线观看| 亚洲国产精品专区欧美| 亚洲在线自拍视频| 久久久久久国产a免费观看| 尾随美女入室| 毛片一级片免费看久久久久| 日本黄大片高清| 久久精品夜色国产| 婷婷色综合大香蕉| 久久久a久久爽久久v久久| 亚洲伊人久久精品综合| 久久久久久久大尺度免费视频| 2022亚洲国产成人精品| 插阴视频在线观看视频| kizo精华| 亚洲内射少妇av| 亚洲不卡免费看| 欧美日韩在线观看h| 麻豆精品久久久久久蜜桃| 亚洲av国产av综合av卡| 亚洲四区av| 天天一区二区日本电影三级| or卡值多少钱| 97超碰精品成人国产| 淫秽高清视频在线观看| 国产午夜精品一二区理论片| 久久久国产一区二区| 高清在线视频一区二区三区| 日本色播在线视频| 熟女人妻精品中文字幕| 午夜福利在线观看免费完整高清在| 午夜爱爱视频在线播放| 久久国产乱子免费精品| 在线观看人妻少妇| 色哟哟·www| 中国美白少妇内射xxxbb| 色播亚洲综合网| 亚洲欧美精品自产自拍| 天堂俺去俺来也www色官网 | 日韩一区二区三区影片| 久久这里只有精品中国| 午夜精品一区二区三区免费看| or卡值多少钱| 黄色日韩在线| 国产成人aa在线观看| 97人妻精品一区二区三区麻豆| av黄色大香蕉| 综合色av麻豆| 亚洲av二区三区四区| 成年人午夜在线观看视频 | 日韩 亚洲 欧美在线| 亚洲熟女精品中文字幕| 国产色爽女视频免费观看| 激情五月婷婷亚洲| 看免费成人av毛片| 啦啦啦中文免费视频观看日本| 在线观看一区二区三区| 国产69精品久久久久777片| 日本wwww免费看| 免费观看精品视频网站| 超碰av人人做人人爽久久| 老司机影院毛片| 午夜爱爱视频在线播放| eeuss影院久久| 国产成人精品一,二区| 听说在线观看完整版免费高清| 国产黄片美女视频| 能在线免费观看的黄片| 日韩中字成人| 大又大粗又爽又黄少妇毛片口| 美女xxoo啪啪120秒动态图| 老司机影院毛片| 国产成人aa在线观看| 五月伊人婷婷丁香| 国产色婷婷99| 久久精品久久精品一区二区三区| 亚洲欧美一区二区三区黑人 | 亚洲av男天堂| 免费av观看视频| 久久6这里有精品| 亚洲av日韩在线播放| 成人鲁丝片一二三区免费| 久久久久久伊人网av| 波多野结衣巨乳人妻| 在线观看免费高清a一片| 纵有疾风起免费观看全集完整版 | 亚洲av中文av极速乱| 内地一区二区视频在线| 久久精品综合一区二区三区| 国产成人精品福利久久| 综合色丁香网| 肉色欧美久久久久久久蜜桃 | 日韩电影二区| 成人特级av手机在线观看| 在线播放无遮挡| 综合色av麻豆| 国产午夜精品一二区理论片| 一级a做视频免费观看| 午夜免费男女啪啪视频观看| 欧美性猛交╳xxx乱大交人| 久久久久久久久久成人| 一区二区三区四区激情视频| 欧美人与善性xxx| 日韩av在线大香蕉| 国内精品一区二区在线观看| 色播亚洲综合网| 欧美区成人在线视频| 2021少妇久久久久久久久久久| 国产免费一级a男人的天堂| 国产精品一区www在线观看| 午夜日本视频在线| 欧美丝袜亚洲另类| 久久99精品国语久久久| 成人二区视频| 免费播放大片免费观看视频在线观看| 女的被弄到高潮叫床怎么办| 日韩不卡一区二区三区视频在线| 中文字幕免费在线视频6| 搡女人真爽免费视频火全软件| 日韩欧美国产在线观看| 欧美精品国产亚洲| 国产午夜精品一二区理论片| 成人漫画全彩无遮挡| 国产一区二区三区综合在线观看 | 国产精品.久久久| 水蜜桃什么品种好| 日韩制服骚丝袜av| 中文资源天堂在线| 亚洲四区av| 国产一区二区在线观看日韩| 亚洲av中文字字幕乱码综合| 卡戴珊不雅视频在线播放| av.在线天堂| 2021少妇久久久久久久久久久| 精品人妻偷拍中文字幕| 欧美3d第一页| 亚洲欧洲国产日韩| 亚洲在线自拍视频| av一本久久久久| 少妇熟女aⅴ在线视频| videossex国产| 寂寞人妻少妇视频99o| 久久久久精品性色| 国产av在哪里看| 99久国产av精品| 久久久久久伊人网av| 黄色欧美视频在线观看| 一级毛片电影观看| 十八禁网站网址无遮挡 | 九九爱精品视频在线观看| 七月丁香在线播放| 秋霞伦理黄片| 超碰av人人做人人爽久久| 伦理电影大哥的女人| 亚洲精品国产av成人精品| 91久久精品电影网| 成年免费大片在线观看| 精华霜和精华液先用哪个| 熟妇人妻不卡中文字幕| 99久久中文字幕三级久久日本| 亚洲三级黄色毛片| 噜噜噜噜噜久久久久久91| 伦精品一区二区三区| 一区二区三区四区激情视频| 国产老妇女一区| 午夜爱爱视频在线播放| 黄片无遮挡物在线观看| 蜜桃久久精品国产亚洲av| 亚洲自偷自拍三级| 欧美高清成人免费视频www| 视频中文字幕在线观看| av福利片在线观看| 成人性生交大片免费视频hd| 亚洲欧美精品专区久久| 日本三级黄在线观看| 久久人人爽人人爽人人片va| 免费av观看视频| 少妇猛男粗大的猛烈进出视频 | freevideosex欧美| 精品99又大又爽又粗少妇毛片| 在线免费观看的www视频| 中国美白少妇内射xxxbb| 成年人午夜在线观看视频 | 18禁动态无遮挡网站| 国产亚洲最大av| 国产免费福利视频在线观看| 国产久久久一区二区三区| 视频中文字幕在线观看| 久久久亚洲精品成人影院|