CN103324983B - A kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm - Google Patents

A kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm Download PDF

Info

Publication number
CN103324983B
CN103324983B CN201310228196.5A CN201310228196A CN103324983B CN 103324983 B CN103324983 B CN 103324983B CN 201310228196 A CN201310228196 A CN 201310228196A CN 103324983 B CN103324983 B CN 103324983B
Authority
CN
China
Prior art keywords
individual
optimum
individuality
colony
fitness
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.)
Expired - Fee Related
Application number
CN201310228196.5A
Other languages
Chinese (zh)
Other versions
CN103324983A (en
Inventor
杨平
石顺义
刘东静
李霞龙
赵艳芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangsu University
Original Assignee
Jiangsu University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu University filed Critical Jiangsu University
Priority to CN201310228196.5A priority Critical patent/CN103324983B/en
Publication of CN103324983A publication Critical patent/CN103324983A/en
Application granted granted Critical
Publication of CN103324983B publication Critical patent/CN103324983B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
  • Peptides Or Proteins (AREA)

Abstract

The invention discloses a kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm, comprise the steps: 1) topological diagram of its correspondence is formed according to the structure of mechanism kinematics chain;2) genetic coding of adjacency matrix is set according to topological diagram;3) two topological diagram adjacency matrix to be determined are taken, it is determined that object function;4) initial antibodies population is randomly generated;Calculate the fitness of each individuality, and utilize fitness size that it is ranked up;5) it is carried out optimum maintaining strategy;6) concentration of each individuality, total individual number in individual bulk concentration=close individual sum/colony are calculated;7) antibody population obtained in step 5) is carried out optimum Crossover Strategy;8) carry out mutation operation, adopt inverse operators to realize;9) local searching operator is combined with genetic algorithm;10) form immune genetic hybrid algorithm, carry out the Isomorphism Identification computing of mechanism kinematics chain.Instant invention overcomes the shortcoming that genetic algorithm easily converges to local optimum.

Description

A kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm
Technical field
The present invention relates to isomorphism identification method of mechanism kinematics chain, refer in particular to a kind of for isomorphism identification method of mechanism kinematics chain in Creative Mechanism Design and intelligent CAD.
Background technology
Owing to, in Creative Mechanism Design and intelligent CAD, the research of mechanism structure is very important, it may be assumed that how to select required type in the kinematic chain being made up of multiple components and kinematic pair.And this problem is exactly MECHANISM KINEMATICS CHAIN ISOMORPHISM IDENTIFICATION at all, along with the proposition of this proposition, the expert of countries in the world this respect and scholar propose a series of solution.As in 2000, Feng Chun etc. are published in the paper of " Machine Design and research " supplementary issue P40~41 page in 2000 to propose genetic algorithm is applied to MECHANISM KINEMATICS CHAIN ISOMORPHISM IDENTIFICATION one straightforward procedure, but the method makes its population diversity existed in calculating process be deteriorated due to oversimplification, difficulty converges to global extremum, slack-off being even difficult to of algorithm computing late convergence is restrained;Chinese patent ZL200910184138.0 " a kind of isomorphism identification method of mechanism kinematics chain based on pseudo-hybridization hybrid genetic algorithm ", although improve genetic algorithm in the slack-off problem of computing late convergence, but also without the problem that the population diversity well solving exist in calculating process is deteriorated.
Summary of the invention
It is an object of the invention to the deficiency for overcoming existing method for designing to exist, better solve that existing Creative Mechanism Design cycle length, cost be high and also inefficient problem.Creative Mechanism Design is promoted to develop to efficient and low cost direction and provide a kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm.
The technical scheme is that and adopt following steps successively:
1) topological diagram of its correspondence is formed according to the structure of mechanism kinematics chain;
2) genetic coding of adjacency matrix is set according to topological diagram;
3) two topological diagram adjacency matrix to be determined are taken, it is determined that object function;
4) randomly generating initial antibodies population, Population Size is M;Calculate the fitness of each individuality, and utilize fitness size that it is ranked up;If just then preserving the individuality of relevance grade optimum for colony in specific one private variable, forwarding step 6) to, otherwise, continuing;
5) antibody population of step 4) is carried out optimum maintaining strategy, its method is as follows: judge t for the relevance grade of optimum individual in colony P (t) whether more than the relevance grade of optimum individual in private variable, if more than, optimum individual in private variable is replicated portion in colony P (t), deleting worst individuality therein and preserve colony P(t simultaneously) in optimum individual in above-mentioned private variable, substitute original numerical value;If less than, only optimum individual in private variable is replicated portion and arrive in colony P (t), delete worst individuality therein simultaneously;
6) calculating the concentration of each individuality, total individual number in individual bulk concentration=close individual sum/colony, individual bulk concentration is big, then suppress this antibody, allow for close number of individuals and tail off, thus maintain the multiformity of individual in population;Select the individuality that fitness is high and concentration is relatively low and carry out clone's duplication, and the individuality that fitness low concentration is higher is suppressed;Antagonist group performs to select and replicate operation;
7) antibody population obtained in step 5) is carried out optimum Crossover Strategy: selecting t for the 10%*M(M at the most end of the fitness in colony P (t) is population number) individual antibody and the optimum individual P being saved in variable*T () is intersected respectively, the method for intersection is: the individual P selectedi(t) and P*T () is put in a copulation pond, i=1~0.1M, according to selected Crossover Strategy (i.e. pseudo-crossover process), two individualities selected is carried out intersection operation, obtains a pair offspring individual Pi(t) and P*T (), then with individual PiT () substitutes the individual P in colonyi(t) and by individuality P*T () is given up to fall;
8) in order to avoid being absorbed in local optimum, adopt inverse operators that antibody each in above-mentioned antibody population is carried out mutation operation;
9) local searching operator is combined with genetic algorithm;
10) form optimal choice criss-cross inheritance hybrid algorithm, carry out the Isomorphism Identification computing of mechanism kinematics chain.
The invention has the beneficial effects as follows:
1, in the present invention, Immunity Operator and genetic algorithm is utilized to combine, it is ensured that algorithm is in the multiformity of running colony.Overcome the shortcoming that genetic algorithm easily converges to local optimum;
2, in the present invention, utilize the feature of optimum maintaining strategy that the optimum individual occurred so far not only will not be made to be chosen, intersect and mutation operation is lost and destroys, it is ensured that algorithmic statement is to the probability of globally optimal solution.Overcome the shortcoming that genetic algorithm is difficult to restrain in the algorithm later stage;
3, in the present invention, the feature utilizing optimum Crossover Strategy is so that evolving in the individual quick direction higher to fitness that relevance grade is relatively low, accelerates the algorithm convergence rate to global optimum.Overcome the shortcoming that genetic algorithm convergence rate is slow.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the present invention.
Detailed description of the invention
The idiographic flow of inventive algorithm describes as follows:
Step 1: by formula C=TATT, it is the initial population P(0 of M that conversion T randomly generates a scale);
Step 2: calculate target function value and the fitness of each individuality, stores individuality maximum for fitness in locally optimal solution to private variable;
Step 3: if just for population, then forwarding step 5) to;Otherwise, continue;
Step 4: determine the relevance grade of each individuality, and press the sequence of fitness size;Carry out optimum maintaining strategy, its method is as follows: judge P(t) in the relevance grade of optimum individual whether more than the relevance grade of optimum individual in private variable, if more than, optimum individual in private variable is replicated portion in P (t), deleting worst individuality therein and preserve P(t simultaneously) in optimum individual in above-mentioned private variable, substitute original numerical value;If less than, only optimum individual in private variable is replicated portion and arrive in P (t), delete worst individuality therein simultaneously;
Step 5: define according to the similarity of antibody and concentration, calculate the concentration of each antibody;
Step 6: select the individuality that fitness is high and concentration is relatively low and carry out clone's duplication, and the individuality that fitness low concentration is higher is suppressed;Antagonist group performs to select and replicate operation;
Step 7: antagonist group carries out optimum Crossover Strategy: selecting t for the 10%*M(M at the most end of the fitness in colony P (t) is population number) individual antibody and the optimum individual P being saved in variable*T () is intersected respectively, the method for intersection is: the individual P selectedi(t) and P*T () is put in a copulation pond, i=1~0.1M, according to selected Crossover Strategy (, i.e. pseudo-crossover process), two individualities selected carries out intersecting operation, obtains a pair offspring individual Pi(t) and P*T (), then with individual PiT () substitutes the individual P in colonyi(t) and by individuality P*T () is given up to fall;
Step 8: in order to avoid being absorbed in local optimum, adopts inverse operators that antibody each in above-mentioned antibody population is carried out mutation operation;
Step 9: antagonist group performs Local Search;
Step 10: find optimal solution, two kinematic chain isomorphisms;Otherwise, tripe systems.
Wherein in step 2, object function and fitness function computational methods are as follows:
Taking the adjacency matrix that A, B are two topological diagram G1, G2 to be determined, object function takes:
min f ( x ) = Σ i , j = 0 d - 1 | C ( i , j ) - B ( i , j ) | ,
C=TATT,
In formula, d is number of vertex, is also individual lengths, T=P1P2……Pn, wherein PiFor matrix
A is the adjacency matrix of mechanism topological diagram G1,
B is the adjacency matrix of mechanism topological diagram G2,
The matrix that C adjacency matrix A obtains after carrying out line translation and rank transformation through transformation matrix T.
Affinity judgement schematics is as follows:
F(x)=1/f(x)
Can extrapolate the value of F (x) (0,1] between, it is seen that f (x) is more little, and F (x) is more big
The wherein Local Search in step 9, for individual x=(x1,x2,…,xd) algorithm is as follows:
Step 1: randomly select down two elements being designated as i, j;
Step 2: it is exactly an original individual neighbours x' that exchange the two element obtains new individuality;
Step 3: compare the fitness F of x and x', if F (x') > F (x), then x=x';Otherwise give up x'.

