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

    Capurability analysis for arbitrarily high-speed maneuvering targets

    2023-11-10 02:16:24ShooWANGYngGUOShichengWANGLixinWANGYnhuTAOZhengfeiPENG
    CHINESE JOURNAL OF AERONAUTICS 2023年10期

    Shoo WANG, Yng GUO,*, Shicheng WANG, Lixin WANG,Ynhu TAO, Zhengfei PENG

    a Xi’an Research Institution of Hi-Technology, Xi’an 710025, China

    b Xi’an Satellite Control Center, Xi’an 710043, China

    KEYWORDS

    Abstract The capturability of an arbitrarily maneuvering target featuring speed superiority over an interceptor is analyzed for Augmented Pure Proportional Navigation (APPN) and Retro-Augmented Proportional Navigation (RAPN) guidance.This paper focuses on intercepting arbitrary maneuvers to study more general interception problems.A comparative analysis of the capture region between head-on interception related to APPN and head-pursuit interception related to RAPN is proposed.The results indicate that RAPN performs better than APPN in capturability.It is concluded that increasing the target velocity, which increases the velocity ratio, significantly weakens the capturability of the interceptor, and the average acceleration and relative distance affect the location of the capture region but not its size.The analysis is based on prior knowledge of the target maneuver, which inevitably leads to deviations from actual maneuvers in practical engagement, so a deviation analysis is implemented.The effective capture region shrinks as the absolute value of acceleration deviation increases, and the RAPN has a better deviation fault tolerance compared with the APPN.The results reveal that a larger relative distance can weaken the deviation fault tolerance, and the target velocity has opposite effects on head-on and head-pursuit interception.

    1.Introduction

    Pure Proportional Navigation (PPN)1–3guidance law, as a major variant of the Proportional Navigation (PN) guidance law, has been analyzed extensively in the literature.The True Proportional Navigation (TPN) guidance laws,4,5with lateral acceleration normal to the Line-of-Sight (LOS), are another major variant of the PN guidance laws.Unlike TPN,PPN lateral acceleration is applied normal to the interceptor’s velocity vector.The LOS rate item in PPN’s analytical acceleration expression is the only quantity that needs to be measured in the implementation process and is readily available from the interceptor’s Infrared (IR) sensor.Because of its robustness and practicality, PPN outperforms TPN guidance laws in applications.6The class of Generalized Proportional Navigation (GPN) laws was firstly proposed in Ref.7.The results manifested that GPN was superior to TPN with a larger capture region and a shorter interception time.The only work on the derivation of closed-form analytical solutions for PPN was presented by Becker,2and the closed-form analytical solution on TPN, realistic TPN, and Generalization of TPN(GTPN) were presented in Refs.4, 8, and 9 respectively.Although a modification to GTPN 9 was conducted to achieve being almost comparable to PPN in capture zone, in general,the capturability of PPN was found to be better than that of all forms of TPN.10A general study of several variants of PN guidance was done by Yang CD and Yang CC,11and then its extension to 3D PN guidance law was analyzed in Ref.12

    Guelman first studied PPN against a non-maneuvering target,1and then presented sufficient conditions on the capturability against a target executing constant maneuver, and at the same time,some results on the interceptor lateral acceleration requirement were obtained.13The capturability for a maneuvering target executing bounded piecewise continuous time-varying maneuver was analyzed by Ghawghawe and Ghose14using Guelman’s approach.The capturability of an arbitrary target with time-varying normal acceleration using three-dimensional (3-D) PPN was studied by introducing a Lyapunov-like approach.15,16Shin and Li17improved 3-D PPN performance by properly selecting the guidance command direction.The Ideal Proportional Navigation (IPN)which used the relative velocity as the referenced direction of interceptor was studied by Yuan and Chern.18According to the closed-form solutions of IPN, the capture criterion of IPN was only related to the navigation constant.IPN was found to be superior to any other PN laws on capturability in the premise of more energy consumption.Using a line-ofsight fixed natural coordinate and three modified polar variables, Tyan19obtained capture zone of General Ideal Proportional Navigation (GIPN) guidance law.However, these studies generally assume a low-speed target or that the interceptor speed is higher than the target speed.The interception of a higher-speed target using PPN was studied by Tyan and Shen.20,21

    The two collision conditions for interception of a nonmaneuvering target with higher velocity using PPN are Head-On (HO) and Head-Pursuit (HP) interception.Most existing guidance laws focus on HO interception,which is generally designed for the interception of lower-speed targets.When this guidance law is used to engage a higher-speed target, there is a larger closing velocity between the interceptor and target, and the capturability deteriorates significantly.

    Compared with guidance laws based on HO interception,those based on HP interception are designed to intercept high-speed targets.As presented by Golan and Shima,22HP interception places the interceptor ahead of the target on its flight trajectory, so the closing velocity is reduced.Ge et al.23extended the HP interception to a 3-D nonlinear variable space using Lyapunov stability theory.Prasanna and Ghose6described an HP interception guidance law named Retro Proportional Navigation (RPN) that uses a negative navigation constant.Analysis of the RPN against a non-maneuvering target with higher constant speed than interceptor manifests that,even though more time is consumed by RPN to intercept the target successfully than by PN, the interception performance of the RPN is better than that of the PN, in terms of capturability, lateral acceleration, and closing velocity.However, the RPN capturability analysis was only implemented for nonmaneuvering targets.Ghosh et al.24proposed a Composite PN (CPN) guidance law that both PPN and RPN guidance laws are used in a 3-D engagement scenario against higherspeed non-maneuvering targets based on the initial situation and terminal impact angle.Ghosh and Ghose augmented RPN in 3-D engagement and presented a nontrivial RPN target capture zone.The concept of Augmented Proportional Navigation(APN)was proposed to solve the problem of intercepting a maneuvering target in Refs.23–26.

    Augmented Pure Proportional Navigation (APPN), as a variation of the PPN guidance law, was presented to account for target maneuvers and the capturability analysis for it was given.Sufficient conditions for speed ratio, navigation gain,and constraint of lateral acceleration were obtained for targets executing piecewise continuous maneuvers in time.27However,only the situation in which the interceptor’s speed is higher than the target’s speed was considered.The capturability research for high-speed maneuvering targets is still lacking.

    Motivated by the RPN which has good performance in intercepting high-speed targets, and the APPN which is suitable for intercepting a maneuvering target, Retro-Augmented Proportional Navigation (RAPN) comprising these two guidance laws was presented by Sun et al.28Although the capturability analysis of this guidance law is available for targets with active maneuverability, a priori knowledge of the target maneuvering model which is related to the relative distance is needed.The restrictions on target maneuvering model inevitably lead to deviations during realistic maneuvers.

    The main works of this paper are presented as follows:Inspired by Refs.28,29, the capturability analysis of the RAPN for arbitrarily high-speed maneuvering targets is conducted.Different with the definition for the capture region in Refs.28,29, this paper defines the capture region in the initial heading angle space based on Ref.6,because capture region in that space is important for the mid-end guidance handover and the analysis is conducted in that space.The analysis of this paper is mainly focused on the parameters and variables in the guidance process, such as target velocity, navigation gain,relative distance, and acceleration deviation.Based on the works proposed by Prasanna and Ghose6who analyzed the capturability of RPN against non-maneuvering target, this paper studies the capturability for arbitrarily high-speed maneuvering targets.The conclusions of this paper are more general than those of Ref.6.

    The organization of the remainder of this paper is given as follows.Section 2 presents the intended research problem and the zero effort interception condition.Section 3 gives the capture region for an arbitrarily high-speed maneuvering targets for both the APN and the RAPN.Section 4 performs the numerical simulations to analyze the capturability using APPN and RAPN for an arbitrarily maneuvering target.Section 5 draws some interesting conclusions for intercepting arbitrarily high-speed maneuvering targets.

    2.Problem formulation

    The initial condition of terminal guidance reached between the high-speed maneuvering targets and the interceptor after the mid-end guidance process can determine the interceptor’s guidance performance against the target.In many existing studies, the initial value of the heading angle is thought to be directly related to the successful interception of the high-speed maneuvering targets.Because this paper focuses on the interception guidance scenario, the high-speed maneuvering target is the intercepted target in this paper.

    2.1.Planar pursuit preliminaries

    Consider a planar pursuit scenario between an arbitrarily maneuvering target T and an APPN-guided interceptor M,as shown in Fig.1.The assumptions are first introduced into a planar engagement scenario for the interceptor and target.

    Assumption 1.The interceptor and target have constant speeds of VMand VT(with a non-dimensional velocity ratio ρ=VT/VM?1), respectively.

    Assumption 2.The effect of the command delay from the overload acceleration to the actual control variable is neglected.

    Assumption 3.The maximum detection range of the interceptor seeker is Rmax.If the relative range between the interceptor and target is less than the maximum detection range Rmax, the target is considered to be successfully detected and captured in the interceptor’s field of view.

    The engagement scenario shown in Fig.1 is established in the coordinate system XI-OI-YI,where the normal acceleration, LOS, range-to-go, and flight-path angle are denoted as a,q,R,and θ,respectively.The target and interceptor heading angles are denoted as η=θT-qMTand γ=qMT-θM,respectively.

    The nonlinear dynamics of the relative motion between the interceptor and target can be formulated as follows:where ˙qMTis the LOS rate between the interceptor and target;and ˙R is the closing velocity, ˙R>0 means that the relative range increases, and ˙R<0 means that the relative range decreases in the guidance process.In this paper’s interception guidance scenario, ˙R<0 is maintained over the entire guidance process.

    Fig.1 Planar engagement geometry.

    The PN guidance law is widely used in the engagement problem because of its simple expression,where the command acceleration is directly proportional to the LOS rate and can be written as aM=N′VMT˙qMT,where N′is the navigation coefficient,and VMTis the closing velocity.Many modern guidance laws have been developed based on the PN guidance law.Although the forms of modern guidance laws are complicated,they are essentially equivalent to proportional guidance with correction terms.Taking the APN guidance law as an example,its general form can be expressed aswhich is the augmented version of TPN.However,the expressions for both the TPN and APN guidance laws contain closing velocity terms VMT, which are difficult to obtain in practical situations.

    From a practical implementation point of view, the lateral acceleration is applied normal to the interceptor velocity vector under the assumption of a small angle of attack.Therefore,the PPN guidance law is applied by N′VTM=NVM, which is more widely used than TPN.Therefore, augmentation of PPN guidance law would be much more realistic and relevant.The lateral acceleration of the APPN guidance law can be expressed by

    where K is the correction coefficient,K=1 makes the terminal LOS rate converge to zero, and K=3 makes the required terminal acceleration converge to zero,which is the optimal guidance law neglecting the inertial link of the flight control system.

    The acceleration of the target aTis arbitrary,bounded,and piecewise continuous in time throughout the engagement.Therefore, by integrating both sides of Eq.(5), we can obtain

    2.2.Zero effort interception condition

    The interception in the case where the angle between the interceptor velocity vector and the target velocity vector is less than π/2 is named as HO interception, and that in the case where the same angle is greater than π/2 is named as HP interception.

    As shown in Fig.2, take the end O′of the target velocity vector VTas the center of the circle,and use the modulus value of the interceptor velocity vector VMas the radius to make a circle, which intersects with the LOS at point A and point B.The interceptor velocity vector is.Ifis established, the condition for HO zero effort interception is satisfied.Ifis established, the condition for HP zero effort interception is satisfied.The HO zero effort interception means that the APN guidance law is applied by interceptor with positive navigation constant, that is, N>0.The RAPN guidance law is applied by interceptor in the HP zero effort interception, where the navigation constant is a negative constant, that is, N<0.

    Fig.2 Schematic diagram of zero effort interception.

    Lemma 1.For a given value ρ>1, if a point in space (γ,η)reaches the manifold defined by

    in finite time, then successful interception is guaranteed.Proof.Substituting Eq.(9) into Eq.(2), we can obtain

    It can be further determined from Eq.(9) that

    Substituting Eq.(11) into Eq.(1), we find

    The result ˙R<0 can be obtained by substituting ρ>1 into Eq.(12), implying that the interceptor and target are close to each other.Based on the result expressed by Eq.(10) that the LOS rate is equal to zero and the result that the interceptor and target are close to each other,a successful interception can be realized.It is noted that another set of cos γ values from Eq.(2)are given by cosfor which ˙qMT=0 but ˙R>0.This completes the proof.

    The heading angle of the interceptor can be obtained by Eq.(9) as

    With the geometric relationship shown in Fig.2, there are only two values of the interceptor’s lead angle when considering the case γ ?[-π,π]:

    Because γ2=π-sin-1(ρ sin η) for ˙R>0, the result γ=sin-1(ρ sin η) is what this paper intends to analyze.

    The zero effort interception condition boundary is obtained when the LOS is tangent to the interceptor’s velocity vector circle.At this time, the heading angle of the interceptor is γ=±π/2.Thus, we can find the heading angle range for the target:

    where ηmax=arcsin(1/ρ).Eq.(15) can also be written as

    It can be determined from Eq.(16) that the target heading angle range is closely related to the velocity ratio of the target and interceptor.As the target-to-interceptor velocity ratio increases, the range of heading angles allowed by the target gradually decreases.

    Lemma 2.The trajectory guided by the APN or the RAPN guidance law in the space(γ,η)can be represented by a straight line for all values of ρ>1, expressed as

    Proof.Substituting ˙θTfrom Eq.(5) after taking the time derivative of Eq.(4), we have

    Taking the time derivative of Eq.(6) and substituting Eqs.(3) and (7) into it, we can find that

    Substituting Eq.(18) into Eq.(19) yields

    Integrating both sides of Eq.(19) yields

    where γ0and η0are the initial heading angles of the interceptor and target, respectively.

    Then,

    The point (γ0,η0+η′) can be defined as the equivalent initial heading angle in the (γ,η) space.If the target implements no maneuvers in the terminal guidance, then η′=0 and the equivalent initial heading angle is equal to the initial heading angle (γ0,η0).

    3.Capture region for an arbitrarily high-speed maneuvering target

    This section presents the capture region analysis for an arbitrarily high-speed maneuvering target in the initial heading angle space (γ0,η0).The complete capture region results against the maneuvering target for both the APN and the RAPN are given.Compared to the capture region for a nonmaneuvering target,the capture region for a maneuvering target is additionally affected by the initial relative distance.The analyzed capture region results provide a heading angle reference for successfully intercepting maneuvering targets.

    3.1.Capture region for HO interception

    Lemma 3.Assume that the equivalent initial heading angle(γ0,η0+η′)is used as the initial value of the heading angle for the terminal guidance.For ρ>1 and N>0 (HO interception with APN guidance law), the region

    defines the capture region of equivalent initial heading angle in (γ,η) space, where

    where

    and

    and then,

    Remark 1.Although Lemma 3 gives the capture region for the interceptor in (γ,η) space, the capture region only defines the equivalent initial heading angle.The effects generated by the average turn rate of the target ^anTand the homing phase time t in η′of Eq.(17) can be neglected in (γ,η) space, where the nonlinearity from the initial heading angles γ0and η0to the manifold f(γ,η)=0 is not reflected.The transformation from the equivalent interceptor and target initial heading angles γ0and η0+η′to the manifold f(γ,η)=0 obeys a linear law where the trajectory in (γ,η) space is a straight line with slopeThus, the capture region defined in (γ0,η0) space needs to be given for the initial heading angle.

    It is noted that the effects generated by the average turn rate of the target ^anTand the homing phase time t in η′of Eq.(17) cannot be neglected when obtaining the capture region of the initial heading angle.The capture region for the equivalent initial heading angle in (γ,η) space can be refined into an initial heading angle in(γ0,η0)space by considering the effects of the nonlinear terms η′.

    Theorem 1.For ρ>1 and N>0 (HO interception with APN guidance law), the region

    defines the capture region of the initial heading angle in(γ0,η0) space, where

    Proof.The target can only be captured by the interceptor when the equivalent initial heading angle γ0, η0+η′of the interceptor and the target lie in the capture region RAPNof the (γ,η) space from Lemma 3.Firstly, the situation that the point of the equivalent initial heading angle (γ0,η0+η′) lying in the capture region A1∩B1is analyzed.The following inequality reflects the relationship between the point and the lines L1(γ,η)=0,and that between the point and the manifold f(γ,η)=0 can be given as

    Next,the situation that the point(γ0,η0+η′)lies in the capture region A2∩B2is also analyzed.The following inequality reflects the relationship between the point and the lines L2(γ,η)=0, and that between the point and the manifold f(γ,η)=0 can also be given as

    The extreme value of the heading angle appears at the terminal guidance end time, which can be estimated by

    where R0and V0are the initial values of the relative distance and the closing velocity, respectively.

    From the constraint |η|≤arcsin(1/ρ)given by Eq.(16),the target heading angle value is regarded as small enough with the velocity ratio ρ ?1.Thus, cos η0≈η0, and the equation can be written as

    Substituting Eq.(48) into Eqs.(43)–(46), we have

    Therefore, the capture region in (γ0,η0) space is

    This completes the proof.

    Fig.3 gives the capture region for HO interception given byin (γ0,η0) space.Compared with the capture region presented by Prasanna and Ghose in Ref.6, the boundaries L1(γ0,η0)=0 and L2(γ0,η0)=0 shown in Fig.3 behave as a nonlinear form,the reasons for which is that the boundary suffers from the nonlinear term.It can be seen from Fig.3 that the capture region values on the y-axis are all negative, meaning that the interceptor cannot intercept the target when the target heading angle η0is positive because the capture region suffers from target turning effects.

    This section is focused on the HO interception condition,where the APN guidance law with N>0 is used to intercept an arbitrarily maneuvering target.As discussed earlier, HP interception is another condition for intercepting the target and is generated by the RAPN guidance law with N<0.The following section discusses the capture region analysis for the RAPN guidance law.

    3.2.Capture region for HP interception

    Lemma 4.Assume that the equivalent initial heading angles(γ0,η0+η′) are used as the initial heading angles for the terminal guidance.For ρ>1 and N<0 (HP interception with RAPN guidance law), the region

    defines the capture region of equivalent initial heading angle in (γ,η) space, where

    where

    and

    Fig.3 Capture region for HO interception given by .

    Proof.There exist the same arguments as given by the proof of Lemma 3, and more details can be obtained from it.

    Theorem 2.For ρ>1 and N<0(HP interception with RAPN guidance law), the region of

    defines the capture region of the initial heading angle in(γ0,η0) space, where

    where

    3.3.Acceleration deviation impact on capture region

    The concept of the average turn rate ^anT(t)= ^aT/VTis introduced in Section 1.1,1.2,2.1 and 2.2.It is important for calculating the capture region where the zero effort interception manifold f(γ0,η0) and the boundary of the capture region L1(γ0,η0), L2(γ0,η0), L3(γ0,η0), and L4(γ0,η0) mainly depend on the average turn rate ^anT(t)= ^aT/VT.Because of the constant velocity of the target,the average turn rate^anT(t)is determined by ^aT, the average lateral acceleration of the target in terminal guidance, denoted asin Section 2.1.Therefore, the average lateral acceleration value is essential in determining the capture region in (γ0,η0) space.

    Remark 2.Although it is said that target maneuvering model is known as a prior knowledge,there still exist great limitations in actual capturability analysis in Ref.28.In target maneuver model assumption of Ref.28, the target maneuvers using an acceleration aTwhich is applied normal to the LOS, and is given by

    where b is a nonnegative constant gain and ρ is the relative distance between target and missile.

    Eq.(74) gives the prior knowledge of target maneuvering model mentioned in Ref.28, which indicates that acceleration aTis related to the relative distance.Ref.28 intends to simplify the dynamic equations in modified polar variables for capturability analysis by associating acceleration aTwith the reciprocal of the relative distance 1/ρ.Although the dynamic equations have been simplified by the assumption expressed by Eq.(74), it still brings great limitations in actual situation.Thus,capturability analysis presented by this paper is based on the prior knowledge that average accelerationis known in advance.

    Remark 3.Because the average lateral acceleration is used in the capture region analysis,Theorems 1 and 2 presented in this paper can be applied to arbitrarily maneuvering targets given prior knowledge about the target’s maneuver.Therefore,this is the point of departure of the analysis presented in this paper from that in Ref.6,29.Apart from measuring the relative motion information at every instant, the implementation of this theorem requires only the a priori knowledge of the target’s maneuvering capability instead of estimating the target lateral acceleration at every instant.Thus, prior knowledge of the average lateral acceleration makes it possible to analyze the capturability so that the heading angle lies in the suitable region to intercept the target successfully.

    Fig.4 Capture region for HP interception given by .

    The above discussion gives the results that the average lateral acceleration of the target is directly related to the accurate solution of capture region and it can be achieved by a priori knowledge about target maneuver.Based on the results, the accurate capture region in the (γ0,η0) space can be given by Theorems 1 and 2 in theory.However, there exists a controversy that the prior knowledge about the target depends on the properties of the target itself which is often not reachable by interceptors, making it difficult to guarantee the accuracy of the average lateral acceleration.The deviation σ is inevitably caused since the average lateral acceleration in the realistic engagement is not equal to that given by priori knowledge.There exists a certain deviation between the capture region calculated by average lateral acceleration in the presence of a priori knowledge and that in realistic engagement, which is caused by the deviation σ.

    Fig.5 gives the different capture regions for HO interception with ^aT=40 m/s2, ^aT=30 m/s2, and ^aT=25 m/s2, where^aT=30 m/s2is the true average lateral acceleration, ^aT=40 m/s2is the value from prior knowledge with deviation σ=10 m/s2compared to the true situation, and ^aT=25 m/s2is the value given from prior knowledge with deviation σ=-5 m/s2compared to the true situation.It can be seen from Fig.5 that the intersection is generated between the realistic capture region and the capture region with deviation σ=10 m/s2,which means that the capture region in intersection still enables the zero effort interception to be satisfied when the average lateral acceleration given by prior knowledge has the deviation σ=-5 m/s2with the realistic value.This indicates that the deviation σ=-5 m/s2is a suitable value given by the prior knowledge within an acceptable margin of deviations.But there is not intersection between the realistic capture region and the capture region with deviation σ=10 m/s2,which means that no intersection is generated by two capture regions to ensure the zero effort interception to be satisfied when the average lateral acceleration given by prior knowledge has the deviation σ=10 m/s2with the realistic value.This indicates that the deviation σ=10 m/s2exceeds the acceptable margin of deviations so that the zero effort interception cannot be guaranteed whatever initial heading angles are.

    The above discussion shows that an acceleration deviation from mismatched prior knowledge impacts capturability.Therefore,we analyze the impact of the acceleration deviation on the capture region as follows.

    Theorem 3.For ρ>1 and N>0 (HO interception with APN guidance law), the region

    Fig.5 Different capture regions for HO interception with^aT =40 m/s2(acceleration deviation σ=10 m/s2),^aT =30 m/s2(true value),and ^aT =25 m/s2(acceleration deviation σ=-5 m/s2).

    defines the sets of acceptable average lateral acceleration deviations σ to intercept the target successfully for an initial interceptor heading angle γ0, where

    where

    Case 1.For σ>0, the capture region calculated using+σ gradually moves down as the deviation increases in (γ0,η0)space from Fig.5.When the deviation σ belongs to an acceptable deviation margin, an intersection between the true capture region and the capture region with deviation σ is guaranteed.

    The boundary of the capture region with acceleration deviation σ can be represented as

    It can be seen from Fig.5 that the intersection between the boundaryof the capture region with deviation and the boundary L2(γ0,η0) of the true capture region can be found.Thus, the inequality-L2(γ0,η0)>0 can be satisfied, yielding

    Then,the lower boundary of the acceleration deviation can be obtained:

    Case 2.For σ<0, the capture region calculated using ^aT+σ gradually moves up as the deviation decreases in(γ0,η0)space from Fig.5.When the deviation σ belongs to an acceptable deviation margin, an intersection between the true capture region and the capture region with acceleration deviation σ is guaranteed.

    It can be seen from Fig.5 that the intersection between the boundary L1(γ0,η0) of the capture region with deviation and the boundary Lσ2(γ0,η0) of the true capture region can be found.Thus, the inequality Lσ2(γ0,η0)-L1(γ0,η0)<0 can be satisfied, yielding

    Then,the upper boundary of the acceleration deviation can be obtained:

    The lower boundary D1(γ0,σ)=0 and the upper boundary D2(γ0,σ)=0 of the acceleration deviation σ have now been obtained, completing the proof.

    Theorem 4.For ρ>1 and N<0(HP interception with RAPN guidance law), the region

    defines the sets of acceptable average lateral acceleration deviations σ to intercept the target successfully for initial interceptor heading angle γ0, where

    where

    Table 1 Initial parameters for interception engagement simulation.

    Fig.6 Capture intervals for HO (top image) and HP (bottom image)interception with η0 =-60° for different navigation gains.

    Proof.The arguments used to prove Theorem 3 can also be used to prove Theorem 4.Please see the proof of Theorem 3 for details.

    Fig.7 Capture interval lengths of γ0 for HO(top image)and HP(bottom image) interception with different values of η0.

    4.Simulation analyses

    This section performs numerical simulations to analyze the capturability using APN and RAPN for an arbitrarily maneuvering target.First, the initial conditions are given for the terminal guidance between the interceptor and the target.Subsequently, capturability analyses for the maneuvering target for HO interception with APN and HP interception with RAPN are implemented.The results are used to show that increasing the target’s velocity can weaken the capturability and that the acceleration and relative distance can change the capture interval’s direction in polar coordinates but not its length.Finally, the acceleration deviation analysis for the average acceleration is implemented,with the significant result that HP interception has a better deviation fault tolerance compared with HO interception and the deviation’s fault tolerance can weaken as the relative distance increases.

    4.1.Capturability analysis for maneuvering target

    The initial conditions for the capturability simulation of the interceptor using HO and HP interception are set in Table 1.The velocity ratio ρ=VT/VM=2>1 can be obtained from Table 1, and the models used are simple kinematic models,as given in Section 2.1.

    Fig.6 presents the possible capture interval of both the HO interception and HP interception with η0=-60°for different values of N.The capture interval corresponding to the initial flight directions of the interceptor is defined in polar coordinates.It can be seen from Fig.6 that the HP interception capture interval is larger than the HO interception capture interval at the same absolute value of N.The capture interval shrinks as the absolute value of N decreases.

    Fig.7 shows the capture interval length of γ0for different values of η0for HO and HP interception.As soon as the target enters the terminal guidance so that the initial parameters are determined, the capture region of the interceptor can be ensured for the different initial target heading angles η0.It can be observed from Fig.7 that the target is capturable by the interceptor only for η0in the interval [-101.22°,-25.24°]for HO interception and only in the interval[-91.29°,-4.55°] for HP interception.Successful interception cannot be realized when the initial value η0is outside these intervals.The maximum value of the capture interval length occurs at η0=-36.75°for HO interception and η0=-63.24°for HP interception.The maximum value of capture interval length for HP interception is larger than that for HO interception, indicating that HP interception, where RAPN is used to intercept the target,has a larger interception range than HO interception.

    Fig.8 shows the possible capture intervals for both HO and HP interception with η0=-60°at different target velocities.The capture interval shrinks as the target velocity increases because the velocity ratio ρ=VT/VMincreases with increasing target velocity.This result verifies that the high-speed maneuvering targets are more difficult to capture than low-velocity targets, and the capturability of high-speed targets decreases as the target velocity increases.Thus, increasing the target velocity, which increases the velocity ratio ρ, weakens the interceptor’s capturing ability.

    As shown in Fig.5, the capture region in (γ0,η0) space moves in the negative direction of the y-axis as the average target acceleration increases.The target cannot be intercepted for positive values of η0because the calculated capture region is completely on the negative part of the y-axis.Because the average acceleration’s impact on the capture region is reflected mainly on the y-axis, a suitable γ0is selected to facilitate the analysis of the acceleration effect on the capture region.Fig.9 presents the possible capture intervals for HO interception with γ0=50°and HP interception with γ0=200°for different average accelerations.It can be seen from Fig.9 that the acceleration can only change the capture interval’s location in polar coordinates but not its length.

    Similarly,the impact of the relative distance on the capture region is mainly reflected on the y-axis.Fig.10 presents the possible capture intervals for HO interception with γ0=50°and HP interception with γ0=200°for different relative distances.These results also indicate that the relative distance can only change the capture interval’s location in polar coordinates but not its length.The important conclusion is that the average target acceleration and the initial relative distance only affect the capture region’s location but not its size in (γ0,η0)space.

    Fig.11 shows the capture region for HO interception guided by PN and HP interception guided by RPN against non-maneuvering target, which is a special case of this paper proposed by Prasanna and Ghose6where the target’s acceleration is equal to zero.Fig.12 shows the capture region for HO interception guided by PN and HP interception guided by RPN against maneuvering target, which is another special case where the correction coefficient K of command acceleration shown by Eq.(7) is equal to zero.Because the above two cases are two special cases of this paper where one case is that the acceleration of target is equal to zero and another case is that the correction coefficient K of command acceleration is equal to zero, the conclusion given by the analysis in this paper can be applied to both two cases as well.

    4.2.Acceleration deviation analysis

    Figs.13 and 14 are effective capture region for HO and HP interception with inaccurate target maneuver caused by the mismatched acceleration deviation σ respectively.As can be seen from Figs.13 and 14, the capture region marked by dark shades is the realistic one which is provided by the real target maneuver=15 m/s2,and the captureregionmarked by light shadesis the inaccurate one which is provided by the inaccurate target maneuver+σ generated by mismatched prior knowledge.In Figs.13 and 14, the values of deviation σ are taken as |σ|=5,10,14 m/s2so that the variation of effective capture region can be displayed.The effective capture region,the intersection of the dark shades and the light shades in Figs.13 and 14,is the region that the successful interception can still be achieved while the inaccurate target maneuver=+σ generated by mismatched prior knowledge cannot be avoided.It is manifested that the initial heading angle of interceptor and target γ0,η0lying in the effective capture region can guarantee the successful interception.The effective capture region shrinks as the absolute value of acceleration deviation |σ|increases for HO and HP interception from Figs.13 and 14, that is, successful interception will become increasingly difficult to achieve.Fig.15 shows the capture interval lengths of η0for HO (top image) interception at γ0=35°and HP(bottom image)interception at γ0=220°with different values of σ.Compared with HO interception, HP interception can produce a larger effective capture region under the same deviation condition from Figs.13-15, which manifests that the HP interception has a better deviation fault tolerance.

    Fig.12 Capture region for HO interception guided by PN (left image) and HP interception guided by RPN (right image) against maneuvering target.

    Fig.13 Effective capture region for HO interception with inaccurate target maneuver caused by mismatched acceleration deviation σ.

    Fig.16 presents the acceptable acceleration deviations for HO and HP interception.The acceptable acceleration deviation formulation is presented in Section 3.3,which means that the acceleration deviation caused by the mismatch between the acceleration values in presence of the prior knowledge and the realistic values causes the deviation of the capture region.If the deviation is within a reasonable range, as shown in the shaded areas of Fig.16, the capture region in intersection between the capture region obtained by prior knowledge and that obtained by realistic value still enables the zero effort interception to be satisfied.If the deviation is not within a reasonable range, the intersection cannot guarantee that a zero effort interception occurs.It can be seen from Fig.16 that the acceptable deviation region for HP interception is larger than that for HO interception and the region is not affected by the values of average acceleration.

    Fig.17 shows the acceptable acceleration deviation ranges for HO and HP interception with different navigation gains.The acceptable acceleration deviation range increases as the absolute value of N increases, and the navigation gains have a greater impact on the HO interception acceptable acceleration deviation range than on HP interception.Larger absolute values of N values make the deviation fault tolerance higher.

    Fig.16 Acceptable acceleration deviation ranges for HO (top image) and HP (bottom image) interception.

    Fig.17 Acceptable acceleration deviation ranges for HO (top image) and HP (bottom image) interception with different navigation gains.

    Fig.18 shows the acceptable acceleration deviation ranges for HO and HP interception at different target velocities.Notably, the acceptable acceleration deviation range shrinks as the target velocity increases for HO interception, but the opposite trend occurs for HP interception.Fig.19 shows the acceptable acceleration deviation ranges for HO and HP interception with different relative distances.The acceptable acceleration deviation range shrinks as the relative distance between the interceptor and target increases, indicating that larger relative distances weaken the fault tolerance.

    The above analysis of the guidance parameters for deviations provides a reference for target interception.Suitable initial conditions for the terminal guidance can improve the probability of successful interception.

    5.Conclusions

    Capturability analysis of APPN and RAPN for arbitrarily high-speed maneuvering targets are conducted in this paper.The analysis is available for any target with arbitrarily maneuvering form since the concept of average acceleration is introduced, which can be known by a prior knowledge about the target maneuver.Different from the previous studies, the acceleration deviation analyses are implemented, where the deviation stems from the mismatch between the prior knowledge and the realistic maneuver.The studies proposed by Prasanna and Ghose6can be seen as a special case of this paper when the acceleration of the target is equal to 0 m/s2.

    Fig.18 Acceptable acceleration deviation ranges for HO (top image) and HP (bottom image) interception at different target velocities.

    Some interesting conclusions can be drawn via numerical simulation.The capturability of interceptor for the highspeed targets decreases as the velocity of that target increases,which verifies that the high-speed targets are more difficult to capture than low-velocity targets.The average acceleration and the initial relative distance only have the effect on the location of the capture region but not the size of it in the (γ0,η0)space.In the analysis of the acceleration deviation, the results manifest that the larger relative distance of the target can weaken the fault tolerance more for the deviation and the velocity has true opposite effect on the HO interception and the HP interception.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (No.62073335), the Science Fund for Distinguished Young People in Shaanxi Province, China(No.2022JC-42),and the China Postdoctoral Science Foundation (Nos.2017M613201, 2019T120944 and 2020M683737).

    Fig.19 Acceptable acceleration deviation ranges for HO (top image)and HP(bottom image)interception with different relative distances.

    √禁漫天堂资源中文www| 91国产中文字幕| 久久精品久久久久久噜噜老黄| xxx大片免费视频| av卡一久久| 秋霞在线观看毛片| 日韩欧美一区视频在线观看| 久久久久国产网址| 3wmmmm亚洲av在线观看| 久久久精品94久久精品| 老司机亚洲免费影院| 久久精品国产鲁丝片午夜精品| 99re6热这里在线精品视频| 男女啪啪激烈高潮av片| 国产国拍精品亚洲av在线观看| 人妻人人澡人人爽人人| 老女人水多毛片| 女人久久www免费人成看片| 国产精品无大码| 国产亚洲最大av| av电影中文网址| 久久精品国产亚洲av天美| 亚洲精品美女久久av网站| 国产精品国产三级专区第一集| 欧美人与性动交α欧美精品济南到 | 午夜影院在线不卡| 不卡视频在线观看欧美| 中文字幕久久专区| 国产探花极品一区二区| 亚洲第一区二区三区不卡| 欧美日韩精品成人综合77777| 国产亚洲午夜精品一区二区久久| 天美传媒精品一区二区| 午夜福利视频精品| 啦啦啦中文免费视频观看日本| av女优亚洲男人天堂| 国产69精品久久久久777片| 亚洲精品一区蜜桃| 2021少妇久久久久久久久久久| 亚洲图色成人| 人人妻人人爽人人添夜夜欢视频| 日日摸夜夜添夜夜添av毛片| 欧美精品一区二区大全| 七月丁香在线播放| 校园人妻丝袜中文字幕| 老司机亚洲免费影院| 成年av动漫网址| 在线观看免费视频网站a站| 天天影视国产精品| 久久久久久久久久久免费av| 精品人妻偷拍中文字幕| 99久久中文字幕三级久久日本| 91久久精品国产一区二区成人| 久久精品久久久久久久性| 亚洲人成网站在线观看播放| 久久99一区二区三区| 亚洲av.av天堂| 午夜免费鲁丝| 午夜福利影视在线免费观看| 午夜激情av网站| 久久精品国产亚洲av天美| 免费看不卡的av| 在线观看一区二区三区激情| 国产精品三级大全| 97超视频在线观看视频| 久久久久视频综合| 国产成人免费观看mmmm| 亚洲精品国产色婷婷电影| 最近中文字幕高清免费大全6| a级毛片免费高清观看在线播放| 少妇的逼好多水| 亚洲怡红院男人天堂| 亚洲色图 男人天堂 中文字幕 | 成人综合一区亚洲| 国产伦理片在线播放av一区| 色94色欧美一区二区| 22中文网久久字幕| 午夜福利网站1000一区二区三区| av在线观看视频网站免费| 亚洲欧美精品自产自拍| 精品少妇黑人巨大在线播放| 午夜影院在线不卡| 内地一区二区视频在线| 久久久国产欧美日韩av| 男女无遮挡免费网站观看| 91精品三级在线观看| 国产男女超爽视频在线观看| 久久久久久久久久久久大奶| 观看av在线不卡| 一边亲一边摸免费视频| 日韩成人av中文字幕在线观看| 免费日韩欧美在线观看| 国产一级毛片在线| 久久国产精品男人的天堂亚洲 | .国产精品久久| 亚洲激情五月婷婷啪啪| 桃花免费在线播放| 9色porny在线观看| 国产精品嫩草影院av在线观看| 晚上一个人看的免费电影| 免费人成在线观看视频色| 亚洲经典国产精华液单| 亚洲av中文av极速乱| 欧美国产精品一级二级三级| 国产老妇伦熟女老妇高清| videos熟女内射| 日韩,欧美,国产一区二区三区| 天堂中文最新版在线下载| 国产 一区精品| 亚洲综合色网址| 久久 成人 亚洲| 国产精品一二三区在线看| 91成人精品电影| 99re6热这里在线精品视频| 成人毛片60女人毛片免费| 国产熟女午夜一区二区三区 | 精品视频人人做人人爽| 久久久国产精品麻豆| 国产黄频视频在线观看| 日韩中文字幕视频在线看片| 日日摸夜夜添夜夜爱| 国产一级毛片在线| 国产av码专区亚洲av| 国产男人的电影天堂91| 亚洲精品456在线播放app| 成人免费观看视频高清| 九色成人免费人妻av| 免费黄色在线免费观看| 女性生殖器流出的白浆| 亚洲内射少妇av| 秋霞在线观看毛片| 看免费成人av毛片| 国产精品人妻久久久影院| 简卡轻食公司| 国产成人一区二区在线| 最近最新中文字幕免费大全7| 老司机影院成人| 考比视频在线观看| 国产免费现黄频在线看| 成人国语在线视频| 一边摸一边做爽爽视频免费| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 又黄又爽又刺激的免费视频.| 成人毛片a级毛片在线播放| 日韩强制内射视频| 超碰97精品在线观看| 伦理电影大哥的女人| 七月丁香在线播放| 欧美丝袜亚洲另类| 日本免费在线观看一区| 九色成人免费人妻av| 热re99久久精品国产66热6| 伦理电影大哥的女人| 日韩中文字幕视频在线看片| 飞空精品影院首页| 精品人妻偷拍中文字幕| 美女内射精品一级片tv| 97超碰精品成人国产| 久久精品久久久久久噜噜老黄| 国产精品久久久久成人av| av专区在线播放| 欧美精品一区二区免费开放| 国产亚洲精品久久久com| 国产av一区二区精品久久| 欧美日韩在线观看h| 91精品伊人久久大香线蕉| tube8黄色片| 91午夜精品亚洲一区二区三区| 日日爽夜夜爽网站| 久久这里有精品视频免费| 哪个播放器可以免费观看大片| 丝袜美足系列| 青春草视频在线免费观看| 只有这里有精品99| 亚洲av免费高清在线观看| 欧美激情 高清一区二区三区| 又大又黄又爽视频免费| 亚洲五月色婷婷综合| 欧美激情极品国产一区二区三区 | 久久国内精品自在自线图片| 日韩,欧美,国产一区二区三区| 亚洲精品乱久久久久久| 亚洲精品第二区| 黄色配什么色好看| 国产亚洲最大av| 欧美成人午夜免费资源| 婷婷色麻豆天堂久久| 国产成人免费无遮挡视频| 精品一区二区三卡| 搡女人真爽免费视频火全软件| 免费日韩欧美在线观看| 成人免费观看视频高清| 精品久久国产蜜桃| 老女人水多毛片| 免费人妻精品一区二区三区视频| 中文字幕精品免费在线观看视频 | 91久久精品电影网| 色哟哟·www| 午夜激情福利司机影院| 欧美老熟妇乱子伦牲交| 免费高清在线观看视频在线观看| 欧美另类一区| 国产乱来视频区| 久久av网站| 日本av免费视频播放| 国国产精品蜜臀av免费| 五月天丁香电影| 久久狼人影院| 在线看a的网站| 五月玫瑰六月丁香| 午夜福利影视在线免费观看| 欧美一级a爱片免费观看看| 一级,二级,三级黄色视频| 大码成人一级视频| 大香蕉97超碰在线| 午夜影院在线不卡| 在线看a的网站| 国产成人精品婷婷| 黄色视频在线播放观看不卡| 日韩电影二区| 国产成人精品无人区| 国产国拍精品亚洲av在线观看| 亚洲高清免费不卡视频| 国产精品久久久久久精品电影小说| 一区二区日韩欧美中文字幕 | 亚洲av日韩在线播放| 九色成人免费人妻av| 韩国高清视频一区二区三区| www.av在线官网国产| 久久久久久久久大av| 少妇熟女欧美另类| 寂寞人妻少妇视频99o| 观看美女的网站| 亚洲av成人精品一区久久| 日本与韩国留学比较| 大香蕉97超碰在线| 久久国产精品大桥未久av| 国产精品一区二区在线不卡| 国产淫语在线视频| 18禁裸乳无遮挡动漫免费视频| 国产69精品久久久久777片| 亚洲欧美精品自产自拍| 午夜91福利影院| 久久av网站| 国产一区亚洲一区在线观看| 亚洲不卡免费看| 亚洲av中文av极速乱| 一本久久精品| 国产精品一区二区在线观看99| 91久久精品国产一区二区三区| 久久精品国产鲁丝片午夜精品| 免费人妻精品一区二区三区视频| 九色亚洲精品在线播放| 中文字幕人妻丝袜制服| 99久久综合免费| 国产不卡av网站在线观看| 亚洲综合色惰| 香蕉精品网在线| 日韩av不卡免费在线播放| 国产免费一区二区三区四区乱码| 欧美一级a爱片免费观看看| 亚洲一区二区三区欧美精品| 美女xxoo啪啪120秒动态图| 国产69精品久久久久777片| 久久久久国产网址| 少妇人妻 视频| 国产免费视频播放在线视频| 精品亚洲乱码少妇综合久久| 国产欧美亚洲国产| 精品国产一区二区久久| 美女cb高潮喷水在线观看| 一级毛片aaaaaa免费看小| 国产黄色免费在线视频| 国产精品一区二区在线不卡| 丝袜在线中文字幕| 色94色欧美一区二区| 亚洲一级一片aⅴ在线观看| 欧美 亚洲 国产 日韩一| 亚洲丝袜综合中文字幕| 精品久久蜜臀av无| 日本91视频免费播放| 日本色播在线视频| 国产69精品久久久久777片| 精品人妻在线不人妻| 国产极品粉嫩免费观看在线 | 亚洲av二区三区四区| 久久人人爽av亚洲精品天堂| 赤兔流量卡办理| 亚洲一级一片aⅴ在线观看| 纯流量卡能插随身wifi吗| 超色免费av| 亚洲国产精品成人久久小说| 免费观看在线日韩| 国语对白做爰xxxⅹ性视频网站| 亚洲国产日韩一区二区| 日韩成人av中文字幕在线观看| 人妻 亚洲 视频| 一本色道久久久久久精品综合| 18禁动态无遮挡网站| 久久99精品国语久久久| 少妇 在线观看| 男女免费视频国产| 我要看黄色一级片免费的| 成人免费观看视频高清| 黄色毛片三级朝国网站| 国产一区二区三区综合在线观看 | 日本黄色片子视频| 18在线观看网站| 国产精品嫩草影院av在线观看| 成人亚洲欧美一区二区av| 欧美日韩视频高清一区二区三区二| 成年女人在线观看亚洲视频| 自拍欧美九色日韩亚洲蝌蚪91| 内地一区二区视频在线| 亚洲欧洲国产日韩| 午夜精品国产一区二区电影| 看非洲黑人一级黄片| 99热全是精品| 亚洲精品乱码久久久久久按摩| 欧美 亚洲 国产 日韩一| 亚洲成人av在线免费| 黄色毛片三级朝国网站| 国产精品.久久久| 视频中文字幕在线观看| 在线精品无人区一区二区三| 欧美性感艳星| 欧美老熟妇乱子伦牲交| 18禁在线播放成人免费| 欧美人与善性xxx| 丰满乱子伦码专区| 久热这里只有精品99| 啦啦啦在线观看免费高清www| 丝袜喷水一区| 成人国语在线视频| 亚洲精品中文字幕在线视频| 精品熟女少妇av免费看| 黄色怎么调成土黄色| 免费黄网站久久成人精品| 亚洲国产精品一区二区三区在线| 国产日韩一区二区三区精品不卡 | videosex国产| 中文精品一卡2卡3卡4更新| 另类精品久久| 国产av码专区亚洲av| 国产亚洲最大av| 乱人伦中国视频| 久久免费观看电影| 亚洲精品久久午夜乱码| 少妇精品久久久久久久| 国产精品国产三级专区第一集| 18禁在线无遮挡免费观看视频| av专区在线播放| 99热国产这里只有精品6| 精品人妻熟女av久视频| 成年人免费黄色播放视频| 一本色道久久久久久精品综合| 丰满饥渴人妻一区二区三| 国产免费现黄频在线看| √禁漫天堂资源中文www| 晚上一个人看的免费电影| 制服丝袜香蕉在线| 久久99热这里只频精品6学生| 国产日韩欧美在线精品| 亚洲av欧美aⅴ国产| 日韩强制内射视频| 涩涩av久久男人的天堂| 18在线观看网站| 美女内射精品一级片tv| 成人午夜精彩视频在线观看| 国产一区二区在线观看日韩| 五月天丁香电影| 人妻制服诱惑在线中文字幕| 一边亲一边摸免费视频| 欧美 日韩 精品 国产| 永久免费av网站大全| 日日爽夜夜爽网站| av专区在线播放| 欧美日韩精品成人综合77777| 亚洲精品日本国产第一区| 岛国毛片在线播放| 少妇被粗大猛烈的视频| 亚洲色图综合在线观看| 精品视频人人做人人爽| 久久精品久久久久久噜噜老黄| 国产精品99久久99久久久不卡 | 最新中文字幕久久久久| 天堂俺去俺来也www色官网| 91国产中文字幕| 超碰97精品在线观看| 成人黄色视频免费在线看| 欧美丝袜亚洲另类| 免费播放大片免费观看视频在线观看| 内地一区二区视频在线| 免费不卡的大黄色大毛片视频在线观看| 免费黄频网站在线观看国产| 亚洲av福利一区| 日韩精品免费视频一区二区三区 | av国产精品久久久久影院| www.av在线官网国产| 男男h啪啪无遮挡| 久久久久久久久久人人人人人人| 最近手机中文字幕大全| 亚洲av欧美aⅴ国产| 狂野欧美白嫩少妇大欣赏| 亚洲激情五月婷婷啪啪| 蜜桃国产av成人99| 亚洲怡红院男人天堂| 乱码一卡2卡4卡精品| 一本色道久久久久久精品综合| 精品久久久噜噜| 日韩av免费高清视频| 少妇的逼水好多| 国产日韩一区二区三区精品不卡 | 国产片内射在线| 国产午夜精品一二区理论片| 各种免费的搞黄视频| 香蕉精品网在线| 亚洲精品,欧美精品| 日韩 亚洲 欧美在线| 亚洲一区二区三区欧美精品| 999精品在线视频| 又大又黄又爽视频免费| 国产精品麻豆人妻色哟哟久久| 午夜激情久久久久久久| 热99久久久久精品小说推荐| 久久青草综合色| 伦理电影大哥的女人| 最后的刺客免费高清国语| 肉色欧美久久久久久久蜜桃| 最近2019中文字幕mv第一页| 亚洲精品日韩av片在线观看| 精品人妻一区二区三区麻豆| 午夜久久久在线观看| h视频一区二区三区| 成人国产av品久久久| 午夜免费鲁丝| h视频一区二区三区| 国产精品秋霞免费鲁丝片| 国产亚洲一区二区精品| 亚洲欧美一区二区三区黑人 | 最新的欧美精品一区二区| 一二三四中文在线观看免费高清| 在线播放无遮挡| 久久人人爽人人爽人人片va| av线在线观看网站| 日本色播在线视频| 亚洲精品日本国产第一区| 夫妻性生交免费视频一级片| 亚洲av电影在线观看一区二区三区| 天堂8中文在线网| 水蜜桃什么品种好| 亚洲人成77777在线视频| 久久 成人 亚洲| 麻豆成人av视频| 成人漫画全彩无遮挡| 人妻 亚洲 视频| 少妇的逼好多水| 国产色爽女视频免费观看| 国产成人91sexporn| 熟妇人妻不卡中文字幕| 亚洲欧洲国产日韩| 国产爽快片一区二区三区| 成人手机av| 狂野欧美激情性bbbbbb| 免费人妻精品一区二区三区视频| 九九爱精品视频在线观看| 日韩三级伦理在线观看| 国产又色又爽无遮挡免| 国产av国产精品国产| 高清午夜精品一区二区三区| 少妇 在线观看| 婷婷色综合大香蕉| 丝袜喷水一区| 制服人妻中文乱码| .国产精品久久| 美女国产高潮福利片在线看| 街头女战士在线观看网站| 国产深夜福利视频在线观看| 一边亲一边摸免费视频| 亚洲在久久综合| 一级毛片我不卡| 天天操日日干夜夜撸| 欧美3d第一页| 亚洲av男天堂| 国产av一区二区精品久久| 久久99热这里只频精品6学生| 人体艺术视频欧美日本| 亚洲五月色婷婷综合| 69精品国产乱码久久久| 天天影视国产精品| 青青草视频在线视频观看| 黑丝袜美女国产一区| 欧美精品亚洲一区二区| 3wmmmm亚洲av在线观看| 国产免费又黄又爽又色| 午夜免费鲁丝| 免费观看在线日韩| 在线观看www视频免费| 精品久久久久久电影网| 赤兔流量卡办理| 免费观看无遮挡的男女| 两个人的视频大全免费| 蜜桃国产av成人99| 日韩三级伦理在线观看| 欧美人与善性xxx| 性色avwww在线观看| 亚洲av中文av极速乱| 蜜桃久久精品国产亚洲av| 国产亚洲欧美精品永久| 国产一区二区三区综合在线观看 | 又黄又爽又刺激的免费视频.| 久久久久人妻精品一区果冻| 成人亚洲欧美一区二区av| 亚洲人与动物交配视频| 少妇人妻 视频| 精品一品国产午夜福利视频| 黄色毛片三级朝国网站| 欧美精品一区二区大全| 午夜精品国产一区二区电影| 亚洲内射少妇av| 国产欧美亚洲国产| av女优亚洲男人天堂| 一区在线观看完整版| 国产一区有黄有色的免费视频| 成人漫画全彩无遮挡| 99国产综合亚洲精品| 国产日韩欧美视频二区| 女人精品久久久久毛片| 亚洲无线观看免费| 高清欧美精品videossex| 在线观看免费日韩欧美大片 | 日韩强制内射视频| 黑丝袜美女国产一区| 韩国av在线不卡| 亚洲国产av新网站| 国产不卡av网站在线观看| 精品国产国语对白av| 2022亚洲国产成人精品| 国产日韩欧美亚洲二区| av天堂久久9| 精品国产一区二区三区久久久樱花| 51国产日韩欧美| 亚洲怡红院男人天堂| 午夜免费男女啪啪视频观看| av卡一久久| 亚洲精品亚洲一区二区| 日日啪夜夜爽| 欧美性感艳星| 国产成人精品久久久久久| 丝袜在线中文字幕| 九色成人免费人妻av| 国产精品.久久久| 狂野欧美白嫩少妇大欣赏| 九九久久精品国产亚洲av麻豆| 热99久久久久精品小说推荐| 亚洲国产最新在线播放| 欧美日韩成人在线一区二区| 丝袜喷水一区| 久久青草综合色| 交换朋友夫妻互换小说| 黑人巨大精品欧美一区二区蜜桃 | 国产日韩欧美视频二区| 中国三级夫妇交换| 国产精品99久久久久久久久| 亚洲精品日本国产第一区| 日韩伦理黄色片| 欧美变态另类bdsm刘玥| 性色avwww在线观看| 91aial.com中文字幕在线观看| 少妇精品久久久久久久| av国产精品久久久久影院| 91久久精品电影网| 内地一区二区视频在线| av天堂久久9| av.在线天堂| 飞空精品影院首页| 国产精品欧美亚洲77777| 久久久久久久久久人人人人人人| 性色av一级| 亚洲综合色网址| 久久久久久久久久久免费av| 日韩一区二区三区影片| 国产一区亚洲一区在线观看| 国产成人精品福利久久| 在线免费观看不下载黄p国产| 我的女老师完整版在线观看| 免费看光身美女| av女优亚洲男人天堂| 国产伦精品一区二区三区视频9| 少妇的逼水好多| 婷婷色综合www| 3wmmmm亚洲av在线观看| 午夜免费男女啪啪视频观看| 免费观看的影片在线观看| 3wmmmm亚洲av在线观看| 久久毛片免费看一区二区三区| 日韩欧美精品免费久久| 日本色播在线视频| 国产av一区二区精品久久| 插阴视频在线观看视频| 成人黄色视频免费在线看| 日韩欧美一区视频在线观看| 性色avwww在线观看| 免费观看a级毛片全部| 欧美变态另类bdsm刘玥| 我要看黄色一级片免费的| 国产欧美日韩综合在线一区二区| 国产精品人妻久久久久久| 热99国产精品久久久久久7| 婷婷色综合大香蕉| 国产成人午夜福利电影在线观看| 亚洲四区av| 一区二区av电影网| 国产精品99久久久久久久久| 日韩 亚洲 欧美在线|