CN107703493A - Sea clutter optimal soft survey instrument and method based on adaptive drosophila optimized algorithm Optimized Least Square Support Vector - Google Patents

Sea clutter optimal soft survey instrument and method based on adaptive drosophila optimized algorithm Optimized Least Square Support Vector Download PDF

Info

Publication number
CN107703493A
CN107703493A CN201711117137.5A CN201711117137A CN107703493A CN 107703493 A CN107703493 A CN 107703493A CN 201711117137 A CN201711117137 A CN 201711117137A CN 107703493 A CN107703493 A CN 107703493A
Authority
CN
China
Prior art keywords
support vector
vector machine
optimization
function
axism
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201711117137.5A
Other languages
Chinese (zh)
Inventor
刘兴高
王文川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201711117137.5A priority Critical patent/CN107703493A/en
Publication of CN107703493A publication Critical patent/CN107703493A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications

Abstract

A most young waiter in a wineshop or an inn is optimized into the sea clutter optimal soft survey instrument and method of SVMs based on TSP question drosophila the invention discloses a kind of, including radar, field intelligent instrument, control station, the spot database for depositing data, the optimal hard measurement host computer based on improved drosophila optimized algorithm Optimization of Wavelet neutral net and forecast hard measurement value display instrument, field intelligent instrument and control station are connected with radar, are connected with spot database;Optimal hard measurement host computer is connected with spot database and hard measurement value display instrument.Described optimizes a most young waiter in a wineshop or an inn into the optimal hard measurement host computer of SVMs, including data preprocessing module, wavelet neural network module, model modification module based on TSP question drosophila.The present invention realizes the online optimal hard measurement of sea clutter, and randomness caused by overcoming human factor influences, and improves the stability of model prediction, reduces the possibility that model prediction is absorbed in local optimum.

Description