Claims (1)

1., based on an isomorphism identification method of mechanism kinematics chain for immune genetic hybrid algorithm, it is characterized in that adopting successively following steps:
1) topological diagram of its correspondence is formed according to the structure of mechanism kinematics chain;
2) genetic coding of adjacency matrix is set according to topological diagram;
3) two topological diagram adjacency matrix to be determined are taken, it is determined that object function;
4) randomly generating initial antibodies population, Population Size is M;Calculate the fitness of each individuality, and utilize fitness size that it is ranked up;If just then preserving the individuality of relevance grade optimum for colony in specific one private variable, forward step 6 to), otherwise, continue;
5) to step 4) antibody population carry out optimum maintaining strategy, its method is as follows: if t=1, namely it is initial population, as long as individuality optimum in initial population is saved in private variable, if t > 1, then must judge t for the relevance grade of optimum individual in colony P (t) whether more than the relevance grade of optimum individual in private variable, if more than, optimum individual in private variable is replicated portion and arrives in colony P (t), delete optimum individual in worst individuality therein and preservation colony P (t) simultaneously and substitute original numerical value in above-mentioned private variable;If less than, only optimum individual in private variable is replicated portion and arrive in colony P (t), delete worst individuality therein simultaneously;
6) calculating the concentration of each individuality, total individual number in individual bulk concentration=close individual sum/colony, individual bulk concentration is big, then suppress this antibody, allow for close number of individuals and tail off, thus maintain the multiformity of individual in population;Select the individuality that fitness is high and concentration is relatively low and carry out clone's duplication, and the individuality that fitness low concentration is higher is suppressed;Antagonist group performs to select and replicate operation;
7) to step 5) in the antibody population that obtains carry out optimum Crossover Strategy: select t for minimum 10%*M the antibody of the fitness in colony P (t) and the optimum individual P being saved in variable*T () is intersected respectively, the method for intersection is: the individual P selectedi(t) and P*T () is put in a copulation pond, i=1~0.1M, according to selected Crossover Strategy, i.e. and pseudo-crossover process, two individualities selected carry out intersecting operating, obtains a pair offspring individual Pi(t) and P*T (), then with individual PiT () substitutes the individual P in colonyi(t) and by individuality P*T () is given up to fall;
8) in order to avoid being absorbed in local optimum, adopt inverse operators that antibody each in above-mentioned antibody population is carried out mutation operation;
9) local searching operator is combined with genetic algorithm;Local Search is for individual x=(x1,x2,…,xd) step be:
Step 1: randomly select down two elements being designated as i, j;
Step 2: it is exactly an original individual neighbours x ' that exchange the two element obtains new individuality;
Step 3: compare the fitness F of x and x ', if F (x ') > F (x), then x=x ';Otherwise give up x ';
10) form optimal choice criss-cross inheritance hybrid algorithm, carry out the Isomorphism Identification computing of mechanism kinematics chain.
CN201310228196.5A 2013-06-08 2013-06-08 A kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm Expired - Fee Related CN103324983B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310228196.5A CN103324983B (en) 2013-06-08 2013-06-08 A kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310228196.5A CN103324983B (en) 2013-06-08 2013-06-08 A kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm

Publications (2)

Publication Number Publication Date
CN103324983A CN103324983A (en) 2013-09-25
CN103324983B true CN103324983B (en) 2016-07-20

Family

ID=49193709

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310228196.5A Expired - Fee Related CN103324983B (en) 2013-06-08 2013-06-08 A kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm

Country Status (1)

