CN108983178A - A kind of Intelligent radar sea target detection system that agility is adaptive - Google Patents
A kind of Intelligent radar sea target detection system that agility is adaptive Download PDFInfo
- Publication number
- CN108983178A CN108983178A CN201810689688.7A CN201810689688A CN108983178A CN 108983178 A CN108983178 A CN 108983178A CN 201810689688 A CN201810689688 A CN 201810689688A CN 108983178 A CN108983178 A CN 108983178A
- Authority
- CN
- China
- Prior art keywords
- radar
- module
- value
- layer
- rule
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/417—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section involving the use of neural networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
Abstract
The invention discloses a kind of Intelligent radar sea target detection systems that agility is adaptive, including radar, database and host computer;Radar, database and host computer are sequentially connected, and radar is irradiated detected sea area, and by radar sea clutter data storage into database, and host computer models the sea clutter data in database and detects target;The host computer includes data preprocessing module, dynamic fuzzy neural network modeling module, differential evolution algorithm optimization module, module of target detection, model modification module and result display module.And propose a kind of Intelligent radar sea target detection method based on differential evolution algorithm optimization dynamic fuzzy neural network.The present invention provides a kind of radar marine target detection system and method for realizing that agility is adaptive.
Description
Technical field
The present invention relates to radar data process fields, particularly, are related to a kind of Intelligent radar sea mesh that agility is adaptive
Mark detection system.
Background technique
Sea clutter, i.e. the radar raster-displaying echo from sea.In recent decades, with the depth recognized sea clutter
Entering, the country such as Germany, Norway attempts to obtain radar wave image using radar observation sea clutter in succession come inverting Wave Information, with
The real time information about sea state, such as wave height, direction and the period of wave are obtained, thus further to marine small objects
It is detected, this has a very important significance offshore activities.
Naval target detection technique has consequence, and providing accurate target decision is to the important of extra large radar work
One of task.Radar automatic checkout system makes judgement under given detection threshold value according to decision rule, and strong sea clutter is past
Toward the main interference for becoming weak target signal.Detection of the radar under marine environment will be directly influenced by how handling sea clutter
Ability: 1) ice of navigation by recognition buoy, small pieces swims in the greasy dirt on sea, these may carry out potential crisis to navigation band;
2) monitoring illegal fishing is an important task of environmental monitoring.
In traditional target detection, sea clutter is considered as that a kind of noise of interference navigation is removed.However, in radar
When to extra large observed object, faint Moving Target Return is usually buried in sea clutter, and signal to noise ratio is lower, and radar is not easy to detect
Target, while a large amount of spikes of sea clutter also will cause serious false-alarm, produce bigger effect to the detection performance of radar.For each
For kind for sea police's ring and early warning radar, the main target of research is to improve the detectability of target under sea clutter background.Therefore,
Not only there is important theory significance and practical significance, but also be also the difficult point and hot spot of domestic and international naval target detection.
Summary of the invention
In order to overcome the shortcomings of that existing radar method for detecting targets at sea detection efficiency is not high, adaptive ability is poor, this
Invention provides a kind of Intelligent radar sea target detection system that agility is adaptive.
The technical solution adopted by the present invention to solve the technical problems is: a kind of Intelligent radar sea mesh that agility is adaptive
Mark detection system, including radar, database and host computer;The host computer includes data preprocessing module, dynamic fuzzy
Neural net model establishing module, differential evolution algorithm optimization module, module of target detection, model modification module and mould as the result is shown
Block, in which:
Data preprocessing module: pre-processing the radar sea clutter data of database input, complete using following process
At:
(1) N number of radar sea clutter echo-signal amplitude x is acquired from databaseiAs training sample, i=1,2 ..., N;
(2) training sample is normalized, obtains normalization amplitude
Wherein, min x indicates the minimum value in training sample, and max x indicates the maximum value in training sample;
(3) training sample after normalization is reconstructed, respectively obtains input matrix X and corresponding output square Y:
Wherein, D indicates reconstruct dimension, and D is natural number, and the value range of D < N, D are 50-70.
Dynamic fuzzy neural network modeling module: it to establish forecasting model, is completed using following process:
(1) network structure: p-th of training sample X is setp=[xp1,xp2,…,xpn], wherein n is the number of input variable.
The structure of dynamic fuzzy neural network is formed by 5 layers, and details are as follows:
1st layer is input layer, and each node respectively indicates the linguistic variable of an input.P-th of training sample known to then
XpThere is the linguistic variable x of n inputp1,xp2,…,xpn。
2nd layer is membership function layer, and each node respectively represents a membership function, which is with following height
This function representation:
Wherein, μijIt is xpiJ-th of membership function, cijIt is xpiJ-th of Gauss member function center, σjIt is xpi's
The width of j-th of Gauss member function, n are the numbers of input variable, and u is the quantity of membership function, also represent the total rule of system
Then count.
3rd layer is T- norm layer, and each node respectively represents the part IF- in a possible fuzzy rule.Therefore, should
Node layer reflects number of fuzzy rules;Meanwhile this layer of each node also represents a RBF unit, so number of fuzzy rules with
RBF number of nodes is equal.J-th of rule RjOutput be
Wherein,It is the center of j-th of RBF unit.
4th layer is normalization layer, and the node of this layer is known as N node.N number of nodes is equal with fuzzy rule number of nodes.J-th
Node NjOutput be
5th layer is output layer, and each node in this layer respectively indicates an output variable, which is all input letters
Number superposition:
Wherein, y (Xp) it is the corresponding reality output of p-th of training sample, ωkIt is the part THEN- or k-th of rule
Connection weight.
ωk=αk0+αk1xp1+…+αknxpnK=1,2 ..., u (7)
Wherein, αk0、αk1、…、αknIt is the consequent parameter of network.
Wushu (4), formula (5), formula (7) substitute into formula (6), then obtain:
(2) dynamic fuzzy neural network learning algorithm: for every group of observation data (Xp,tp), wherein XpIt is p-th of training
Sample, tpIt is the corresponding desired output of p-th of training sample, it is online to carry out system structure and parameter identification, implement step
It is as follows:
(2.1) determination can accommodate the effective radius k on boundaryd, according to desired precision previously selected keIt is worth, according to rule
Then importance and preset constant kerr。
(2.2) as the first group of observation data (X that arrive1,t1) when generate the first rule, and determine premise parameter and conclusion
Parameter.
(2.3) since second group of observation data, often carry out one group of observation data, calculate dmin.If dmin>kd, then into
Enter step 4;Otherwise, into step 5.
dmin=argmin (dp(j))=argmin (| | Xp-Cj| |) j=1,2 ..., u (9)
Wherein, XpIt is p-th of training sample, CjThe center of existing RBF unit, u be existing number of fuzzy rules or
The quantity of RBF unit.
(2.4) calculate | | ep||.If | | ep||>ke, then new rule is generated, otherwise adjust automatically well-regulated ginseng
Number.
||ep| |=| | tp-yp|| (10)
Wherein, tpIt is the observation desired output of data of pth group, ypIt is the observation actual output of data of pth group.
(2.5) calculate | | ep||.If | | ep||>ke, then the Gaussian function width of rule is adjusted, then adjust consequent parameter,
Otherwise the conclusion ginseng of adjust automatically rule.
(2.6) η is calculatedj, ηjIt is the importance of j-th strip rule.If ηj>kerr, then j-th strip rule is deleted, otherwise automatically
Adjust the consequent parameter of rule.
(2.7) judge whether to terminate, if returning to step 3 without if, otherwise terminate entire learning process.
Differential evolution algorithm optimization module: for using differential evolution algorithm holding to dynamic fuzzy neural network module
Receive the effective radius k on boundarydValue, the previously selected k according to desired precisioneValue, the preset constant according to regular importance
kerrValue optimizes, the specific steps are as follows:
(1) algorithm initialization initializes population in the solution space of problemIt is a
Body Xi(k)=[xi,1(k),xi,2(k),…,xi,D(k)], i=1,2 ..., NpThe solution of characterization problems, D are the dimension of solution space, k
Indicate the number of iterations, NpFor population invariable number.Initialization population should cover entire search space as much as possible, and defined
Individual is unanimously randomized in the search space of minimum and maximum bound of parameter limitation, setting minimum and maximum boundary difference
For Xmin={ xmin,1,…,xmin,DAnd Xmax={ xmax,1,…,xmax,D}.In the number of iterations k=1, j-th of i-th of individual
The initial value of component generates as the following formula:
xi,j(0)=xmin,j+rand(0,1)·(xmax,j-xmin,j) (11)
Wherein rand (0,1) is equally distributed random number between 0 to 1.
(2) mutation operation.After being initialized to population, for each individual Xi(k) change can correspondingly be generated
Incorgruous amount Vi(k)=[vi,1(k),vi,2(k),…,vi,D(k)], individual Xi(k) object vector, the expression of mutation operation are also referred to as
Formula is as follows:
Vi(k)=Xbest(k)+F·(Xr1(k)-Xr2(k)),1≤r1≠r2≠i≤Np (12)
Wherein, XbestIt (k) is the individual vector in population in current kth time iteration with optimal adaptation angle value, ratio
Factor F is the positive control parameter for scaling difference vector.
(3) crossover operation.For each pair of object vector Xi(k) and its accordingly make a variation vector Vi(k) crossover operation is carried out,
To generate trial vector Ui(k)=[ui,1(k),ui,2(k),…,ui,D(k)].The binomial interleaved scheme is carried out according to the following formula:
Wherein, CR is a customized crossing-over rate, and for value range usually between 0 to 1, it controls the multiplicity of population
Property, and algorithm is avoided to fall into local optimum.If the value of certain parameters of newly-generated trial vector has been more than above and below corresponding
Limit, then reinitialize it in predefined scope uniformly at random.Then, the fitness of all trial vectors is assessed
Value.
(4) selection operation.The fitness of trial vector is compared by the operation with corresponding object vector, and selection is wherein
Preferably solution.For minimization problem, the expression formula of selection operation is as follows:
Wherein, f (x) is fitness value.In the evolutionary process of every generation, each individual vector is as target individual one
Secondary, algorithm retains defect individual by constantly iterating to calculate, and eliminates worst individual, and guiding search process is to globally optimal solution
It approaches.
Module of target detection: it to carry out target detection, is completed using following process:
(1) D sea clutter echo-signal amplitude is acquired in sampling instant t obtain TX=[xt-D+1,…,xt], xt-D+1It indicates
The sea clutter echo-signal amplitude of t-D+1 sampling instant, xtIndicate the sea clutter echo-signal amplitude of t sampling instant;
(2) it is normalized:
(3) the sea clutter predicted value that sampling instant (t+1) is calculated in Fuzzy Wavelet Network modeling module is substituted into.
(4) difference e for calculating sea clutter predicted value and radar return measured value calculates control limit Qα:
Wherein, α is confidence level, θ1,θ2,θ3,h0It is intermediate variable,Indicate the i of j-th of characteristic value of covariance matrix
Power, k are sample dimension, CαIt is the statistics that normal distribution confidence level is α;
(5) it carries out detection judgement: working as e2Difference is greater than control limit QαWhen, there are targets for the point, otherwise without target.
The adaptive radar marine target detection system of the agility, the host computer further include: model modification module is used
To acquire data by the sampling time interval of setting, by obtained measured data compared with model prediction value, if relative error
Greater than 10%, then training sample data are added in new data, update forecasting model.And result display module, to by target
The testing result of detection module is shown in host computer.
Beneficial effects of the present invention are mainly manifested in: the present invention detects naval target, overcomes in existing Radar Sea
The deficiency that object detection method detection efficiency is not high, adaptive ability is poor, the chaotic characteristic foundation for radar sea clutter are dynamic
Morphotype pastes neural Network Model for Forecasting, differential evolution optimization algorithm is further introduced into, thus the agility under realizing sea clutter background
Adaptive targets detection.
Detailed description of the invention
Fig. 1 is the hardware structure diagram of system proposed by the invention;
Fig. 2 is the functional block diagram of host computer proposed by the invention.
Specific embodiment
The present invention is illustrated below according to attached drawing.
Referring to Fig.1, Fig. 2, a kind of Intelligent radar sea target detection system that agility is adaptive, including radar 1, database
2 and host computer 3, radar 1, database 2 and host computer 3 be sequentially connected, 1 pair of detected sea area of the radar is irradiated, and will
To the database 2, the host computer 3 includes: radar sea clutter data storage
Data preprocessing module 4: pre-processing the radar sea clutter data of database input, complete using following process
At:
(1) N number of radar sea clutter echo-signal amplitude x is acquired from databaseiAs training sample, i=1,2 ..., N;
(2) training sample is normalized, obtains normalization amplitude
Wherein, min x indicates the minimum value in training sample, and max x indicates the maximum value in training sample;
(3) training sample after normalization is reconstructed, respectively obtains input matrix X and corresponding output square Y:
Wherein, D indicates reconstruct dimension, and D is natural number, and the value range of D < N, D are 50-70.
Dynamic fuzzy neural network modeling module 5: it to establish forecasting model, is completed using following process:
(1) network structure: p-th of training sample X is setp=[xp1,xp2,…,xpn], wherein n is the number of input variable.
The structure of dynamic fuzzy neural network is formed by 5 layers, and details are as follows:
1st layer is input layer, and each node respectively indicates the linguistic variable of an input.P-th of training sample known to then
XpThere is the linguistic variable x of n inputp1,xp2,…,xpn。
2nd layer is membership function layer, and each node respectively represents a membership function, which is with following height
This function representation:
Wherein, μijIt is xpiJ-th of membership function, cijIt is xpiJ-th of Gauss member function center, σjIt is xpi's
The width of j-th of Gauss member function, n are the numbers of input variable, and u is the quantity of membership function, also represent the total rule of system
Then count.
3rd layer is T- norm layer, and each node respectively represents the part IF- in a possible fuzzy rule.Therefore, should
Node layer reflects number of fuzzy rules;Meanwhile this layer of each node also represents a RBF unit, so number of fuzzy rules with
RBF number of nodes is equal.J-th of rule RjOutput be
Wherein,It is the center of j-th of RBF unit.
4th layer is normalization layer, and the node of this layer is known as N node.N number of nodes is equal with fuzzy rule number of nodes.J-th
Node NjOutput be
5th layer is output layer, and each node in this layer respectively indicates an output variable, which is all input letters
Number superposition:
Wherein, y (Xp) it is the corresponding reality output of p-th of training sample, ωkIt is the part THEN- or k-th of rule
Connection weight.
ωk=αk0+αk1xp1+…+αknxpnK=1,2 ..., u (7)
Wherein, αk0、αk1、…、αknIt is the consequent parameter of network.
Wushu (4), formula (5), formula (7) substitute into formula (6), then obtain:
(2) dynamic fuzzy neural network learning algorithm: for every group of observation data (Xp,tp), wherein XpIt is p-th of training
Sample, tpIt is the corresponding desired output of p-th of training sample, it is online to carry out system structure and parameter identification, implement step
It is as follows:
(2.1) determination can accommodate the effective radius k on boundaryd, according to desired precision previously selected keIt is worth, according to rule
Then importance and preset constant kerr。
(2.2) as the first group of observation data (X that arrive1,t1) when generate the first rule, and determine premise parameter and conclusion
Parameter.
(2.3) since second group of observation data, often carry out one group of observation data, calculate dmin.If dmin>kd, then into
Enter step 4;Otherwise, into step 5.
dmin=argmin (dp(j))=argmin (| | Xp-Cj| |) j=1,2 ..., u (9)
Wherein, XpIt is p-th of training sample, CjThe center of existing RBF unit, u be existing number of fuzzy rules or
The quantity of RBF unit.
(2.4) calculate | | ep||.If | | ep||>ke, then new rule is generated, otherwise adjust automatically well-regulated ginseng
Number.
||ep| |=| | tp-yp|| (10)
Wherein, tpIt is the observation desired output of data of pth group, ypIt is the observation actual output of data of pth group.
(2.5) calculate | | ep||.If | | ep||>ke, then the Gaussian function width of rule is adjusted, then adjust consequent parameter,
Otherwise the conclusion ginseng of adjust automatically rule.
(2.6) η is calculatedj, ηjIt is the importance of j-th strip rule.If ηj>kerr, then j-th strip rule is deleted, otherwise automatically
Adjust the consequent parameter of rule.
(2.7) judge whether to terminate, if returning to step 3 without if, otherwise terminate entire learning process.
Differential evolution algorithm optimization module 6: for using differential evolution algorithm to dynamic fuzzy neural network module can
Accommodate the effective radius k on boundarydValue, the previously selected k according to desired precisioneIt is worth, is preset normal according to regular importance
Number kerrValue optimizes, the specific steps are as follows:
(1) algorithm initialization initializes population in the solution space of problemIt is a
Body Xi(k)=[xi,1(k),xi,2(k),…,xi,D(k)], i=1,2 ..., NpThe solution of characterization problems, D are the dimension of solution space, k
Indicate the number of iterations, NpFor population invariable number.Initialization population should cover entire search space as much as possible, and defined
Individual is unanimously randomized in the search space of minimum and maximum bound of parameter limitation, setting minimum and maximum boundary difference
For Xmin={ xmin,1,…,xmin,DAnd Xmax={ xmax,1,…,xmax,D}.In the number of iterations k=1, j-th of i-th of individual
The initial value of component generates as the following formula:
xi,j(0)=xmin,j+rand(0,1)·(xmax,j-xmin,j) (11)
Wherein rand (0,1) is equally distributed random number between 0 to 1.
(2) mutation operation.After being initialized to population, for each individual Xi(k) change can correspondingly be generated
Incorgruous amount Vi(k)=[vi,1(k),vi,2(k),…,vi,D(k)], individual Xi(k) object vector, the expression of mutation operation are also referred to as
Formula is as follows:
Vi(k)=Xbest(k)+F·(Xr1(k)-Xr2(k)),1≤r1≠r2≠i≤Np (12)
Wherein, XbestIt (k) is the individual vector in population in current kth time iteration with optimal adaptation angle value, ratio
Factor F is the positive control parameter for scaling difference vector.
(3) crossover operation.For each pair of object vector Xi(k) and its accordingly make a variation vector Vi(k) crossover operation is carried out,
To generate trial vector Ui(k)=[ui,1(k),ui,2(k),…,ui,D(k)].The binomial interleaved scheme is carried out according to the following formula:
Wherein, CR is a customized crossing-over rate, and for value range usually between 0 to 1, it controls the multiplicity of population
Property, and algorithm is avoided to fall into local optimum.If the value of certain parameters of newly-generated trial vector has been more than above and below corresponding
Limit, then reinitialize it in predefined scope uniformly at random.Then, the fitness of all trial vectors is assessed
Value.
(4) selection operation.The fitness of trial vector is compared by the operation with corresponding object vector, and selection is wherein
Preferably solution.For minimization problem, the expression formula of selection operation is as follows:
Wherein, f (x) is fitness value.In the evolutionary process of every generation, each individual vector is as target individual one
Secondary, algorithm retains defect individual by constantly iterating to calculate, and eliminates worst individual, and guiding search process is to globally optimal solution
It approaches.
Module of target detection 7: it to carry out target detection, is completed using following process:
(1) D sea clutter echo-signal amplitude is acquired in sampling instant t obtain TX=[xt-D+1,…,xt], xt-D+1It indicates
The sea clutter echo-signal amplitude of t-D+1 sampling instant, xtIndicate the sea clutter echo-signal amplitude of t sampling instant;
(2) it is normalized:
(3) the sea clutter predicted value that sampling instant (t+1) is calculated in Fuzzy Wavelet Network modeling module is substituted into.
(4) difference e for calculating sea clutter predicted value and radar return measured value calculates control limit Qα:
Wherein, α is confidence level, θ1,θ2,θ3,h0It is intermediate variable,Indicate i times of j-th of characteristic value of covariance matrix
Side, k is sample dimension, CαIt is the statistics that normal distribution confidence level is α;
(5) it carries out detection judgement: working as e2Difference is greater than control limit QαWhen, there are targets for the point, otherwise without target.
Model modification module 8 acquires data by the sampling time interval of setting, by obtained measured data and model prediction
Value compares, if relative error is greater than 10%, training sample data is added in new data, update forecasting model.
Result display module 9 shows the testing result of module of target detection in host computer.
The hardware components of the host computer 3 include: I/O element, for the acquisition of data and the transmitting of information;Data storage
Device, data sample and operating parameter etc. needed for storage operation;The software program of functional module is realized in program storage, storage;
Arithmetic unit executes program, realizes specified function;Display module shows the parameter and testing result of setting.
The embodiment of the present invention is used to illustrate the present invention, rather than limits the invention, in spirit of the invention
In scope of protection of the claims, to any modifications and changes that the present invention makes, protection scope of the present invention is both fallen within.
Claims (5)
1. a kind of radar marine target detection system that agility is adaptive, including radar, database and host computer;Radar, number
It is sequentially connected according to library and host computer, radar is irradiated detected sea area, and by radar sea clutter data storage to database
In, host computer models the sea clutter data in database and detects target;The host computer includes data prediction
Module, dynamic fuzzy neural network modeling module, differential evolution algorithm optimization module, module of target detection, model modification module
And result display module.
2. quick adaptive radar marine target detection system according to claim 1, which is characterized in that the data are pre-
Processing module pre-processes the radar sea clutter data that database inputs, and is completed using following process:
(1) N number of radar sea clutter echo-signal amplitude x is acquired from databaseiAs training sample, i=1,2 ..., N;
(2) training sample is normalized, obtains normalization amplitude
Wherein, min x indicates the minimum value in training sample, and max x indicates the maximum value in training sample;
(3) training sample after normalization is reconstructed, respectively obtains input matrix X and corresponding output square Y:
Wherein, D indicates reconstruct dimension, and D is natural number, and the value range of D < N, D are 50-70.
3. quick adaptive radar marine target detection system according to claim 1, which is characterized in that the dynamic analog
Neural net model establishing module is pasted, to establish forecasting model, is completed using following process:
(1) network structure: p-th of training sample X is setp=[xp1,xp2,…,xpn], wherein n is the number of input variable.
The structure of dynamic fuzzy neural network is formed by 5 layers, and details are as follows:
1st layer is input layer, and each node respectively indicates the linguistic variable of an input.P-th of training sample X known to thenpThere is n
The linguistic variable x of a inputp1,xp2,…,xpn。
2nd layer is membership function layer, and each node respectively represents a membership function, which is with following Gaussian function
What number indicated:
Wherein, μijIt is xpiJ-th of membership function, cijIt is xpiJ-th of Gauss member function center, σjIt is xpiJ-th
The width of Gauss member function, n are the numbers of input variable, and u is the quantity of membership function, also represent the total regular number of system.
3rd layer is T- norm layer, and each node respectively represents the part IF- in a possible fuzzy rule.Therefore, the layer section
Point reflects number of fuzzy rules;Meanwhile this layer of each node also represents a RBF unit, so number of fuzzy rules and RBF are saved
It counts equal.J-th of rule RjOutput be
Wherein,It is the center of j-th of RBF unit.
4th layer is normalization layer, and the node of this layer is known as N node.N number of nodes is equal with fuzzy rule number of nodes.J-th of node
NjOutput be
5th layer is output layer, and each node in this layer respectively indicates an output variable, which is all input signals
Superposition:
Wherein, y (Xp) it is the corresponding reality output of p-th of training sample, ωkIt is the connection of the part THEN- or k-th of rule
Power.
ωk=αk0+αk1xp1+…+αknxpnK=1,2 ..., u (7)
Wherein, αk0、αk1、…、αknIt is the consequent parameter of network.
Wushu (4), formula (5), formula (7) substitute into formula (6), then obtain:
(2) dynamic fuzzy neural network learning algorithm: for every group of observation data (Xp,tp), wherein XpIt is p-th of trained sample
This, tpIt is the corresponding desired output of p-th of training sample, online to carry out system structure and parameter identification, specific implementation step is such as
Under:
(2.1) determination can accommodate the effective radius k on boundaryd, according to desired precision previously selected keIt is worth, according to rule weight
The property wanted and preset constant kerr。
(2.2) as the first group of observation data (X that arrive1,t1) when generate the first rule, and determine premise parameter and consequent parameter.
(2.3) since second group of observation data, often carry out one group of observation data, calculate dmin.If dmin>kd, then enter the 4th
Step;Otherwise, into step 5.
dmin=arg min (dp(j))=arg min (| | Xp-Cj| |) j=1,2 ..., u (9)
Wherein, XpIt is p-th of training sample, CjIt is the center of existing RBF unit, u is existing number of fuzzy rules or RBF
The quantity of unit.
(2.4) calculate | | ep||.If | | ep||>ke, then new rule is generated, otherwise adjust automatically well-regulated parameter.
||ep| |=| | tp-yp|| (10)
Wherein, tpIt is the observation desired output of data of pth group, ypIt is the observation actual output of data of pth group.
(2.5) calculate | | ep||.If | | ep||>ke, then the Gaussian function width of rule is adjusted, then adjust consequent parameter, otherwise
The conclusion of adjust automatically rule is joined.
(2.6) η is calculatedj, ηjIt is the importance of j-th strip rule.If ηj>kerr, then j-th strip rule is deleted, otherwise adjust automatically
The consequent parameter of rule.
(2.7) judge whether to terminate, if returning to step 3 without if, otherwise terminate entire learning process.
4. quick adaptive radar marine target detection system according to claim 1, which is characterized in that the difference into
Change algorithm optimization module, for accommodate boundary effectively half using differential evolution algorithm to dynamic fuzzy neural network module
Diameter kdValue, the previously selected k according to desired precisioneValue, the preset constant k according to regular importanceerrValue optimizes,
Specific step is as follows:
(1) algorithm initialization carries out initialization X={ X to population in the solution space of problem1,X2,…,XNp, individual Xi(k)
=[xi,1(k),xi,2(k),…,xi,D(k)], i=1,2 ..., NpThe solution of characterization problems, D are the dimension of solution space, and k expression changes
Generation number, NpFor population invariable number.Initialization population should cover entire search space as much as possible, and defined minimum and
Individual is unanimously randomized in the search space of maximum bound of parameter limitation, setting minimum and maximum boundary is respectively Xmin
={ xmin,1,…,xmin,DAnd Xmax={ xmax,1,…,xmax,D}.In the number of iterations k=1, j-th of component of i-th of individual
Initial value generate as the following formula:
xi,j(0)=xmin,j+rand(0,1)·(xmax,j-xmin,j) (11)
Wherein rand (0,1) is equally distributed random number between 0 to 1.
(2) mutation operation.After being initialized to population, for each individual Xi(k) a variation vector can correspondingly be generated
Vi(k)=[vi,1(k),vi,2(k),…,vi,D(k)], individual Xi(k) it is also referred to as object vector, the expression formula of mutation operation is as follows:
Vi(k)=Xbest(k)+F·(Xr1(k)-Xr2(k)),1≤r1≠r2≠i≤Np (12)
Wherein, XbestIt (k) is the individual vector in population in current kth time iteration with optimal adaptation angle value, scale factor F
It is the positive control parameter for scaling difference vector.
(3) crossover operation.For each pair of object vector Xi(k) and its accordingly make a variation vector Vi(k) crossover operation is carried out, with life
At trial vector Ui(k)=[ui,1(k),ui,2(k),…,ui,D(k)].The binomial interleaved scheme is carried out according to the following formula:
Wherein, CR is a customized crossing-over rate, and for value range usually between 0 to 1, it controls the diversity of population,
And algorithm is avoided to fall into local optimum.If the value of certain parameters of newly-generated trial vector has been more than corresponding bound,
Then it is reinitialized uniformly at random in predefined scope.Then, the fitness value of all trial vectors is assessed.
(4) selection operation.The fitness of trial vector is compared by the operation with corresponding object vector, is selected wherein more preferable
Solution.For minimization problem, the expression formula of selection operation is as follows:
Wherein, f (x) is fitness value.In the evolutionary process of every generation, each individual vector is primary as target individual, calculates
Method retains defect individual, eliminates worst individual, guiding search process is approached to globally optimal solution by constantly iterating to calculate.
5. quick adaptive radar marine target detection system according to claim 1, which is characterized in that the target inspection
Module is surveyed, to carry out target detection, is completed using following process:
(1) D sea clutter echo-signal amplitude is acquired in sampling instant t obtain TX=[xt-D+1,…,xt], xt-D+1Indicate t-D
The sea clutter echo-signal amplitude of+1 sampling instant, xtIndicate the sea clutter echo-signal amplitude of t sampling instant;
(2) it is normalized:
(3) the sea clutter predicted value that sampling instant (t+1) is calculated in Fuzzy Wavelet Network modeling module is substituted into.
(4) difference e for calculating sea clutter predicted value and radar return measured value calculates control limit Qα:
Wherein, α is confidence level, θ1,θ2,θ3,h0It is intermediate variable,Indicate the i power of j-th of characteristic value of covariance matrix, k
It is sample dimension, CαIt is the statistics that normal distribution confidence level is α;
(5) it carries out detection judgement: working as e2Difference is greater than control limit QαWhen, there are targets for the point, otherwise without target.
The adaptive radar marine target detection system of the agility, the host computer further include: model modification module, to press
The sampling time interval of setting acquires data, by obtained measured data compared with model prediction value, if relative error is greater than
10%, then training sample data are added in new data, update forecasting model.And result display module, to by target detection
The testing result of module is shown in host computer.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810689688.7A CN108983178A (en) | 2018-06-28 | 2018-06-28 | A kind of Intelligent radar sea target detection system that agility is adaptive |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810689688.7A CN108983178A (en) | 2018-06-28 | 2018-06-28 | A kind of Intelligent radar sea target detection system that agility is adaptive |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108983178A true CN108983178A (en) | 2018-12-11 |
Family
ID=64539489
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810689688.7A Pending CN108983178A (en) | 2018-06-28 | 2018-06-28 | A kind of Intelligent radar sea target detection system that agility is adaptive |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108983178A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110736980A (en) * | 2019-11-04 | 2020-01-31 | 广东博智林机器人有限公司 | method, device and system for monitoring climbing rack state information |
CN110852440A (en) * | 2019-12-02 | 2020-02-28 | 中国人民解放军国防科技大学 | Ocean front detection method based on dynamic fuzzy neural network |
CN112364678A (en) * | 2019-09-30 | 2021-02-12 | 山东省科学院海洋仪器仪表研究所 | Buoy identification and positioning method based on NPU board card and shipborne device thereof |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7372393B2 (en) * | 2006-07-07 | 2008-05-13 | Mitsubishi Electric Research Laboratories, Inc. | Method and system for determining unwrapped phases from noisy two-dimensional wrapped-phase images |
CN101718870A (en) * | 2009-11-13 | 2010-06-02 | 西安电子科技大学 | High-speed weak target flight path detection method of image field |
CN102147463A (en) * | 2011-03-03 | 2011-08-10 | 浙江大学 | System and method for forecasting Qunzhi radar sea clutters |
CN102147465A (en) * | 2011-03-03 | 2011-08-10 | 浙江大学 | System and method for detecting sea target by chaos optimizing radar |
CN102183749A (en) * | 2011-03-03 | 2011-09-14 | 浙江大学 | Sea target detecting system of adaptive radar and method thereof |
-
2018
- 2018-06-28 CN CN201810689688.7A patent/CN108983178A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7372393B2 (en) * | 2006-07-07 | 2008-05-13 | Mitsubishi Electric Research Laboratories, Inc. | Method and system for determining unwrapped phases from noisy two-dimensional wrapped-phase images |
CN101718870A (en) * | 2009-11-13 | 2010-06-02 | 西安电子科技大学 | High-speed weak target flight path detection method of image field |
CN102147463A (en) * | 2011-03-03 | 2011-08-10 | 浙江大学 | System and method for forecasting Qunzhi radar sea clutters |
CN102147465A (en) * | 2011-03-03 | 2011-08-10 | 浙江大学 | System and method for detecting sea target by chaos optimizing radar |
CN102183749A (en) * | 2011-03-03 | 2011-09-14 | 浙江大学 | Sea target detecting system of adaptive radar and method thereof |
Non-Patent Citations (2)
Title |
---|
吴燕玲: "遗传规划及其在数据驱动软测量建模中的应用", 《中国博士学位论文全文数据库 信息科技辑》 * |
袁国平: "航天器姿态系统的自适应鲁棒控制", 《中国博士学位论文全文数据库 工程科技I辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112364678A (en) * | 2019-09-30 | 2021-02-12 | 山东省科学院海洋仪器仪表研究所 | Buoy identification and positioning method based on NPU board card and shipborne device thereof |
CN112364678B (en) * | 2019-09-30 | 2023-04-07 | 山东省科学院海洋仪器仪表研究所 | Buoy identification and positioning method based on NPU board card and shipborne device thereof |
CN110736980A (en) * | 2019-11-04 | 2020-01-31 | 广东博智林机器人有限公司 | method, device and system for monitoring climbing rack state information |
CN110736980B (en) * | 2019-11-04 | 2021-09-21 | 广东博智林机器人有限公司 | Method, device and system for monitoring state information of climbing frame |
CN110852440A (en) * | 2019-12-02 | 2020-02-28 | 中国人民解放军国防科技大学 | Ocean front detection method based on dynamic fuzzy neural network |
CN110852440B (en) * | 2019-12-02 | 2022-07-26 | 中国人民解放军国防科技大学 | Ocean front prediction method based on dynamic fuzzy neural network |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102183749B (en) | Sea target detecting system of adaptive radar and method thereof | |
CN102147465B (en) | System and method for detecting sea target by chaos optimizing radar | |
Jain et al. | Neural networks in ocean engineering | |
CN108230302B (en) | Detection and disposal method for marine organism invading from cold source sea area of nuclear power plant | |
Cheng et al. | Modeling and analysis of motion data from dynamically positioned vessels for sea state estimation | |
CN108983178A (en) | A kind of Intelligent radar sea target detection system that agility is adaptive | |
CN108983177A (en) | A kind of radar sea clutter forecast system that parameter is optimal and method | |
CN107942312A (en) | A kind of Intelligent radar sea target detection system and method based on differential evolution invasive weed optimization algorithm | |
CN108983179A (en) | A kind of radar marine target detection system of colony intelligence agility | |
CN107656250A (en) | A kind of Intelligent radar sea target detection system and method based on artificial bee colony algorithm | |
CN107703491A (en) | Sea clutter optimal soft survey instrument and method based on improved drosophila optimized algorithm optimization RBF neural | |
CN101893711B (en) | Ocean distributed target detection method based on collective evolution of agents | |
CN108983181A (en) | A kind of radar marine target detection system of gunz optimizing | |
CN116502478B (en) | Submarine topography monitoring-based pile-off auxiliary decision-making method for self-lifting platform | |
CN102183751B (en) | Intelligent radar sea target detection system and method | |
CN107942300A (en) | A kind of Intelligent radar sea target detection system and method based on improvement artificial bee colony algorithm | |
CN108983185A (en) | A kind of radar marine target detection system and method for intelligent adaptive | |
CN102183754B (en) | System and method for detecting sea target by using robust intelligent radar | |
CN102183746A (en) | Radar marine target detection system and method | |
CN108983183A (en) | A kind of adaptive radar sea clutter forecast system | |
CN107942303A (en) | A kind of Intelligent radar sea clutter forecast system and method based on improvement artificial bee colony algorithm | |
Kerbaa et al. | CNN-LSTM based approach for parameter estimation of k-clutter plus noise | |
CN108983184A (en) | A kind of adaptive high-precision Intelligent radar sea target detection system | |
CN107976662A (en) | A kind of Intelligent radar sea target detection system and method based on invasive weed optimization algorithm | |
CN102183744B (en) | Swarm-intelligence radar sea target detecting system and method |
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: 20181211 |
|
RJ01 | Rejection of invention patent application after publication |