CN108154233A - The complex network community detection method of global convergence performance can be improved - Google Patents
The complex network community detection method of global convergence performance can be improved Download PDFInfo
- Publication number
- CN108154233A CN108154233A CN201810035790.5A CN201810035790A CN108154233A CN 108154233 A CN108154233 A CN 108154233A CN 201810035790 A CN201810035790 A CN 201810035790A CN 108154233 A CN108154233 A CN 108154233A
- Authority
- CN
- China
- Prior art keywords
- individual
- population
- pop
- community
- value
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 26
- 238000005457 optimization Methods 0.000 claims abstract description 33
- 230000035772 mutation Effects 0.000 claims abstract description 22
- 230000003044 adaptive effect Effects 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims description 34
- 238000012360 testing method Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 16
- 230000006978 adaptation Effects 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 10
- 230000002708 enhancing effect Effects 0.000 claims description 8
- 238000002474 experimental method Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 7
- 239000004615 ingredient Substances 0.000 claims description 5
- 230000001149 cognitive effect Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 108090000623 proteins and genes Proteins 0.000 claims description 3
- 238000005728 strengthening Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 230000006870 function Effects 0.000 description 14
- 230000010429 evolutionary process Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 6
- 244000257790 Brassica carinata Species 0.000 description 4
- 235000005156 Brassica carinata Nutrition 0.000 description 4
- 241001481833 Coryphaena hippurus Species 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000002028 premature Effects 0.000 description 3
- 201000004569 Blindness Diseases 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000019771 cognition Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004540 process dynamic Methods 0.000 description 2
- 230000001737 promoting effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 240000000233 Melia azedarach Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013442 quality metrics Methods 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/01—Social networking
Abstract
The invention discloses a kind of complex network community detection methods that can improve global convergence performance, specifically include：The step of improving the global convergence performance of differential evolution algorithm, which is specifically：(1) classification adaptive differential class Mutation Strategy；(2) dynamic self-adapting parameter adjustment；(3) the carry out difference selection operation based on historical information.TSP question based on classification will act on all individuals in every generation population to be terminated, therefore the variation of each individual can be accessed and targetedly be adjusted until evolving.On the one hand, the exploring ability of excellent individual can be strengthened, to increase the possibility that global optimum is found in its neighborhood；On the other hand, the producing capacity of poor individual can be strengthened, to accelerate its search speed to global optimization.
Description
Technical field
The present invention relates to a kind of community detection method, specifically a kind of complex network that can improve global convergence performance
Community detection method.
Background technology
It is proposed in succession there are many community detection method in past few years, wherein most widely used is based on modularity
Optimal method.However, modularity optimization is substantially a typical np hard problem, traditional deterministic method,
Such as mathematical programming approach, greedy algorithm, spectrum analysis method and extremal optimization algorithm, it will usually there is Premature Convergence or convergence to stagnate now
As.In addition, with the enhancing of real world network size and structural fuzzy, the extreme value degenerate problem during optimizing becomes
More serious, this is meant that in numerous locally optimal solutions with exponential increase, and finding community of global optimum and dividing becomes more
Add difficulty, therefore be severely impacted the Stability and veracity of community structure obtained by detection.
Modularity optimization method based on EA is due to powerful global optimization's ability, in a variety of test problems upper tables
Reveal notable superiority.In addition, it is contemplated that the acquisition of the prior information in real world network is more difficult, such algorithm is not required to
Want any prior information (such as community's number) and specific mathematical model.However, although the modularity optimization method based on EA is more
Achieved on kind of Web Community test problems it is satisfactory as a result, but Premature Convergence and extreme value the problem of degenerating do not obtain
Adequately solve.
Quality is divided in order to overcome the above problem and improve optimal community, the modularity optimization based on EA should be further improved
Convergence energy.Previous experiments the result shows that, the constringency performance of the modularity optimization algorithm based on EA depends primarily on two
A key factor, primary factor are also most important factor is how to improve the global convergence abilities of EA in itself, and another factor is
How to efficiently use network topological information and reduce search space huge in modularity optimization process.However as far as we know, it is existing
Have in algorithm and basic EAs is usually ignored into its convergence capabilities directly as optimisation strategy, so as to cause the Premature Convergence of EAs,
It is also poor that the optimal community obtained divides quality.At the same time, although existing some algorithm carries out the evolutional operation in EAs
It improves, community's detection demand is met, but the inappropriate use of topology information destroys the overall situation most by converged network topology information
The search space that excellent community divides.
Invention content
In view of the deficiencies of the prior art, the present invention proposes the complex network community detection that can improve global convergence performance
Method has redesigned three main evolutional operations, including the TSP question strategy based on classification, dynamic self-adapting parameter
Adjustable strategies and the selection operation based on historical information.On the other hand, preferably to utilize network topological information, it is proposed that a kind of
Improved community's adjustable strategies based on neighborhood information are drawn using ensureing while DE search spaces are reduced as community of global optimum
Divide and enough search spaces are provided.Finally, the new modularity optimization algorithm CDEMO based on DE is proposed.
To achieve the above object, the present invention provides a kind of complex network community detections that can improve global convergence performance
Method specifically includes：The step of improving the global convergence performance of differential evolution algorithm, which is specifically：
(1) classification adaptive differential class Mutation Strategy；
(2) dynamic self-adapting parameter adjustment；
(3) the carry out difference selection operation based on historical information.
Further, classification adaptive differential class Mutation Strategy, concrete operations are as follows：
For each target individual Xi,tIf its ideal adaptation angle value fiMore than current entire population at individual fitness
The average of value, then be classified as excellent individual, and globally optimal solution is closer in the position of search space；Therefore, in Xi,t
In good gene be reserved for strengthening local search around individual, the corresponding vector V that makes a variationi,tGenerating mode is as follows：
Vi,t=Fi,t.Xpbesti,t+Wi,t.(Xr2,t-Xr3,t) (1)
Wherein, Xpbesti,tRepresent individual Xi,tIn the history optimal solution in preceding t generations, for enhancing individual exploring ability；Xr2,tWith
Xr3,tIt is randomly selected two Different Individuals from population, and meets condition r2 ≠ r3 ≠ i；Fi,tAnd Wi,tIt is XiControl
Parameter, numerical value is according to evolutionary generation and Xi,tIdeal adaptation angle value dynamic adjust；
For each target individual Xi,tIf its ideal adaptation angle value fiLess than current entire population at individual fitness
The average of value is then classified as poor individual, in the position of search space and globally optimal solution farther out；Therefore, strengthen it
The exchanging to promote global search between excellent individual in population, the corresponding vector V that makes a variationi,tGenerating mode is as follows：
Vi,t=Wi,t.Xr1,t+Ki,t.(Xgbest,t-Xi,t) (2)
Wherein Xr1,tIt is the randomly selected individual from population, and meets condition r1 ≠ i；Xgbest,tRepresent current iteration kind
Optimal solution in group, for enhancing Xi,tExploring ability；Wi,tAnd Ki,tIt is XiControl parameter, numerical value is according to evolutionary generation
And Xi,tIdeal adaptation angle value carry out dynamic adjustment.
Further, dynamic self-adapting parameter adjustment：Three control parameters W, K, F, it is respectively random in mutation process
Ingredient, social ingredient and cognitive component；In addition, also there are one crucial control parameter CR in crossover operation, it is each for determining
Test individual ui,tIn from variation individual Vi,tThe percentage of middle succession；Adjustment process is specific as follows：
1. parameter adaptive adjustment is carried out according to the fitness value of individual：To poor individual, strengthen variation and the journey intersected
Degree, to introduce more directivity informations during evolution.Therefore, the random element in mutation process, social ingredient with
And the succession in crossover process all enhances, the CR values in W and K and intersection in corresponding formula (2) are larger；On the contrary, for
For excellent individual, strengthen the cognition part in mutation process, parameter adjustment should defer to opposite principle, corresponding to formula (1)
In larger F values and smaller W values.
2. according to evolution iterations dynamic self-adapting：In early stage of evolving, strengthen the exploring ability of individual, with
Ensure fully to be searched in each individual neighborhood.On the contrary, in the later stage of evolution stage, strengthen the producing capacity of individual, strengthen
The convergence of entire group is accelerated in exchange between individual.According to this principle, F in evolutionary process, the value of W, CR is gradually reduced,
And K values gradually increase.
Based on mentioned above principle, parameter value can be adaptively adjusted, and each individual can obtain in evolutionary process
Dynamic control.Specific operation process is as follows：
Further, the step of carry out difference selection operation based on historical information, specifically：
The outstanding solution generated in entire evolutionary process will be saved in historical information, and for follow-up evolutional operation.
To realize this target, special population pbest_pop is introduced, by the history optimal solution X of individual each in populationpbesti,tForm kind
Group pbest_pop, and generated in initial phase, it is updated after each evolutional operation；To individual X each in populationi,t,
If its fitness value is improved during a certain evolutional operation, then newly-generated individual will be used as Xi,tCurrent go through
History optimal solution, and be saved in pbest_pop；It is all individual by substitute species in pbest_pop after every generation evolutional operation
All individuals in group pop, and current optimal solution X is selected from pbest_popgbest,t。
Further, the step of this method is further included using improved progress community's amendment based on neighborhood information, specifically：
If a node meets community's correction conditions, then the node will likely be placed in again in its all affiliated community of neighborhood node,
And the probability of merging is directly proportional to the scale of neighborhood community.
Further, this method further includes the modularity optimization algorithm based on classification differential evolution algorithm, specifically：
S1：Initialization of population；
S1.1 sets network parameter, including number of nodes n, adjacency matrix adj, community correction threshold δ；DE algorithms ginseng is set
Number, including individual dimension D, Population Size NP, population iterations t and maximum iteration tmax；
S1.2 is with the individual representation random initializtion population pop of community's label；
S2：It identifies and records optimal solution
S2.1 is identified and is recorded t for the optimum individual X in population popgbest,t；
S2.2 is identified and is recorded t for individual X each in population popi,tHistory optimal solution Xpbesti,t；By all populations
The X of individualpbesti,tBuild initial population pbest_pop；
S3：When population iterations are less than population maximum iteration, population iterations are unsatisfactory for condition from adding one
Then terminate the cycle of S3.1-S3.5；
S3.1 passes through adaptive classification differential variation construction of strategy variation population mutation_pop；
When the value of i is arrived for 1 in Population Size numberical range, step a) is to cycle e) for progress, if the value of i is not arrived 1
In Population Size numberical range, then step a) is jumped out to e), end loop；
A) 3 different individual X are randomly selected from population popr1,t, Xr2,t, Xr3,t；
B) dynamic adjustment Mutation parameter Fi,t、wi,t、Ki,t；
C) according to fitness value Q to Xi,tClassify；
D) according to adaptive classification differential variation strategy generating variation individual Vi,t；
E) V is calculatedi,tModule angle value and and Xi,tIndividual is made comparisons, and more excellent individual is stored in pbest_pop；
If i is more than NP, leapfrog goes out a) suddenly to cycle e)；
S3.2 is based on neighborhood information and carries out community's amendment；
S3.3 is according to variation population mutation_pop and population pop structure cross-species crossover_pop；
When the value of i is arrived for 1 in Population Size numberical range, step a) is to cycle d) for progress, if the value of i is not arrived 1
In Population Size numberical range, then step a) is jumped out to d), end loop；
A) i-th of individual u in cross-species is initializedi,t=xi,t；
B) dynamic adjustment cross parameter CRi,t；
C) by from variation individual Vi,tIt inherits community information and carrys out Adjustment Tests individual ui,t；
D) u is calculatedi,tModule angle value and be compared with i-th of individual in pbest_pop, retain compared with the figure of merit extremely
pbest_pop；
S3.4 is based on neighborhood information and carries out community's amendment；
S3.5 is by replacing all individual update pop in pbest_pop；
S4：Export the X in popgbest,tIt is divided as final optimal community, otherwise returns to S3 steps.
The present invention due to using the technology described above, can obtain following technique effect：
TSP question based on classification will act on all individuals in every generation population to be terminated, therefore each until evolving
The variation of individual, which can access, targetedly to be adjusted.On the one hand, the exploring ability of excellent individual can be strengthened, to increase
Its neighborhood finds the possibility of global optimum；On the other hand, the producing capacity of poor individual can be strengthened, to accelerate it to the overall situation
The search speed of optimization.In short, the evolution demand with different fitness characteristic individuals, can be obtained by new Mutation Strategy
To better meeting.Under the guiding of directivity information, the blindness in search process can efficiently reduce, and offspring individual
It can also be improved with the quality of optimal solution.And the degree of variation of each individual of dynamic self-adapting during evolution.
Historical information will be saved as, and for follow-up evolutional operation by also achieving the outstanding solution generated in entire evolutionary process.
New correction strategy can effectively reduce search space, additionally it is possible to relax limitation when community is corrected, be global optimum
Solution provides sufficient search space, so as to preferably using topology information known to network, and promote the convergence of CDEMO algorithms.
CDEMO algorithms can efficiently identify the community structure of complex network, improve accuracy that optimal community divides, steady
Qualitative and scalability has the complex network of very fuzzy community structure including those.
Description of the drawings
Fig. 1 is the adaptive differential evolution algorithm flow chart based on classification；
Fig. 2 is the modularity optimization algorithm CDEMO flow charts based on differential evolution；
Fig. 3 is the CDEMO that the different zout of GN networks is worth to and the average NMI values figure of other algorithms；
Fig. 4 is the CDEMO that the different μ of LFR networks is worth to and the average NMI values figure of other algorithms；
Fig. 5 is that community structure of the CDEMO algorithms on Karate networks divides identification figure；
Fig. 6 is that community structure of the CDEMO algorithms on Dolphin networks divides identification figure；
Fig. 7 is that community structure of the CDEMO algorithms on Polbooks networks divides identification figure；
Fig. 8 is that community structure of the CDEMO algorithms on Football networks divides identification figure.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, it is right in the following with reference to the drawings and specific embodiments
The present invention is described in detail.
Embodiment 1
A kind of network community detection method that can improve global convergence performance is present embodiments provided, is specifically included：
First, the global convergence performance for raising DE algorithms, has redesigned three main evolutional operations：
(1) classification adaptive differential Mutation Strategy
Corrective measure mainly includes following several respects：
1. utilize current population optimal solution Xgbest,tWith the history optimal solution X of each individualpbesti,tChange randomly selected
Body guiding variation direction；
2. propose and balanced using a kind of new adaptive classification mechanism the spy of the individual with different compliance characteristics
Rope and producing capacity；
3. the degree of variation of each individual carries out dynamic self-adapting by parameter in evolutionary process.
New Mutation Strategy concrete operations are described as follows：
For each target individual Xi,tIf its ideal adaptation angle value fiMore than current entire population at individual fitness
The average of value, then be classified as excellent individual, and globally optimal solution is closer in the position of search space.Therefore, in Xi,t
In good gene should be reserved for strengthening the local search around individual, the corresponding vector V that makes a variationi,tGenerating mode is as follows：
Vi,t=Fi,t.Xpbesti,t+Wi,t.(Xr2,t-Xr3,t) (1)
Wherein Xpbesti,tRepresent individual Xi,tIn the history optimal solution in preceding t generations, for enhancing individual exploring ability.Xr2,tWith
Xr3,tIt is randomly selected two Different Individuals from population, and meets condition r2 ≠ r3 ≠ i.Fi,tAnd Wi,tIt is XiControl
Parameter, numerical value is according to evolutionary generation and Xi,tIdeal adaptation angle value dynamic adjust.
For each target individual Xi,tIf its ideal adaptation angle value fiLess than current entire population at individual fitness
The average of value is then classified as poor individual, in the position of search space and globally optimal solution farther out.Therefore, should strengthen
Its exchanging to promote global search between excellent individual in population, the corresponding vector V that makes a variationi,tGenerating mode is as follows：
Vi,t=Wi,t.Xr1,t+Ki,t.(Xgbest,t-Xi,t) (2)
Wherein Xr1,tIt is the randomly selected individual from population, and meets condition r1 ≠ i.Xgbest,tRepresent current iteration kind
Optimal solution in group, for enhancing Xi,tExploring ability.Wi,tAnd Ki,tIt is XiControl parameter, numerical value is according to evolutionary generation
With this Xi,tIdeal adaptation angle value dynamic adjust.
The advantageous effect of generation：It is straight that the above-mentioned TSP question based on classification will act on all individuals in every generation population
Terminate to evolution, therefore the variation of each individual can be accessed and targetedly be adjusted.On the one hand, excellent individual can be strengthened
Exploring ability, with increase its neighborhood find global optimum possibility；On the other hand, the exploitation of poor individual can be strengthened
Ability, to accelerate its search speed to global optimization.In short, the evolution demand with different fitness characteristic individuals, it can
To be better met by new Mutation Strategy.Under the guiding of directivity information, blindness in search process can be with
It efficiently reduces, and the quality of offspring individual and optimal solution can also be improved.
(2) dynamic self-adapting parameter adjustment
Three control parameters W, K, F correspond respectively to random element, social ingredient and cognitive component in mutation process.
In addition, also there are one crucial control parameter CR in crossover operation, for determining each experiment individual ui,tIn from variation individual
Vi,tThe percentage of middle succession.
1. parameter adaptive adjustment is carried out according to the fitness value of individual.To poor individual, it should strengthen variation and intersect
Degree, to introduce more directivity informations during evolution.Therefore, the random element in mutation process, society into
Point and crossover process in succession should all enhance, W and K in corresponding formula (2) and the CR values in intersecting are larger.
On the contrary, for excellent individual, it should strengthen the cognition part in mutation process, parameter adjustment should defer to opposite principle,
Corresponding to larger F values and smaller W values in formula (1).
2. according to evolution iterations dynamic self-adapting.In early stage of evolving, it should strengthen the exploration energy of individual
Power, to ensure fully to be searched in each individual neighborhood.On the contrary, in the later stage of evolution stage, it should strengthen the exploitation of individual
Ability strengthens the exchange between individual, accelerates the convergence of entire group.According to this principle, F in evolutionary process, the value of W, CR
It is gradually reduced, and K values gradually increase.
Based on mentioned above principle, parameter value can be adaptively adjusted, and each individual can obtain in evolutionary process
Dynamic control.Specific operation process is as follows：
The advantageous effect of generation：The degree of variation of each individual of dynamic self-adapting during evolution.
(3) the difference selection operation based on historical information
The outstanding solution generated in entire evolutionary process will be saved as historical information, and for follow-up evolutional operation.For
It realizes this target, introduces special population pbest_pop, by the history optimal solution X of individual each in populationpbesti,tForm population
Pbest_pop, and generated, and be updated after each evolutional operation in initial phase.To individual each in population
Xi,tIf its fitness value is improved during a certain evolutional operation, then newly-generated individual will be used as Xi,tWork as
Preceding history optimal solution, and be saved in pbest_pop.After every generation evolutional operation, all individuals will replace in pbest_pop
For all individuals in population pop, and current optimal solution X is selected from pbest_popgbest,t。
The advantageous effect of generation：Realize that the outstanding solution generated in entire evolutionary process will be saved as historical information, and
For follow-up evolutional operation.
Improved differential evolution Algorithm Convergence is tested：
Above-mentioned three corrective measures are provided to improve the global convergence of DE algorithms, algorithm flow chart such as Fig. 1 after improvement
It is shown.
Different from standard DE algorithms, more directivity informations are combined in new mutation operation, therefore individual can more have
Pointedly into row variation.Performed in addition, selection operation is not placed on after crossover operation, but by each evolutional operation it
Population Regeneration pbest_pop is selected and is retained outstanding solution afterwards.
In order to verify the above-mentioned corrective measure for DE, carried out using 18 standard Benchmark function pairs innovatory algorithms
Test, wherein f1-f5 is single mode state function, and f6-f14 is basic multi-modal function, and f15-f16 is spread function, and f17-f18
It is composite function.Table 1 provides the details of standard Benchmark functions.
Improved DE algorithms and 4 efficient and widely used DE algorithm patterns carry out performance comparison, including DE/
Rand/2/dir, DE/rand/1/bin, DE/current-to-best/2/bin and DE/best/1/bin.For ease of comparing,
DE_version1 is named as using the algorithm of new Mutation Strategy and parameter adaptive adjustable strategies, and uses all three improvement
The DE algorithms of measure are named as DE_version2.
During the experiment, all algorithms use same initial population scale NP=100 on each test problem, together
The variable dimension D=30 of sample and same stop criterion Max_FEs=5.0e+0.5.In addition, in the DE algorithms of all patterns
Parameter F and CR be adjusted all in accordance with the adaptive mode shown in formula (5) and (6).The value range of relevant parameter is W ∈
[0.1,0.9], K ∈ [0.3,0.9], F ∈ [0.3,0.9], CR ∈ [0.1,0.9].
The DE algorithms of six kinds of patterns carry out performance comparison in terms of optimal solution accuracy and robustness.Experimental result such as table 2
It is shown, acquire the average value of optimal solution and standard deviation including 30 independent operatings on each test function (in bracket).Each is surveyed
Optimal solution on trial function is shown with runic.It will be seen that DE_version1 and DE_version2 almost exist from table 2
It is better than other 4 kinds of algorithms on all test functions.DE_version2 is successfully converged to really in 50.0% test function
Globally optimal solution, and show on 88.9% test function optimal.The above results prove, the TSP question based on classification
Strategy can effectively improve the accuracy of offspring individual quality and optimal solution.In addition, compared with DE_version1, DE_
Version2 is significantly improved in terms of accuracy, illustrates that the new selection operation based on historical information can effectively improve DE algorithms
Global convergence ability.
Above-mentioned experimental result illustrates that improved method proposed in this paper is successfully effective, can effectively improve primary standard DE
The global convergence performance of algorithm, the modularity optimization problem in being detected for complex network community provide a kind of effective overall situation most
Optimization method.
2nd, preferably to utilize network topological information, it is proposed that a kind of improved community's adjustment plan based on neighborhood information
Slightly, enough search spaces are provided to ensure to divide as community of global optimum while DE search spaces are reduced.
For DE algorithms after improvement, preferably to utilize network topological information, it is proposed that a kind of improved based on neighborhood letter
Community's correction strategy of breath：
In order to avoid this inappropriate use to topology information, a kind of new community's correction strategy is proposed in CDEMO.
If a node meets community's correction conditions, then the node will likely be placed in again in its all affiliated community of neighborhood node,
And the probability of merging is directly proportional to the scale of neighborhood community.
The advantageous effect of generation：It can effectively reduce search space new correction strategy and former strategy, and heavier be
Limitation when community is corrected can be relaxed, sufficient search space is provided for globally optimal solution, so as to preferably utilize network
Know topology information, and promote the convergence of CDEMO algorithms.
3rd, the new modularity optimization algorithm CDEMO based on DE
(1) CDEMO algorithms, algorithm flow chart are as shown in Figure 2：
1：Initialization of population；
1.1 setting network parameters, including number of nodes n, adjacency matrix adj, community correction threshold δ.DE algorithm parameters are set,
Including individual dimension D, Population Size NP, population iterations t and maximum iteration tmax；
1.2 with the individual representation random initializtion population pop of community's label；
2：It identifies and records optimal solution
2.1 identify and record t for the optimum individual X in population popgbest,t；
2.2 identify and record t for individual X each in population popi,tHistory optimal solution Xpbesti,t.By all populations
The X of bodypbesti,tBuild initial population pbest_pop；
3：When population iterations are less than population maximum iteration, population iterations are unsatisfactory for condition from adding one
Then terminate the cycle of S3.1-S3.5；
3.1 pass through adaptive classification differential variation construction of strategy variation population mutation_pop；
When the value of i is arrived for 1 in Population Size numberical range, step a) is to cycle e) for progress, if the value of i is not arrived 1
In Population Size numberical range, then step a) is jumped out to cycle e), end loop；
A) 3 different individual X are randomly selected from population popr1,t, Xr2,t, Xr3,t；
B) dynamic adjustment Mutation parameter Fi,t、wi,t、Ki,t；
C) according to fitness value Q to Xi,tClassify；
D) according to adaptive classification differential variation strategy generating variation individual Vi,t；
E) V is calculatedi,tModule angle value and and Xi,tIndividual is made comparisons, and more excellent individual is stored in pbest_pop；
If i is more than NP, leapfrog goes out a) suddenly to cycle e)；
3.2 carry out community's amendment based on neighborhood information；
3.3 according to variation population mutation_pop and population pop structure cross-species crossover_pop；
When the value of i is arrived for 1 in Population Size numberical range, step a) is to cycle d) for progress, if the value of i is not arrived 1
In Population Size numberical range, then step a) is jumped out to cycle d), end loop；
A) i-th of individual u in cross-species is initializedi,t=xi,t；
B) dynamic adjustment cross parameter CRi,t；
C) by from variation individual Vi,tIt inherits community information and carrys out Adjustment Tests individual ui,t；
D) u is calculatedi,tModule angle value and be compared with i-th of individual in pbest_pop, retain compared with the figure of merit extremely
pbest_pop；
3.4 carry out community's amendment based on neighborhood information；
3.5 by replacing all individual update pop in pbest_pop.
4：Stop algorithm if stopping criterion is met, export the X in popgbest,tIt is divided as final optimal community,
Otherwise the 3rd step is returned.
The advantageous effect of generation：CDEMO algorithms can efficiently identify the community structure of complex network, improve optimal community
Accuracy, stability and the scalability of division have the complex network of very fuzzy community structure including those.
CDEMO algorithm performances are tested：
Strategy validity will be changed by experimental verification new communities, and verify whether DE algorithmic statements performance is promoted advantageous
In its application in modularity optimization.6 kinds of modularity optimization algorithms based on DE are built, are named as DEMO1-6.These algorithms
Using the optimisation strategy of different DE algorithms (there is different experiment individual generation strategies) modularity as an optimization.DEMO1-4
It is middle to apply different DE algorithms respectively, including DE/rand/2/dir, DE/rand/1/bin, DE/current-to-best/2/
Bin, DE/best/1/bin.DEMO5 employs a kind of widely used random variation strategy, i.e. node community ownership is with complete
Random manner is adjusted.DEMO6 is tactful as an optimization by improved DE_version2.It is combined on the basis of DEMO6
The improvement of preceding proposition is to improve the global convergence of algorithm and reduce algorithm search space on this basis and ensure global optimum
CDEMO is constructed in new communities' modification operation of solution search space.
All algorithms are tested on 4 real world social networks, as shown in table 3, including karate club net
Network, dolphin network, American politics books network and American university rugby network.Experimental result is as shown in table 4, including each calculation
The average value and standard deviation of method gained modularity Q values after 30 independent operatings on each test network.
From table 4 we it can be clearly seen that due to different convergences, the DE algorithms of different mode optimize in modularity
Performance has larger difference in problem.Compared with DEMO3 and DEMO4, in DEMO1-2 and DEMO5 Mutation Strategies it is random into
Divide and make it have stronger exploring ability, therefore better Qavg and Qstd values can be obtained.In addition, divided in optimal community
In terms of precision and stability, DEMO6 performances are better than DEMO1-5, it was demonstrated that DE convergences, which can be promoted, contributes to it in mould
Application in lumpiness optimization.Compared with DEMO6, CDEMO is obtained on Karate networks, Dolphin networks and PolBooks networks
Better Qavg and Qstd values are obtained, the accuracy of detected community is further promoted.
Based on above-mentioned test result, we may safely draw the conclusion, believes from promoting DE algorithm global convergence abilities and promoting topology
Two aspect enhancing algorithmic statement performance of service efficiency is ceased, contributes positively to improve what optimal community in modularity optimization problem divided
Quality.
4th, community's detection performance is tested
1. experimental setup
Performance Evaluation is carried out to CDEMO algorithms on artificial synthesized network and real world social networks.CDEMO algorithms exist
7.0 software programmings of MATLAB are realized, and in 7 systems of Windows for using Pentium Dual Core 2.5GHz processors and 2.0GB memories
On tested.Parameter setting in CDEMO is as follows：Population scale NP values 100, maximum iteration tmax values 200, control
The value range of parameter processed is set as, W ∈ [0.1,0.9], K ∈ [0.3,0.9], F ∈ [0.3,0.9], CR ∈ [0.1,0.9].
2. Performance evaluation criterion
(1) modularity Q：For the real world network of unknown community structure, usually refer to by the use of modularity function as performance
Mark weighs the significance degree of detection gained community structure.Modularity is defined as follows：
Wherein, M is the total number of edges of network；A=(aij) n*n is network adjacent matrix；Ki and kj represents node i and j's respectively
Degree；δ (i, j) represents community's attaching relation of node i and node j, if it is 1 that the two, which belongs to same community value, otherwise value
It is 0.It being represented in network there are community structure when Q values are more than 0, more than 0.3 when represents that the community structure of network is more apparent,
Q values are bigger to illustrate that community structure is more notable.It is that current use is most extensive although there are resolution ratio restricted problems for modularity
Community divide quality metric.
(2) normalized mutual information NMI：For the artificial synthesized network of known community structure, usually refer to by the use of NMI as performance
Mark weighs detection gained community and divides the approximation ratio divided with community content, and calculation formula is such as shown in (8).Assuming that A is network
Community content divide, B be detection gained community divide, define hybrid matrix C, wherein row represent A in community divide, list
Show that the community detected in B divides.Elements C ij represents the section identical with dividing j-th of community in B of i-th of community in division A
It counts out.According to the definition of C, evaluation criterion NMI is defined as follows：
Wherein, N represents the interstitial content in network；CA and CB represents to divide community's number in A and B respectively；Ci is mixed
The sum of i-th row element in Matrix C of confusing represents and divides i-th of community's interstitial content in A；Cj is jth column element in confusion matrix C
The sum of, it represents and divides j-th of community's interstitial content in B.If A is identical with B, NMI gets maximum value 1, if on the contrary,
A and B are entirely different, and NMI values are 0.
3. the experimental result of artificial synthesized network
Community's detection of verification CDEMO algorithms on the extension GN Benchmark networks of the propositions such as Lancichinetti
Performance.Comprising 128 nodes in each GN networks, it is divided into 4 communities, each community includes 32 nodes.Each node and society
Other nodes connect number of edges mesh for Zin inside area, and connect number of edges purpose for Zout with community's external node, and sum of the two is equal to node
The 16 of degree.Zout values are bigger, and the Lian Bianyue of node and community's external node is more, and community structure gets over unobvious, to detection
The detection performance requirement of algorithm is higher.
CDEMO algorithms are tested in Zout values gradually incremental 9 difference GN networks, only according to algorithm on each network
The accuracy and stability of the average value measure algorithm of vertical 30 gained NMI of operation, and calculated with 10 kinds of typical modularity optimizations
Method is compared (including CNM, GN, GATHB, ECGA, LGA, MA, UMDA, MOEA/D-Net, DECD and IDDE)), experimental result
As shown in Figure 3.
From figure 3, it can be seen that all algorithms can obtain optimal N MI values as Zout≤3, that is, it successfully is detected GN networks
Community structure.However, being gradually incremented by with Zout, the community structure of network, which becomes more to obscure, to be also more difficult to identify, is owned
NMI values obtained by algorithm all continuously decrease.It is worth noting that, CDEMO algorithms testing result is better than other 10 kinds of algorithms always, especially
It is to work as Zout>Afterwards, this illustrates that CDEMO algorithms are more accurate in community's detection of computer synthesis network and stablize.
It is further to test the Scalable Performance of CDEMO algorithms, it is gradually increased more massive in hybrid parameter μ
Test experiments are carried out on LFR Benchmark networks.The node degree of LFR networks is distributed as power-law distribution and community's scale can
Become, therefore closer real world network characteristic.Mixture of networks parameter μ is determined between community's interior nodes and other community's nodes altogether
The quantity on side is enjoyed, numerical value is bigger, and corresponding network community structure is fuzzyyer.0.7 interval is increased to from 0 using μ values in experiment
0.1 8 LFR networks, each LFR networks include 1000 nodes, community's scale value range be [10,50], each node
Average degree for 20, maximal degree 50.On each LFR networks, CDEMO algorithms independent operating 30 times, same to CNM, GATHB,
8 kinds of algorithms of MOGA-Net, MPSOA, ECGA, UMDA, MOEA/D-Net and DECD are compared, and detection gained is measured using NMI
The accuracy and stability of community structure, experimental result are as shown in Figure 4.
From Fig. 4 we may notice that being compared with other modularity optimization algorithms, CDEMO algorithms can be in 8 LFR nets
Optimal NMI values are obtained on network.Work as μ<The performance advantage of CDEMO algorithms is not obvious when 0.2, and being incremented by with μ values,
Advantage in the accuracy and stability of CDEMO algorithms gradually highlights.It is above-mentioned the experimental results showed that, CDEMO is in artificial synthesized net
There is preferable accuracy, stability and scalability on community's test problems of network.
4. real world network experiment result
CDEMO algorithm performances are verified on the real world social networks shown in table 3, and use 16 kinds of module recognizers
Performance comparison is carried out with CDEMO.Comparison algorithm is divided into three groups：First group comprising 6 kinds of tradition, the optimization of Qualitative module degree is calculated really
Method, including Fast Nm, CNM, GN, BGLL, MSFCM, FMM/H1；Second group comprising 4 kinds the modularity optimization algorithm based on GA,
Including GATHB, MOGA-Net, ECGA, and MOEA/D-Net；The last one group is excellent comprising 5 kinds of modularities based on PSO and DE
Change algorithm, including Meme-Net, MODPSO, DECD, CCDECD and IDDE.All algorithms independent operating on each test network
30 times, and measure optimal community using modularity Q and divide quality, table 5-7 record CDEMO and optimal Q obtained by other comparison algorithms
Value.
It is in table 5-7 the experimental results showed that, although all algorithms can identify the community structure in live network, and
Algorithm performance in first category has no larger difference, but the modularity optimization algorithm based on EA is compared to traditional certainty
Modularity optimization algorithm has apparent superiority.In the algorithm based on EA, DECD, CCDECD, IDDE and CDEMO are obtained
Q values it is relatively high, it was demonstrated that superiority based on DE optimisation strategies.Although the modularity optimization algorithm (Meme- based on PSO
Net and MODPSO) it can detect the optimal communities of Karate networks, but their performances on other networks are not to the utmost such as people
Meaning.Compared with DECD, CCDECD and IDDE, only CDEMO algorithms can always obtain optimal Q values, especially in Dolphin and
On Polbooks networks.
Fig. 5-8 shows the optimal community division result that CDEMO algorithms detect in 4 real world social networks.
The experimental results showed that other than artificial synthesized network, CDEMO algorithms can also efficiently identify community's knot of true social networks
Structure, it is more accurate, more stable compared with the effective modularity optimization algorithm in a variety of forward positions, thus also further demonstrate that algorithm entirety
The validity and advance that constringency performance improves.
7 tables involved in the present embodiment are introduced below：
1 Benchmark function details of table
1 Benchmark function details (Continued) of table
Table 2.DE algorithmic statement performances compare
Table 3.Benchmark real world network characteristics
Table 4.DE algoritic modules degree optimization performance compares
The optimal Q values of 5 Fast Nm CNM of table, GN, BGLL, MSFCM and FMM/H1 on real world network
The optimal Q values of table 6 GATHB, MOGA-Net, ECGA and MOEA/D-Net on real world network
The optimal Q values of table 7 Meme-Net, MODPSO, DECD, CCDECD and IDDE on real world network
The foregoing is only a preferred embodiment of the present invention, but protection scope of the present invention be not limited thereto,
Any one skilled in the art in the technical scope of present disclosure, according to the technique and scheme of the present invention and its
Inventive concept is subject to equivalent substitution or change, should be covered by the protection scope of the present invention.
Claims (6)
1. a kind of complex network community detection method that can improve global convergence performance, which is characterized in that specifically include：It improves
The step of global convergence performance of differential evolution algorithm, the step are specifically：
(1) classification adaptive differential class Mutation Strategy；
(2) dynamic self-adapting parameter adjustment；
(3) the carry out difference selection operation based on historical information.
2. a kind of complex network community detection method that can improve global convergence performance according to claim 1, feature
It is, adaptive differential class Mutation Strategy of classifying, concrete operations are as follows：
For each target individual Xi,tIf its ideal adaptation angle value fiIt is flat more than current entire population at individual fitness value
Mean is then classified as excellent individual, and globally optimal solution is closer in the position of search space；Therefore, in Xi,tIn it is good
Gene is reserved for strengthening the local search around individual, the corresponding vector V that makes a variationi,tGenerating mode is as follows：
Vi,t=Fi,t.Xpbesti,t+Wi,t.(Xr2,t-Xr3,t) (1)
Wherein, Xpbesti,tRepresent individual Xi,tIn the history optimal solution in preceding t generations, for enhancing individual exploring ability；Xr2,tAnd Xr3,t
It is randomly selected two Different Individuals from population, and meets condition r2 ≠ r3 ≠ i；Fi,tAnd Wi,tIt is XiControl ginseng
Number, numerical value is according to evolutionary generation and Xi,tIdeal adaptation angle value dynamic adjust；
For each target individual Xi,tIf its ideal adaptation angle value fiIt is flat less than current entire population at individual fitness value
Mean is then classified as poor individual, in the position of search space and globally optimal solution farther out；Therefore, strengthen it in population
In exchanging to promote global search between excellent individual, the corresponding vector V that makes a variationi,tGenerating mode is as follows：
Vi,t=Wi,t.Xr1,t+Ki,t.(Xgbest,t-Xi,t) (2)
Wherein Xr1,tIt is the randomly selected individual from population, and meets condition r1 ≠ i；Xgbest,tIt represents in current iteration population
Optimal solution, for enhancing Xi,tExploring ability；Wi,tAnd Ki,tIt is XiControl parameter, numerical value is according to evolutionary generation and Xi,t
Ideal adaptation angle value carry out dynamic adjustment.
3. a kind of complex network community detection method that can improve global convergence performance according to claim 1, feature
It is, dynamic self-adapting parameter adjustment：Three control parameters W, K, F, random element, social ingredient respectively in mutation process
And cognitive component；In addition, also there are one crucial control parameter CR in crossover operation, for determining each experiment individual ui,tIn
From variation individual Vi,tThe percentage of middle succession；Adjustment process is specific as follows：
4. a kind of complex network community detection method that can improve global convergence performance according to claim 1, feature
It is, the step of carry out difference selection operation based on historical information, specifically：By the history optimal solution of individual each in population
Xpbesti,tPopulation pbest_pop is formed, and is generated in initial phase, is updated after each evolutional operation；To in population
Each individual Xi,tIf its fitness value is improved during a certain evolutional operation, then newly-generated individual using as
Xi,tCurrent history optimal solution, and be saved in pbest_pop；After every generation evolutional operation, own in pbest_pop
Individual will substitute all individuals in population pop, and current optimal solution X is selected from pbest_popgbest,t。
5. a kind of complex network community detection method that can improve global convergence performance according to claim 1, feature
It is, this method is further included is based on the step of neighborhood information carries out community's amendment using improved, specifically：If a node meets
Community's correction conditions, then the node will be placed in again in its all affiliated community of neighborhood node, and the probability being placed in and neighbour
The scale of domain community is directly proportional.
6. a kind of complex network community detection method that can improve global convergence performance according to claim 1, feature
It is, this method further includes the modularity optimization algorithm based on classification differential evolution algorithm, specifically：
S1：Initialization of population；
S1.1 sets network parameter, including number of nodes n, adjacency matrix adj, community correction threshold δ；DE algorithm parameters, packet are set
Include individual dimension D, Population Size NP, population iterations t and maximum iteration tmax；
S1.2 is with the individual representation random initializtion population pop of community's label；
S2：It identifies and records optimal solution
S2.1 is identified and is recorded t for the optimum individual X in population popgbest,t；
S2.2 is identified and is recorded t for individual X each in population popi,tHistory optimal solution Xpbesti,t；By all population at individual
Xpbesti,tBuild initial population pbest_pop；
S3：When population iterations are less than population maximum iteration, population iterations are unsatisfactory for condition and then tie from adding one
The cycle of beam S3.1-S3.5；
S3.1 passes through adaptive classification differential variation construction of strategy variation population mutation_pop；
When i value for 1 to step a) in Population Size numberical range, is carried out to cycle e), if the value of i is not 1 to population
In the range of magnitude numerical value, then step a) is jumped out to e), end loop；
A) 3 different individual X are randomly selected from population popr1,t, Xr2,t, Xr3,t；
B) dynamic adjustment Mutation parameter Fi,t、wi,t、Ki,t；
C) according to fitness value Q to Xi,tClassify；
D) according to adaptive classification differential variation strategy generating variation individual Vi,t；
E) V is calculatedi,tModule angle value and and Xi,tIndividual is made comparisons, and more excellent individual is stored in pbest_pop；
If i is more than NP, leapfrog goes out a) suddenly to cycle e)；
S3.2 is based on neighborhood information and carries out community's amendment；
S3.3 is according to variation population mutation_pop and population pop structure cross-species crossover_pop；
When i value for 1 to step a) in Population Size numberical range, is carried out to cycle d), if the value of i is not 1 to population
In the range of magnitude numerical value, then step a) is jumped out to d), end loop；
A) i-th of individual u in cross-species is initializedi,t=xi,t；
B) dynamic adjustment cross parameter CRi,t；
C) by from variation individual Vi,tIt inherits community information and carrys out Adjustment Tests individual ui,t；
D) u is calculatedi,tModule angle value and be compared with i-th of individual in pbest_pop, retain compared with the figure of merit to pbest_
pop；
S3.4 is based on neighborhood information and carries out community's amendment；
S3.5 is by replacing all individual update pop in pbest_pop；
S4：Export the X in popgbest,tIt is divided as final optimal community, otherwise returns to S3 steps.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810035790.5A CN108154233A (en) | 2018-01-15 | 2018-01-15 | The complex network community detection method of global convergence performance can be improved |
| PCT/CN2018/086543 WO2019136893A1 (en) | 2018-01-15 | 2018-05-11 | Complex network community detection method capable of improving global convergence performance |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810035790.5A CN108154233A (en) | 2018-01-15 | 2018-01-15 | The complex network community detection method of global convergence performance can be improved |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108154233A true CN108154233A (en) | 2018-06-12 |
Family
ID=62461321
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810035790.5A Pending CN108154233A (en) | 2018-01-15 | 2018-01-15 | The complex network community detection method of global convergence performance can be improved |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN108154233A (en) |
| WO (1) | WO2019136893A1 (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102521655A (en) * | 2012-01-05 | 2012-06-27 | 西安电子科技大学 | Method for detecting dynamic network community on basis of non-dominated neighbor immune algorithm |
| CN103457800A (en) * | 2013-09-08 | 2013-12-18 | 西安电子科技大学 | Network community detection method based on M elite coevolution strategy |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106022464A (en) * | 2016-05-17 | 2016-10-12 | 大连理工大学 | Cluster-based high-dimensional multi-objective evolution method |
| CN106355193A (en) * | 2016-08-18 | 2017-01-25 | 大连理工大学 | Adaptive mutation particle swarm optimization algorithm |
-
2018
- 2018-01-15 CN CN201810035790.5A patent/CN108154233A/en active Pending
- 2018-05-11 WO PCT/CN2018/086543 patent/WO2019136893A1/en active Application Filing
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102521655A (en) * | 2012-01-05 | 2012-06-27 | 西安电子科技大学 | Method for detecting dynamic network community on basis of non-dominated neighbor immune algorithm |
| CN103457800A (en) * | 2013-09-08 | 2013-12-18 | 西安电子科技大学 | Network community detection method based on M elite coevolution strategy |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2019136893A1 (en) | 2019-07-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108133272A (en) | A kind of method of complex network community detection | |
| Tian et al. | Diversity assessment of multi-objective evolutionary algorithms: Performance metric and benchmark problems [research frontier] | |
| CN110519128B (en) | Random forest based operating system identification method | |
| Ekárt et al. | Maintaining the diversity of genetic programs | |
| CN110401564B (en) | Method for constructing command control hyper-network model based on relative hybrid preference | |
| CN103888541B (en) | Method and system for discovering cells fused with topology potential and spectral clustering | |
| CN108171331A (en) | A kind of modularity optimization method based on differential evolution algorithm | |
| CN107169871B (en) | Multi-relationship community discovery method based on relationship combination optimization and seed expansion | |
| Melamed et al. | Community structure in multi-mode networks: Applying an eigenspectrum approach | |
| Mishra et al. | Kohonen self organizing map with modified k-means clustering for high dimensional data set | |
| Lozano et al. | Modified fuzzy C-means algorithm for cellular manufacturing | |
| Ding et al. | Ceo: identifying overlapping communities via construction, expansion and optimization | |
| CN107016080A (en) | A kind of high-efficiency network packet classification method | |
| CN108154233A (en) | The complex network community detection method of global convergence performance can be improved | |
| CN109858533A (en) | A kind of RBF-MLP mixed structure neural network classifier of prepositus nucleus | |
| CN105488247A (en) | K-mean community structure mining method and apparatus | |
| CN108073748A (en) | A kind of guided missile damage effects group decision appraisal procedure based on intuitionistic fuzzy entropy weight method | |
| Deart et al. | Agglomerative Clustering of Network Traffic Based on Various Approaches to Determining the Distance Matrix | |
| Yuan et al. | A multi-granularity backbone network extraction method based on the topology potential | |
| CN114298414A (en) | Equipment system efficiency prediction and index optimization method | |
| CN113408602A (en) | Tree process neural network initialization method | |
| CN106055883A (en) | Transient stability assessment input characteristic validity analysis method based on improved Sammon mapping | |
| Jiang et al. | Evolving hard and easy traveling salesman problem instances: a multi-objective approach | |
| CN110049039A (en) | A kind of information centre's network-caching contamination detection method based on GBDT | |
| CN109753986A (en) | A kind of clustering method and device of the index information based on data block |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180612 |