• <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.

    亚洲国产毛片av蜜桃av| 777米奇影视久久| 国产99久久九九免费精品| 国产精品免费视频内射| 成人午夜精彩视频在线观看| 日韩av在线免费看完整版不卡| 少妇的丰满在线观看| 国产一区二区 视频在线| 精品高清国产在线一区| 成人亚洲精品一区在线观看| 国产成人一区二区在线| 亚洲精品久久成人aⅴ小说| 国产欧美日韩一区二区三 | 777久久人妻少妇嫩草av网站| 又紧又爽又黄一区二区| 飞空精品影院首页| 国产在线一区二区三区精| 巨乳人妻的诱惑在线观看| 国产又色又爽无遮挡免| 久久精品人人爽人人爽视色| 久久99精品国语久久久| 精品一区在线观看国产| 妹子高潮喷水视频| 日本黄色日本黄色录像| 亚洲av成人精品一二三区| 我的亚洲天堂| 亚洲专区中文字幕在线| 99国产综合亚洲精品| 一区二区av电影网| 老司机深夜福利视频在线观看 | 91精品三级在线观看| 亚洲五月婷婷丁香| 国产1区2区3区精品| 女人被躁到高潮嗷嗷叫费观| 天天添夜夜摸| 国产一区亚洲一区在线观看| 精品福利观看| 深夜精品福利| 女警被强在线播放| 99九九在线精品视频| 亚洲精品国产av蜜桃| 亚洲成av片中文字幕在线观看| 青草久久国产| 国产成人av激情在线播放| 国产日韩欧美视频二区| 亚洲成av片中文字幕在线观看| 中国美女看黄片| 久久久久精品国产欧美久久久 | 国产男女超爽视频在线观看| 免费在线观看黄色视频的| 脱女人内裤的视频| 亚洲国产精品国产精品| 高清黄色对白视频在线免费看| 成年女人毛片免费观看观看9 | 亚洲精品美女久久久久99蜜臀 | 少妇精品久久久久久久| 在线看a的网站| 狂野欧美激情性xxxx| 又黄又粗又硬又大视频| 午夜福利乱码中文字幕| 这个男人来自地球电影免费观看| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品国产精品| 精品久久久久久久毛片微露脸 | 天天操日日干夜夜撸| 国产精品一区二区在线观看99| 人人妻人人爽人人添夜夜欢视频| 亚洲精品美女久久久久99蜜臀 | 一级毛片我不卡| 一边摸一边做爽爽视频免费| avwww免费| 久久人妻熟女aⅴ| 十八禁人妻一区二区| 女人精品久久久久毛片| 亚洲伊人久久精品综合| 啦啦啦在线免费观看视频4| 蜜桃国产av成人99| 国产淫语在线视频| 久久鲁丝午夜福利片| 久久精品国产亚洲av高清一级| 国产精品偷伦视频观看了| 在线av久久热| 校园人妻丝袜中文字幕| 夜夜骑夜夜射夜夜干| 国产无遮挡羞羞视频在线观看| 欧美97在线视频| 国产精品.久久久| 日韩大码丰满熟妇| 日韩大片免费观看网站| 精品少妇久久久久久888优播| www日本在线高清视频| 无遮挡黄片免费观看| 亚洲欧美精品自产自拍| 欧美久久黑人一区二区| 久久精品国产a三级三级三级| 叶爱在线成人免费视频播放| 精品一区在线观看国产| 在线亚洲精品国产二区图片欧美| 免费观看人在逋| av网站在线播放免费| 亚洲专区国产一区二区| 性高湖久久久久久久久免费观看| 精品欧美一区二区三区在线| 一级黄色大片毛片| 国产亚洲欧美在线一区二区| 777米奇影视久久| 国产精品免费视频内射| 美女高潮到喷水免费观看| 99久久精品国产亚洲精品| 99国产精品一区二区三区| 精品一品国产午夜福利视频| 在线观看免费日韩欧美大片| 国产免费视频播放在线视频| 国产精品香港三级国产av潘金莲 | 国产成人欧美| 男女之事视频高清在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 美女大奶头黄色视频| 久久午夜综合久久蜜桃| av天堂在线播放| 欧美xxⅹ黑人| 亚洲欧美日韩高清在线视频 | 999精品在线视频| 蜜桃在线观看..| 国产在线免费精品| 久热爱精品视频在线9| 在线亚洲精品国产二区图片欧美| 国产免费又黄又爽又色| 亚洲成色77777| 国产一区二区三区av在线| 欧美中文综合在线视频| 欧美 亚洲 国产 日韩一| 欧美+亚洲+日韩+国产| 男女床上黄色一级片免费看| 美女大奶头黄色视频| 久久久久久亚洲精品国产蜜桃av| 亚洲,一卡二卡三卡| 亚洲国产精品成人久久小说| 男女午夜视频在线观看| 欧美激情极品国产一区二区三区| 一级毛片 在线播放| 亚洲精品日本国产第一区| 18禁国产床啪视频网站| 成人国产av品久久久| 亚洲av成人精品一二三区| 高清欧美精品videossex| 久久亚洲精品不卡| 久久这里只有精品19| 国产老妇伦熟女老妇高清| 熟女av电影| 国产成人av教育| 女性被躁到高潮视频| 欧美变态另类bdsm刘玥| 国产97色在线日韩免费| 天堂俺去俺来也www色官网| 国产三级黄色录像| 久久久久国产一级毛片高清牌| 九色亚洲精品在线播放| av有码第一页| 欧美精品av麻豆av| av在线app专区| av网站免费在线观看视频| 青青草视频在线视频观看| 嫩草影视91久久| 永久免费av网站大全| 人人澡人人妻人| 涩涩av久久男人的天堂| 欧美日韩视频精品一区| 水蜜桃什么品种好| 大陆偷拍与自拍| 男女国产视频网站| 日本91视频免费播放| 国产精品一区二区免费欧美 | 久久精品人人爽人人爽视色| 久久午夜综合久久蜜桃| 国产成人系列免费观看| 青草久久国产| 一区二区三区乱码不卡18| 少妇被粗大的猛进出69影院| 99热网站在线观看| 国产日韩欧美视频二区| 日本a在线网址| 天堂8中文在线网| 91成人精品电影| 婷婷色麻豆天堂久久| 中文字幕高清在线视频| 人妻一区二区av| 久久久久久亚洲精品国产蜜桃av| 一边摸一边做爽爽视频免费| 最近中文字幕2019免费版| 国产片内射在线| av线在线观看网站| 国产一区二区三区av在线| 最近手机中文字幕大全| 欧美黄色片欧美黄色片| 亚洲综合色网址| 国产无遮挡羞羞视频在线观看| 免费观看人在逋| 黑人欧美特级aaaaaa片| 精品一区二区三卡| 美女主播在线视频| 777久久人妻少妇嫩草av网站| 国产在线免费精品| 中国国产av一级| 国产亚洲精品第一综合不卡| 精品熟女少妇八av免费久了| 成年人午夜在线观看视频| 欧美日本中文国产一区发布| av天堂久久9| 国产日韩欧美亚洲二区| 亚洲人成电影免费在线| 丰满饥渴人妻一区二区三| 国产精品av久久久久免费| 国产日韩一区二区三区精品不卡| avwww免费| 日本一区二区免费在线视频| 久久精品久久久久久久性| 夫妻性生交免费视频一级片| av天堂久久9| 国产亚洲精品第一综合不卡| 99热国产这里只有精品6| 欧美精品亚洲一区二区| 久久人妻熟女aⅴ| 女人精品久久久久毛片| 妹子高潮喷水视频| 欧美在线一区亚洲| 免费在线观看日本一区| 欧美av亚洲av综合av国产av| 成人手机av| 老司机亚洲免费影院| 十八禁网站网址无遮挡| 男女国产视频网站| 男人添女人高潮全过程视频| 黄色片一级片一级黄色片| 国产精品免费视频内射| 国产成人精品久久二区二区91| 国产成人免费无遮挡视频| 亚洲精品久久成人aⅴ小说| 国产成人一区二区三区免费视频网站 | 中文字幕制服av| 老鸭窝网址在线观看| 9191精品国产免费久久| 深夜精品福利| 国产高清国产精品国产三级| 天天躁狠狠躁夜夜躁狠狠躁| 日韩一卡2卡3卡4卡2021年| 热re99久久精品国产66热6| 亚洲一码二码三码区别大吗| 又黄又粗又硬又大视频| 天天添夜夜摸| 中文字幕高清在线视频| 看免费成人av毛片| 国产一区二区 视频在线| 欧美97在线视频| 亚洲av成人精品一二三区| 午夜两性在线视频| 大话2 男鬼变身卡| 又黄又粗又硬又大视频| 午夜福利免费观看在线| 天堂8中文在线网| 一区福利在线观看| 国产97色在线日韩免费| 午夜免费成人在线视频| 中文字幕人妻熟女乱码| 国产一级毛片在线| 中文欧美无线码| 亚洲午夜精品一区,二区,三区| 亚洲av电影在线进入| 精品一区二区三区av网在线观看 | 国产片特级美女逼逼视频| 亚洲av日韩精品久久久久久密 | 国产精品香港三级国产av潘金莲 | 午夜视频精品福利| 你懂的网址亚洲精品在线观看| 久久久久视频综合| 亚洲 欧美一区二区三区| 国产成人欧美| 久久鲁丝午夜福利片| 在线观看免费午夜福利视频| 91九色精品人成在线观看| 波多野结衣av一区二区av| 国产精品偷伦视频观看了| 黄色毛片三级朝国网站| av线在线观看网站| 国产伦人伦偷精品视频| 国产精品三级大全| 久久久精品区二区三区| 男女床上黄色一级片免费看| 在线观看免费午夜福利视频| 欧美 亚洲 国产 日韩一| 国产精品国产av在线观看| 日本av免费视频播放| 丰满少妇做爰视频| 五月开心婷婷网| 欧美国产精品一级二级三级| 欧美97在线视频| 日本vs欧美在线观看视频| 女人被躁到高潮嗷嗷叫费观| 香蕉国产在线看| 亚洲欧美精品综合一区二区三区| 永久免费av网站大全| 午夜福利影视在线免费观看| 人体艺术视频欧美日本| 99久久99久久久精品蜜桃| 国产午夜精品一二区理论片| 男女边吃奶边做爰视频| 久久久精品免费免费高清| 久久久久视频综合| 亚洲成人免费av在线播放| bbb黄色大片| 老汉色∧v一级毛片| 免费看十八禁软件| 久久人人爽人人片av| 亚洲av日韩在线播放| 精品国产国语对白av| 欧美黑人欧美精品刺激| 制服人妻中文乱码| 两个人免费观看高清视频| 狂野欧美激情性xxxx| 极品少妇高潮喷水抽搐| 成人国产av品久久久| 丁香六月欧美| 国产精品.久久久| 脱女人内裤的视频| 国产精品九九99| 50天的宝宝边吃奶边哭怎么回事| 一二三四社区在线视频社区8| 制服人妻中文乱码| 午夜福利,免费看| 国产精品一区二区在线不卡| 国产伦理片在线播放av一区| 日本av免费视频播放| 又紧又爽又黄一区二区| 十八禁人妻一区二区| 日韩 亚洲 欧美在线| 在线看a的网站| 国产精品久久久久成人av| 国产精品一区二区在线观看99| 男女高潮啪啪啪动态图| 少妇被粗大的猛进出69影院| 国产高清国产精品国产三级| 黄片小视频在线播放| 国产高清videossex| 精品一品国产午夜福利视频| 国产成人啪精品午夜网站| 国产淫语在线视频| 亚洲综合色网址| 免费黄频网站在线观看国产| 日韩中文字幕视频在线看片| 国产淫语在线视频| 黑人猛操日本美女一级片| 国产主播在线观看一区二区 | 男女床上黄色一级片免费看| 精品少妇久久久久久888优播| 国产91精品成人一区二区三区 | 国产成人欧美在线观看 | 欧美久久黑人一区二区| 99re6热这里在线精品视频| 国产高清国产精品国产三级| 一本一本久久a久久精品综合妖精| 亚洲精品国产色婷婷电影| 国产av一区二区精品久久| svipshipincom国产片| 久久久精品免费免费高清| 国产亚洲欧美精品永久| 国产午夜精品一二区理论片| 9色porny在线观看| 亚洲国产日韩一区二区| 高清不卡的av网站| 夫妻性生交免费视频一级片| 又黄又粗又硬又大视频| 男女下面插进去视频免费观看| 日本欧美国产在线视频| 欧美精品一区二区大全| 久久人人爽人人片av| 日韩大码丰满熟妇| 亚洲精品中文字幕在线视频| 婷婷色综合大香蕉| 亚洲图色成人| 中文字幕色久视频| 日韩 欧美 亚洲 中文字幕| 又粗又硬又长又爽又黄的视频| 欧美在线黄色| 伦理电影免费视频| 91精品三级在线观看| 看免费av毛片| 精品久久久久久久毛片微露脸 | 曰老女人黄片| 色精品久久人妻99蜜桃| 超碰97精品在线观看| 一区在线观看完整版| 午夜精品国产一区二区电影| 日本wwww免费看| 国产日韩欧美亚洲二区| 久久久国产一区二区| 久久亚洲精品不卡| 两个人看的免费小视频| 国产在线免费精品| 一本综合久久免费| 满18在线观看网站| 少妇 在线观看| 婷婷色麻豆天堂久久| 精品少妇一区二区三区视频日本电影| 如日韩欧美国产精品一区二区三区| 老汉色∧v一级毛片| 精品熟女少妇八av免费久了| videos熟女内射| 久久久久国产一级毛片高清牌| 日韩精品免费视频一区二区三区| 精品一区在线观看国产| 十八禁网站网址无遮挡| netflix在线观看网站| 免费女性裸体啪啪无遮挡网站| a级毛片黄视频| 日韩人妻精品一区2区三区| 欧美精品亚洲一区二区| 亚洲美女黄色视频免费看| 老司机靠b影院| 黄色片一级片一级黄色片| 国产日韩欧美亚洲二区| 2018国产大陆天天弄谢| 亚洲成人国产一区在线观看 | 国产av精品麻豆| 热99国产精品久久久久久7| 电影成人av| 亚洲欧美激情在线| 观看av在线不卡| 成年av动漫网址| 国产成人免费观看mmmm| 搡老岳熟女国产| 亚洲国产欧美在线一区| 1024香蕉在线观看| 久久 成人 亚洲| 免费高清在线观看日韩| 在线 av 中文字幕| 只有这里有精品99| 日韩制服丝袜自拍偷拍| 又粗又硬又长又爽又黄的视频| 另类亚洲欧美激情| 啦啦啦在线免费观看视频4| 大陆偷拍与自拍| 男人爽女人下面视频在线观看| 亚洲av国产av综合av卡| 纯流量卡能插随身wifi吗| e午夜精品久久久久久久| 亚洲国产看品久久| 久久人人爽av亚洲精品天堂| 精品久久久久久久毛片微露脸 | 亚洲三区欧美一区| 少妇人妻 视频| 国产精品久久久av美女十八| 久久午夜综合久久蜜桃| 视频区欧美日本亚洲| 超碰成人久久| 只有这里有精品99| 好男人电影高清在线观看| 欧美大码av| 国产视频首页在线观看| av在线老鸭窝| 性少妇av在线| 国产爽快片一区二区三区| 一区二区三区精品91| 亚洲专区国产一区二区| 欧美精品一区二区大全| 久久精品亚洲熟妇少妇任你| 看十八女毛片水多多多| 亚洲色图综合在线观看| 欧美黑人欧美精品刺激| 搡老乐熟女国产| 国产精品国产av在线观看| 多毛熟女@视频| 精品一区二区三卡| 精品国产一区二区久久| 黄网站色视频无遮挡免费观看| 美女大奶头黄色视频| 亚洲欧美清纯卡通| 天天添夜夜摸| 久久精品国产综合久久久| 久久99一区二区三区| 老司机影院成人| 中国美女看黄片| 狠狠精品人妻久久久久久综合| 国产精品亚洲av一区麻豆| 成人午夜精彩视频在线观看| 黄色a级毛片大全视频| 五月开心婷婷网| av电影中文网址| 国产女主播在线喷水免费视频网站| 麻豆乱淫一区二区| 久久精品久久久久久噜噜老黄| 丝袜喷水一区| 在线观看国产h片| www.999成人在线观看| 国产免费一区二区三区四区乱码| 纵有疾风起免费观看全集完整版| 国产无遮挡羞羞视频在线观看| 久久这里只有精品19| 美女午夜性视频免费| 久久国产精品人妻蜜桃| av有码第一页| 国产老妇伦熟女老妇高清| 这个男人来自地球电影免费观看| 国产成人欧美| 精品久久久久久电影网| 欧美av亚洲av综合av国产av| 国产日韩一区二区三区精品不卡| 久久女婷五月综合色啪小说| 精品福利观看| 亚洲人成电影免费在线| 天天影视国产精品| 欧美大码av| 成人亚洲精品一区在线观看| 欧美日韩av久久| 亚洲天堂av无毛| 丝袜喷水一区| 青春草亚洲视频在线观看| 国产成人精品在线电影| 国产精品一区二区精品视频观看| 日本色播在线视频| 在线观看一区二区三区激情| 久久久久久久久免费视频了| 视频区欧美日本亚洲| 美国免费a级毛片| 老鸭窝网址在线观看| 国产高清视频在线播放一区 | 免费观看人在逋| 秋霞在线观看毛片| 亚洲国产av影院在线观看| 午夜福利乱码中文字幕| 一区二区三区四区激情视频| 久久精品国产亚洲av高清一级| 啦啦啦啦在线视频资源| 一本—道久久a久久精品蜜桃钙片| 亚洲视频免费观看视频| 亚洲成人手机| 最新的欧美精品一区二区| 乱人伦中国视频| 男女高潮啪啪啪动态图| 久久精品人人爽人人爽视色| 国产黄色视频一区二区在线观看| 亚洲欧洲国产日韩| 人妻 亚洲 视频| 国产免费视频播放在线视频| 伊人亚洲综合成人网| 午夜精品国产一区二区电影| 一边亲一边摸免费视频| 黄色视频不卡| a 毛片基地| 女警被强在线播放| 精品一品国产午夜福利视频| 国产成人免费无遮挡视频| 精品福利永久在线观看| 少妇粗大呻吟视频| xxx大片免费视频| 国产欧美日韩精品亚洲av| 久久精品亚洲熟妇少妇任你| 一区在线观看完整版| 国产深夜福利视频在线观看| 男女高潮啪啪啪动态图| 免费高清在线观看日韩| 亚洲自偷自拍图片 自拍| 国产精品成人在线| 国产免费又黄又爽又色| 久久青草综合色| 伊人亚洲综合成人网| 老熟女久久久| 黑人欧美特级aaaaaa片| 久久99一区二区三区| 亚洲 国产 在线| 国产高清国产精品国产三级| 免费观看av网站的网址| 久久久久久久大尺度免费视频| 高清av免费在线| 亚洲成色77777| 亚洲精品成人av观看孕妇| 成人国产av品久久久| 精品人妻1区二区| 黑人欧美特级aaaaaa片| 免费观看av网站的网址| 在线观看www视频免费| 啦啦啦在线观看免费高清www| 亚洲中文av在线| 麻豆乱淫一区二区| 亚洲欧美成人综合另类久久久| 亚洲第一av免费看| 制服人妻中文乱码| 男女之事视频高清在线观看 | 大陆偷拍与自拍| 老司机亚洲免费影院| 波多野结衣一区麻豆| 日韩欧美一区视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 香蕉国产在线看| 久久亚洲精品不卡| 国语对白做爰xxxⅹ性视频网站| 欧美成狂野欧美在线观看| 国产亚洲av片在线观看秒播厂| 久久99热这里只频精品6学生| 大型av网站在线播放| 亚洲精品一区蜜桃| 人人妻人人添人人爽欧美一区卜| 纵有疾风起免费观看全集完整版| 亚洲伊人色综图| 国产日韩一区二区三区精品不卡| 性色av乱码一区二区三区2| 丝袜在线中文字幕| √禁漫天堂资源中文www| 欧美激情极品国产一区二区三区| 亚洲欧美中文字幕日韩二区| 久久久精品国产亚洲av高清涩受| 亚洲第一青青草原| 高清欧美精品videossex| 男人操女人黄网站| 国产又爽黄色视频| 丁香六月欧美|