Country Link
CN (1) CN103324983B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104021425B (en) * 2014-05-19 2017-04-19 中国人民解放军国防科学技术大学 Meme evolutionary algorithm for solving advancing-delay scheduling problem
CN104268628B (en) * 2014-07-23 2017-02-22 江苏大学 Mechanism kinematic link isomorphism recognition design method based on fish swarm algorithm
CN104268629B (en) * 2014-09-15 2017-02-15 西安电子科技大学 Complex network community detecting method based on prior information and network inherent information
CN105447277B (en) * 2015-12-28 2018-11-09 泉州装备制造研究所 A kind of isomorphism identification method containing multiple hinge kinematic chain based on topological characteristic loop code
CN107694094A (en) * 2017-10-10 2018-02-16 电子科技大学 A kind of bridge gambling party sample generating method
CN109408504B (en) * 2018-08-13 2021-07-13 浙江理工大学 Method for screening separable topological graph based on adjacency matrix and total loop number
CN110851926B (en) * 2019-10-31 2023-09-01 武汉科技大学 Planetary gear train isomorphism judging method, system and medium based on traversal loop
CN110852022B (en) * 2019-10-31 2023-05-23 武汉科技大学 Planetary gear train isomorphism judging method, system and medium based on circuit model
CN110851925B (en) * 2019-10-31 2024-02-20 武汉科技大学 Planetary gear train isomorphism judging method, system and medium based on improved adjacency matrix
CN110826214B (en) * 2019-10-31 2023-08-11 武汉科技大学 Prime number layering-based planetary gear train isomorphism judging method, system and medium
CN110909424B (en) * 2019-10-31 2023-08-15 武汉科技大学 Planetary gear train isomorphism judging method, system and medium based on adjacency matrix

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101655928A (en) * 2009-08-25 2010-02-24 江苏大学 Isomorphism identification method of mechanism kinematics chain based on pseudo-hybridization hybrid genetic algorithm
JP2012008983A (en) * 2010-06-28 2012-01-12 Hiroshima Univ Graph isomorphism determination device, graph isomorphism determination method, program and evaluation system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101655928A (en) * 2009-08-25 2010-02-24 江苏大学 Isomorphism identification method of mechanism kinematics chain based on pseudo-hybridization hybrid genetic algorithm
JP2012008983A (en) * 2010-06-28 2012-01-12 Hiroshima Univ Graph isomorphism determination device, graph isomorphism determination method, program and evaluation system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
An improved artificial immune algorithm for mechanism kinematic chain isomorphism identification;Kehan Zeng et al;《International Journal of Materials and Structural Integrity》;20081231;第2卷(第4期);第383-395页 *
基于混合人工免疫优化算法的机构运动链同构识别研究;曾科翰;《中国优秀硕士学位论文全文数据库(工程科技II辑)》;20091031;C029-5 *
改进的免疫遗传算法在地产决策支持系统中的研究及应用;史特特;《中国优秀硕士学位论文全文数据库(信息科技辑)》;20120229;I138-1445 *

Also Published As

Publication number Publication date
CN103324983A (en) 2013-09-25

Similar Documents

Publication Publication Date Title
CN103324983B (en) A kind of isomorphism identification method of mechanism kinematics chain based on immune genetic hybrid algorithm
CN104764980B (en) A kind of distribution line failure Section Location based on BPSO and GA
CN104810829B (en) A kind of network structure simplified process method of power distribution network reconfiguration
CN108921294A (en) A kind of gradual piece of knowledge distillating method accelerated for neural network
CN104332995A (en) Improved particle swarm optimization based power distribution reconstruction optimization method
Wilck IV et al. A genetic algorithm for the split delivery vehicle routing problem
CN112671029A (en) Multi-stage fault recovery method for distribution network with distributed power supply
CN104268628B (en) Mechanism kinematic link isomorphism recognition design method based on fish swarm algorithm
CN103593287A (en) Genetic-algorithm-based method for automatically generating data stream test cases
CN107704959B (en) Electric power system grid reconstruction method for achieving dimensionality reduction based on ant colony algorithm
CN108413963A (en) Bar-type machine people's paths planning method based on self study ant group algorithm
CN104268629A (en) Complex network community detecting method based on prior information and network inherent information
CN106611378A (en) Priority encoding-based hybrid genetic algorithm for solving job-shop scheduling problem
CN102196527A (en) Route recovery method and recovery protocol for mobile Sink wireless sensor network
CN104463328A (en) Order crossover multi-filial-generation genetic algorithm for solving traveling salesman problem
CN105512755A (en) Decomposition-based multi-objective distribution estimation optimization method
CN102915407B (en) Prediction method for three-dimensional structure of protein based on chaos ant colony algorithm
CN110687397B (en) Active power distribution network fault positioning method based on improved artificial fish swarm algorithm
CN108537370A (en) Especially big basin water station group Optimization Scheduling based on hybrid intelligent dimension-reduction algorithm
CN107579517A (en) The whether feasible determination methods of gained solution in the quantum telepotation reconstruct of Complicated Distribution Network
CN104217252A (en) Power transmission network power flow chart automatic layout optimization algorithm and system based on genetic algorithm
CN103426127A (en) Urban power network planning method based on pseudo-crossover taboo hybrid genetic algorithm
CN107368454A (en) A kind of GPU of the sparse lower trigonometric equation group of a large amount of isomorphisms pushes away method before accelerating
CN104392317A (en) Project scheduling method based on genetic culture gene algorithm
CN103679270A (en) Heuristic adaptive immune clonal method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160720

Termination date: 20200608