Sea clutter optimal soft measuring instrument and method for optimizing least square support vector machine based on self-adaptive fruit fly optimization algorithm
Technical Field
The invention relates to the field of optimal soft measuring instruments and methods, in particular to a sea clutter optimal soft measuring instrument and method based on a self-adaptive drosophila optimization algorithm optimization least square support vector machine.
Background
In the field of radar, echo signals reflected from the surface of seawater are called sea clutter, and the sea clutter is related to various factors such as sea conditions, tides, radar parameters and the like. For coastal warning radars, ship-borne radars and other radars working in marine environments, the detection and tracking performance of sea surface targets are influenced by serious sea surface reflection echoes, the properties of sea clutter are mastered, and the establishment of an accurate sea clutter model is a precondition for analyzing and improving the radar performance. The statistical properties of the sea clutter include amplitude properties and correlation properties. The correlation properties of the sea clutter include temporal correlation and spatial correlation. The time correlation is also called as inter-pulse correlation, reflects the fluctuation of the amplitude of the sea clutter with time, and can be equivalently represented by a power spectrum. The spatial correlation of sea clutter is divided into azimuth correlation and distance correlation. The wave energy of the sea helps us to study the dynamic characteristics of the sea, but if targets are to be detected from the background of the sea clutter, such as floating ice, ships, etc., it becomes a great obstacle and must be suppressed as much as possible to reduce or eliminate these disturbances. The main purpose of the research on sea clutter is as follows: on one hand, the natural mechanism of the sea clutter is explained, and a reasonable model is provided; on the other hand, the method is to reduce the interference of the sea clutter on the detection target and find out how to extract the target signal submerged in the Jiang Hai clutter background. The establishment of an accurate sea clutter model is the key to achieving the above objectives
Most of the research work on modeling and forecasting of the sea clutter in recent years is focused on the artificial neural network, and good effects are achieved. However, artificial neural networks also have their own drawbacks, such as overfitting, the number of nodes in the hidden layer, and poor parameter determination. Secondly, the data collected in the observation field has certain uncertain errors due to noise, manual operation errors and the like, so that the forecasting model using the artificial neural network with strong certainty generally has weak popularization capability. Support vector machines, introduced by Vapnik in 1998, are widely used in pattern recognition, fitting and classification problems due to their good generalization ability. Since the standard support vector machine is sensitive to isolated points and noise points, a least squares support vector machine was proposed later. The least squares support vector machine is better able to process noisy sample data than the standard support vector machine, and is used here for modeling. The Fruit Fly Optimization Algorithm, namely, the Fruit Fly Optimization Algorithm, is a biological intelligent Optimization Algorithm, called FOA for short, which is proposed by professor Wen-Tsao Pan in Taiwan and is deduced based on foraging behavior of Fruit flies. In the whole optimizing process of the traditional fruit fly optimizing algorithm, once the optimal individual of the iteration is found, the fruit fly population can be gathered to the individual. If the found optimal fruit fly individual is not the global optimal point, the optimal fruit fly individual is likely to fall into local optimal, the convergence speed and the convergence precision are reduced, and the premature convergence problem is caused. In order to overcome the defects, a self-adaptive variant drosophila optimization algorithm is provided, and the algorithm can ensure that local optima are jumped out, and other spaces are searched continuously until a global optima point is found. The adaptive variant drosophila optimization algorithm is used for searching the optimal parameter combination of the least square support vector machine so as to achieve the purpose of optimizing the model.
Disclosure of Invention
In order to overcome the defects of low measurement precision, low noise sensitivity and poor popularization performance of the conventional radar, the invention provides the optimal sea clutter soft measurement instrument and method based on the adaptive drosophila optimization algorithm optimization least square support vector machine, which have the advantages of online measurement, high calculation speed, automatic model updating, strong noise resistance and good popularization performance.
The purpose of the invention is realized by the following technical scheme: a sea clutter optimal soft measuring instrument based on a self-adaptive drosophila optimization algorithm optimization least square support vector machine comprises a radar, a field intelligent instrument for measuring easily-measured variables, a control station for measuring operation variables, a field database for storing data and a sea clutter forecast soft measurement value display instrument; the field intelligent instrument and the control station are connected with the radar, the field intelligent instrument and the control station are connected with a field database, the soft measuring instrument further comprises an optimal soft measuring upper computer based on an adaptive fruit fly optimization algorithm optimization least square support vector machine, the field database is connected with the input end of the optimal soft measuring upper computer based on the adaptive fruit fly optimization algorithm optimization least square support vector machine, and the output end of the optimal soft measuring upper computer based on the adaptive fruit fly optimization algorithm optimization least square support vector machine is connected with a sea clutter soft measurement value display instrument; the optimal soft measurement upper computer based on the least square support vector machine optimized by the adaptive fruit fly optimization algorithm comprises:
the data preprocessing module is used for preprocessing the model training samples input from the field database, centralizing the training samples, namely subtracting the average value of the samples, and then normalizing the training samples:
calculating an average value:
calculating the variance:
and (3) standardization:
wherein TX is a training sample, N is the number of training samples,is the mean of the training samples, X is the normalized training sample, σ x To calculate the variance.
And the least square support vector machine module adopts a least square support vector machine for modeling. Model ith normalized training sample X i Target output of (2) is O i The least squares support vector machine equates the fitting problem to the following quadratic programming problem by transformation:
where R (w, ξ) is the objective function of the optimization problem,a non-linear mapping function, N is the number of training samples, ξ = { ξ = 1 ,…,ξ N Is the relaxation variable, w is the normal vector of the hyperplane of the least squares support vector machine, b is the corresponding offset, and ω is i I =1, …, N and γ are the weight and penalty factors, respectively, of the least squares support vector machine, whereIs the relaxation variable xi i Estimate of standard deviation, constant c 1 ,c 2 Is usually taken as c 1 =2.5,c 2 =3, from which a normalized training sample X can be obtained i The output of (c) is:
wherein, K<·&gt is the kernel function of a least squares support vector machine, where K<·&Taking a linear kernel function; alpha is alpha m M =1, …, N is the corresponding lagrange multiplier.
The adaptive fruit fly optimization algorithm module is used for optimizing punishment factors and error tolerance values of the least square support vector machine by adopting an adaptive fruit fly optimization algorithm, and comprises the following specific steps of:
(1) determining optimization parameters of the self-adaptive drosophila optimization algorithm as a penalty factor and an error tolerance value of a least square support vector machine module, the individual number of particle swarm popsize and the maximum cyclic optimization number iter max Variance threshold delta of group fitness 1 Theoretical bestFigure of merit delta 2 And the initial position regions X _ axis, Y _ axis of the p-th particle.
(2) Setting an optimization objective function, converting the optimization objective function into fitness, calculating the fitness function through a corresponding error function, considering that the fitness of the particles with large errors is small, and expressing the fitness function f of the particles p as:
f p =1/(E p +1) (8)
in the formula, E p Is the error function of the least squares support vector machine model, expressed as:
in the formula (I), the compound is shown in the specification,is the predicted output of the least squares support vector machine model, O i Outputting a target of a least squares support vector machine model; n is the number of training samples;
(3) according to the following formula, the particles are searched,
in the formula, randomValue is the search distance;
(4) for the particle p, the distance Dist from the origin is estimated in advance, and the taste concentration determination value S is calculated as the reciprocal distance:
S i =1/Dist i (12)
(5) judging the taste concentration value S i A substitute taste concentration judgment function (or fitness function) for determining the taste concentration Smell of the individual positions of the drosophila i :
(6) The mean taste of the Drosophila population Smell was calculated according to equation (14) ave Then, the fruit fly population fitness variance tau is calculated according to the formula (15) 2
(7) If τ is 2 ≤δ 1 And Smellbest>δ 2 Or distributed in [0,1]Random number r between&P, then, firstly, the optimal drosophila individuals (X _ axism) are replicated according to the formula (10) j ,Y_axism j ) (j =1,2, …, M); secondly, the copied optimal fruit fly individual is mutated according to the formula (11), and the position of the copied optimal fruit fly individual is updated to be a new position (X _ axism) j ,Y_axism j )(j=1,2,…, M):
(8) The new position (X _ axism) is estimated again first according to the following equation j ,Y_axism j ) And the distance Dist 'from the original point, and then calculating a new position taste concentration judgment value S' according to a formula:
(9) the taste concentration determination value S 'calculated again is substituted into the taste concentration determination function to calculate the taste concentration Smell' at the new position.
Smell i '=Function(S i ') (19)
If the R satisfies Smell i '&l < Smellbest > then Smellbest = Smell' j ,X_axis=X_axism' j , Y_axis=Y_axism' j ,(j=1,2,…,M).
Judging whether the performance requirements are met, if so, finishing the optimization to obtain a group of optimized parameters of the least square support vector machine; otherwise, returning to the step (5), and continuing the iteration optimization until the maximum iteration number iter is reached max
As a preferred scheme, the optimal soft measurement upper computer for optimizing the least square support vector machine based on the adaptive drosophila optimization algorithm further comprises: and the model updating module is used for updating the model on line, inputting offline verification data into a training set regularly and updating the least square support vector machine model.
A sea clutter optimal soft measurement method based on an adaptive drosophila optimization algorithm optimization least square support vector machine comprises the following steps:
1) Selecting an operation variable and an easily-measured variable as the input of a model for the radar object according to characteristic analysis and climate analysis, wherein the operation variable and the easily-measured variable are obtained by a field database;
2) Preprocessing a model training sample input from a field database, centralizing the training sample, namely subtracting the average value of the sample, and then normalizing the training sample so that the average value is 0 and the variance is 1. The processing is accomplished using the following mathematical process:
2.1 Calculate the mean value:
2.2 Calculate variance:
2.3 Normalization:
wherein TX is a training sample, N is the number of training samples,is the mean of the training samples, X is the normalized training sample, σ x To calculate the variance.
3) And modeling the training sample transmitted from the data preprocessing module by adopting a least square support vector machine. Training sample X after model standardization i Target output of O i The least squares support vector machine equates the fitting problem to the following quadratic programming problem by transformation:
where R (w, ξ) is the objective function of the optimization problem,a non-linear mapping function, N is the number of training samples, ξ = { ξ = 1 ,…,ξ N Is the relaxation variable, w is the normal vector of the hyperplane of the least squares support vector machine, b is the corresponding offset, and ω is i I =1, …, N and γ are the weight and penalty factors, respectively, of the least squares support vector machine, whereIs the relaxation variable xi i Estimate of standard deviation, constant c 1 ,c 2 Is usually taken as c 1 =2.5,c 2 =3, from which a normalized training sample X can be obtained i The output of (c) is:
wherein, K<·&gt is the kernel function of a least squares support vector machine, where K<·&Taking a linear kernel function; alpha is alpha m M =1, …, N is the corresponding lagrange multiplier.
4) The method adopts a self-adaptive drosophila optimization algorithm to optimize the punishment factor and the error tolerance value of the least square support vector machine, and comprises the following specific steps:
(1) determining optimization parameters of the self-adaptive drosophila optimization algorithm as a penalty factor and an error tolerance value of a least square support vector machine module, the individual number of particle swarm popsize and the maximum cyclic optimization number iter max Variance threshold delta of group fitness 1 Theoretical optimum value delta 2 And the initial position regions X _ axis, Y _ axis of the p-th particle.
(2) Setting an optimization objective function, converting the optimization objective function into fitness, calculating the fitness function through a corresponding error function, considering that the fitness of the particles with large errors is small, and expressing the fitness function f of the particles p as:
f p =1/(E p +1) (8)
in the formula, E p Is the error function of the least squares support vector machine model, expressed as:
in the formula (I), the compound is shown in the specification,is the predicted output of the least squares support vector machine model, O i For least-squares support vector machine modelsOutputting a target; n is the number of training samples;
(3) according to the following formula, the particles are searched,
in the formula, randomValue is the search distance;
(4) for the particle p, the distance Dist from the origin is estimated in advance, and the taste concentration determination value S is calculated as the reciprocal distance:
Dist i =(X i 2 +Y i 2 ) 1/2 (11)
S i =1/Dist i (12)
(5) judging the taste concentration value S i A substitute taste concentration judgment function (or fitness function) for determining the taste concentration Smell of the individual positions of the drosophila i :
(6) The mean taste of the Drosophila population Smell was calculated according to equation (14) ave Then calculating the fruit fly population fitness variance tau according to the formula (15) 2
(7) If τ is 2 ≤δ 1 And Smellbest>δ 2 Or distributed in [0,1]Random number r between&P, then, firstly, the optimal fruit fly individual (X _ axism) is duplicated according to the formula (10) j ,Y_axism j ) (j =1,2, …, M); according to the formula (11 Carrying out mutation on the copied optimal fruit fly individual, and updating the position of the copied optimal fruit fly individual to a new position (X _ axism) j ,Y_axism j )(j=1,2,…, M):
(8) The new position (X _ axism) is estimated again first according to the following equation j ,Y_axism j ) And the distance Dist 'from the original point, and then calculating a new position taste concentration judgment value S' according to a formula:
(9) the taste concentration determination value S 'calculated again is substituted into the taste concentration determination function to calculate the taste concentration Smell' at the new position. Smell i '=Function(S i ')(19)
If the R satisfies Smell i '&l < Smellbest > then Smellbest = Smell' j ,X_axis=X_axism' j , Y_axis=Y_axism' j ,(j=1,2,…,M).
Judging whether the performance requirements are met, if so, ending the optimization, and obtaining a group of optimized parameters of the least square support vector machine;
otherwise, returning to the step (5), and continuing the iteration optimization until the maximum iteration number iter is reached max
As a preferred solution: the soft measurement method further comprises the following steps: 5) And inputting the offline experimental data into a training set regularly, and updating the least square support vector machine model.
The technical conception of the invention is as follows: the method has the advantages that online optimal soft measurement is carried out on the sea clutter, the defects that an existing sea clutter measuring instrument is low in measuring precision, low in noise sensitivity and poor in popularization performance are overcome, the adaptive drosophila optimization algorithm is introduced to carry out automatic optimization on a least square support vector machine model, and parameters of the least square support vector machine are adjusted for multiple times without human experience, so that an optimal soft measurement result is obtained. Compared with the existing sea clutter soft measurement model, the model has the following advantages: modeling is carried out through a least square support vector machine model, and high forecasting precision is achieved; the existing model parameters are generally determined through experience of operators, the existing model parameters have certain limitation and uncertainty, once the values are not properly taken, oscillation of model prediction output and larger prediction errors can be caused, and the model automatically optimizes the parameters of the model through a self-adaptive variation drosophila optimization algorithm to obtain an optimal soft measurement model.
The invention has the following beneficial effects: the online soft measurement of the sea clutter is realized, and the online soft measurement has the advantages of strong noise resistance, high interference resistance, high precision and strong popularization capability.
Drawings
FIG. 1 is a schematic diagram of the basic structure of an optimal soft measurement instrument and method for a sea clutter modeling process based on an adaptive drosophila optimization algorithm optimization least square support vector machine;
FIG. 2 is a schematic diagram of an optimal soft measurement upper computer structure based on an adaptive drosophila optimization algorithm optimization least square support vector machine.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The examples are intended to illustrate the invention, but not to limit the invention, and any modifications and variations of the invention within the spirit and scope of the claims are intended to fall within the scope of the invention.
Example 1
Referring to fig. 1 and 2, the sea clutter optimal soft measuring instrument based on the adaptive drosophila optimization algorithm optimization least square support vector machine comprises a radar 1, a field intelligent instrument 2 used for measuring easily-measured variables, a control station 3 used for measuring operation variables, a field database 4 used for storing data and a sea clutter soft measurement value display instrument 6, wherein the field intelligent instrument 2 and the control station 3 are connected with the radar 1, the field intelligent instrument 2 and the control station 3 are connected with the field database 4, the soft measuring instrument further comprises an optimal soft measurement upper computer 5 of the adaptive drosophila optimization algorithm optimization least square support vector machine, the field database 4 is connected with the input end of the optimal soft measurement upper computer 5 based on the adaptive drosophila optimization algorithm optimization support vector machine, and the output end of the optimal soft measurement upper computer 5 of the adaptive drosophila optimization least square support vector machine is connected with the sea clutter soft measurement value display instrument 6. The field intelligent instrument 2 measures the easily measurable variable of the radar object and transmits the easily measurable variable to the field database 4; the control station 3 controls the manipulated variables of the radar object, and transmits the manipulated variables to the field database 4. The variable data recorded in the field database 4 is used as the input of the optimal soft measurement upper computer 5 based on the adaptive fruit fly optimization algorithm optimization least square support vector machine, and the sea clutter soft measurement value display instrument 6 is used for displaying the output, namely the soft measurement value, of the optimal soft measurement upper computer 5 based on the adaptive fruit fly optimization algorithm optimization least square support vector machine. The optimal soft measurement upper computer 5 based on the least square support vector machine optimized by the adaptive fruit fly optimization algorithm comprises the following 4 parts:
a data preprocessing module 7, configured to preprocess model training samples input from the field database, centralize the training samples, that is, subtract an average value of the samples, and then normalize the training samples:
calculating an average value:
calculating the variance:
and (3) standardization:
wherein TX is a training sample, N is the number of training samples,is the mean of the training samples, X is the normalized training sample, σ x To calculate the variance.
And the least square support vector machine module 8 adopts a least square support vector machine to carry out modeling. Model ith normalized training sample X i Target output of O i The least squares support vector machine equates the fitting problem to the following quadratic programming problem by transformation:
where R (w, ξ) is the objective function of the optimization problem,a non-linear mapping function, N is the number of training samples, xi = { xi = 1 ,…,ξ N Is the relaxation variable, w is the normal vector of the hyperplane of the least squares support vector machine, b is the corresponding offset, and ω is i I =1, …, N and γ are the weight and penalty factors, respectively, of the least squares support vector machine, whereIs the relaxation variable xi i The estimation of the standard deviation is carried out,constant c 1 ,c 2 Is usually taken as c 1 =2.5,c 2 =3, from which a normalized training sample X can be obtained i The output of (c) is:
wherein, K<·&gt is the kernel function of a least squares support vector machine, where K<·&Taking a linear kernel function; alpha is alpha m M =1, …, N is the corresponding lagrange multiplier.
The adaptive fruit fly optimization algorithm module 9 is used for optimizing the punishment factor and the error tolerance value of the least square support vector machine by adopting an adaptive fruit fly optimization algorithm, and comprises the following specific steps:
(1) determining optimization parameters of the self-adaptive drosophila optimization algorithm as a penalty factor and an error tolerance value of a least square support vector machine module, the individual number of particle swarm popsize and the maximum cyclic optimization number iter max Variance threshold delta of group fitness 1 Theoretical optimum value delta 2 And the initial position regions X _ axis, Y _ axis of the p-th particle.
(2) Setting an optimization objective function, converting the optimization objective function into fitness, calculating the fitness function through a corresponding error function, considering that the fitness of the particles with large errors is small, and expressing the fitness function f of the particles p as:
f p =1/(E p +1) (8)
in the formula, E p Is an error function of a least squares support vector machine model, expressed as:
in the formula (I), the compound is shown in the specification,is the predicted output of the least squares support vector machine model, O i For least squares support vector machine modelA target output of the model; n is the number of training samples;
(3) according to the following formula, the particles are searched,
in the formula, randomValue is the search distance;
(4) for the particle p, the distance Dist from the origin is estimated in advance, and the taste concentration determination value S is calculated as the reciprocal distance:
Dist i =(X i 2 +Y i 2 ) 1/2 (11)
S i =1/Dist i (12)
(5) judging the taste concentration value S i A substitute taste concentration judgment function (or fitness function) for calculating the taste concentration Smell of the individual positions of the fruit flies i :
(6) The mean taste of the Drosophila population Smell was calculated according to equation (14) ave Then, the fruit fly population fitness variance tau is calculated according to the formula (15) 2
(7) If τ is 2 ≤δ 1 And Smellbest>δ 2 Or distributed in [0,1]Random number r between&P, then, firstly, the optimal fruit fly individual (X _ axism) is duplicated according to the formula (10) j ,Y_axism j ) (j =1,2, …, M); followed byThe copy optimal fruit fly individual is mutated in the formula (11), and the position of the copy optimal fruit fly individual is updated to a new position (X _ axism) j ,Y_axism j )(j=1,2,…, M):
(8) The new position (X _ axism) is estimated again first according to the following equation j ,Y_axism j ) And the distance Dist 'from the original point, and then calculating a new position taste concentration judgment value S' according to a formula:
(9) the taste concentration determination value S 'calculated again is substituted into the taste concentration determination function to calculate the taste concentration Smell' at the new position.
Smell i '=Function(S i ') (19)
If the R satisfies Smell i '&l < Smellbest > then Smellbest = Smell' j ,X_axis=X_axism' j , Y_axis=Y_axism' j ,(j=1,2,…,M).
Judging whether the performance requirements are met, if so, finishing the optimization to obtain a group of optimized parameters of the least square support vector machine;
otherwise, returning to the step (5), and continuing the iteration optimization until the maximum iteration number iter is reached max
As a preferred scheme, the optimal soft measurement upper computer for optimizing the least square support vector machine based on the adaptive drosophila optimization algorithm further comprises: and the model updating module 10 is used for updating the model on line, inputting offline experimental data into a training set regularly, and updating the least square support vector machine model.
Example 2
Referring to fig. 1 and 2, a sea clutter optimal soft measurement method based on an adaptive drosophila optimization algorithm optimization least square support vector machine comprises the following steps:
1) Selecting an operation variable and an easily-measured variable as input of a model for the radar object according to process analysis and operation analysis, wherein the operation variable and the easily-measured variable are obtained by a field database;
2) Preprocessing a model training sample input from a field database, centralizing the training sample, namely subtracting the average value of the sample, and then normalizing the training sample so that the average value is 0 and the variance is 1. The processing is accomplished using the following mathematical process:
2.1 Calculate the mean value:
2.2 Calculate variance:
2.3 Normalization:
wherein TX is a training sample, N is the number of training samples,is the mean of the training samples, X is the normalized training sample, σ x To calculate the variance.
3) And modeling the training sample transmitted from the data preprocessing module by adopting a least square support vector machine. Training sample X after model standardization i Target output of O i The least squares support vector machine equates the fitting problem to the following quadratic programming problem by transformation:
where R (w, ξ) is the objective function of the optimization problem,a non-linear mapping function, N is the number of training samples, ξ = { ξ = 1 ,…,ξ N Is the relaxation variable, w is the normal vector of the hyperplane of the least squares support vector machine, b is the corresponding offset, and ω is i I =1, …, N and γ are the weight and penalty factors, respectively, of the least squares support vector machine, whereIs the relaxation variable xi i Estimate of standard deviation, constant c 1 ,c 2 Is usually taken as c 1 =2.5,c 2 =3, from which normalized training sample X can be obtained i The output of (c) is:
wherein, K<·&gt is the kernel function of a least squares support vector machine, where K<·&Taking a linear kernel function; alpha (alpha) ("alpha") m M =1, …, N is the corresponding lagrange multiplier.
4) The method adopts a self-adaptive drosophila optimization algorithm to optimize the punishment factor and the error tolerance value of the least square support vector machine, and comprises the following specific steps:
(1) determining an adaptive fruit fly optimization algorithmThe optimization parameters are a penalty factor and an error tolerance value of a least square support vector machine module, the individual number of particle swarms popsize and the maximum cyclic optimization number iter max Variance threshold delta of group fitness 1 Theoretical optimum value delta 2 And the initial position regions X _ axis, Y _ axis of the p-th particle.
(2) Setting an optimization objective function, converting the optimization objective function into fitness, calculating the fitness function through a corresponding error function, considering that the fitness of the particles with large errors is small, and expressing the fitness function f of the particles p as:
f p =1/(E p +1) (8)
in the formula, E p Is the error function of the least squares support vector machine model, expressed as:
in the formula (I), the compound is shown in the specification,is the predicted output of the least squares support vector machine model, O i Outputting a target of a least squares support vector machine model; n is the number of training samples;
(3) according to the following formula, the particles are searched,
in the formula, randomValue is the search distance;
(4) for the particle p, the distance Dist from the origin is estimated in advance, and the taste concentration determination value S is calculated as the reciprocal distance:
Dist i =(X i 2 +Y i 2 ) 1/2 (11)
S i =1/Dist i (12)
(5) judging the taste concentration value S i Concentration of taste of the personA determination function (or fitness function) for determining the taste concentration Smell of the individual positions of the fruit flies i :
(6) The mean taste of the Drosophila population Smell was calculated according to equation (14) ave Then, the fruit fly population fitness variance tau is calculated according to the formula (15) 2
(7) If τ is 2 ≤δ 1 And Smellbest>δ 2 Or distributed in [0,1]Random number r between&P, then, firstly, the optimal fruit fly individual (X _ axism) is duplicated according to the formula (10) j ,Y_axism j ) (j =1,2, …, M); secondly, the copied optimal fruit fly individual is mutated according to the formula (11), and the position of the copied optimal fruit fly individual is updated to be a new position (X _ axism) j ,Y_axism j )(j=1,2,…, M):
(8) The new position (X _ axism) is estimated again first according to the following equation j ,Y_axism j ) And the distance Dist 'from the original point, and then calculating a new position taste concentration judgment value S' according to a formula:
(9) the taste concentration determination value S 'calculated again is substituted into the taste concentration determination function to calculate the taste concentration Smell' at the new position.
Smell i '=Function(S i ') (19)
If the R satisfies Smell i '&l < Smellbest > then Smellbest = Smell' j ,X_axis=X_axism' j ,
Y_axis=Y_axism' j ,(j=1,2,…,M).
Judging whether the performance requirements are met, if so, finishing the optimization to obtain a group of optimized parameters of the least square support vector machine;
otherwise, returning to the step (5), and continuing the iteration optimization until the maximum iteration number iter is reached max
As a preferred solution: the soft measurement method further comprises the following steps: 4) And inputting the offline experimental data into a training set regularly, and updating the least square support vector machine model.
The method of the embodiment comprises the following specific implementation steps:
step 1: for the radar object 1, the manipulated variables and the easily measurable variables are selected as the inputs of the model based on the characteristic analysis and the climate analysis. The manipulated variables and easily measurable variables are obtained from the field database 4.
Step 2: and sample data is preprocessed and completed by a data preprocessing module 7.
And step 3: and establishing an initial least square support vector machine model 8 based on model training sample data. Input data is obtained as described in step 2 and output data is obtained from an off-line assay.
And 4, step 4: the parameters of the initial least squares support vector machine model are optimized by the adaptive drosophila optimization algorithm module 9.
And 5: the model updating module 10 periodically inputs offline experimental data into a training set, updates the least square support vector machine model, and optimizes the optimal soft measurement upper computer 5 of the least square support vector machine model based on the adaptive drosophila optimization algorithm to complete establishment.
And 6: and the sea clutter soft measurement value display instrument 6 displays the output of the optimal soft measurement upper computer 5 based on the self-adaptive drosophila optimization algorithm optimized least square support vector machine model, and the optimal soft measurement display of the sea clutter is completed.
The optimal soft measurement upper computer based on the adaptive variation fruit fly optimization minimum two-component support vector machine comprises a data preprocessing module, a wavelet neural network module and a model updating module, and provides a soft measurement method realized by a soft measurement instrument. The method realizes the online optimal soft measurement of the sea clutter, overcomes the random influence caused by human factors, improves the stability of model prediction, and reduces the possibility that the model prediction falls into local optimization.

