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

    Close-in weapon system planning based on multi-living agent theory

    2022-07-27 06:46:18TangTangYueWangLijuanJiaJinHuChengMa
    Defence Technology 2022年7期

    Tang Tang, Yue Wang, Li-juan Jia, Jin Hu, Cheng Ma

    School of Information and Electronic, Beijing Institute of Technology, Beijing,100081, China

    Keywords:Multi-living agent theory Close-in weapon system Effectiveness evaluation Livelihood degree Queueing system Missile-gun integrated weapon system

    ABSTRACT The close-in weapon system(CIWS)is a combat system that faces a complex environment full of dynamic and unknown challenges, whose construction and planning require a systematic design method. Multiliving agent (MLA) theory is a methodology for the combat system design, which uses the livelihood degree to evaluate the multi-dimensional long-term operational effectiveness of the system; whereas,there is still no uniform quantization framework for the livelihood degree,and the adjustment methods of livelihood degree need to be further improved. In this paper, we propose the uniform quantization framework for the livelihood degree and detailed discuss the methods of livelihood adjustment.Based on the MLA theory, the multi-dimensional operational effectiveness of the missile-gun integrated weapon system (MGIWS) is analyzed, and the long-term combat effectiveness against the saturation attack is assessed. Furthermore, the planning problem of the equipment deployment and configuration is investigated. Two objectives, including the overall livelihood degree and cost-effectiveness (CE), are proposed, and the optimization method based on genetic algorithm (GA) is studied for the planning problem.

    1. Introduction

    In modern warfare,high-tech equipment has been put into use,which has gradually changed the form of war into joint operations based on information systems.Combat effectiveness is frequently a topic of significant interest in military combat operations. The methods of assessing equipment’s effectiveness are mainly based on the concept of the effectiveness index to represent the effect of weapons on a wide variety of targets[1].Currently,there are many combat effectiveness evaluation methods for specific equipment or systems, such as detection equipment [2-4], strike equipment[5-7], communication system [8-10], and so on. However, it is difficult to configure and deploy the system only depending on the effectiveness indices of single equipment or system to enable the system to carry out combat tasks efficiently and reliably in a complex and dynamic operational environment. In order to meet the combat requirements in a specific operational scenario, it is necessary to conduct a systematic effectiveness evaluation for the system in a specific environment. At present, there are some categories for evaluating system combat effectiveness, including expert evaluation [11], analytical methods [12-14] and combat simulation methods [15-18].

    A close-in weapon system (CIWS) is a point-defense combat system to detect and destroy short-range incoming missile and enemy aircraft which have penetrated the outer defenses [19].There are some researches that study the combat effectiveness of CIWS.Yang et al.[20]analyzed the combat effectiveness of airborne tactical ballistical anti-ship missile against the CIWS,which mainly focused on the relation between the speed of missile and the penetration capability. Zhang et al. [21] used Monte Carlo simulation to analyze the interception effectiveness of CIWS on different prototypes of hypersonic missile(HM)and explored the feasibility of using CIWS to intercept HM and the improvement measures. In the context of distributed combat[22],saturation attack formed by a large number of unmanned aerial vehicles (UAVs) with strike capability might be a major challenge for the CIWS, which has the advantages of strong autonomy,low price,high degree of freedom and rapid upgrade[23].Considering the interception effectiveness against the saturation attack,queueing theory can be used to model the combat process of CIWS.Liu et al.[24]synthetically researched the calculating problem about the penetration probability of antiship missile with the queueing theory under the condition that surface ships formation adopt various means to help to oppugn mutually.Zheng et al.[25]analyzed the penetration effectiveness of anti-ship missile and stealth airplane to the “Phalanx” system and the influence of flying parameter of aircraft. Zhao et al. [26]calculated the combat effectiveness of the “Sidewinder” and “Tor-M1” defense missile weapon system with the queueing theory.Previous effectiveness evaluations of CIWS against the saturation attack mostly used the single-layer queuing model to analyze the interception probability of surface-to-air missile (SAM) or antiaircraft gun (AAG) systems, rather than that of the missile-gun integrated weapon system (MGIWS), which should be further explored.Besides,the studies did not take the states and reliability of the system under long-term combat situation into consideration,which is also important effectiveness should be considered in the combat. For example, in the Gulf War, ”Patriot” failed to intercept the”Scud”missile because of the failure caused by the long hours of continuous operation [27].

    For the design and planning of the combat system under the strongly constrained and hostile environment, Wang [28,29] proposed the Multi-living agent (MLA) system theory. Livelihood degree was studied to evaluate the multi-dimensional continuous effectiveness of the system.It is a new method for research,design,and analysis of a combat system. MLA method was established based on the dissipative self-organization theory, and the humandominated heter-organization was emphasized to maintain the effectiveness of the system[30].Based on the MLA system theory,a series of studies have been carried out.Wang et al.[31]studied the distributed intrusion detection system based on MLA with the queueing theory. It was shown that the security and reliability of the proposed system are effectively enhanced, and the agent can normally work in a strongly constrained and hostile environment.Liu et al.[32]proposed an analytic hierarchy process based method with a dependent pairwise comparison matrix for the evaluation of the livelihood degree of the MLA system. An application to the netted radar system was illustrated to support the proposed method. Tang et al. [33] put forward the overall design method of the high-performance artificial system based on MLA, which took the functional livelihood degree as the basis for research and development. Zhao et al. [34] studied the architecture adjustment method based on MLA theory. By replacing the principal-agent or adding a new living agent, the system architecture is adjusted to maintain the system’s functional livelihood to meet the task requirements. Although many studies have been proposed on the MLA theory, there is still no uniform standard for the quantitative characterization of the livelihood degree. The mechanism of livelihood adjustment needs to be summarized and explained.

    In this paper, we propose a unified framework for the assessment of the livelihood degree and discuss the methods for adjusting livelihood. On this basis, we use the two-layer queuing model to evaluate the interception probabilities of a CIWS consists of multiple joint operated MGIWSs. Considering the equipment failure and maintenance, we further analyze the long-term operational livelihood degree in different configuration strategies. In a multi-point deployment scenario, we propose two objectives to plan different configuration schemes: combat effectiveness preferred and cost-effectiveness (CE) preferred. The optimization method based on the genetic algorithm (GA) is proposed to solve the planning problem.

    The rest of this paper is arranged as follows: In Section 2, the uniform quantization framework for the livelihood degree and the adjustment methods of livelihood degree are proposed. In Section 3, we propose the long-term operational livelihood degree assessment method of MGIWS under different configuration strategies. In Section 4, for the planning of the battlefield equipment deployment and configuration, we propose two strategies with different optimization objectives and use GA to solve them. In Section 5,numerous simulations are carried out to show the result of the proposed method. The conclusion in Section 6 presents the main innovation of the proposed method and some further work.

    2. MLA system theory

    2.1. The define of living agent and livelihood degree

    In the define of MLA, the term”agent”derives from the field of computer software design, which indicates an autonomous entity that carries out tasks on behalf of users. ”Livelihood” means the ability to maintain specific service functions. The ”living agent”combines the human intentions and the structure of systems in terms of livelihood to achieve a particular functional purpose. The being of a”living agent”is reflected in its continuous effectiveness in the system. The being has two forms: the tangible form of the structure which is constructed by hardware and software, and the intangible form which is composed of man’s related concepts and the mutual relations they formed. The livelihood can be quantitatively characterized by the livelihood degree F [30], that is, the possibility to complete a specified task.

    According to Ref. [29], the six-element relationship group of a living agent is shown as follows to characterize the dynamic characteristics:

    where P represents the set of service probability in all the possible cases.The MLA theory was proposed to design and analyze systems in a complex and hostile environment.Therefore,in the evaluation of livelihood degree,the service probability in the ideal state is not appropriate to be used as an evaluation standard, because the antagonistic factors might be designed to counter the weakness of the system. Thus, we use the long-term service probability in the rigorous case to evaluate the livelihood degree.

    2.2. The structure of MLA system

    Fig.1. The structure of the MLA system.

    The structure of the MLA system is shown in Fig.1. In order to meet the real-time scheduling requirements of the system, the living agents should have a unified interface and should be universal and tailorable to support efficient resource interaction and collaboration. The MLA system consists of three layers, where the underlying layer consists of the multiple functional living agents,the middle layer is the input and output interface layer,and the high layer is the management layer backstage.The living agents at the underlying layer perform perception, processing, state adjustment and output control tasks according to the selforganizing mechanism of the system. Living agents cooperate through interactions to achieve overall requirements. The inputoutput agent is responsible for the interaction between the MLA system and other systems and environment, including the input and output of information, substances, and energy. The management and control (MC) agent is responsible for MLA’s information collection, information processing, and overall control.

    In the case that the self-organization mechanism of the system is able to meet the functional requirement, the MC agent is responsible for the negotiation and coordination among the living agents in the system, which is based on the logic and the expert systems. When the self-organization mechanism of the system cannot meet the functional requirements, the heter-organization should be implemented under the guidance of humans, such as adjusting the system decision-making strategy,adding or replacing agents, modifying the system architecture. In a complex environment full of strong constraint and confrontation, to meet the dynamic and complex demands, the MLA system should be able to maintain and adjust the living agents through the self-organization of the MC agent and the high-level heter-organization.

    2.3. The livelihood adjustment methods

    As shown in Fig.2,there are two levels for the adjustment of an MLA system:agent level adjustment and system level adjustment.

    2.3.1. Agent level adjustment

    When a weak link occurs and leads to the decline of the overall performance,the agent level adjustment is often first considered to maintain the overall livelihood of the system.

    Replacing agent is the most straightforward approach to make up for the weak link;however,this method has a high requirement on the capability of the new agent, which might be difficult to achieve due to technical constraints.

    Except for replacing the agent directly, we can consider adding agents to improve the livelihood. The advantage of adding living agents is that it can satisfy the service requirement by a combination of agents available, while the disadvantage is that it increases the cost and the redundancy of the system.The following types of agent can be added to improve the livelihood of the original agent.

    (1) Serial agent: When a type of task can be accomplished step by step and the original agent can only achieve one of the steps, a serial agent can be added to pre-process or postprocess the task. The service probability of this link is the probability that both the two agents operate successfully.

    (2) Homogeneous parallel agent:When a weak link is due to the high task intensity or low service probability and the original agent has no function dimension loss,a homogeneous agent can be paralleled to increase livelihood. This link is failed only if all the agents failed to provide service.

    (3) Heterogeneous parallel agent: When there is a defect that occurs in the original agent’s function dimensions due to the change of environment or requirement, we can consider adding a heterogeneous parallel agent with complementary function dimensions to enhance the livelihood.

    (4) Livelihood maintenance agent: When an agent will be degraded or broken due to long-term work, a livelihood maintenance agent can be added to monitor and maintain the state of the agent.This method is effective in maintaining the long-term effectiveness of the agent.

    2.3.2. System level adjustment

    In the case that the inadequacy of system livelihood is not caused by a single agent, the system-level adjustment should be considered.

    If the system could not meet the new requirements under the previous system architecture, the method of architecture adjustment can be adopted. By reconfiguring the system structure,operation processes,coordination strategies to enable the system to fulfill the task requirements in the target environment.

    In the case that the system architecture can meet the requirements,but the livelihood is difficult to reach the optimum due to the changeable environment and tasks, a variety of operational modes can be set up to realize the switching between modes based on the dynamic perception and evaluation of the environment and tasks.

    Fig. 2. The livelihood adjustment methods.

    In the next section, the livelihood degree and the adjustment methods of CIWS are investigated on the proposed method.

    3. Livelihood degree assessment of CIWS

    3.1. Background and architecture of the CIWS

    The mission of CIWS is to defend the important strategic location, which is the last line of defense to intercept enemy attacks.The target types are various aircraft and air-to-surface missiles(ASMs). In the case that the enemy adopts a full range of multibatch intensive attack, CIWS needs to perceive the dynamic information and situation of the enemy’s intrusion flow in time and coordinate the subsystems to carry out the effective interception.

    The main requirements of the system can be divided into several sub-requirements, including searching, tracking, destroying the target, computing support, commanding and controlling, humancomputer interaction. The sub-requirements can be achieved through corresponding functional agents and the system architecture is shown in Fig. 3.

    The searching agent is composed of the searching radars,which searches for the aircraft within the detection range and reports target information to the subsequent agents.

    The tracking agent is responsible for continuously tracking after the searching agent detects the target. It complements the effectiveness dimension of radar and optical equipment through the heterogeneous parallel livelihood adjustment method. The radar has a strong ability to track moving targets and is less affected by visibility. Optical detection equipment has the advantages of high concealment, strong passive detection ability and strong tracking ability to the extremely slow-moving and static targets.Besides,the optical equipment has the imaging ability,which makes up for the shortcomings of the radar in target identification and greatly improves the reconnaissance and monitoring effects.

    The firing agent receives the firing command of the MC agent and strikes the enemy’s invading targets.This system considers the parallel architecture of SAM and AAG,which is regarded as MGIWS.AAG has the advantages of high initial velocity, high rate of fire,high close-range hit accuracy, high probability of damage, short response time,strong anti-interference ability and low combat cost,which is responsible for the super-close-range interception mission. SAM has a high single-shot kill probability and a large strike range, and the kill probability is basically unchanged within the range.It can be responsible for the interception mission outside the range of AAG. As the strike range shown in Fig. 4, through the heterogeneous parallel livelihood adjustment method of the SAM agent and AAG agent,the combat effectiveness of the firing agent in the space dimension can be supplemented.

    Fig. 3. The architecture of CIWS.

    Fig. 4. The strike range of MGIWS.

    The computing agent and the MC agent are responsible for the collection, processing and distribution of information in the system. They support the commander’s unified management and control of the system.

    Considering the high effectiveness requirement of CIWS, multiple MGIWSs can be combined for joint operation,which enhances the livelihood degree through the homogeneous parallel method.

    3.2. Interception probability assessment

    For the CIWS described in Section 3.1, the interception probability can be evaluated as p(q,s,c), where sis the input information of the enemy’s targets that might be mixed with interference,c is the combat environment of CIWS,q is the operation state of CIWS. Besides, sis the interception flow output by the system. The prerequisite for the successful destruction of the target is that each functional sub-agent of CIWS successfully completes the task. Therefore, considering the process of the interception mission, the interception probability is

    where pis the searching probability of a single searching radar.

    pis the probability that the system successfully track the detected target and meet the fire control command requirements.

    pis the probability that one of the SAM or AAG strike and destroy the tracked target, that is,

    pis the probability that the SAM agent fires and destroys the target in the SAM strike range,that is,

    where pis the probability that AAG agent can fire at a tracked target that enters the strike range of AAG agent,and pis the kill probability of a round of AAG fire.

    Among the above parameters, most of them can be calculated according to the equipment parameters; whereas, the fire probabilities of SAM and AAG need to be evaluated combining with the specific attack flow intensity and the number of parallel systems.In section 3.3, the destruction probability is evaluated with the queueing model.

    3.3. Evaluation of the destruction probability

    Due to the uncertainty of military confrontation, we consider the situation that we cannot predict the follow-up situation based on the current situation and take corresponding countermeasures.Therefore, a model without aftereffect is used to simulate the interception process of CIWS. Construct a two-layer impatient M/M/n queuing model as shown in Fig. 5 and make the following assumptions:

    First,consider the first layer of the queuing model,which is the combat process of SAM. The arrival rate of targets is the rate that targets are detected and tracked, that is,

    According to Ref.[33],with the arrival rate ω,the service rate μ,the leaving rate υ and the number of servers n, the service probability of the impatient M/M/n queueing model can be calculated as

    Fig. 5. The two-layer M/M/n queueing model.

    Similarly, substituting the parameters ω, μ, υand n into Eq. (12), the fire probability of AAG can be obtained as p.

    As the fire probabilities pand pcan be obtained, the destruction probability pcan be calculated according to Eqs.(7)-(9). Further, the interception probability can be obtained via Eq. (5).

    3.4. Long-term livelihood degree assessment

    The previous sections studied the interception probability p(q,s,c)of the joint operated MGIWSs in the face of the input sunder a specific environment c. To further consider the long-term operational livelihood degree of the system, it is also necessary to consider the change of the system state q,including the equipment failure and maintenance.Consider that the system consists of N sets of MGIWS for the strike and M sets of maintenance equipment for the maintenance,which adjusts the long-term livelihood degree of the system by means of adding livelihood maintenance agents.

    3.5. Make the following assumptions

    (1) Each operating MGIWS is subjected to deterioration over time, and the time to failure follows an exponential distribution with the failure rate η.

    (2) Each maintenance equipment can repair an MGIWS in failure, and multiple maintenance systems cannot carry out coordinated maintenance on a single failed MGIWS. The maintenance time follows an exponential distribution with the maintenance rate γ.

    (3) MGIWS can stay in the sleep state, that is, it does not participate in the combat missions and will not turns into a failure state.

    Define qto represent the state of CIWS that there are n MGIWSs in normal state,where 0 ≤n ≤N,and the state set of the system is Q = {q,q,…,q,q}.

    According to the principle of diminishing marginal utility [35],when the number of paralleled MGIWSs exceeds a certain number,more working agents might not be able to greatly improve combat effectiveness but increase unnecessary losses. We can consider an operational strategy to avoid the unnecessary wastage of the system during long-term combat.When the number of MGIWSs in the normal state is greater than K,where 1 ≤K ≤N,the extra MGIWSs can stay in the sleep state, and only K MGIWSs will participate in the combat missions.

    Therefore, we can define two parameters to represent the numbers of operating MGIWSs and maintenance equipment at state qas in Eqs. (15) and (16), respectively.

    Consequently, the forward Chapman-Kolmogorov (C-K) equations can be written as follows

    where

    Fig. 6. Transition diagram of CIWS states.

    Furthermore, taking the limit of Eq. (27) as t→∞, the steadystate livelihood degree of CIWS can be obtained as F(N,M,K).

    This section considers the long-term combat livelihood degree assessment method of a CIWS under a certain configuration strategy, and the next section will study the overall deployment and configuration strategy planning method of CIWSs on the battlefield.

    4. Planning method of the deployment and configuration strategy

    4.1. Planning problem description

    On the battlefield,CIWSs are needed to be deployed at multiple strategic locations.With the limit of the amount of equipment,it is an important issue to optimize the overall combat effectiveness of the battlefield through the planning method. Assume that the numbers of MGIWSs and maintenance equipment that can be deployed on the battlefield are Nand M,respectively.CIWSs must be deployed at I strategic locations, and the CIWS at each location are independent of each other.

    The strategic locations can be divided into G important grades.For the i-th strategic location,its importance grade is g?{1,2,…,G},the maximum target intensity is ω,the interception probability requirement is W, and the long-term livelihood degree requirement is D.

    The purpose of the battlefield commander is to make the CIWS of each strategic location meet the interception requirement and optimize the overall combat effectiveness through rational equipment deployment and operational strategy planning under the constraint of equipment amount.In this study,we achieve the goal of overall optimization by adjusting the number of MGIWSs N,the number of maintenance equipment Mand the maximum operational number Kat each location.

    Note that,we add a subscript i to a parameter to indicate that it belongs to the ith location.

    4.2. Optimization objective

    Differ from the general equipment deployment and configuration problems that mainly consider the overall performance of the system in the ideal state[38-40],the optimization objective of this planning problem is to optimize the long-term overall combat effectiveness. We propose two objectives to achieve this goal,including the weighted arithmetic mean of livelihood degree and the CE.

    For the first case, we only consider the overall interception effectiveness and use the weighted sum function of interception effectiveness livelihood degree as the objective,which is shown as follows

    In addition to the combat effectiveness,the cost may also be an important dimension that needs to be considered in the decision.We can consider the overall combat CE as the objective rather than merely combat effectiveness.

    Assume that the operating costs of the operating MGIWS and the maintenance equipment are cand c, respectively. The expected operation cost of location i is

    4.3. Solution to the planning problem

    To plan the deployment and operation strategy to implement overall optimization, the planning problem is given by

    Constraint(34)determines that the long-term livelihood of each CIWS should satisfy the livelihood degree requirement of each location. Constraints (35) and (36) indicate that the deployment must be executed with limited MGIWSs and maintenance equipment.Constraint(37)determines that the number of maintenance equipment should not exceed that of MGIWSs. The maximum operating number should be no more than the number of MGIWSs,which is expressed in constraint (38).

    In this article,the planning problem is solved with GA,which is widely used in operations research to deal with combinatorial optimization problems and generate global high-quality solutions.GA has been demonstrated to converge to the optimal solution for many diverse difficult problems, although optimality cannot be guaranteed.

    GA starts by randomly initializing a population of chromosomes in the first generation,and each chromosome represents a possible solution to the optimization problem.As shown in Fig.7,we express a chromosome as a row vector X. Each chromosome consists of 3I genes.The first I genes represent the numbers of MGIWSs deployed at different locations. The genes from Xto Xrepresent the numbers of deployed maintenance equipment,and the genes from Xto Xrepresent the maximum numbers of operating MGIWSs. The pseudocode of the proposed optimization algorithm is shown on the bottom of this page. The fitness function, i.e. the objective function in Eq.(29)or(32),is used to evaluate the quality of individuals. After calculating the fitness values of the chromosomes, existing individuals are recombined using cross-over, mutation and reproduction procedures,to obtain the next generation.GA terminates when a given number of generations Nis attained or the improvement of fitness value is not greater than a threshold ε.

    5. Simulation and analysis

    Several numerical simulations are conducted in this section to assess the performance of the proposed livelihood degree evaluation method and planning method.

    Simulation 1. Consider a combat scenario where a CIWS composed of MGIWSs faces the saturation attack.The specifications of a type of MGIWS are shown in Table 1. Countering the invading targets with maximum speed 500 m/s,the interception probability of a CIWS consists of n operating MGIWSs can be calculated via Eq.(5),and the required parameters can be derived from the indicators in Table 1.

    We assume that the SAM is responsible for interception tasks within the SAM strike range but outside the AAG strike range,and AAG fires after the target enters the AAG strike range.AAG fires at the target for 1 s each time, that is, 10 rounds per fire. Then the diversion fire rate of AAG is

    Algorithm 1 GA based on livelihood degree Input:f ?{f1,f2}: the objective function NP: the size of initial population NG: the max number of generations pC: the probability of crossover operation pM: the probability of mutation operation pR: the probability of recombination operation ε: the accuracy threshold Output: The value of X that optimizes f Method:1: g←0 2: Initializing Pg to NP randomly generated individuals 3: Fitness evaluation for each individual in Pg with f 4: while termination condition not met do 5: Select individuals from Pg 6: Recombine individuals 7: Mutate individuals 8: Fitness evaluation for each modified individual 9: Pg+1←new created individuals 10: g←g+ 1 11: return The best individual in Pg

    Based on the specifications in Table 1, the parameters required by the proposed model can be obtained as shown in Table 2.

    According to the proposed 2-layer M/M/n queueing model in Section 3, the interception probabilities of CIWS consists of n operating MGIWSs against different target intensities ωcan be calculated via Eq. (5), and the results are shown in Table 3 on the top of the next page. On this basis, the curves of the interception probability in different conditions can be obtained as shown in Fig. 8.

    It can be concluded that the increase of the intensity of the targets will result in a decrease in the fire probabilities of the firing agents and further reduce the interception probability of CIWS.Besides, the method of paralleling MGIWSs can be adopted to enhance the interception probability, which verifies the effectiveness of the homogeneous parallel livelihood adjustment method.

    Fig. 7. Chromosome representation of GA.

    Table 1 Specifications of an MGIWS.

    Table 2 Parameters of the proposed model.

    Simulation 2. On the basis of simulation 1, we consider the long-term interception mission livelihood degree of CIWS consists of the above-mentioned MGIWSs under different configuration strategies.

    In the situation that the intensity of the intrusion targets is ω=0.3 and the interception requirement of CIWS is W =0.85,the livelihood degree of the CIWS consists of N MGIWS and M maintenance equipment,with the maximum number of operating agents K, can be calculated via Eq. (27).

    There are numerous configuration strategies that can be adopted,and we choose the five strategies as shown in Table 4 to analyze the influence of different parameters on the livelihood degree.

    The curves of the livelihood degree in different configuration strategies over time are shown in Fig. 9 and the steady-state livelihood degree are shown in Table 4. It can be inferred that the livelihood degree gradually decreases over time and finally convergences to the steady-state. In the initial stage, the livelihood decreased because the failure rate of MGIWSs is greater than the maintenance rate.With the failure of MGIWSs and the operation ofmaintenance equipment, the failure and maintenance of MGIWS gradually reach a dynamic balance, and the system also tends to a steady-state.

    Table 3 Interception probabilities in different conditions.

    Comparing strategies 1 and 2, the increase of the paralleled MGIWSs enhances the livelihood degree of the system, which indicates the effectiveness of the homogeneous parallel method to improve the livelihood degree.

    Comparing strategies 2,3 and 5, it can be inferred that the increase of maintenance equipment can improve the livelihood degree, which implies that the method of adding livelihood maintenance agent can effectively adjust the livelihood degree of the system. Besides, according to strategy 5, the system will eventually fail without maintenance.

    Strategy 4 represents the situation where all normal MGIWSs operate without considering the sleep strategy. Comparing strategies 2 and 4, it can be inferred that rational sleep strategy can improve the livelihood degree of the system, which verifies the availability of operational mode adjustment method.

    Simulation 3. Consider the situation that there are I =6 strategic locations must be defensed by CIWSs. The maximum target intensity ω, the importance level g?{1,2,3}, the interception probability requirement W, and the livelihood degree requirement Dof each strategic location are shown in Table 5.

    The costs of the operation of a MGIWS and a maintenance equipment are 1 per hour and 3 per hour, respectively.

    To optimize the deployment and configuration strategy of the overall battlefield,the planning problem proposed in this paper can be adopted. The solution to the deployment and configuration problem can be represented by a chromosome X consists of 18 genes. The initial setting values of the parameters used in GA are:N= 200, N= 1000,P= 0.8, P= 0.01, P= 0.05, ε = 0.001.

    Four solutions are given in Table 6.

    Fig. 8. Interception probabilities.

    Table 4 Livelihood degree in different strategies.

    Fig. 9. Livelihood degree in different strategies.

    Table 5 Parameters of different strategic locations.

    The first solution is obtained via the proposed planning method based on the overall livelihood degree f, and the second one is obtained based on the overall CE f. For comparison, solution 3 is the configuration method that does not consider the sleep strategy on the basis of solution 1.Solution 4 is the solution to optimize the global ideal state performance, which can be expressed by the weighted arithmetic mean of the interception probabilities of multiple locations, that is,

    where pis the performance when all the MGIWSs deployed at the strategic locations i are operating normally,which can be calculated via Eq.(5).At the same time,without considering the sleep strategy, making as many MGIWSs as possible work can make the system in the ideal state perform optimally. Then allocating the maintenance equipment according to the number of MGIWSs at each location, solution 4 can be obtained.

    As shown in Fig. 10, the comparison of the overall livelihood degree,the overall cost,the overall CE,the number of locations that satisfied the livelihood degree requirements, and the ideal performance are presented.

    Comparing solution 1 and solution 2,it can be inferred that the overall livelihood of solution 1 is higher than that of solution 2.However, due to the higher operating cost, its CE is lower. These two solutions match their optimization objectives, respectively.Because the planning processes of the two strategies are carried out under the constraints of livelihood degree, their deployment and configuration strategies at all the 6 strategic locations can meet the long-term combat requirements.

    Compared with soulution1, solution 3 does not consider the sleep strategy proposed in this paper. It can be concluded that the operating cost of solution 1 is much lower than that of solution 3,and the livelihood degree is higher;whereas,its initial performance is lower than that of solution 3. Through the sleep strategy proposed in this paper, solution 1 keeps working with the minimum MGIWSs while satisfying the interception requirements, which reduces the operating cost of MGIWSs. Besides, the sleep strategy also reduces the maintenance pressure of maintenance equipment,which enables deteriorated MGIWS to be repaired in time. Due to the high operating intensity of solution 3, there are only two strategic locations that can meet the long-term operational requirements.

    Comparing solution 1 and solution 4,the difference between the proposed method and the traditional performance-oriented method can be concluded. Solution 4 can achieve optimal performance when the system is in perfect condition,while solution 1 can ensure optimal combat effectiveness when the system is in longterm operation. In long-term operations, only 2 of the 6 strategic locations of solution 4 can meet the long-term operational requirements.

    To analyze the computational effectiveness of the proposed method, mean square error (MSE) [41] is obtained through N=1000 independent GA optimization iterations. MSE can be calculated as

    where Xis the ideal solution, Xis the best solution of the g-th generation in the k-th iteration. Using GA to optimize the deployment and configuration strategy X with the global livelihood fas the objective. Solution 1 in Table 6 is the optimal solution X. The curve of MSE with the change of generation is shown in Fig.11.TheMSE of proposed method gradually converges to a low level with the increase of generation.

    Table 6 Solutions and performance of different strategies.

    Fig.10. Comparison of parameters of different strategies.

    Based on the above results, we can conclude that the proposed planning method can effectively optimize the deployment and configuration strategy of the equipment on the battlefield according to the long-term combat requirements.

    6. Conclusion

    In this paper, we propose a unified quantitative framework for the livelihood degree of the MLA system theory and the basic modes of livelihood adjustment. Based on the MLA theory, the effectiveness evaluation of CIWS and the optimization of battlefield equipment deployment and configuration strategies are carried out. This paper adopts a two-layer queuing model to evaluate the interception probability of MGIWSs in the joint operation.Considering the equipment failure and maintenance,the long-term combat effectiveness of CIWS under different configurations is evaluated with the livelihood degree.For the planning problem of equipment deployment and configuration on the battlefield, we propose two optimization objectives as overall livelihood and CE.Moreover,the optimization method based on GA is investigated to solve this planning problem.Finally,simulations are carried out to demonstrate the effectiveness of the proposed methods.

    Fig.11. MSE with the change of generation.

    On the basis of this article,there are some further studies can be done.In this paper,the interception probability is analyzed with the queueing model,and the service rate of the firing agent is obtained via the diversion fire time,which is the total time of the process.To consider the details of the process for more detailed evaluation,the method of modeling and simulation can be adopted. In the investigation of livelihood degree,the binary-state as normal and failure of the MGIWS is considered.In future work,the multi-state,such as the gradual degradation of equipment, can be studied to describe the transition of system states in more detail. Besides, this article mainly considers the deployment and configuration of equipment through prior information. In combat, the adjustment strategies based on the real-time changes of the situation should be further studied on the basis. Some parameters in this article are obtained through theoretical analysis, such as the theoretical searching probability of the paralleled radars,etc.In practical situations,these parameters should be analyzed in detail based on specific equipment and application environments.

    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.

    We would like to express our sincere thanks to the editors and reviewers for their constructive comments and suggestions which improved and enriched the presentation of the paper. Further, we are grateful to the experts who offered help to us.

    This work is supported by the Beijing Natural Science Foundation under contract number L191004 and the National Natural Science Foundation of China under contract number U1833203.

    在线观看免费午夜福利视频| 亚洲熟女毛片儿| 亚洲第一电影网av| 亚洲一区高清亚洲精品| 久久性视频一级片| 亚洲精品中文字幕一二三四区| 精品国产三级普通话版| 老司机在亚洲福利影院| 99精品在免费线老司机午夜| 在线看三级毛片| 久久久久久国产a免费观看| 999久久久国产精品视频| 免费在线观看视频国产中文字幕亚洲| 国产成人欧美在线观看| 桃红色精品国产亚洲av| 深夜精品福利| 熟女电影av网| 欧美日韩精品网址| 亚洲专区字幕在线| xxx96com| 午夜激情福利司机影院| 亚洲人成电影免费在线| 国产精品女同一区二区软件 | 老熟妇乱子伦视频在线观看| 真实男女啪啪啪动态图| 非洲黑人性xxxx精品又粗又长| 亚洲精品中文字幕一二三四区| 亚洲av日韩精品久久久久久密| 国产成人精品久久二区二区91| 免费看十八禁软件| 国产精品久久久久久人妻精品电影| 九九热线精品视视频播放| 成年女人永久免费观看视频| 精品国产三级普通话版| 黄色 视频免费看| 欧美精品啪啪一区二区三区| 国产午夜精品论理片| 久久久久精品国产欧美久久久| 国产真实乱freesex| 国产黄色小视频在线观看| 久久久水蜜桃国产精品网| 一a级毛片在线观看| 亚洲av成人不卡在线观看播放网| 亚洲av五月六月丁香网| 在线观看舔阴道视频| 老司机午夜十八禁免费视频| 91久久精品国产一区二区成人 | 日韩欧美在线二视频| 午夜免费激情av| 国产黄色小视频在线观看| 舔av片在线| 国产欧美日韩一区二区三| 午夜免费成人在线视频| 国产三级中文精品| 亚洲 欧美 日韩 在线 免费| 岛国在线观看网站| 国产1区2区3区精品| 伦理电影免费视频| 日本三级黄在线观看| 国产精品美女特级片免费视频播放器 | 国产精品1区2区在线观看.| 99久久精品一区二区三区| 久久伊人香网站| 久久精品国产99精品国产亚洲性色| 操出白浆在线播放| 亚洲精品久久国产高清桃花| 国产精品免费一区二区三区在线| 午夜a级毛片| 老汉色∧v一级毛片| 黄色日韩在线| 不卡av一区二区三区| 国产精品亚洲美女久久久| 精品久久久久久久久久免费视频| 成熟少妇高潮喷水视频| 欧美另类亚洲清纯唯美| 18禁黄网站禁片午夜丰满| 午夜精品久久久久久毛片777| 亚洲精品久久国产高清桃花| 国产黄色小视频在线观看| 久久中文字幕一级| 免费搜索国产男女视频| 99国产综合亚洲精品| 免费大片18禁| 久久99热这里只有精品18| 亚洲激情在线av| 国产精品免费一区二区三区在线| 久久欧美精品欧美久久欧美| 久久久久亚洲av毛片大全| 精品国产超薄肉色丝袜足j| 午夜两性在线视频| 少妇熟女aⅴ在线视频| 亚洲av成人精品一区久久| 国产av一区在线观看免费| 狂野欧美白嫩少妇大欣赏| 后天国语完整版免费观看| 久久精品91无色码中文字幕| 日日干狠狠操夜夜爽| av视频在线观看入口| 日韩欧美免费精品| 性色av乱码一区二区三区2| 欧美一区二区国产精品久久精品| 在线永久观看黄色视频| 美女免费视频网站| 亚洲成人中文字幕在线播放| 亚洲av成人av| 国产三级黄色录像| 国产精品自产拍在线观看55亚洲| 免费观看人在逋| 国产成人福利小说| 极品教师在线免费播放| 69av精品久久久久久| 亚洲性夜色夜夜综合| 亚洲九九香蕉| 免费看美女性在线毛片视频| 91在线精品国自产拍蜜月 | 国产精品亚洲一级av第二区| 国产亚洲av高清不卡| 午夜福利在线观看免费完整高清在 | 午夜精品在线福利| 青草久久国产| 亚洲国产色片| xxxwww97欧美| 亚洲欧美精品综合久久99| 天天一区二区日本电影三级| 成年人黄色毛片网站| 三级毛片av免费| 黄片大片在线免费观看| 国产淫片久久久久久久久 | 亚洲国产精品999在线| 日本免费一区二区三区高清不卡| 老汉色∧v一级毛片| 亚洲专区中文字幕在线| 亚洲片人在线观看| 九九久久精品国产亚洲av麻豆 | 免费av毛片视频| 黄片小视频在线播放| 国产精品久久久久久精品电影| 国产欧美日韩一区二区三| 婷婷六月久久综合丁香| 日韩欧美精品v在线| av中文乱码字幕在线| 大型黄色视频在线免费观看| 又粗又爽又猛毛片免费看| 大型黄色视频在线免费观看| 欧美日韩黄片免| 国内精品美女久久久久久| 久久人人精品亚洲av| av天堂在线播放| 在线观看一区二区三区| 亚洲av美国av| 99热只有精品国产| 欧美午夜高清在线| 美女黄网站色视频| 偷拍熟女少妇极品色| 成在线人永久免费视频| 99国产精品一区二区蜜桃av| 久久九九热精品免费| 特大巨黑吊av在线直播| 国产1区2区3区精品| 亚洲专区国产一区二区| 免费在线观看影片大全网站| 一级黄色大片毛片| 在线a可以看的网站| 日韩人妻高清精品专区| 亚洲精品在线观看二区| 日韩成人在线观看一区二区三区| 美女高潮喷水抽搐中文字幕| 国产精品久久视频播放| 精品午夜福利视频在线观看一区| 日韩大尺度精品在线看网址| 婷婷精品国产亚洲av| 亚洲国产高清在线一区二区三| 中文亚洲av片在线观看爽| 欧美不卡视频在线免费观看| 亚洲一区二区三区色噜噜| 男人舔奶头视频| 精品久久久久久,| 欧美在线黄色| 99热只有精品国产| 香蕉丝袜av| 亚洲欧美日韩高清专用| 女人被狂操c到高潮| 亚洲第一电影网av| 欧美绝顶高潮抽搐喷水| 亚洲熟女毛片儿| 亚洲精品一区av在线观看| 男女床上黄色一级片免费看| 一个人免费在线观看电影 | 法律面前人人平等表现在哪些方面| 丰满人妻熟妇乱又伦精品不卡| 色综合婷婷激情| 成人国产综合亚洲| 欧美又色又爽又黄视频| 18禁国产床啪视频网站| 淫秽高清视频在线观看| 亚洲七黄色美女视频| 久久久久性生活片| 国产精品一区二区免费欧美| 小蜜桃在线观看免费完整版高清| 中文亚洲av片在线观看爽| 国产三级中文精品| 97人妻精品一区二区三区麻豆| 色噜噜av男人的天堂激情| 哪里可以看免费的av片| 88av欧美| 琪琪午夜伦伦电影理论片6080| 十八禁人妻一区二区| 午夜福利视频1000在线观看| 无人区码免费观看不卡| 成熟少妇高潮喷水视频| 真实男女啪啪啪动态图| 国产精品99久久久久久久久| 国产精品亚洲一级av第二区| 成人国产一区最新在线观看| 国产精品久久久久久精品电影| 国产亚洲欧美在线一区二区| 制服丝袜大香蕉在线| 一本久久中文字幕| 免费在线观看亚洲国产| 国产 一区 欧美 日韩| 亚洲国产欧美网| 国内精品一区二区在线观看| 999久久久精品免费观看国产| 久久这里只有精品中国| 男女床上黄色一级片免费看| 叶爱在线成人免费视频播放| 欧美成狂野欧美在线观看| 欧美色欧美亚洲另类二区| svipshipincom国产片| 又黄又粗又硬又大视频| 亚洲午夜精品一区,二区,三区| 12—13女人毛片做爰片一| 午夜a级毛片| 亚洲第一欧美日韩一区二区三区| 国产人伦9x9x在线观看| 久久久久久久久中文| 午夜视频精品福利| 精品久久久久久,| 视频区欧美日本亚洲| 禁无遮挡网站| 精品一区二区三区视频在线 | 亚洲精华国产精华精| 少妇的丰满在线观看| 国产av麻豆久久久久久久| 亚洲自偷自拍图片 自拍| 此物有八面人人有两片| 亚洲国产欧洲综合997久久,| 中亚洲国语对白在线视频| 在线免费观看不下载黄p国产 | 国产成人aa在线观看| 亚洲第一电影网av| 天堂av国产一区二区熟女人妻| av福利片在线观看| 一进一出抽搐gif免费好疼| 啦啦啦免费观看视频1| 色吧在线观看| 国产高清有码在线观看视频| 欧美色欧美亚洲另类二区| 亚洲av成人精品一区久久| 99国产精品一区二区三区| 999久久久国产精品视频| 久久这里只有精品中国| 最近视频中文字幕2019在线8| 亚洲五月婷婷丁香| www国产在线视频色| 成人18禁在线播放| 好看av亚洲va欧美ⅴa在| 亚洲成人免费电影在线观看| 午夜免费成人在线视频| 老汉色av国产亚洲站长工具| a级毛片在线看网站| 亚洲人与动物交配视频| 校园春色视频在线观看| 又黄又粗又硬又大视频| 欧美三级亚洲精品| 男女下面进入的视频免费午夜| 听说在线观看完整版免费高清| 禁无遮挡网站| 婷婷丁香在线五月| 在线国产一区二区在线| 亚洲成人免费电影在线观看| 日韩欧美 国产精品| 国产精品久久久人人做人人爽| 日本免费a在线| 99久久国产精品久久久| 亚洲中文日韩欧美视频| 999久久久国产精品视频| 国产91精品成人一区二区三区| 亚洲欧美日韩卡通动漫| 非洲黑人性xxxx精品又粗又长| 亚洲欧洲精品一区二区精品久久久| 国产免费男女视频| 午夜视频精品福利| 国产伦精品一区二区三区视频9 | 欧美日本视频| 日本黄色视频三级网站网址| 国产又色又爽无遮挡免费看| 舔av片在线| 美女cb高潮喷水在线观看 | 国产乱人视频| 视频区欧美日本亚洲| 亚洲五月天丁香| 亚洲avbb在线观看| 亚洲专区字幕在线| 亚洲 欧美 日韩 在线 免费| 欧美在线黄色| 免费观看人在逋| 欧美极品一区二区三区四区| 在线十欧美十亚洲十日本专区| 1024香蕉在线观看| 黄片大片在线免费观看| 国产精品野战在线观看| 国产精品久久久久久亚洲av鲁大| 一个人免费在线观看的高清视频| 色视频www国产| 成年版毛片免费区| 亚洲美女视频黄频| 舔av片在线| 色综合欧美亚洲国产小说| 国产黄片美女视频| 久久这里只有精品中国| 18美女黄网站色大片免费观看| 中文在线观看免费www的网站| 国产精品亚洲一级av第二区| 天堂√8在线中文| 三级毛片av免费| 欧美日韩国产亚洲二区| 成人三级黄色视频| 母亲3免费完整高清在线观看| 午夜精品一区二区三区免费看| 国产爱豆传媒在线观看| 亚洲中文av在线| 非洲黑人性xxxx精品又粗又长| 国产精品99久久久久久久久| 久久亚洲精品不卡| 一级作爱视频免费观看| 国产亚洲精品综合一区在线观看| 久久精品人妻少妇| 天堂动漫精品| 亚洲欧美日韩高清专用| 波多野结衣巨乳人妻| 亚洲精品在线观看二区| 少妇裸体淫交视频免费看高清| av女优亚洲男人天堂 | 19禁男女啪啪无遮挡网站| 免费在线观看亚洲国产| 亚洲专区国产一区二区| 国产精品精品国产色婷婷| 久久草成人影院| 免费看光身美女| 亚洲片人在线观看| 亚洲国产精品成人综合色| 久久这里只有精品19| 国产伦在线观看视频一区| 亚洲自拍偷在线| 免费人成视频x8x8入口观看| 亚洲国产中文字幕在线视频| av福利片在线观看| 99热只有精品国产| 99精品久久久久人妻精品| 成人18禁在线播放| 成人特级黄色片久久久久久久| 在线永久观看黄色视频| 毛片女人毛片| 国产欧美日韩精品一区二区| 好男人电影高清在线观看| 日日干狠狠操夜夜爽| 亚洲av成人精品一区久久| 黑人巨大精品欧美一区二区mp4| 97超视频在线观看视频| 可以在线观看的亚洲视频| 一个人看视频在线观看www免费 | 亚洲精品在线美女| 99热只有精品国产| 久久性视频一级片| 国产伦精品一区二区三区四那| 2021天堂中文幕一二区在线观| 国产黄色小视频在线观看| 法律面前人人平等表现在哪些方面| 在线观看舔阴道视频| av欧美777| 久久久久九九精品影院| 精品国内亚洲2022精品成人| 欧美一区二区精品小视频在线| 久久久水蜜桃国产精品网| 午夜激情福利司机影院| 欧美绝顶高潮抽搐喷水| 国产高清激情床上av| 最近最新免费中文字幕在线| 国产精品一区二区免费欧美| 黄色成人免费大全| 在线永久观看黄色视频| 老鸭窝网址在线观看| 国产精品久久电影中文字幕| 九九在线视频观看精品| 久久久成人免费电影| 欧美黑人巨大hd| 高潮久久久久久久久久久不卡| 亚洲精品一区av在线观看| 最近在线观看免费完整版| 中出人妻视频一区二区| 久久精品91无色码中文字幕| 婷婷亚洲欧美| 岛国在线免费视频观看| 午夜日韩欧美国产| 宅男免费午夜| 欧美3d第一页| av中文乱码字幕在线| 亚洲人成电影免费在线| 精品福利观看| cao死你这个sao货| 欧美乱色亚洲激情| 90打野战视频偷拍视频| www.熟女人妻精品国产| 村上凉子中文字幕在线| 床上黄色一级片| 亚洲九九香蕉| av黄色大香蕉| 久久久久久人人人人人| 日本免费a在线| 免费在线观看成人毛片| 精品久久久久久久人妻蜜臀av| 日韩欧美 国产精品| 欧美午夜高清在线| xxxwww97欧美| 午夜两性在线视频| 亚洲中文av在线| 久久伊人香网站| 欧美日韩综合久久久久久 | 国产在线精品亚洲第一网站| 999久久久国产精品视频| 欧美黑人巨大hd| 国产精品永久免费网站| 91字幕亚洲| 色播亚洲综合网| 9191精品国产免费久久| 午夜免费观看网址| 久久久色成人| 网址你懂的国产日韩在线| 十八禁网站免费在线| 亚洲精品456在线播放app | 黑人巨大精品欧美一区二区mp4| 长腿黑丝高跟| 国产视频内射| 熟女电影av网| 免费看美女性在线毛片视频| 校园春色视频在线观看| 久久久久久九九精品二区国产| 国内毛片毛片毛片毛片毛片| 色尼玛亚洲综合影院| 中文亚洲av片在线观看爽| 黄色 视频免费看| 日本一本二区三区精品| 男人舔女人的私密视频| 国产欧美日韩一区二区精品| 老司机在亚洲福利影院| 九九热线精品视视频播放| 男女那种视频在线观看| 草草在线视频免费看| 久久这里只有精品19| 国内少妇人妻偷人精品xxx网站 | 99久久精品一区二区三区| 97超视频在线观看视频| 中出人妻视频一区二区| 人妻夜夜爽99麻豆av| 岛国在线免费视频观看| 99热精品在线国产| 免费电影在线观看免费观看| 国产成人欧美在线观看| 在线免费观看的www视频| 国产麻豆成人av免费视频| 亚洲专区国产一区二区| 18禁美女被吸乳视频| 欧美日本视频| 国产精品日韩av在线免费观看| 男女那种视频在线观看| 国产精品一及| 国产精品 欧美亚洲| 18禁黄网站禁片免费观看直播| 国产成人一区二区三区免费视频网站| 国产乱人伦免费视频| a级毛片a级免费在线| 精品国产乱子伦一区二区三区| 午夜福利欧美成人| www国产在线视频色| 大型黄色视频在线免费观看| 午夜福利在线观看吧| 国产午夜精品论理片| 亚洲精品美女久久久久99蜜臀| 精品日产1卡2卡| 亚洲真实伦在线观看| 成人无遮挡网站| 老司机深夜福利视频在线观看| 综合色av麻豆| 怎么达到女性高潮| 国产精品一区二区三区四区久久| 国内精品久久久久久久电影| 熟女电影av网| 亚洲av熟女| 亚洲国产欧美网| 88av欧美| 国产人伦9x9x在线观看| 男女下面进入的视频免费午夜| 欧美丝袜亚洲另类 | 亚洲精品一卡2卡三卡4卡5卡| 国产伦人伦偷精品视频| 看黄色毛片网站| 热99在线观看视频| 很黄的视频免费| 亚洲 欧美 日韩 在线 免费| 欧美黑人巨大hd| 搞女人的毛片| 精品国产乱码久久久久久男人| 国产精品久久电影中文字幕| 男女午夜视频在线观看| 在线免费观看不下载黄p国产 | 91九色精品人成在线观看| 亚洲激情在线av| 校园春色视频在线观看| 欧美日本亚洲视频在线播放| 亚洲人成伊人成综合网2020| 精品国产超薄肉色丝袜足j| 欧美在线一区亚洲| 亚洲最大成人中文| 国产一区二区在线av高清观看| 久久久久精品国产欧美久久久| 久久人人精品亚洲av| 午夜激情福利司机影院| 麻豆久久精品国产亚洲av| 国模一区二区三区四区视频 | 老汉色av国产亚洲站长工具| 精品一区二区三区av网在线观看| 黄色女人牲交| 久久草成人影院| 99国产极品粉嫩在线观看| 国产真实乱freesex| 黄色成人免费大全| 欧美激情在线99| 岛国在线观看网站| 亚洲无线在线观看| 中文字幕久久专区| 深夜精品福利| 精品一区二区三区av网在线观看| 97超视频在线观看视频| 真实男女啪啪啪动态图| 五月伊人婷婷丁香| 99久久国产精品久久久| 我的老师免费观看完整版| 香蕉国产在线看| 悠悠久久av| 在线观看午夜福利视频| 网址你懂的国产日韩在线| 亚洲黑人精品在线| 超碰成人久久| 18禁美女被吸乳视频| 一级毛片高清免费大全| 19禁男女啪啪无遮挡网站| 高潮久久久久久久久久久不卡| 九色成人免费人妻av| 99精品在免费线老司机午夜| 波多野结衣高清作品| 亚洲国产中文字幕在线视频| 亚洲国产精品sss在线观看| 欧美3d第一页| 色综合亚洲欧美另类图片| 国产又黄又爽又无遮挡在线| 国产91精品成人一区二区三区| 欧美大码av| 18禁国产床啪视频网站| 欧美成人免费av一区二区三区| 美女午夜性视频免费| 欧美一级毛片孕妇| 亚洲自拍偷在线| 欧美3d第一页| 欧美av亚洲av综合av国产av| 在线看三级毛片| 日本五十路高清| 首页视频小说图片口味搜索| 村上凉子中文字幕在线| 性欧美人与动物交配| 亚洲国产精品999在线| 欧美中文综合在线视频| 男女下面进入的视频免费午夜| 色av中文字幕| 国产不卡一卡二| 欧美中文综合在线视频| 美女大奶头视频| 久久久久国产一级毛片高清牌| www日本黄色视频网| 日本在线视频免费播放| 少妇丰满av| 国产久久久一区二区三区| 丰满的人妻完整版| xxx96com| 变态另类丝袜制服| 精品国产乱子伦一区二区三区| 日本精品一区二区三区蜜桃| 怎么达到女性高潮| 久久香蕉国产精品| 精品久久久久久久人妻蜜臀av| 久久久精品欧美日韩精品| 超碰成人久久| 国产精品 国内视频| 一区二区三区激情视频| 给我免费播放毛片高清在线观看| 999精品在线视频| 亚洲成人久久爱视频| 欧美又色又爽又黄视频| 桃色一区二区三区在线观看| av国产免费在线观看| 人妻久久中文字幕网| 精品国产三级普通话版| av黄色大香蕉| 亚洲av成人av| 无限看片的www在线观看|