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

    Automatic landing system using neural networks and radio-technical subsystems

    2017-11-21 12:54:31RomulusLunguMihaiLungu
    CHINESE JOURNAL OF AERONAUTICS 2017年1期

    Romulus Lungu,Mihai Lungu

    University of Craiova,Faculty of Electrical Engineering,107 Decebal Blvd.,Craiova,Romania

    Automatic landing system using neural networks and radio-technical subsystems

    Romulus Lungu,Mihai Lungu*

    University of Craiova,Faculty of Electrical Engineering,107 Decebal Blvd.,Craiova,Romania

    Adaptive control;Automatic landing;Neural network;PID;Reference model

    The paper focuses on the design of a new automatic landing system(ALS)in longitudinal plane;the new ALS controls the aircraft trajectory and longitudinal velocity.Aircraft control is achieved by means of a proportional-integral(PI)controller and the instrumental landing system–the first phase of landing(the glide slope)and a proportional-integral-derivative(PID)controller together with a radio-altimeter–the second phase of landing(the flare);both controllers modify the reference model associated with aircraft pitch angle.The control of the pitch angle and longitudinal velocity is performed by a neural network adaptive control system,based on the dynamic inversion concept,having the following as components:a linear dynamic compensator,a linear observer,reference models,and a Pseudo control hedging(PCH)block.The theoretical results are software implemented and validated by complex numerical simulations;compared with other ALSs having the same radio-technical subsystems but with conventional or fuzzy controllers for the control of aircraft pitch angle and longitudinal velocity,the architecture designed in this paper is characterized by much smaller overshoots and stationary errors.

    1.Introduction

    Most aircraft have automatic landing systems(ALSs)based on the instrumental landing system1,2and different conventional controllers(proportional-derivative–PD,proportionalintegral – PI,proportional-integral-derivative – PID),1–4for aircraft trajectory’s control during landing,and the use of optimal control laws(H2,H∞,H2/H∞),together with full-or reduced-order observers,provides good results.3,5Because the atmospheric conditions and the dynamics of aircraft are drastically changing during flight and,of course,during landing,it is difficult to land safely by using conventional controllers.To design perfect conventional controllers,one has to know the precise mathematical model of the system to be controlled.Furthermore,the aircraft dynamics may vary with respect to the altitude and the flight conditions.Therefore,the adaptive controllers are better choices.

    The presence of unknown or partially known nonlinearities in aircraft dynamics leads to the necessity of using evolved adaptive control architectures in various stages of the flight and especially during landing.Also,the actuators have strongly nonlinear behavior because of the saturation of theirdisplacements and/or mobile elements’velocities.In these cases,in the context of using the inversion of the nonlinearities’dynamics,the adaptive control based on dynamic inversion and neural networks theories is a very good choice.6–9The adapting and the train of the neural networks(NNs)are based on the signals provided by the observers which receive information relative to the error of the automatic control system.The adaptive component of the automatic control law(provided by the NNs)must compensate the inversion(approximation)error of the aircraft dynamics’nonlinear subsystem.For a safety landing,the required information in longitudinal plane is obtained by means of gyro transducers(for aircraft pitch angle–θ and pitch rate–q=˙θ),accelerometers(providing the acceleration signal and,by integration,the aircraft longitudinal velocity),or aerodynamic transducers(for aircraft attack angle10)whether or not the landing control architecture includes an observer.

    One of the most commonly employed nonlinear control method is the feedback linearization;in the automatic landing systems’design it is used in some papers,11but the paper presents limited insight into the performance of simulations of this controller and no tests are performed outside of these simulations;the main disadvantage of the feedback linearization method is that all parametric plant uncertainties must appear in the same equation of the state-space representation as the control.

    Feed-forward neural networks based on the back propagation learning algorithm have also been used12;the main disadvantage is that the neural networks require a priori training on normal and faulty operating data.Other approaches involves the usage of the time delay neural networks;a controller based on this type of neural networks has been designed,13but its main drawback is related to the flight path track accuracy and the fact that it is enable only under limited conditions.Several neural network control approaches have been proposed based on Lyapunov stability theory.13,14The main advantage of these control schemes is that the adaptive laws were obtained from the Lyapunov synthesis and,therefore,guarantee the system’s stability;the disadvantage is that some conditions should be assumed;these requirements are not easy to be satisfied in practical control application.15Juang designed a new learning technique using a time delay network or networks with backpropagation through time algorithms to control the landing16;the main drawbacks are:(1)the number of hidden units was determined by trial;(2)the convergence time is high.Seven different neural network structures(including critic or Radial Basis Function Neural Networks)have been used for obtaining intelligent auto-landing controllers by means of linearized inverse dynamic model17,18;also,the fuzzy logic technique was used to design controllers that track a pre-determined flight path trajectory for safe landing.19In the research area of optimal synthesis,Ref.20have developed a mixed technique for theH2/H∞control of landing,while Ref.21have used theH∞control technique to design an approach for aircraft automatic and landing.In these papers,the authors did not analyze the robustness of the designed controllers in the presence of sensor errors and external disturbances.22

    The auto-landing systems designed in the above mentioned works are characterized by insufficient generality or accuracy;the neural network and dynamic inversion based control approaches could bring improvements.Thus,the paper presents a new adaptive landing architecture for aircraft control in longitudinal plane.According to the authors of this paper,little progress has been reported for the landing flight control systems in longitudinal plane by using neural networks,dynamic inversion concept,linear dynamic compensator,state observer,and PCH block;this motivates the present study.Also,it is interesting to see if the aircraft’s trajectory during landing in longitudinal plane can be tracked with high accuracy by a neural network based controller which uses both the dynamic inversion technique and PCH blocks.The main advantages of the dynamic inversion are:(1)the plant nonlinearities are canceled;(2)the closed loop plant behaves like stable linear system;(3)simplicity in the control structure,ease of implementation,global exponential stability of the tracking error,etc.11On the other hand,the strong point of the neural networks is their approximation ability,these being capable to approximate an unknown system dynamics through learning.A PCH block eliminates the NNs’adapting difficulties.Having in mind the advantages of the NNs,dynamic inversion approach,PCH blocks and the combining of these elements with linear dynamic compensators and state observers,the present paper brings absolute novelty in the search area of ALSs’design.

    The paper is organized as follows:the structure of the ALS is given in the second section;the design of the adaptive system for the control of the pitch angle and the longitudinal component of the flight velocity is presented in the third section;in the next section,complex simulations to validate the new designed ALS have been performed and the obtained results are analyzed;finally,some conclusions are shared in the fifth section of the paper.

    2.Structure of the new automatic landing system

    The automatic control of aircraft during landing(longitudinal plane)is achieved by means of two systems:an automatic system for the trajectory’s control and an automatic system for the control of flight velocity.The automatic system for the control of the flight trajectory in longitudinal plane has two subsystems:(1)the first one is for the aircraft control during the glide slope phase(control of the angular deviation Γ = γ- γc; γ and γcare the real and the calculated slope angles of the aircraft trajectory,respectively),by using an ILS system for the determination of the angle Γ at altitudesH≥H0;(2)the second one is for the aircraft control during the flare phase by means of a radio-altimeter(control of the altitudeH,H<H0,H0– the altitude at which the glide slope phase ends and the second landing phase begins).There are many papers in the literature which threat the control of aircraft in lateral-directional plane in the presence of crosswind.These systems cancel the angular deviation of aircraft with respect to the runway direction,the deviation of the flight direction relative to the runway,and aircraft lateral velocity.Concluding,the control of aircraft in the lateral-directional plane(one landing phase:initial approach)can be achieved by means of other automatic control systems,while the control in longitudinal plane(two landing phases:glide slope and flare)can be achieved by using the neural networks based adaptive system presented in this paper.

    The controllers for the glide slope and flare act toward changing the calculated pitch angle θc;the pitch angle θ may be controlled,for example,by using a PD controller.This contains controllers for the glide slope,flare,pitch angle,and flight velocity11;there,the wind shears and the errors of the sensors have been taken into consideration.

    To compensate the influence of the unknown or partially known nonlinearities from the dynamics of aircraft and actuators,the PD controllers for aircraft pitch angle and longitudinal velocity should be replaced with adaptive architectures based on dynamic inversion and neural networks.The structure of the new automatic system for aircraft control in longitudinal plane is given in Fig.1;U0is the velocity from the moment when aircraft begins the descent along the glide slope,R–the distance between aircraft and the runway threshold,τ–the time constant associated to the flare curve,α–the attack angle,U– aircraft longitudinal velocity,y=[θ,U]T– the controlled output vector,yc=[θc,Uc]T– the calculated output vector,X–horizontal distance covered during landing,x=[Vx,α,ωy,θ,H]T– the state vector of the longitudinal dynamics, δeand δT– the states associated to the actuators(δe– the elevator deflection and δT– the throttle command),u=[δe,δT]T–the command vector,δecand δTcare the commands applied to elevator and to engines,respectively;the significance of the other variables in Fig.1 are to be presented later in this paper.The actuators are described by the equation:

    TeandTTare the time constants associated to the elevator and engine,respectively.The linearized dynamics of aircraft longitudinal motion(with respect to the landing trajectory)is4:

    Above,V0is the nominal value of aircraft speed.

    The linearization of an aircraft’s nonlinear dynamics is generally based on the small disturbances’method with respect to an equilibrium trajectory.Here,the linear model of the aircraft during landing consists of Eqs.(1)–(3).

    3.Design of the new automatic landing system

    The obtaining of the aircraft desired landing trajectory mainly involves two variables’control:the forward speed(Vx)and the pitch angle(θ).According to the landing requirements,all types of aircraft must firstly descend from the cruising altitude to a lower one;meanwhile,the aircraft speed is also reduced from the cruising value to an approach value and,after that,it remains constant.So,when we choose the desired trajectory,we firstly choose the desired forward speed.

    The structure of the adaptive system for the control of the output vector y=[θ,U]T– is presented in Fig.3(Fig.3(a) –complete form,(b)–simplified form),where the block A(having the input u and the output y)represents the longitudinal dynamics of aircraft together with the actuators.This system may be described by the general equations:

    with x=[U,α,q,θ,H,δe,δT]T,u=[δec,δTc]T,and y=[θ,U]T;fandhare nonlinear functions,generally unknown or linear functions obtain in the linearization process with respect to the landing trajectory in longitudinal plane;in this paper,we assume thath(x)is a uniquely invertible function.The aircraft linear dynamics in longitudinal plane is obtained(deduced)from the aircraft nonlinear dynamics by linearization and by using the exact total differential method.4The linearization has been made around a speci fic operating point(regime)–flight at constant altitude,the initial angular rates(roll,pitch,and yaw angular rates)being null.

    Let us considerri–the relative degrees of the system described by Eq. (4) with respect to the variablesyi,i=1,2;y=[y1,y2]T=[θ,U]T.The system described by Eq.(4)may be described by means of the equations:

    one has separated the derivatives that do not depend onufrom the componentshridepending on u.The vectors λiand Yihave the following forms:

    The system described by Eq.(5)is modeled by the block diagram in Fig.4,where

    The system described by Eq.(5)may be put under the vectorial form:

    with

    otherwise,we have:

    ε=[ε1,ε2]Tis the approximation(inversion)error of the function hr,which acts as a disturbance that must be compensated by an adaptive component of the control law.

    This component is calculated with respect to the signals provided by the reference model(command filter).

    The control law^v may be chosen of the following form6:

    with the notations

    Thus,the relative degrees are:r1=3 andr2=2,r=5.Therefore,the transfer functionsHdi(s)have the forms:

    Putting together the two Eq.(17),one gets:

    From Eqs.(21)and(22),by identification,it results^v and ε.Thus,according to Eq.(10),

    Now,replacing Eq.(25)in Eq.(14),one obtains the equation:

    The main purpose of the Pseudo Control Hedging block is to prevent the adaptive element of a control system from trying to adapt to a class of system input characteristics(characteristics of the plant or of the controller).To achieve this issue,the control law is prevented from ‘seeing”the system characteristic;preventing the adaptive element from ‘knowing” about a system’s characteristics means to prevent that adaptive element from seeing the system characteristic as model tracking error.On the other hand,the adaptive controllers are sensitive to actuator nonlinearity.That is why,in these controllers’architecture,one may introduce a block which limits the pseudocontrol by means of a component representing an estimation of the execution element’s dynamics.In this paper,we introduce a PCH block which limits the signalvwith a component representing the actuator dynamics’estimation.Thus,PCH‘moves back” the reference model,introducing correction responses of the reference models with respect to the estimation of the execution element’s position.The signal provided by PCH(vh1)is a reference model’s additional input.23

    modeled by the block diagram in Fig.5(b).

    The states Yi,i=1,2 having the form in Eq.(6)i.e.the vectors associated to the reduced-order subsystem having the inputviand the outputyi,with the transfer functionHdi(s)verify the equation:

    From the state Eq.(32),the vectorial state equation(the equation of the reduced-order subsystem in Fig.3(b),having the input v,the output y,and the transfer function Hd(s))is obtained;it is:

    The two components of the output and of the output vector of the linear dynamic compensator are expressed as:

    Replacing Eq.(42)in Eq.(40),one obtains the state equation:

    and which is modeled by the state equation5:

    where L is the gain matrix of the observer;Eq.(46)is equivalent with the following one:

    The gain L is calculated such that the matrix~A is asymptotically stable.

    The adaptive control law is calculated by using the formula7:

    where W and V–the weights of a neural network(NNc)are the solutions of the equations’system24:

    withn1– the number of neurons in the neural network’s input layer,d– the sample step,whileIi,i=1,n1are the outputs of the neural network.We considered that NNcis a feed-forward type neural network having 3 layers:an input layer(withn1neurons),a hidden layer(withn2neurons),and an output layer(withn3neurons);the weights of the neural network have the forms5:

    vi,jare the weights associated to the connections between the input layer and the hidden layer,whilewj,kare the weights associated to the connections between the hidden and the output layer.The relationship between the functions σj(for the hidden layer)and the outputs of the NNc(vak,k=1,n3)are25:

    bvandbware the biases of the neural network.The sigmoid function σ(z)= σ(VTη)has the form:

    ajbeing the activation potentials whose values are different for each neurons;

    Now,we make a brief stability analysis of the system in Fig.1:

    (2)The stability of the two outer loops(one for the first landing phase and one for the second one):imposing negative roots for the afferent characteristic equations,one obtains the coef ficients of the two controllers.

    Concluding,the synthesis(the design)of the system in Fig.1 involves the calculation of the controllers’coefficients in conditions of stability.

    4.Numerical simulation results

    To study the functionality of the new designed ALS,we consider a light aircraft(Charlie type)having the matrices from Eq.(3)of the following forms2:

    To calculate the parameters of the two controllers(glide slope and flare),we impose negative roots for the linear subsystems which have on feedback ways the ILS system and the radio-altimeter,respectively;these two subsystems have,on their direct ways,the glide slope/flare controllers in series with the linear automatic subsystem for the adaptive control of the pitch angle;this adaptive subsystem (closed by unitary negativefeedbacks)hason itsdirectwaystheblocks(kp1+kd1s) andHd1(s)=1/s2(s+1). The following values resulted:kc=0.2°/V,T1=0.4s,T2=0.04s,Tc=30s,k′c=-1.5(°)/(m ·s-1),Ti=7.5 × 103s,Td=0.9 s.

    The new architecture of ALS(Fig.1)is software implemented in Matlab/Simulink environment and we obtained the time characteristics for the glide slope landing phase(Fig.6)and flare landing phase(Fig.7);in Fig.8,the characteristicsH(X)is plotted;Xrepresents the horizontal distance covered by aircraft during landing.If the actuator associated to the elevator deflection is nonlinear,its models is the one presented in Fig.2;in this case,it is good to use a PCH block(vh1≠0)because it allows the system to work in the linear zones of the nonlinearities.As we already stated above,the signal provided by the PCH block is a reference model’s additional input.Thus,in Figs.6–8,we analyze the influence of Pseudo Control Hedging usage on the variables describing the new designed ALS,taking or not into account the signal from the PCH block.All the characteristics prove the new ALS’s stability and its small overshoots;on the other hand,the differences between the two cases(with and without PCH block)are visible.The usage of a PCH block(when the actuator is nonlinear)does not modify the final values of the variables,but it decrease the transient regime period.The landing begins at a longitudinal speed initially exceeding the nominal speed by 3 m/s(first mini-graphic in Fig.6).The speed should be reduced to the normal speed(67 m/s)and then kept at this value;the landing process begins at 100 m(seventh mini-graphic in Figs.6 and 8),while the flare phase begins at 3.25 m(seventh mini-graphic in Figs.7 and 8).From Figs.6 and 7,we can see that the final error between the desired path and the actual path(altitude error)is less than 0.5 m during the glide slope phase and 0 m for the flare.These errors are very good if the Federal Aviation Administration(FAA)accuracy requirements for Category III(the best category)are analyzed28;according to FAA Category III accuracy requirements,the vertical error(altitude deviation with respect to its nominal value)must be less than 0.5 m(first landing phase)and 0 m(second landing phase).From the forward speed’s point of view,the error is less than 0.3 m/s(first mini-graphics in Figs.6 and 7);also,the aircraft slope angle tends to its desired value: γ→ γc=-2.5°,for the first landing phase,and γ→ γc=0°for the second phase.

    The problem of landing has been discussed in other papers,different types of ALSs being designed.For the same aircraft type,we can make a brief comparison between the obtained characteristics with the ones in other paper1;this ALS consists of:an ILS system,a radio-altimeter,and two conventional controllers(a proportional-derivative controller for the pitch angle and a proportional controller for the flight velocity).We can remark that in the case of the ALS designed in the present paper,the dynamic regimes are slightly slower than the onesforconventionalcontrollersand similarto fuzzy controllers,but the overshoots and stationary errors are much better,this leading to the conclusion that the neural networksbased adaptive controllers are more efficient than the conventional and fuzzy ones for aircraft landing(see Table 1).Our new ALS uses modern design technique based on the usage of a feed-forward neural network and of the dynamic inversion concept,this bringing advantage over classical and fuzzy control techniques:our new adaptive control method has applicability to problems involving multivariate systems with cross-coupling between channels;moreover,by using PCH blocks,the non-linear constraints are generally well-handled.

    Table 1 Comparison between overshoot of the neural networks-based adaptive controllers,conventional and fuzzy ones.

    5.Conclusions

    The ALS designed in this paper controls the landing trajectory and the longitudinal velocity by means of two controllers(a glide slope controller and a flare one-P.and PID forms,respectively)and an adaptive control system with neural network,based on the dynamic inversion concept,having as components:a linear dynamic compensator(PD type),a linear observer,two reference models,and a PCH block for the delay of the signal provided by the pitch angle’s reference model(to eliminate the neural network adapting difficulties which appear because of the actuators’nonlinearities).In the design process of the adaptive control system,we have split the aircraft longitudinal dynamics into two vectorial components:Hd(s)(the matrix consisting of the transfer functions of the linear reduced-order subsystems associated to aircraft pitch angle and longitudinal velocity)and hr(y,u).This decomposition is based on the two differential equations associated to aircraft dynamics;these equations express the relative degrees with respect to the main two variables(aircraft pitch angle and its longitudinal velocity)and provide the vectorial functions hr(x,u)and ε– the inversion error(error of the hrfunction’s approximation).Then,we chose the reference model and we designed the control law,the PCH block,the linear dynamic compensator,the linear observer,and the neural network.The aircraft dynamic characteristics(obtained in Matlab/Simulink through numerical simulation for the new ALS–Fig.1)are characterized by much smaller overshoots and stationary errors than the ones associated to the ALSs based on conventional controllers or fuzzy controllers.

    Acknowledgements

    This work was supported by the Grant No.89/1.10.2015(Modern architectures for the control of aircraft landing)of the Romanian National Authority for Scientific Research and Innovation,CNCS–UEFISCDI,project code PN-IIRU-TE-2014-4-0849.

    1.Lungu R,Lungu M,Grigorie TL.ALSs with conventional and fuzzy controllers considering wind shears and gyro errors.J Aerospace Eng2000;26(4):794–813.

    2.McLean D.Automatic flight control systems.New York:Prentice Hall Publisher;1990.

    3.Lungu R,Lungu M,Grigorie TL.Automatic control of aircraft in longitudinal plane during landing.IEEE Trans Aerospace Electron Syst2013;49(2):1338–50.

    4.Lungu M.Sisteme de conducere a zborului.Craiova:Sitech Publisher;2008.

    5.Lungu M,Lungu R.Estimarea sta?rii aparatelor de zbor(State estimation of the flying objects).Craiova:Sitech Publisher;2014.

    6.Calise AJ,Hovakymyan N,Idan M.Adaptive output feedback control of nonlinear systems using neural networks.Automatica2001;37(8):1201–11.

    7.Calise AJ,Lee S,Sharma M.Direct adaptive reconfigurable control of a tailless fighter aircraft.AIAA guidance,navigation,and control conference and exhibit.Reston:AIAA;1998.

    8.Calise AJ,Lee H,Kim N.High bandwidth adaptive flight control.AIAA guidance,navigation and control conference;2000 August 14–17;Montreal,Canada.Reston:AIAA;2000.

    9.Singh S,Padhi R.Landing of unmanned aerial vehicles using dynamic inversion.Proceedings of international conference on aerospace science and technology;2008 June 26–28;Bangalore,India;2008.

    10.Aron I,Lungu R,Cismaru C.Sisteme de navigatie aerospatiala.Craiova:Scrisul Romanesc Publisher;1989.

    11.Singh S,Padhi R.Automatic path planning and control design for autonomous landing of UAVs using dynamic inversion.American control conference riverfront.2009 June 10-12;St.Louis,U.S.A.Piscataway(NJ):IEEE Press;2009.

    12.Wagner T,Valasek J.Digital autoland control laws using quantitative feedback theory and direct digital design.J Guid,Control,Dyn2007;30(5):1399–413.

    13.Jang JO,Jeon GJ.A parallel neuro-controller for DC motors containing nonlinear friction.Neurocomputing2000;30(1):233–48.

    14.Seshagiri S,Khalil HK.Output feedback control of nonlinear systems using RBF neural network.IEEE Trans Neural Networks2000;11(1):69–79.

    15.Wai RJ.Development of new training algorithms for neurowavelet systems on the robust control of induction servo motor drive.IEEE Trans Ind Electron2002;49(6):1323–41.

    16.Juang JG,Cheng KC.Application of neural networks to disturbances encountered landing control.IEEE Trans Intell Transp Syst2006;7(4):582–8.

    17.Saini G,Balakrishnan SN.Adaptive critic based neuro-controller for autolanding of aircraft.Proceedings of the American control conference;1997 June 6–6;New Mexico,USA,Piscataway(NJ):IEEE Press;1997.

    18.Juang JG,Chang HH,Cheng KC.Intelligent landing control using linearized inverse aircraft model.Proceedings of American control conference;2002.p.3269–74.

    19.Raj KD,Tattikota G.Design of fuzzy logic controller for auto landing applications.Int J Sci Res Publ2013;3(5):1–9.

    20.Kumar V,Rana KP,Gupt V.Real-time performance evaluation of a fuzzy PI+fuzzy PD controller for liquid-level process.Int J Intell Control Syst2008;13(2):89–96.

    21.Lau K,Lopez R,Onate E.Neural networks for optimal control of aircraft landing systems.Proceedings of the world congress on engineering;2007 July 2–4;London,U.K.;2007.p.1–8.

    22.Lungu R,Lungu M.Automatic control of aircraft in lateraldirectional plane during landing.Asian J Control2016;18(3):1–16.

    23.Johnson EN,Calise AJ.Pseudo-control hedging:a new method for adaptive control.Advances in navigation guidance and control technology workshop;2000 November 1–2;Redstone Arsenal,Alabama,U.S.A.;2000.

    24.Chwa D,Choi J.Adaptive nonlinear guidance law considering control loop dynamics.IEEE Trans Aerospace Electron Syst2003;39(4):1134–43.

    25.Johnson EN,Calise AJ,Rysdyk R,Shirbiny EA.Feedback linearization with neural network augmentation applied to X-33 attitude control.AIAA guidance,navigation and control conference and exhibit;2000 August 14–17;Dever,U.S.A.Reston:AIAA;2000.p.1–11.

    26.Calise A,Johnson EN,Johnson MD,Corban JE.Applications of adaptive neural–networks control to unmanned aerial vehicles.J Harbin Inst Technol2006;38(11):1865–9.

    27.Mori R,Suzuki S.Neural network modelling of lateral pilot landing control.J Aircraft2009;46(5):1721–6.

    28.Braff R,Powell JD,Dorfler J.Applications of GPS to air traffic control.Global positioning system:theory and applications,Reston,vol.II.Progress in astronautics and aeronautics.Reston:AIAA;1996.p.327–74.

    16 February 2016;revised 29 April 2016;accepted 27 May 2016

    Available online 7 January 2017

    ?2017 Chinese Society of Aeronautics and Astronautics.Production and hosting by Elsevier Ltd.This is anopenaccessarticleundertheCCBY-NC-NDlicense(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    *Corresponding author.

    E-mail address:Lma1312@yahoo.com(M.Lungu).

    Peer review under responsibility of Editorial Committee of CJA.

    91狼人影院| 亚洲av成人精品一二三区| 最近视频中文字幕2019在线8| 免费观看性生交大片5| 国产中年淑女户外野战色| 亚洲在线观看片| 九草在线视频观看| av国产免费在线观看| 又粗又硬又长又爽又黄的视频| 性插视频无遮挡在线免费观看| 久久精品国产亚洲网站| 国产伦一二天堂av在线观看| 国产 一区精品| 嫩草影院精品99| 亚洲欧美中文字幕日韩二区| 午夜精品一区二区三区免费看| 亚洲经典国产精华液单| 菩萨蛮人人尽说江南好唐韦庄 | 欧美3d第一页| 成人午夜精彩视频在线观看| 神马国产精品三级电影在线观看| 最近视频中文字幕2019在线8| 欧美xxxx性猛交bbbb| 国产在视频线在精品| 成年版毛片免费区| 99久久精品热视频| 亚洲va在线va天堂va国产| 欧美高清性xxxxhd video| 精品少妇黑人巨大在线播放 | 国产91av在线免费观看| 久久久久久久久久成人| 五月伊人婷婷丁香| 久久午夜福利片| 午夜精品国产一区二区电影 | 如何舔出高潮| 久久久久久久久久成人| 久久亚洲精品不卡| 性插视频无遮挡在线免费观看| 成人亚洲精品av一区二区| 内射极品少妇av片p| 国产三级在线视频| 青青草视频在线视频观看| 精品一区二区三区视频在线| 国产精品一二三区在线看| 天堂影院成人在线观看| 亚洲av二区三区四区| 少妇被粗大猛烈的视频| 国产精品久久久久久av不卡| 国产大屁股一区二区在线视频| 久久久久久久久久久免费av| 高清日韩中文字幕在线| 1024手机看黄色片| 波多野结衣巨乳人妻| 22中文网久久字幕| av在线天堂中文字幕| 99热全是精品| 久久这里有精品视频免费| 一本一本综合久久| 国产v大片淫在线免费观看| 中文资源天堂在线| 亚洲成人精品中文字幕电影| 99国产精品一区二区蜜桃av| 国产亚洲av片在线观看秒播厂 | 国产精品国产高清国产av| 成人三级黄色视频| 少妇裸体淫交视频免费看高清| 成年免费大片在线观看| 大香蕉久久网| 免费播放大片免费观看视频在线观看 | 好男人在线观看高清免费视频| 亚洲av成人精品一区久久| 全区人妻精品视频| 精品久久久久久久久av| 国产人妻一区二区三区在| 一级毛片aaaaaa免费看小| 永久免费av网站大全| 变态另类丝袜制服| 建设人人有责人人尽责人人享有的 | 午夜福利高清视频| 国产一区二区亚洲精品在线观看| 免费黄色在线免费观看| 日韩在线高清观看一区二区三区| 日韩中字成人| 亚洲av熟女| 亚洲人成网站在线播| 免费无遮挡裸体视频| 国产精品av视频在线免费观看| 99久久无色码亚洲精品果冻| 真实男女啪啪啪动态图| 中文字幕av成人在线电影| a级毛片免费高清观看在线播放| 最近2019中文字幕mv第一页| 精品少妇黑人巨大在线播放 | 黄色欧美视频在线观看| 日日啪夜夜撸| 成人漫画全彩无遮挡| 麻豆av噜噜一区二区三区| 欧美+日韩+精品| 国产 一区精品| 少妇猛男粗大的猛烈进出视频 | 成人二区视频| 久久久久国产网址| 免费电影在线观看免费观看| 国产精品久久久久久精品电影小说 | 久久婷婷人人爽人人干人人爱| 亚洲高清免费不卡视频| 午夜精品一区二区三区免费看| 国内少妇人妻偷人精品xxx网站| 国产视频首页在线观看| 91精品国产九色| 欧美不卡视频在线免费观看| 国产精品国产三级国产av玫瑰| 国产亚洲av嫩草精品影院| 国产免费又黄又爽又色| 2021天堂中文幕一二区在线观| 日本av手机在线免费观看| 成人亚洲精品av一区二区| 中文天堂在线官网| 少妇的逼好多水| 一级黄片播放器| 伦精品一区二区三区| 少妇被粗大猛烈的视频| 99在线视频只有这里精品首页| 成人欧美大片| 午夜日本视频在线| 色播亚洲综合网| 日韩欧美国产在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产国拍精品亚洲av在线观看| 国产精品美女特级片免费视频播放器| 亚洲在线观看片| 亚洲成人av在线免费| 26uuu在线亚洲综合色| av又黄又爽大尺度在线免费看 | 欧美区成人在线视频| 国产欧美另类精品又又久久亚洲欧美| 成年免费大片在线观看| av.在线天堂| 床上黄色一级片| 国产精品久久久久久精品电影小说 | 久久久久久久久久黄片| a级毛色黄片| 欧美激情久久久久久爽电影| 亚洲国产欧洲综合997久久,| 国产一级毛片在线| 高清午夜精品一区二区三区| 久久久精品欧美日韩精品| 国产精品伦人一区二区| 亚洲电影在线观看av| 男女国产视频网站| 尤物成人国产欧美一区二区三区| 日韩av不卡免费在线播放| 免费看美女性在线毛片视频| 国产在线一区二区三区精 | 在线观看一区二区三区| 日本欧美国产在线视频| 伦精品一区二区三区| 婷婷六月久久综合丁香| 国产片特级美女逼逼视频| 国产精品不卡视频一区二区| 九九久久精品国产亚洲av麻豆| 2021少妇久久久久久久久久久| 麻豆一二三区av精品| 亚洲欧美日韩无卡精品| 一级二级三级毛片免费看| 亚洲国产精品成人久久小说| 久久亚洲国产成人精品v| 3wmmmm亚洲av在线观看| 久久久久网色| 日韩成人伦理影院| 国产精品国产三级国产av玫瑰| 三级毛片av免费| 天天一区二区日本电影三级| 久久人人爽人人爽人人片va| 韩国高清视频一区二区三区| av福利片在线观看| 中国国产av一级| 18禁动态无遮挡网站| 全区人妻精品视频| 69人妻影院| 国产综合懂色| 午夜福利在线在线| 嫩草影院入口| 免费人成在线观看视频色| 久久国产乱子免费精品| 乱系列少妇在线播放| 搞女人的毛片| 久久久久久伊人网av| 乱系列少妇在线播放| 九草在线视频观看| 狂野欧美激情性xxxx在线观看| 18禁动态无遮挡网站| 午夜福利在线观看免费完整高清在| 色综合站精品国产| 午夜久久久久精精品| 亚洲精品日韩在线中文字幕| 午夜免费激情av| 18+在线观看网站| 少妇丰满av| 男人和女人高潮做爰伦理| 国产精品不卡视频一区二区| 欧美精品国产亚洲| 99热这里只有是精品在线观看| 亚洲av二区三区四区| 成年女人永久免费观看视频| 国产乱人偷精品视频| 日韩一本色道免费dvd| av视频在线观看入口| 亚洲在线自拍视频| 男人舔女人下体高潮全视频| 中文字幕精品亚洲无线码一区| 亚洲av二区三区四区| a级毛片免费高清观看在线播放| 九九爱精品视频在线观看| 天天躁夜夜躁狠狠久久av| av卡一久久| 美女被艹到高潮喷水动态| 亚洲精品成人久久久久久| 国产亚洲精品久久久com| 日韩成人伦理影院| 99在线人妻在线中文字幕| 亚洲精品456在线播放app| 国产成人精品久久久久久| 日本爱情动作片www.在线观看| 丝袜喷水一区| 美女大奶头视频| 亚洲四区av| 久久精品国产亚洲av涩爱| 国产欧美另类精品又又久久亚洲欧美| 国产久久久一区二区三区| 久久热精品热| 91精品伊人久久大香线蕉| 国产精品野战在线观看| 国产精品,欧美在线| 青春草国产在线视频| 欧美区成人在线视频| 久久综合国产亚洲精品| 国产激情偷乱视频一区二区| 丰满人妻一区二区三区视频av| 久久99热6这里只有精品| 精品国产露脸久久av麻豆 | 好男人视频免费观看在线| 小蜜桃在线观看免费完整版高清| 亚洲自偷自拍三级| 亚洲图色成人| 婷婷色麻豆天堂久久 | 日本免费一区二区三区高清不卡| 亚洲va在线va天堂va国产| 熟妇人妻久久中文字幕3abv| 国产一区二区在线av高清观看| 99久久人妻综合| 久久欧美精品欧美久久欧美| 好男人视频免费观看在线| 日日摸夜夜添夜夜添av毛片| 国产黄a三级三级三级人| 欧美极品一区二区三区四区| 少妇人妻一区二区三区视频| 人妻少妇偷人精品九色| 全区人妻精品视频| 女人十人毛片免费观看3o分钟| 国产欧美日韩精品一区二区| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| 国产三级中文精品| 好男人在线观看高清免费视频| 国产白丝娇喘喷水9色精品| 国产精品精品国产色婷婷| 我的女老师完整版在线观看| 干丝袜人妻中文字幕| 亚洲欧美精品自产自拍| 国产精品久久久久久av不卡| 黄色一级大片看看| 我要看日韩黄色一级片| 国产在线男女| 国产精品一及| 国产91av在线免费观看| 五月伊人婷婷丁香| 亚洲精品456在线播放app| 91午夜精品亚洲一区二区三区| 亚洲综合精品二区| 亚洲经典国产精华液单| 91久久精品国产一区二区成人| 久久久久久久久中文| 免费不卡的大黄色大毛片视频在线观看 | 内地一区二区视频在线| 黄片无遮挡物在线观看| 99热6这里只有精品| 亚洲国产高清在线一区二区三| 久久99热这里只频精品6学生 | 国产欧美另类精品又又久久亚洲欧美| 高清在线视频一区二区三区 | 国产精品国产三级国产专区5o | 少妇熟女欧美另类| 纵有疾风起免费观看全集完整版 | 简卡轻食公司| 美女脱内裤让男人舔精品视频| 一个人免费在线观看电影| 久久精品熟女亚洲av麻豆精品 | 久久精品综合一区二区三区| 一级毛片aaaaaa免费看小| 最近最新中文字幕免费大全7| 亚洲欧美精品综合久久99| 日韩精品青青久久久久久| 禁无遮挡网站| h日本视频在线播放| 一级二级三级毛片免费看| 可以在线观看毛片的网站| 免费电影在线观看免费观看| 乱码一卡2卡4卡精品| 精品国产三级普通话版| 中文字幕制服av| 国产av在哪里看| 熟女人妻精品中文字幕| 天天躁夜夜躁狠狠久久av| 久久久色成人| 欧美日韩综合久久久久久| 久久99热这里只频精品6学生 | 九色成人免费人妻av| 国产在视频线在精品| 亚洲欧美成人综合另类久久久 | 色吧在线观看| 亚洲精品一区蜜桃| 热99在线观看视频| 亚洲av电影不卡..在线观看| 五月伊人婷婷丁香| 中文字幕亚洲精品专区| 日本色播在线视频| 免费看日本二区| 亚洲精品久久久久久婷婷小说 | 九九在线视频观看精品| 一本久久精品| 黄色日韩在线| 69av精品久久久久久| 麻豆成人av视频| 我要看日韩黄色一级片| 国产伦理片在线播放av一区| 久久久精品欧美日韩精品| 天堂av国产一区二区熟女人妻| 国产成年人精品一区二区| 色尼玛亚洲综合影院| 国内揄拍国产精品人妻在线| 亚洲美女搞黄在线观看| 精品国产三级普通话版| 久久久久久久国产电影| 亚洲怡红院男人天堂| 蜜桃久久精品国产亚洲av| 亚洲国产精品sss在线观看| a级毛色黄片| 波多野结衣高清无吗| 成人性生交大片免费视频hd| 男女啪啪激烈高潮av片| 三级毛片av免费| 日本免费一区二区三区高清不卡| 午夜激情欧美在线| 亚洲精品日韩在线中文字幕| 国产亚洲一区二区精品| 97人妻精品一区二区三区麻豆| 欧美丝袜亚洲另类| 只有这里有精品99| 欧美丝袜亚洲另类| 免费观看a级毛片全部| 免费看日本二区| 少妇人妻一区二区三区视频| 老司机福利观看| 美女大奶头视频| 蜜桃亚洲精品一区二区三区| www日本黄色视频网| 成人二区视频| www日本黄色视频网| 18禁动态无遮挡网站| 精品久久久久久久久亚洲| 超碰97精品在线观看| 欧美日韩在线观看h| 桃色一区二区三区在线观看| 91aial.com中文字幕在线观看| 久久精品国产亚洲网站| 国产 一区精品| 99久国产av精品| 精华霜和精华液先用哪个| 久久久国产成人精品二区| 亚洲精品自拍成人| 少妇熟女欧美另类| 91av网一区二区| 男人狂女人下面高潮的视频| 亚洲精品成人久久久久久| 中文字幕亚洲精品专区| 免费av毛片视频| 国产大屁股一区二区在线视频| 亚洲一级一片aⅴ在线观看| 国产综合懂色| 欧美精品国产亚洲| 看黄色毛片网站| 男女国产视频网站| 熟妇人妻久久中文字幕3abv| 免费搜索国产男女视频| 国产高清不卡午夜福利| 2021少妇久久久久久久久久久| 国产免费福利视频在线观看| 狂野欧美白嫩少妇大欣赏| 汤姆久久久久久久影院中文字幕 | 级片在线观看| 五月玫瑰六月丁香| 少妇的逼水好多| 国产精品精品国产色婷婷| 久久久久久久久久成人| 毛片一级片免费看久久久久| av.在线天堂| 国产老妇伦熟女老妇高清| 欧美一区二区精品小视频在线| 最近的中文字幕免费完整| 日韩精品青青久久久久久| 日本黄色片子视频| 精品久久久久久久久亚洲| 亚洲国产精品久久男人天堂| 亚洲精品色激情综合| 国产精品电影一区二区三区| 日韩强制内射视频| 伊人久久精品亚洲午夜| 毛片一级片免费看久久久久| 人体艺术视频欧美日本| 成人性生交大片免费视频hd| 日本与韩国留学比较| 免费看av在线观看网站| 久久精品久久久久久噜噜老黄 | 九九热线精品视视频播放| 国产又色又爽无遮挡免| 狠狠狠狠99中文字幕| 国产精品国产高清国产av| 97热精品久久久久久| 亚洲第一区二区三区不卡| 日韩亚洲欧美综合| 日产精品乱码卡一卡2卡三| 超碰97精品在线观看| 免费播放大片免费观看视频在线观看 | 网址你懂的国产日韩在线| 性色avwww在线观看| 午夜福利在线观看免费完整高清在| 精品一区二区三区视频在线| 欧美日本亚洲视频在线播放| 国产免费福利视频在线观看| 国产成人aa在线观看| 只有这里有精品99| 日韩av不卡免费在线播放| 欧美xxxx性猛交bbbb| 两个人的视频大全免费| 丰满少妇做爰视频| 午夜免费激情av| 午夜免费男女啪啪视频观看| 精品久久久久久电影网 | 夜夜看夜夜爽夜夜摸| 最近中文字幕高清免费大全6| 看黄色毛片网站| 变态另类丝袜制服| 丝袜喷水一区| 久久久成人免费电影| 如何舔出高潮| 亚洲欧美日韩无卡精品| 青春草国产在线视频| 91午夜精品亚洲一区二区三区| 99在线视频只有这里精品首页| 欧美日韩一区二区视频在线观看视频在线 | 久久精品久久久久久久性| 两个人的视频大全免费| 男人的好看免费观看在线视频| 免费av不卡在线播放| 69av精品久久久久久| 亚洲最大成人中文| 最后的刺客免费高清国语| 麻豆久久精品国产亚洲av| 欧美xxxx性猛交bbbb| 亚洲成人中文字幕在线播放| 日本一本二区三区精品| 淫秽高清视频在线观看| 九草在线视频观看| 国产精品嫩草影院av在线观看| 日韩av在线免费看完整版不卡| 成人鲁丝片一二三区免费| 成人国产麻豆网| 天堂网av新在线| 综合色av麻豆| 黄色一级大片看看| 亚洲欧美日韩高清专用| 少妇人妻一区二区三区视频| 国产色爽女视频免费观看| 久久久久久大精品| 岛国毛片在线播放| 好男人在线观看高清免费视频| 国产三级在线视频| 午夜久久久久精精品| 97在线视频观看| 国产探花极品一区二区| 一本久久精品| 免费看a级黄色片| 国产在线一区二区三区精 | 国产高清三级在线| 99九九线精品视频在线观看视频| 小说图片视频综合网站| 亚洲人与动物交配视频| 午夜a级毛片| 免费人成在线观看视频色| 国产黄色小视频在线观看| 夜夜看夜夜爽夜夜摸| 日韩一区二区三区影片| 草草在线视频免费看| 亚洲三级黄色毛片| 狠狠狠狠99中文字幕| 国产成人精品久久久久久| 啦啦啦啦在线视频资源| 国产精品国产三级国产av玫瑰| 91精品一卡2卡3卡4卡| 免费搜索国产男女视频| 精品人妻偷拍中文字幕| 亚洲av不卡在线观看| 97超视频在线观看视频| 伦精品一区二区三区| 国产在视频线精品| av播播在线观看一区| 最近最新中文字幕大全电影3| 欧美性猛交黑人性爽| 亚洲国产精品合色在线| 久久久色成人| 精品久久久久久久久亚洲| 欧美97在线视频| 亚洲av电影在线观看一区二区三区 | 直男gayav资源| 久久久精品94久久精品| 免费av不卡在线播放| 最近的中文字幕免费完整| 草草在线视频免费看| 赤兔流量卡办理| 色播亚洲综合网| 天天躁夜夜躁狠狠久久av| 内地一区二区视频在线| 九九热线精品视视频播放| av天堂中文字幕网| 少妇人妻精品综合一区二区| 草草在线视频免费看| 中文字幕av在线有码专区| 久久6这里有精品| 欧美一区二区精品小视频在线| 亚洲人成网站在线播| 久久久久久国产a免费观看| 国产黄片美女视频| 国产中年淑女户外野战色| 精品一区二区免费观看| 日韩成人伦理影院| 国产亚洲av嫩草精品影院| 亚洲国产精品sss在线观看| 精品久久久噜噜| 午夜亚洲福利在线播放| 精品熟女少妇av免费看| 女人被狂操c到高潮| 日韩欧美三级三区| 不卡视频在线观看欧美| 简卡轻食公司| 久久这里只有精品中国| 97超碰精品成人国产| 国产在视频线精品| 日本色播在线视频| 国产欧美日韩精品一区二区| 美女内射精品一级片tv| 亚洲真实伦在线观看| 十八禁国产超污无遮挡网站| 免费大片18禁| 我的老师免费观看完整版| 国产一级毛片在线| 精品久久国产蜜桃| 成人亚洲欧美一区二区av| 麻豆av噜噜一区二区三区| 国产高清视频在线观看网站| 黄色配什么色好看| 午夜激情福利司机影院| 在线免费观看的www视频| 国产成人freesex在线| 亚洲欧美日韩高清专用| 久久亚洲精品不卡| 亚洲国产日韩欧美精品在线观看| 亚洲四区av| 少妇丰满av| 看黄色毛片网站| 亚洲欧美清纯卡通| 精品人妻视频免费看| 夜夜爽夜夜爽视频| 欧美成人a在线观看| 自拍偷自拍亚洲精品老妇| 国产免费又黄又爽又色| 少妇猛男粗大的猛烈进出视频 | 22中文网久久字幕| 精品人妻一区二区三区麻豆| 欧美性感艳星| 永久网站在线| 亚洲精品乱码久久久v下载方式| 亚洲av电影在线观看一区二区三区 | 久久久久久久亚洲中文字幕| 99国产精品一区二区蜜桃av| 午夜福利在线在线| 国产免费又黄又爽又色| 亚洲五月天丁香| 欧美日韩在线观看h| av专区在线播放| 18禁在线无遮挡免费观看视频| 午夜福利在线在线| 国产色婷婷99| 成人午夜高清在线视频| 久久久亚洲精品成人影院| 观看免费一级毛片| 69av精品久久久久久| 只有这里有精品99| 国产探花在线观看一区二区| 欧美成人一区二区免费高清观看| 深夜a级毛片| 欧美变态另类bdsm刘玥| 免费观看的影片在线观看| 亚洲欧洲国产日韩| 亚洲国产日韩欧美精品在线观看|