Claims (2)

1. A sea clutter optimal soft measuring instrument based on a self-adaptive drosophila optimization algorithm optimization least square support vector machine comprises a radar, a field intelligent instrument for measuring easily-measured variables, a control station for measuring operation variables, a field database for storing data and a sea clutter soft measurement value display instrument; the field intelligent instrument and the control station are connected with the radar, and the field intelligent instrument and the control station are connected with a field database, and the radar is characterized in that: the soft measuring instrument further comprises an optimal soft measuring upper computer based on an adaptive fruit fly optimization algorithm optimization least square support vector machine, the field database is connected with the input end of the optimal soft measuring upper computer based on the adaptive fruit fly optimization algorithm optimization least square support vector machine, and the output end of the optimal soft measuring upper computer based on the adaptive fruit fly optimization algorithm optimization least square support vector machine is connected with a sea clutter soft measurement value display instrument; the optimal soft measurement upper computer based on the least square support vector machine optimized by the adaptive fruit fly optimization algorithm comprises:
the data preprocessing module is used for preprocessing the model training samples input from the field database, centralizing the training samples, namely subtracting the average value of the samples, and then normalizing the training samples:
calculating an average value:
calculating the variance:
and (3) standardization:
wherein TX is a training sample, N is the number of training samples,is the mean of the training samples, X is the normalized training sample, σ x To calculate the variance.
And the least square support vector machine module adopts a least square support vector machine for modeling. Model ith normalized training sample X i Target output of O i The least squares support vector machine equates the fitting problem to the following quadratic programming problem by transformation:
where R (w, ξ) is the objective function of the optimization problem,a non-linear mapping function, N is the number of training samples, ξ = { ξ = 1 ,…,ξ N Is the relaxation variable, w is the normal vector of the hyperplane of the least squares support vector machine, b is the corresponding offset, and ω is i I =1, …, N and γ are the weight and penalty factors, respectively, of the least squares support vector machine, whereIs the relaxation variable xi i Estimate of standard deviation, constant c 1 ,c 2 Is usually taken as c 1 =2.5,c 2 =3, from which the ith normalized training sample X can be obtained i The output of (c) is:
wherein, K<·&gt is the kernel function of a least squares support vector machine, where K<·&Taking a linear kernel function; alpha is alpha m M =1, …, N is the corresponding lagrange multiplier.
The self-adaptive fruit fly optimization algorithm module is used for optimizing a punishment factor and an error tolerance value of a least square support vector machine by adopting a self-adaptive fruit fly optimization algorithm, and comprises the following specific steps of:
(1) determining optimization parameters of the self-adaptive drosophila optimization algorithm as a penalty factor and an error tolerance value of a least square support vector machine module, the individual number of particle swarm popsize and the maximum cyclic optimization number iter max Variance threshold delta of group fitness 1 Theoretical optimum value delta 2 And the initial position regions X _ axis, Y _ axis of the p-th particle.
(2) Setting an optimization objective function, converting the optimization objective function into fitness, calculating the fitness function through a corresponding error function, considering that the fitness of the particles with large errors is small, and expressing the fitness function f of the particles p as:
f p =1/(E p +1) (8)
in the formula, E p Is a least squares supportError function of the vector machine model, expressed as:
in the formula (I), the compound is shown in the specification,is the predicted output of the least squares support vector machine model, O i Outputting a target of a least squares support vector machine model; n is the number of training samples;
(3) according to the following formula, the particles are searched,
in the formula, randomValue is the search distance;
(4) for the particle p, the distance Dist from the origin is estimated in advance, and the taste concentration determination value S is calculated as the reciprocal distance:
Dist i =(X i 2 +Y i 2 ) 1/2 (11)
S i =1/Dist i (12)
(5) judging the taste concentration value S i A substitute taste concentration judgment function (or fitness function) for calculating the taste concentration Smell of the individual positions of the fruit flies i :
(6) The mean taste of the Drosophila population Smell was calculated according to equation (14) ave Then, the fruit fly population fitness variance tau is calculated according to the formula (15) 2
(7) If τ is 2 ≤δ 1 And Smellbest>δ 2 Or distributed in [0,1]Random number r between&P, then, firstly, the optimal fruit fly individual (X _ axism) is duplicated according to the formula (10) j ,Y_axism j ) (j =1,2, …, M); secondly, the copied optimal fruit fly individual is mutated according to the formula (11), and the position of the copied optimal fruit fly individual is updated to be a new position (X _ axism) j ,Y_axism j )(j=1,2,…,M):
(8) The new position (X _ axism) is estimated again first according to the following equation j ,Y_axism j ) And the distance Dist 'from the original point, and then calculating a new position taste concentration judgment value S' according to a formula:
(9) the taste concentration determination value S 'calculated again is substituted into the taste concentration determination function to calculate the taste concentration Smell' at the new position.
Smell i '=Function(S i ') (19)
If the R satisfies Smell i '&l < Smellbest > then Smellbest = Smell' j ,X_axis=X_axism' j ,Y_axis=Y_axism' j ,(j=1,2,…,M).
Judging whether the performance requirements are met, if so, ending the optimization, and obtaining a group of optimized parameters of the least square support vector machine;
otherwise, returning to the step (5), and continuing to iterate and optimize until reaching the maximum iteration number iter max
And the model updating module is used for updating the model on line, inputting the offline verification data into the training set periodically and updating the least square support vector machine model.
2. The method for realizing the soft measurement of the sea clutter optimal soft measurement instrument based on the adaptive drosophila optimization algorithm optimized least square support vector machine according to the claim 1, is characterized in that: the soft measurement method comprises the following steps:
1) Selecting an operation variable and an easily-measured variable as the input of a model for the radar object according to characteristic analysis and climate analysis, wherein the operation variable and the easily-measured variable are obtained by a field database;
2) Preprocessing a model training sample input from a field database, centralizing the training sample, namely subtracting the average value of the sample, and then normalizing the training sample so that the average value is 0 and the variance is 1. The process is accomplished using the following algorithm:
2.1 Calculate the mean value:
2.2 Calculate variance:
2.3 Normalization:
wherein TX is a training sample, N is the number of training samples,is the mean of the training samples, X is the normalized training sample, σ x To calculate the variance.
3) For pre-processing slave dataAnd modeling the training sample transmitted by the processing module by adopting a least square support vector machine. Let the ith normalized training sample X i Target output of O i The least squares support vector machine equates the modeling problem to the following quadratic programming problem by transformation:
wherein R (w, ξ) is the objective function of the optimization problem,a non-linear mapping function, N is the number of training samples, ξ = { ξ = 1 ,…,ξ N Is the relaxation variable, w is the normal vector of the hyperplane of the least squares support vector machine, b is the corresponding offset, and ω is i I =1, …, N and γ are the weight and penalty factors, respectively, of the least squares support vector machine, whereIs the relaxation variable xi i Estimate of standard deviation, constant c 1 ,c 2 Is usually taken as c 1 =2.5,c 2 =3, from which normalized training sample X can be obtained i The output of (c) is:
wherein, K<·&gt is the kernel function of a least squares support vector machine, where K<·&Taking a linear kernel function; alpha is alpha m M =1, …, N is the corresponding lagrange multiplier.
4) Optimizing the punishment factor and the error tolerance value of the least square support vector machine by adopting a self-adaptive drosophila optimization algorithm, and specifically comprising the following steps of:
(1) determining optimization parameters of the self-adaptive drosophila optimization algorithm as a penalty factor and an error tolerance value of a least square support vector machine module, the individual number of particle swarm popsize and the maximum cyclic optimization number iter max Variance threshold delta of group fitness 1 Theoretical optimum value delta 2 And the initial position regions X _ axis, Y _ axis of the p-th particle.
(2) Setting an optimization objective function, converting the optimization objective function into fitness, calculating the fitness function through a corresponding error function, considering that the fitness of the particles with large errors is small, and expressing the fitness function f of the particles p as:
f p =1/(E p +1) (8)
in the formula, E p Is the error function of the least squares support vector machine model, expressed as:
in the formula (I), the compound is shown in the specification,is the predicted output of the least squares support vector machine model, O i Outputting a target of a least squares support vector machine model; n is the number of training samples;
(3) according to the following formula, the particles are searched,
wherein, random Value is the search distance;
(4) for the particle p, the distance Dist from the origin is estimated in advance, and the taste concentration determination value S is calculated as the reciprocal distance:
Dist i =(X i 2 +Y i 2 ) 1/2 (11)
S i =1/Dist i (12)
(5) judging the taste concentration value S i A substitute taste concentration judgment function (or fitness function) for determining the taste concentration Smell of the individual positions of the drosophila i :
(6) The mean taste of the Drosophila population Smell was calculated according to equation (14) ave Then calculating the fruit fly population fitness variance tau according to the formula (15) 2
(7) If τ is 2 ≤δ 1 And Smellbest>δ 2 Or distributed in [0,1]Random number r between&P, then, firstly, the optimal fruit fly individual (X _ axism) is duplicated according to the formula (10) j ,Y_axism j ) (j =1,2, …, M); secondly, the copied optimal fruit fly individual is mutated according to the formula (11), and the position of the copied optimal fruit fly individual is updated to be a new position (X _ axism) j ,Y_axism j )(j=1,2,…,M):
(8) The new position (X _ axism) is estimated again first according to the following equation j ,Y_axism j ) And the distance Dist 'from the original point, and then calculating a new position taste concentration judgment value S' according to a formula:
(9) the taste concentration determination value S 'calculated again is substituted into the taste concentration determination function to calculate the taste concentration Smell' at the new position.
Smell i '=Function(S i ') (19)
If R satisfies Smell i '&l < Smellbest > then Smellbest = Smell' j ,X_axis=X_axism' j ,Y_axis=Y_axism' j ,(j=1,2,…,M).
Judging whether the performance requirements are met, if so, finishing the optimization to obtain a group of optimized parameters of the least square support vector machine;
otherwise, returning to the step (5), and continuing the iteration optimization until the maximum iteration number iter is reached max
5) And inputting the offline experimental data into a training set regularly, and updating the least square support vector machine model.
CN201711117137.5A 2017-11-13 2017-11-13 Sea clutter optimal soft survey instrument and method based on adaptive drosophila optimized algorithm Optimized Least Square Support Vector Pending CN107703493A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711117137.5A CN107703493A (en) 2017-11-13 2017-11-13 Sea clutter optimal soft survey instrument and method based on adaptive drosophila optimized algorithm Optimized Least Square Support Vector

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711117137.5A CN107703493A (en) 2017-11-13 2017-11-13 Sea clutter optimal soft survey instrument and method based on adaptive drosophila optimized algorithm Optimized Least Square Support Vector

