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

    Fractionally Delayed Kalman Filter

    2020-02-29 14:19:24AbhinoyKumarSingh
    IEEE/CAA Journal of Automatica Sinica 2020年1期

    Abhinoy Kumar Singh

    Abstract—The conventional Kalman filter is based on the assumption of non-delayed measurements. Several modifications appear to address this problem, but they are constrained by two crucial assumptions: 1) the delay is an integer multiple of the sampling interval, and 2) a stochastic model representing the relationship between delayed measurements and a sequence of possible non-delayed measurements is known. Practical problems often fail to satisfy these assumptions, leading to poor estimation accuracy and frequent track-failure. This paper introduces a new variant of the Kalman filter, which is free from the stochastic model requirement and addresses the problem of fractional delay.The proposed algorithm fixes the maximum delay (problem specific), which can be tuned by the practitioners for varying delay possibilities. A sequence of hypothetically defined intermediate instants characterizes fractional delays while maximum likelihood based delay identification could preclude the stochastic model requirement. Fractional delay realization could help in improving estimation accuracy. Moreover, precluding the need of a stochastic model could enhance the practical applicability. A comparative analysis with ordinary Kalman filter shows the high estimation accuracy of the proposed method in the presence of delay.

    I. INTRODUCTION

    THE conventional Kalman filter [1] assumes that measurements are non-delayed in time [2], [3], even though every practical system involves certain delays. If the delay is large,the measurement carries inappropriate information about the state dynamics. Subsequently, the state estimation with delayed measurement suffers from poor estimation accuracy and frequent track-failure.

    In many real-life problems measurements are time stamped and the measurement system has access to a clock, e.g. a network time protocol [4] based time distribution in large and diverse network systems, leading to a known delay. Such a delay can be modelled with a time-shift in measurements and the ordinary Kalman filter can be implemented. In the absence of time stamping and clock access, e.g. linear frequency modulator radar systems [5] where delay is in micro-seconds and a clock has such small resolution that it is infeasible in low budget applications, the delay remains unknown. An unknown and time varying delay is commonly known as random delay [6], [7], however early filtering literature cited such a delay with different names viz. out of sequence measurements (OOSM) [8], [9], filtering with random sample delay [10], [11], filtering with random time delayed measurements [12], etc. The filtering problems with randomly delayed measurements commonly appear in communication[13], control [14], target tracking [15], network systems [16],[17] etc.

    The literature on state estimation with delayed measurements begins with the work of Rayet al. [18] which involved using a reformulated minimum error variance estimator to deal with randomly delayed measurements. The estimation in this approach was based on sensor data statistics and the maximum delay was assumed to be one-step i.e. the measurement at an instant is either one-step delayed or non-delayed. Yaz and Ray modified this approach in [19] by introducing Grammiam assignment for stochastic parameters of the delayed model, and further modified it in [20] by introducing a linear matrix inequality based stochastic parameter estimation technique.Later, a series of developments [21]-[26] appeared in the literature to further improve filtering performance. In [21], an optimal estimation technique was developed for multiple delayed measurements by extending the ordinary Kalman filter to work in delayed environments. In later developments[22]-[24], the robust and adaptive Kalman filtering techniques were extended for delayed measurements to counter constraints such as model and process error uncertainties. Another popular method [25] was introduced by S. Sun, which modelled the unknown delay with a set of Bernoulli random variables and obtained an optimal estimation in terms of the Riccati difference equation and Lyapunov difference equation. In[26], the delayed measurement Kalman filtering algorithm was further extended for the multiplicative noises.

    The discussed filtering techniques for delayed measurements suffer from several drawbacks. To begin with,they are restricted to zero mean, white, Gaussian and uncorrelated noises, similar to the ordinary Kalman filter.Thankfully, some extensions of ordinary Kalman filter, e.g.[27] and [28], can be implemented with to mitigate these restrictions. Another drawback is that the above algorithms are sensitive to model validity, which is minimized up to some extent by implementing robust extension [29] of the ordinary Kalman filter. Another disadvantage is the discrete-time approximation of continuous process model which can be mitigated with continuous-time extension [30] of the ordinary Kalman filter. Apart from these drawbacks, there are certain limitations against which no proven solution is available in the literature. This paper introduces a solution for two of such limitations as mentioned below:

    1) Existing delay filters approximate the delay as an integer multiple of the sampling interval (Fig. 1(a)), though it is ideally fractional (Fig. 1(b)). Thus, poor estimation accuracy and frequent track-failures are observed if the sampling interval is sufficiently large.

    Fig. 1. Delay representation.

    2) Existing delay filters essentially require a stochastic model representing a stochastic relationship between the delayed measurement and a sequence of possible non-delayed measurements. Such models are formulated with a sequence of delay probabilities which are mostly unpredictable in practical problems. Therefore, existing delay filters are inapplicable to many real-life problems (where the probabilities are unknown), and have poor estimation accuracy (where an approximate stochastic model can be formulated with poorly estimated delay probabilities).

    To counter the aforementioned limitations, this paper proposes a novel modification of the Kalman filter which addresses the problem of fractional delay and does not require the aforementioned stochastic model. The proposed algorithm is named as fractionally delayed Kalman filter.

    In the case of fractionally delayed measurements, the ideal measurement time, which is defined as the arrival time for a measurement if there was no delay, can be in intermediate of previous sampling instants. Subsequently, the delay τdcan be a fractional multiple of sampling intervalT, besides being an integer multiple [18], [21], [25]. To characterize such a delay,a series of intermediate instants is defined hypothetically between two sampling instants. At timetk, the delayed instanttk-τdis one of these intermediate instants. Thus, the delay τdcould be estimated by defining a maximum likelihood based cost function and comparing it across all the intermediate instants (in the range of presumed maximum delay). The likelihood based delay identification could preclude the concerning stochastic model. Based on the τdestimate, the ideal measurement time foryk(measurement received attk)could be identified astk-τd. Subsequently,ykis used for state estimation attk-τdinstead oftk. In this regard, the covariance and optimal Kalman gain are re-derived for state estimation attk-τdusing the delayed measurementyk. Note that the filtering objective is to estimate the states at current instanttk(with measurementyk). Therefore, a sequential timeupdate approach, adopted from [31], [32], is implemented fromtk-τdtotk, and the desired estimates are obtained attk.

    The proposed algorithm is tested in an environment of randomly delayed measurements. The test analysis shows an improved estimation accuracy for the proposed method compared to the ordinary Kalman filter.

    II. KALMAN FILTER AND CHALLENGES WITH RANDOMLY DELAYED MEASUREMENTS

    The Kalman filter is implemented over the dynamic state space model which is composed of process and measurement equations, as

    where xk∈Rnandyk∈Rpare state and measurement vectors at thekth instant,k∈{0,1,2,...,N} withNbeing the number of samples,FkandHkare constants with appropriate dimension and,qkandvkare process and measurement noises, respectively. The process noise accounts for modeling error while the measurement noise compensates for the sensor error. The noises,qkandvk, are assumed to be white and Gaussian with zero mean and covarianceQkandRk, respectively.

    The Kalman filter is implemented in two steps: time update and measurement update.

    1). Time Update

    This step predicts the estimate and covariance of states, one step forward in time, as [2], [3]

    2). Measurement Update

    This step determines the posterior estimateand covariance (Pk|k) at thekth instant. Prior to computingand, the estimate of measurement (y?k|k-1) , innovationinnovation covariance (Sk) and an optimal Kalman gain (Kk)are computed as [2], [3]

    A. Challenges With Randomly Delayed Measurements

    As in (4), the state estimation accuracy depends on the accuracy of innovation,Due to delay,ykconsists of the state information associated withtk-τdwhere τd∈[0,tk). On the other hand,which is the predictedykand obtained by propagating xk-1|k-1through the state space model, consists of the state information attk. Thus, the state information comprised byykanddiffer in time by τd,resulting in an inappropriate computation ofThe state estimation with an inappropriateleads to poor estimation accuracy (equation (4)). The existing delay filters compensate for this difference up to some extent by presuming a stochastic model for delayed measurements;wherezkrepresents a possible non-delayed measurement attk, βiis a Bernoulli random variable with probabilityP(βi=1)=pi(also representsi-step delay probability) anddis the maximum possible delay at any instant. Based on this stochastic model,could be obtained by averaging the estimates ofzkat the current and delayed instantsThus, we getwhich provides a better match withykcompared to the traditional estimate

    Stochastic model based filtering has two major drawbacks:1) the stochastic model is based on the integer delay assumption, hence fractional delay is beyond the scope, and 2)the computation of estimateis subject to an accurate knowledge of probabilitiespi?i∈{1,2,...,d} which is very much unexpected in practical problems.

    The problem of fractional delay can be simplified to an integer delay if the delay is known and the sampling interval is adjustable. However the random delay (which is the interest of this paper) refers to an unknown and time-varying delay[6], [7], hence such simplification is overruled. Moreover, the sampling interval is often not under the control of practitioners, e.g., [15], causing another restriction on such simplification.

    Filtering with fractionally delayed measurements is challenging for following reasons: 1) in case of the stochastic model based filtering, the number of probabilities in the stochastic model is very large due to intermediate instants involved, causing higher uncertainty in the model, and 2)ordinary filters, including the existing delay filters, are designed to perform the filtering updates (time update and measurement update) at sampling instants, however filtering with a fractional delay essentially requires the updates at intermediate instants.

    III. FRACTIONALLY DELAYED KALMAN FILTER

    In this section, a modified Kalman filter is introduced to improve the estimation accuracy in filtering with randomly delayed measurements. The proposed algorithm determines the ideal measurement time forykusing a likelihood based approach. The ideal measurement time can be in intermediate of two sampling instants, given that the delay is fractional.Subsequently, the delay could be characterized by definingmintermediate instants between successive sampling instants.The time duration between two intermediate instants is δ=T/m(Tis the sampling interval). The intermediate instants are denoted asthus(asjbecomesm) representstk+1. In the remaining part of this paper, the filter parameters at intermediate instants are represented with following notations:

    with A being a hypothetical parameter.

    Remark 1: A highermoffers better precision in estimating ideal measurement time, but increases the computational burden. Thus,mshould be chosen as large as possible under the preassigned computational budget.

    The proposed algorithm consists of three steps:

    a) time update,

    b) ideal measurement time estimation,

    c) measurement update.

    A. Time Update

    The time update fromtk-1totkis performed through intermediate instants. Subsequently, the estimate and covariance at an intermediate instanttk-1+jδ could be obtained as

    B. Ideal Measurement Time Estimation

    Let us assume, the maximum delay does not exceeddsampling instants (i.e.,d×T). Let us defineDtas a set of time instants that consists of every possible ideal measurement time foryk, thenA simplified notation isδ} wherej∈{1,2,···,m}. Note thatdis a practitioner’s choice, hence the proposed algorithm is not restricted in terms of maximum delay. Let us assumetsis a past sampling instant within the range of the maximum delay, i.e.,ts∈{tk-d,thenandcan be determined for allj∈{1,2,...,m} using (6). Based on these quantities, the outcome of following Lemma offers a mathematical tool for estimating the ideal measurement time foryk.

    Lemma 1:Ifrepresents a negative Gaussian log likelihood thatykcorresponds to a delayed instantts+jδ, then

    Proof: The Kalman filter is based on Gaussian assumption of state vector xk.Thus,a linear state space model (which is the assumption of this development) insures a Gaussian distribution forykas well. Let us assume that the Gaussian distribution parameters, mean and covariance, ofykare μ and Σ, respectively, then

    whereP(·) represents probability density function (pdf) andp is the dimension ofyk.

    Let us denoteas the likelihood thatykcorresponds to a delayed instant (ts+jδ ). Ifholds true, thenykmust follow the estimate and covariance of a hypothetical measurement corresponding to (ts+jδ ), i.e.,andHence,

    The likelihood analysis is typically an arithmetic comparison,hence, a logarithmic simplification of likelihoods does not harm the generality. Thus, a simplified expression forcould be obtained as

    Following the similar argument,could be further simplified by removing the constant additive termand thereafter the constant multiplication factor 1/2, i.e.

    Note that x andyare linearly related (equation (2)), hence

    similar to (3).

    Substituting (9) into (8),

    is the likelihood ofykassociating with a delayed instant (ts+jδ), which can be obtained for every possible element ofDt. The element ofDtfor which the likelihood ofykis the largest is naturally the ideal measurement time.

    Note:represents a negative likelihood function,hence the maximum likelihood is obtained as the minimum of.

    Afterykis received atcould be computed for alls∈{k-d,k-d+1,...,k-1} and for allj∈{1,2,...,m}, i.e. for every element ofDt. Subsequently, if (ts*+j*δ) is the ideal measurement time foryk, then (s*,j*) could be obtained as

    C. Measurement Update

    The objective of this step is to construct the posterior estimate parameters attk, i.e.and Pk|k, as soon asykis received. As discussed above,ykis most likely to be arriving from (ts*+j*δ) (ideal measurement time) instead oftk. Hence,it is appropriate for computingandinstead ofand Pk|k. The proposed algorithm computesand Pk|kin two steps: i) based onyk,the posterior estimation is performed to obtainandat (ts*+j*δ), and ii)andare sequentially updated [31], [32] from (ts*+j*δ) totkto obtainand Pk|k.

    1) Posterior Estimation at(ts*+j*δ): Following the Kalman filtering strategy [2], the posterior estimation at(ts*+j*δ)requires an optimal Kalman gainat (ts*+j*δ). For the Kalman gainto be optimal, it must minimize the trace of error covariance matrix[2 ]. This paper computes the optimal Kalman gainthrough two subsequent lemmas. The first lemma formulates a generalized expression forwhile the second lemma determinesfor which the trace ofis minimum.innovation and innovation covariance at (ts*+j*δ) may be

    Before proceeding forward, the estimate of measurement,computed as

    Lemma2:If (ts*+j*δ) is the ideal measurement time foryki.e.then a generalized expression forcould be given as

    Proof: The posterior covariance at (ts*+j*δ) is

    similar to (4).

    Following the system dynamics,

    Substitutingykandfrom (15) and (11) in (14),

    Simplifying further,

    Asis uncorrelated,

    A further simplification to above equation gives

    Asand

    The representation of the covariance matrix in form of (12)is also known as Joseph’s form [3] of covariance matrix

    Lemma3:The optimal Kalman gainand corresponding error covarianceat intermediate instant (ts*+j*δ) are

    and

    Proof: As the optimal value ofminimizes the trace of[2],

    Expanding and further simplifying theexpression in (12),

    Substitutingfrom (11),

    Thus, the optimal Kalman gain is obtained as

    After the optimal Kalman gainis obtained, the posterior estimate at (ts*+j*δ) may be determined as

    At this st age, the posterior covariance at (ts*+j*δ) can be given as

    2) Computation ofandPk|k: Aswe can denoteasandas. It is to be noted that the objective of filtering is to estimate the states in real-time, i.e. to determineand Pk|kusingyk.Therefore,andare sequentially updated[31], [32] from (ts*+j*δ) totkto obtainand Pk|k.

    Let us assume thatmdis the number of intermediate instants between (ts*+j*δ) andtk; then,md=(k-s*-1)m+(m-j*).Asandcorrespond to (ts*+j*δ),mdtimes sequential update is required for obtainingand Pk|k. For the linear state dynamics,ith update could be obtained as

    This update step is performed recursively untilireachesmd.Oncei=mdis reached, the desired posterior estimate and covariance are obtained as

    A complete block-diagram for implementing the proposed method is provided in Fig. 2.

    The filter design technique introduced in this section does not use a stochastic model and accounts for fractional delay,leading to higher estimation accuracy. On the other hand, it requires the storage of the estimated data from delayed instants, hence requires a larger storage budget compared to the ordinary Kalman filter. The general assumptions involving noise in Kalman filter design, such as white, zero mean Gaussian and uncorrelated noises, are imposed over the proposed algorithm as well.

    IV. SIMULATION

    In this section, the proposed method is implemented for a filtering problem with randomly delayed measurements. The performance of the proposed method is compared with ordinary Kalman filter in terms of root mean square error(RMSE).

    A. Dynamic State-Space Model

    The system under consideration follows a linear state dynamics, given as

    Fig. 2. Flow chart for implementing the proposed algorithm.

    whereqkis process noise. The process noise is assumed to be normally distributed with zero mean and covariance

    whereT=0.5 is sampling interval and c=10×10-6is constant.

    With ideal behaviour (with no delay in measurements), the measurement system observes the first state as a measurement at every sampling interval. Thus, the expected non-delayed measurement could be modelled as

    whereHk=[1 0], andvkis the measurement noise component which is approximated as Gaussian with zero mean and covarianceRk=1.

    B. Delayed Measurement Simulation

    With the lack of real-life data, a simulated dataset has been generated in Matlab to characterize the randomly delayed measurements. The data simulation is based on a presumed stochastic model for delayed measurements, in a dynamic system represented by (24) and (25). The maximum possible delay is assumed to bedsampling intervals. The fractional delay is characterized by assigningmintermediate instants between two sampling instants. Subsequently, the measurementykreceived attkcorresponds to an instant inDtdefined in Section III-B. Thus, the delayed measurementykcould be modelled using a set of Bernoulli random variablesand ?j∈{1,2,...,m} i.e.Note thatindicates thatts+jδ is the ideal measurement time foryk. Thus,ykcould be modelled as

    Note thatis a hypothetical non-delayed measurement atts+jδ.

    During the data simulation, it was assumed that multiple measurements are not received at any instant, hence β(s,j)can be 1 for only one pair ofsandjattk. To be noted, it is possible that?s∈{k-d,k-d+1,...,k-1} and ?j∈{1,2,...,m} attk, demonstrating the case where there is a missing measurement, which is observed if the delay has a sudden increase.

    It is understood in above discussion thatp(k-1,m)is the probability of receiving a non-delayed measurement. To simulate the measurement data,p(k-1,m)was taken as 0.3 and the remaining probability was equally distributed overdsampling intervals (including the intermediate instants).

    It should be noted that the delay measurement model (26) is used only for generating simulated measurement data, with no role in the filtering algorithm.

    C. Implementation and Results

    The proposed algorithm is implemented over 500 time-steps(k∈{1,2,...,500}) to estimate the states of the dynamic system represented by (24) and (25). For comparison purposes, the simulated data for the states were generated using (24) with initial state being x0=[0 0]Twhere T denotes transpose of matrix. The delayed measurement data was simulated using the stochastic model given in (26). For filter implementation, the initial estimatewas assumed to follow the Gaussian distribution with mean x0and initial covariance P0|0=diag([1 1]).

    1) Estimation Accuracy: To analyze the estimation accuracy, the estimateand the true state xk?k∈{1,2,...,500}are plotted for individual states in Fig. 3. An appreciable match between the true state and the estimated state reflects a high estimation accuracy for the proposed method. It also reflects the ability of the proposed method in dealing with fractionally delayed measurements.

    Fig. 3. True and estimated states.

    2)Comparison with Ordinary Kalman Filter: The comparison between the proposed algorithm and the ordinary Kalman filter is based on RSME analysis. The RMSE is computed by implementing several independent Monte-Carlo runs. At thekth step, the RMSE forith state is determined as

    whereM=1000 is the number of Monte-Carlo runs.

    For most of the practical systems, the delay is not muc h larger than the sampling interval [33], hence, the analysis has been limited to the case where the maximum delay does not exceed 5 sampling intervals. For the case of maximum delay not exceeding 1 and 2 sampling intervals, the RMSE for the proposed method and the ordinary Kalman filter are plotted in Figs. 4 and 5, respectively. In these figures, the Kalman filter is abbreviated as KF. For higher delays, the average RMSE(between 100 to 500 steps) is compared in Table I and II. The proposed method has a convergence time of 100 steps, hence the comparative analysis begins after this duration. The tables also analyze the relative reduction in RMSE for the proposed method compared to the ordinary Kalman filter.

    The plots and tables indicate a significantly reduced RMSE for the proposed method compared to the ordinary Kalman filter. It concludes a better estimation accuracy for the proposed algorithm in the presence of delay.

    Fig. 4. RMSE plots: maximum delay being 1 sampling interval.

    Fig. 5. RMSE plots: maximum delay being 2 sampling intervals.

    TABLE IAVERAGE RMSE OF STATE 1

    TABLE II AVERAGE RMSE OF STATE 2

    To analyze the rise in computational time, the execution time was studied for the proposed algorithm as well as the ordinary Kalman filter by implementing them over 500 timesteps. The implementation was performed in the Matlab environment on a personal computer with a 64-bit operating system, 4GB RAM, and 2GHz intel core i3 processor. The execution time (s) for the proposed method was 0.2400,0.4740, 0.6705, 0.8689 and 1.0785 for maximum delay varying from 1 to 5 sampling intervals, respectively.However, the execution time was constant around 0.02 seconds for all delay cases using the ordinary Kalman filter.

    V. DISCUSSION AND CONCLUSION

    Time delay is a common phenomena in physical measurement systems, though it may be sufficiently small to ignore in several cases. The Kalman filter in its traditional form ignores the measurement delay, and suffers from poor estimation accuracy and frequent track-failure if the delay is sufficiently large. The existing variants of the Kalman filter for addressing measurement delay are constrained with two major assumptions: 1) the delay is an integer multiple of sampling interval, and 2) a precise stochastic model is known for delayed measurement (in terms of a sequence of possible non-delayed measurements). The practical problems usually do not fulfill these assumptions, leading to poor estimation accuracy and often track-failure. A new variant of the Kalman filter is proposed in this paper to address the fractional delay and remove the constraints of the stochastic model.Addressing the fractional delay is advantageous in improving estimation accuracy, while not using the stochastic model could enhance practical applicability in addition to improving the estimation accuracy. The simulation results report a reduced RMSE for the proposed method compared to the ordinary Kalman filter, demonstrating an improved estimation accuracy in presence of delay. On the other hand, the computational time for the proposed method is reported to be larger than the ordinary Kalman filter. The fundamental assumptions on noises, including zero mean, white, Gaussian and uncorrelated, when designing an ordinary kalman filter,are imposed with the proposed method as well. The future extension of the proposed algorithm has multiple scopes, e.g.improving estimation accuracy further, reducing computational time, and removing the constraints on noisesetc.

    狠狠精品人妻久久久久久综合| 国产成人91sexporn| 男女边摸边吃奶| 欧美精品人与动牲交sv欧美| 中文天堂在线官网| 午夜福利网站1000一区二区三区| 成人毛片60女人毛片免费| 成人毛片a级毛片在线播放| 交换朋友夫妻互换小说| 2022亚洲国产成人精品| 成年美女黄网站色视频大全免费 | 国产成人免费观看mmmm| 一区二区av电影网| av又黄又爽大尺度在线免费看| 国产成人精品婷婷| 多毛熟女@视频| 精品久久久久久久久av| 日本猛色少妇xxxxx猛交久久| 国产亚洲最大av| 日韩三级伦理在线观看| 七月丁香在线播放| 啦啦啦在线观看免费高清www| 国产精品免费大片| 一级毛片久久久久久久久女| 啦啦啦中文免费视频观看日本| 久久人人爽av亚洲精品天堂 | 麻豆成人av视频| 国产高清不卡午夜福利| 1000部很黄的大片| 一区二区三区精品91| 亚洲av.av天堂| 97在线视频观看| 欧美一区二区亚洲| 亚洲精品久久午夜乱码| 国产成人a∨麻豆精品| 丰满少妇做爰视频| 国内精品宾馆在线| 国产精品久久久久成人av| 日日摸夜夜添夜夜爱| 女性被躁到高潮视频| 久久精品久久久久久久性| 麻豆成人午夜福利视频| 性色avwww在线观看| 丝袜脚勾引网站| 婷婷色麻豆天堂久久| 最黄视频免费看| 日本vs欧美在线观看视频 | 在线观看一区二区三区| 日韩三级伦理在线观看| 97在线视频观看| 国产亚洲欧美精品永久| 18禁在线播放成人免费| 王馨瑶露胸无遮挡在线观看| 欧美精品亚洲一区二区| 欧美xxxx性猛交bbbb| 如何舔出高潮| 免费少妇av软件| 亚洲,欧美,日韩| 人人妻人人澡人人爽人人夜夜| 日韩欧美 国产精品| 国产精品成人在线| 亚洲国产av新网站| 免费久久久久久久精品成人欧美视频 | 久久久久性生活片| 欧美日韩亚洲高清精品| 亚洲国产精品国产精品| 少妇精品久久久久久久| 一级毛片aaaaaa免费看小| 免费看光身美女| 天堂中文最新版在线下载| 久久国产乱子免费精品| 免费久久久久久久精品成人欧美视频 | 婷婷色综合大香蕉| xxx大片免费视频| 一区二区三区四区激情视频| 国精品久久久久久国模美| 国产精品久久久久成人av| 欧美3d第一页| 国产成人精品一,二区| 亚洲国产日韩一区二区| 国产女主播在线喷水免费视频网站| 久久久久精品久久久久真实原创| 一边亲一边摸免费视频| 亚洲高清免费不卡视频| 青春草国产在线视频| 亚洲第一av免费看| 人人妻人人澡人人爽人人夜夜| 亚洲av中文av极速乱| 搡女人真爽免费视频火全软件| 久久久精品94久久精品| 亚洲va在线va天堂va国产| 欧美 日韩 精品 国产| 久久精品熟女亚洲av麻豆精品| 麻豆成人av视频| 国产精品久久久久久久电影| 亚洲精品,欧美精品| 国产精品人妻久久久影院| 天堂8中文在线网| 五月天丁香电影| 亚洲人成网站在线播| 久久久欧美国产精品| 国产熟女欧美一区二区| 国产黄频视频在线观看| 爱豆传媒免费全集在线观看| 免费观看性生交大片5| 久久久久网色| 亚洲人成网站在线播| 亚洲av中文字字幕乱码综合| 黄色一级大片看看| 婷婷色麻豆天堂久久| 一本—道久久a久久精品蜜桃钙片| 中国国产av一级| 久久国内精品自在自线图片| 久久精品久久久久久久性| 久久久a久久爽久久v久久| 久久av网站| 国产爽快片一区二区三区| 99久久人妻综合| 国产精品国产av在线观看| 女性生殖器流出的白浆| 男人舔奶头视频| 亚洲av.av天堂| 亚洲精品久久午夜乱码| a 毛片基地| 97在线人人人人妻| 久久av网站| 我的老师免费观看完整版| 精品国产乱码久久久久久小说| 国国产精品蜜臀av免费| 久久 成人 亚洲| 欧美xxxx性猛交bbbb| av在线老鸭窝| 看免费成人av毛片| 午夜福利网站1000一区二区三区| 街头女战士在线观看网站| 又粗又硬又长又爽又黄的视频| 在线观看免费高清a一片| 亚洲人与动物交配视频| av线在线观看网站| 亚洲性久久影院| 亚洲av日韩在线播放| 成人国产麻豆网| 一级毛片电影观看| 国产在视频线精品| 国产v大片淫在线免费观看| 欧美一级a爱片免费观看看| 精品少妇黑人巨大在线播放| 妹子高潮喷水视频| 亚洲国产毛片av蜜桃av| 日本黄色日本黄色录像| av在线蜜桃| 日本vs欧美在线观看视频 | 久久人人爽人人爽人人片va| 毛片一级片免费看久久久久| 最近的中文字幕免费完整| 最近最新中文字幕大全电影3| 人妻制服诱惑在线中文字幕| 精品人妻偷拍中文字幕| 99热网站在线观看| 女人十人毛片免费观看3o分钟| 国产黄色视频一区二区在线观看| 国产 一区 欧美 日韩| 在线精品无人区一区二区三 | 热99国产精品久久久久久7| 水蜜桃什么品种好| 超碰97精品在线观看| 日韩一区二区视频免费看| 久久人妻熟女aⅴ| 日韩不卡一区二区三区视频在线| 亚洲国产欧美在线一区| 美女福利国产在线 | 91精品一卡2卡3卡4卡| 亚洲第一区二区三区不卡| 26uuu在线亚洲综合色| 亚洲精品国产成人久久av| 国产亚洲精品久久久com| 99久久人妻综合| 伊人久久精品亚洲午夜| 只有这里有精品99| 自拍欧美九色日韩亚洲蝌蚪91 | 高清欧美精品videossex| 久久精品国产a三级三级三级| 日韩国内少妇激情av| 久久久午夜欧美精品| 亚洲av福利一区| 美女内射精品一级片tv| 国产精品久久久久久精品古装| 午夜免费男女啪啪视频观看| 日本免费在线观看一区| 在线观看一区二区三区| 一级毛片久久久久久久久女| 免费在线观看成人毛片| 成人综合一区亚洲| 亚洲电影在线观看av| 草草在线视频免费看| 亚洲国产最新在线播放| 熟女人妻精品中文字幕| 老女人水多毛片| 黄色日韩在线| 亚洲天堂av无毛| 欧美一区二区亚洲| 久久99热6这里只有精品| 日韩欧美一区视频在线观看 | 91久久精品电影网| 五月玫瑰六月丁香| 午夜福利影视在线免费观看| 亚洲人成网站在线播| 高清视频免费观看一区二区| 亚洲av欧美aⅴ国产| 国产一区亚洲一区在线观看| 97热精品久久久久久| 亚洲图色成人| 国产白丝娇喘喷水9色精品| 少妇人妻一区二区三区视频| 欧美性感艳星| a级毛色黄片| 免费看不卡的av| 欧美精品一区二区大全| 国产男人的电影天堂91| 日本一二三区视频观看| 国产91av在线免费观看| 十八禁网站网址无遮挡 | 国产极品天堂在线| 伦理电影免费视频| 水蜜桃什么品种好| 人人妻人人看人人澡| 在线观看av片永久免费下载| 日本猛色少妇xxxxx猛交久久| 久久久久网色| 大又大粗又爽又黄少妇毛片口| 欧美三级亚洲精品| 亚洲综合色惰| 一区二区三区乱码不卡18| 乱码一卡2卡4卡精品| 日本免费在线观看一区| 在线观看三级黄色| 久久久色成人| 日本欧美视频一区| 三级国产精品欧美在线观看| 国产精品免费大片| 国产黄色视频一区二区在线观看| 欧美高清成人免费视频www| 婷婷色麻豆天堂久久| 超碰av人人做人人爽久久| 热re99久久精品国产66热6| 日韩,欧美,国产一区二区三区| 日韩精品有码人妻一区| 在线观看免费高清a一片| 尤物成人国产欧美一区二区三区| av.在线天堂| 22中文网久久字幕| 国产欧美另类精品又又久久亚洲欧美| 爱豆传媒免费全集在线观看| 男女边吃奶边做爰视频| 久久国内精品自在自线图片| 成人综合一区亚洲| 亚洲精品乱码久久久久久按摩| 男女国产视频网站| 日本黄色片子视频| 五月开心婷婷网| 国产伦精品一区二区三区视频9| 少妇猛男粗大的猛烈进出视频| 国产久久久一区二区三区| 亚洲成人av在线免费| 最近最新中文字幕免费大全7| 午夜福利影视在线免费观看| 王馨瑶露胸无遮挡在线观看| 久久精品国产鲁丝片午夜精品| 亚洲电影在线观看av| 国产精品人妻久久久影院| 国产中年淑女户外野战色| 色5月婷婷丁香| 美女主播在线视频| 国产精品精品国产色婷婷| 亚洲精品成人av观看孕妇| 亚洲国产日韩一区二区| 亚洲av.av天堂| 成年人午夜在线观看视频| 午夜视频国产福利| 欧美日韩在线观看h| 国产高清有码在线观看视频| 亚洲国产精品国产精品| 成人一区二区视频在线观看| 国产爱豆传媒在线观看| 国产综合精华液| 久久精品国产亚洲av天美| 高清午夜精品一区二区三区| 国产精品欧美亚洲77777| 亚洲成人一二三区av| 网址你懂的国产日韩在线| 七月丁香在线播放| 97精品久久久久久久久久精品| 亚洲国产精品999| 欧美bdsm另类| 大片免费播放器 马上看| 在线观看一区二区三区激情| 日韩av不卡免费在线播放| 偷拍熟女少妇极品色| 熟妇人妻不卡中文字幕| 国产av码专区亚洲av| 欧美国产精品一级二级三级 | 99久久综合免费| 免费人成在线观看视频色| 亚洲一区二区三区欧美精品| 日本av手机在线免费观看| 蜜桃久久精品国产亚洲av| 天堂中文最新版在线下载| 丰满乱子伦码专区| 99热6这里只有精品| 中文资源天堂在线| 国产午夜精品久久久久久一区二区三区| 少妇裸体淫交视频免费看高清| 久久久久性生活片| 人妻制服诱惑在线中文字幕| 亚洲伊人久久精品综合| 久久久久精品性色| 欧美日本视频| 少妇被粗大猛烈的视频| 最后的刺客免费高清国语| 久久久精品免费免费高清| 亚洲天堂av无毛| 我要看黄色一级片免费的| 国产成人a区在线观看| 午夜福利在线在线| 久久青草综合色| 丝袜喷水一区| 日韩av在线免费看完整版不卡| 天美传媒精品一区二区| 国产深夜福利视频在线观看| 亚洲精品aⅴ在线观看| 久久久久久久精品精品| 极品少妇高潮喷水抽搐| 日韩一本色道免费dvd| 女人十人毛片免费观看3o分钟| 国产成人精品一,二区| 日本与韩国留学比较| 97热精品久久久久久| 美女中出高潮动态图| 哪个播放器可以免费观看大片| 国产免费视频播放在线视频| 大又大粗又爽又黄少妇毛片口| 精品人妻视频免费看| 午夜福利影视在线免费观看| 亚洲精品一二三| 99热网站在线观看| 免费在线观看成人毛片| 国产熟女欧美一区二区| 97在线视频观看| 成年av动漫网址| 国产高潮美女av| 男男h啪啪无遮挡| 免费少妇av软件| 精品人妻视频免费看| 婷婷色综合www| 日日摸夜夜添夜夜添av毛片| 亚洲欧美成人综合另类久久久| a 毛片基地| 亚洲色图av天堂| 色哟哟·www| 中文字幕免费在线视频6| 精品国产露脸久久av麻豆| 菩萨蛮人人尽说江南好唐韦庄| 熟女人妻精品中文字幕| 亚洲av二区三区四区| 黑人猛操日本美女一级片| 老熟女久久久| 午夜日本视频在线| 欧美xxxx黑人xx丫x性爽| 久久这里有精品视频免费| 久久久亚洲精品成人影院| 久久韩国三级中文字幕| 国产伦精品一区二区三区四那| 国产一区二区三区av在线| 91久久精品国产一区二区成人| 午夜免费观看性视频| 亚洲国产av新网站| 人妻系列 视频| 欧美日韩一区二区视频在线观看视频在线| 一级毛片久久久久久久久女| 国产成人精品一,二区| 欧美日韩视频精品一区| 亚洲欧美中文字幕日韩二区| 国产精品.久久久| 精品久久国产蜜桃| 99久久人妻综合| 只有这里有精品99| 99久久精品国产国产毛片| 精品少妇黑人巨大在线播放| 亚洲国产欧美在线一区| 中文字幕制服av| av网站免费在线观看视频| 一区二区三区四区激情视频| 国产精品国产av在线观看| 成人国产av品久久久| 91精品伊人久久大香线蕉| h视频一区二区三区| 精品国产三级普通话版| 国产成人aa在线观看| 久久久久久久久久人人人人人人| 我要看日韩黄色一级片| 老师上课跳d突然被开到最大视频| 亚洲真实伦在线观看| 久久久久网色| 亚洲av二区三区四区| 亚洲av电影在线观看一区二区三区| 日日啪夜夜爽| av天堂中文字幕网| 一本色道久久久久久精品综合| 国产精品久久久久久精品古装| 国产免费一级a男人的天堂| 午夜福利视频精品| 蜜桃在线观看..| 日韩三级伦理在线观看| 纯流量卡能插随身wifi吗| 亚洲av中文av极速乱| 亚洲va在线va天堂va国产| 免费看光身美女| 精品一区在线观看国产| 一个人免费看片子| 欧美日韩精品成人综合77777| 亚洲精品第二区| 一级毛片黄色毛片免费观看视频| 嫩草影院新地址| 少妇的逼水好多| 日韩大片免费观看网站| 中文字幕制服av| 欧美亚洲 丝袜 人妻 在线| 女性被躁到高潮视频| 国产高清不卡午夜福利| 国产精品女同一区二区软件| 欧美97在线视频| 久久99热这里只有精品18| 亚洲精品中文字幕在线视频 | 另类亚洲欧美激情| 免费不卡的大黄色大毛片视频在线观看| 欧美 日韩 精品 国产| 午夜日本视频在线| 99久国产av精品国产电影| 免费av不卡在线播放| 精品少妇久久久久久888优播| 久久毛片免费看一区二区三区| 黑人高潮一二区| 熟妇人妻不卡中文字幕| 日本与韩国留学比较| 国产黄色视频一区二区在线观看| 免费观看在线日韩| 亚洲成人手机| 五月开心婷婷网| 91在线精品国自产拍蜜月| 午夜免费男女啪啪视频观看| 成人特级av手机在线观看| 观看免费一级毛片| 草草在线视频免费看| 香蕉精品网在线| 欧美精品一区二区免费开放| 欧美丝袜亚洲另类| .国产精品久久| 亚洲综合色惰| 中文字幕亚洲精品专区| 成人毛片a级毛片在线播放| 免费观看的影片在线观看| 久久女婷五月综合色啪小说| 最近中文字幕高清免费大全6| 亚洲精品一二三| 日本猛色少妇xxxxx猛交久久| 伊人久久国产一区二区| 看非洲黑人一级黄片| 国产中年淑女户外野战色| 免费观看在线日韩| 日韩欧美精品免费久久| 岛国毛片在线播放| 综合色丁香网| 色5月婷婷丁香| 久久久久人妻精品一区果冻| 午夜免费观看性视频| 日韩av免费高清视频| 欧美一级a爱片免费观看看| 97超视频在线观看视频| 一个人免费看片子| 国产精品一及| 又爽又黄a免费视频| 久久精品久久久久久噜噜老黄| 嘟嘟电影网在线观看| 国产精品福利在线免费观看| 国产视频内射| 国产黄色免费在线视频| 久久久久国产精品人妻一区二区| 两个人的视频大全免费| 国产精品熟女久久久久浪| av.在线天堂| 精品久久久精品久久久| 人人妻人人澡人人爽人人夜夜| 啦啦啦视频在线资源免费观看| 十八禁网站网址无遮挡 | 高清不卡的av网站| 亚洲久久久国产精品| 欧美一区二区亚洲| 日韩成人伦理影院| 狂野欧美激情性xxxx在线观看| 日韩一区二区视频免费看| 麻豆国产97在线/欧美| 观看免费一级毛片| 国产精品一及| 老女人水多毛片| av免费在线看不卡| 国产久久久一区二区三区| 日韩精品有码人妻一区| 亚洲人成网站在线播| 国内揄拍国产精品人妻在线| 91精品国产九色| 婷婷色综合大香蕉| 欧美 日韩 精品 国产| 日本色播在线视频| 国产精品久久久久久精品古装| 欧美日韩在线观看h| 成年免费大片在线观看| 免费人妻精品一区二区三区视频| 久久久久国产精品人妻一区二区| 亚洲av欧美aⅴ国产| av女优亚洲男人天堂| 亚洲av日韩在线播放| 亚洲国产精品成人久久小说| 欧美精品人与动牲交sv欧美| 亚洲欧美一区二区三区国产| 久久久午夜欧美精品| 建设人人有责人人尽责人人享有的 | 日本黄大片高清| 亚洲四区av| 亚洲电影在线观看av| 国产精品久久久久成人av| 午夜福利在线观看免费完整高清在| av不卡在线播放| 身体一侧抽搐| 亚洲精品日本国产第一区| 国产女主播在线喷水免费视频网站| 热99国产精品久久久久久7| 亚洲电影在线观看av| 久久久色成人| 日本黄大片高清| 欧美精品亚洲一区二区| 日韩av不卡免费在线播放| 狂野欧美白嫩少妇大欣赏| 国产极品天堂在线| 欧美3d第一页| 在线观看一区二区三区| 午夜精品国产一区二区电影| 久久国产乱子免费精品| 久久久久久久国产电影| 美女国产视频在线观看| 国产高清三级在线| 亚洲国产欧美在线一区| 国产伦精品一区二区三区四那| av一本久久久久| 小蜜桃在线观看免费完整版高清| 免费av中文字幕在线| 久久亚洲国产成人精品v| 人人妻人人添人人爽欧美一区卜 | 精品久久久久久久久亚洲| 欧美极品一区二区三区四区| 菩萨蛮人人尽说江南好唐韦庄| 在线观看人妻少妇| 亚洲内射少妇av| 欧美高清成人免费视频www| 久久综合国产亚洲精品| 色综合色国产| 97热精品久久久久久| 成年美女黄网站色视频大全免费 | 亚洲欧洲国产日韩| 中文字幕制服av| 亚洲欧美成人综合另类久久久| 日日啪夜夜爽| 大码成人一级视频| 小蜜桃在线观看免费完整版高清| 久久人人爽av亚洲精品天堂 | 综合色丁香网| 国产精品.久久久| 婷婷色综合www| 中文字幕av成人在线电影| 久久人人爽人人片av| 91aial.com中文字幕在线观看| 国产v大片淫在线免费观看| 国产一区二区在线观看日韩| 日日啪夜夜爽| 国产亚洲欧美精品永久| kizo精华| 一二三四中文在线观看免费高清| 精品国产一区二区三区久久久樱花 | 日韩免费高清中文字幕av| 两个人的视频大全免费| 99热6这里只有精品| av免费观看日本| 精品亚洲乱码少妇综合久久| 各种免费的搞黄视频| 午夜免费男女啪啪视频观看| 日本爱情动作片www.在线观看| .国产精品久久| 国产淫片久久久久久久久| 只有这里有精品99| 高清毛片免费看| 大片电影免费在线观看免费| 午夜免费男女啪啪视频观看| 国产av国产精品国产| 免费观看av网站的网址| 久久久久国产精品人妻一区二区| 观看av在线不卡| 国产精品麻豆人妻色哟哟久久| 国产高清国产精品国产三级 | 嫩草影院入口| 久久久久久九九精品二区国产| 久热久热在线精品观看| 成人国产av品久久久| 黄色欧美视频在线观看| 中国国产av一级| 波野结衣二区三区在线| 国产片特级美女逼逼视频|