Publications (1)

Publication Number Publication Date
CN107703493A true CN107703493A (en) 2018-02-16

Family

ID=61178502

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711117137.5A Pending CN107703493A (en) 2017-11-13 2017-11-13 Sea clutter optimal soft survey instrument and method based on adaptive drosophila optimized algorithm Optimized Least Square Support Vector

Country Status (1)

Country Link
CN (1) CN107703493A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108536955A (en) * 2018-04-08 2018-09-14 西安邮电大学 A kind of hybrid forecasting method of underwater sound signal
CN108681249A (en) * 2018-05-14 2018-10-19 浙江大学 A kind of probabilistic type that parameter independently optimizes output fault diagnosis system
CN108681250A (en) * 2018-05-14 2018-10-19 浙江大学 A kind of improvement machine learning fault diagnosis system based on colony intelligence optimization
CN108681248A (en) * 2018-05-14 2018-10-19 浙江大学 A kind of autonomous learning fault diagnosis system that parameter is optimal
CN109002918A (en) * 2018-07-16 2018-12-14 国网浙江省电力有限公司经济技术研究院 Based on drosophila optimization algorithm-support vector machines electricity sales amount prediction technique
CN109598083A (en) * 2018-12-18 2019-04-09 江南大学 A kind of glutamic acid fermentation process soft-measuring modeling method based on the twin support vector regression of drosophila algorithm optimization
CN110031847A (en) * 2018-09-29 2019-07-19 浙江师范大学 The dynamic of wavelet transformation and support vector machines combination radar reflectivity is short to face Quantitative Precipitation estimating and measuring method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080007450A1 (en) * 2006-07-07 2008-01-10 Jonathan Yedidia Method and system for determining unwrapped phases from noisy two-dimensional wrapped-phase images
CN102147463A (en) * 2011-03-03 2011-08-10 浙江大学 System and method for forecasting Qunzhi radar sea clutters
CN102147464A (en) * 2011-03-03 2011-08-10 浙江大学 Intelligent system and method for forecasting robust radar sea clutter
CN102147465A (en) * 2011-03-03 2011-08-10 浙江大学 System and method for detecting sea target by chaos optimizing radar
CN103675006A (en) * 2013-09-22 2014-03-26 浙江大学 Least-squares-based industrial melt index soft measuring meter and method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080007450A1 (en) * 2006-07-07 2008-01-10 Jonathan Yedidia Method and system for determining unwrapped phases from noisy two-dimensional wrapped-phase images
CN102147463A (en) * 2011-03-03 2011-08-10 浙江大学 System and method for forecasting Qunzhi radar sea clutters
CN102147464A (en) * 2011-03-03 2011-08-10 浙江大学 Intelligent system and method for forecasting robust radar sea clutter
CN102147465A (en) * 2011-03-03 2011-08-10 浙江大学 System and method for detecting sea target by chaos optimizing radar
CN103675006A (en) * 2013-09-22 2014-03-26 浙江大学 Least-squares-based industrial melt index soft measuring meter and method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王文川: "基于人工智能算法的聚丙烯熔融指数预报建模与优化研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108536955A (en) * 2018-04-08 2018-09-14 西安邮电大学 A kind of hybrid forecasting method of underwater sound signal
CN108681249A (en) * 2018-05-14 2018-10-19 浙江大学 A kind of probabilistic type that parameter independently optimizes output fault diagnosis system
CN108681250A (en) * 2018-05-14 2018-10-19 浙江大学 A kind of improvement machine learning fault diagnosis system based on colony intelligence optimization
CN108681248A (en) * 2018-05-14 2018-10-19 浙江大学 A kind of autonomous learning fault diagnosis system that parameter is optimal
CN109002918A (en) * 2018-07-16 2018-12-14 国网浙江省电力有限公司经济技术研究院 Based on drosophila optimization algorithm-support vector machines electricity sales amount prediction technique
CN110031847A (en) * 2018-09-29 2019-07-19 浙江师范大学 The dynamic of wavelet transformation and support vector machines combination radar reflectivity is short to face Quantitative Precipitation estimating and measuring method
CN110031847B (en) * 2018-09-29 2022-08-23 浙江师范大学 Dynamic short-term quantitative rainfall estimation method combining wavelet transformation and support vector machine with radar reflectivity
CN109598083A (en) * 2018-12-18 2019-04-09 江南大学 A kind of glutamic acid fermentation process soft-measuring modeling method based on the twin support vector regression of drosophila algorithm optimization

Similar Documents

Publication Publication Date Title
CN107703493A (en) Sea clutter optimal soft survey instrument and method based on adaptive drosophila optimized algorithm Optimized Least Square Support Vector
Valavanis et al. Modelling of essential fish habitat based on remote sensing, spatial analysis and GIS
CN112001270B (en) Ground radar automatic target classification and identification method based on one-dimensional convolutional neural network
CN110689179A (en) Water bloom prediction method based on space-time sequence mixed model
CN107703491A (en) Sea clutter optimal soft survey instrument and method based on improved drosophila optimized algorithm optimization RBF neural
Waters et al. Data assimilation of partitioned HF radar wave data into Wavewatch III
CN108613645B (en) A kind of Pb-Zn deposits absorbing well, absorption well surveying on sludge thickness method based on parameter Estimation
CN107907872A (en) Sea clutter optimal soft survey instrument and method based on TSP question drosophila optimization algorithm optimization RBF neural
Zhao et al. XGB model: Research on evaporation duct height prediction based on XGBoost algorithm
Dorn Fishing behavior of factory trawlers: a hierarchical model of information processing and decision-making
CN108983179A (en) A kind of radar marine target detection system of colony intelligence agility
CN110501683A (en) A kind of extra large land Clutter Classification method based on 4 D data feature
Sattari et al. Estimation of monthly and seasonal precipitation: A comparative study using data-driven methods versus hybrid approach
CN108896996A (en) A kind of Pb-Zn deposits absorbing well, absorption well water sludge interface ultrasonic echo signal classification method based on random forest
CN115906667B (en) Ocean environment parameter inversion model construction method and device
CN107818224A (en) Sea clutter optimal soft survey instrument and method based on drosophila optimized algorithm Optimization of Wavelet neutral net
CN108983178A (en) A kind of Intelligent radar sea target detection system that agility is adaptive
CN107942304A (en) Sea clutter optimal soft survey instrument and method based on drosophila optimization algorithm optimization least square method supporting vector machine
CN107942313A (en) Sea clutter optimal soft survey instrument and method based on TSP question drosophila optimization algorithm optimization wavelet neural network
Liu et al. Prediction of the Antarctic marginal ICE zone extent based upon its multifractal property
CN107843879A (en) Sea clutter optimal soft survey instrument and method based on free searching algorithm Optimization of Wavelet neutral net
CN107832831A (en) Sea clutter optimal soft survey instrument and method based on free searching algorithm optimization RBF neural
CN107831482A (en) Sea clutter optimal soft survey instrument and method based on improved free searching algorithm optimization RBF neural
Mathew et al. Inverse estimation of wind from the waves measured by high-frequency radar
De et al. Preference of echo features for classification of seafloor sediments using neural networks

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20180216

RJ01 Rejection of invention patent application after publication