CN111104563B - Method for determining isomorphism of kinematic chain based on prime number asymmetric adjacency matrix - Google Patents

Method for determining isomorphism of kinematic chain based on prime number asymmetric adjacency matrix Download PDF

Info

Publication number
CN111104563B
CN111104563B CN201911094538.2A CN201911094538A CN111104563B CN 111104563 B CN111104563 B CN 111104563B CN 201911094538 A CN201911094538 A CN 201911094538A CN 111104563 B CN111104563 B CN 111104563B
Authority
CN
China
Prior art keywords
chain
matrix
motion
discrimination
kinematic
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.)
Active
Application number
CN201911094538.2A
Other languages
Chinese (zh)
Other versions
CN111104563A (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.)
Wuhan University of Science and Engineering WUSE
Original Assignee
Wuhan University of Science and Engineering WUSE
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 Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN201911094538.2A priority Critical patent/CN111104563B/en
Publication of CN111104563A publication Critical patent/CN111104563A/en
Application granted granted Critical
Publication of CN111104563B publication Critical patent/CN111104563B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90Details of database functions independent of the retrieved data types
    • G06F16/901Indexing; Data structures therefor; Storage structures
    • G06F16/9024Graphs; Linked lists
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Software Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Complex Calculations (AREA)

Abstract

The invention relates to a kinematic chain isomorphism judging method based on a prime asymmetric adjacency matrix, which comprises the following steps: numbering each component in the motion chain, assigning corresponding prime numbers to the components according to the number of the motion pairs of the components to construct a feature code M of the motion chain, comparing whether the feature codes of the motion chains are identical, if not, isomerism is realized, and if so, performing the next step; constructing an asymmetric adjacency matrix A for the motion chains with the same feature codes; constructing a discrimination matrix D for each matrix A, calculating the sum array of the discrimination matrix D, comparing whether the sum array of the discrimination matrix D of each motion chain is identical or not, if not, carrying out isomerism, and if so, carrying out the next step; and calculating the eigenvalue and the eigenvector of the discrimination matrix D of the motion chains which are the same as the array, comparing whether the eigenvalue and the eigenvector of each motion chain are the same or not, if so, carrying out isomerism, and if so, carrying out isomorphism. The invention has the advantages of visual expression and discrimination, and is very simple and efficient.

Description

Method for determining isomorphism of kinematic chain based on prime number asymmetric adjacency matrix
Technical Field
The invention relates to the technical field of mechanical kinematic chains, in particular to a kinematic chain isomorphism judging method based on a prime asymmetric adjacency matrix.
Background
In 1964, graph theory was introduced into the field of kinematic chain topology structure research, and the relationship between a mechanism diagram and a topology diagram is established by using a vertex to represent a connecting rod and a side to represent a joint. Because of the convenience of the computer in matrix calculation, graph theory provides a powerful mathematical tool in the research and development of institutions. In the motion chain type synthesis, whatever the method, it is very important to solve the uniqueness and comprehensiveness of the motion chain description, especially the motion chain with double twisting, which is a bottleneck problem in the field of mechanism topology research. To describe the compound hinges in the kinematic chain, students propose a method of bi-color topology, matrix identification and conversion. The size of the matrix constructed in this way changes.
In the kinematic chain analysis, the description of the kinematic chain plays an important role. The information contained in the adjacency matrix describing the motion chain construction is limited for conventional methods and may fail for some highly symmetric motion chains. Other methods may have the disadvantages of application range, insufficient visual expression, complex judgment method, and the like. Therefore, a motion chain isomorphism determination method with simpler and more efficient determination needs to be found.
Disclosure of Invention
The invention aims to provide a method for determining isomorphism of a moving chain based on a prime number asymmetric adjacency matrix, which can uniquely describe the structure of the moving chain and determine whether the moving chain is isomorphic or not, and is simpler and more efficient.
The invention solves the technical problems by adopting the following scheme:
a kinematic chain isomorphism judging method based on prime asymmetric adjacency matrix comprises the following steps:
s1: numbering each component in each moving chain, assigning corresponding prime numbers to each component in each moving chain according to the number of kinematic pairs of each component, and constructing feature codes M, M= [ M ] of each moving chain 1 ,m 2 ,m 3 …m n ]Wherein m is n Comparing whether the feature codes of all the motion chains are identical or not for the prime numbers corresponding to the number of the motion pairs of the component with the number of n, if not, isomerism is achieved, and if so, performing the next step;
s2: constructing an asymmetric adjacency matrix A for the kinematic chains with the same feature codes in the step S1, assigning corresponding prime numbers to the kinematic pairs of each component in each kinematic chain according to the number of the kinematic pairs, and determining corresponding element values in the matrix according to whether the two components are adjacent to each other or not to obtain an n multiplied by n asymmetric adjacency matrix of each kinematic chain:
wherein i, j represents the number of members in the kinematic chain, n represents the number of members in the kinematic chain, and the element a of the diagonal of the matrix i,j (i=j) is 0; other elements a of the matrix i,j (i+.j, i= … n, j= … n) is the prime number corresponding to the number of j kinematic pairs constructed when the component i and the construction j are contiguous; when building block i and build j are not contiguous then a i,j The value of (i+.j, i= … n, j= … n) is zero;
s3: constructing a corresponding discrimination matrix D for each motion chain asymmetric adjacent matrix A constructed in the step S2, wherein D=A' and calculating the sum array of the discrimination matrices D of each motion chain, and then comparing whether the sum array of the discrimination matrices D of each motion chain is identical or not, if not, the discrimination matrices D of each motion chain are heterogeneous, and if so, the next step is carried out;
s4: and (3) calculating a discrimination matrix D of the motion chains with the same discrimination matrix D as the array in the step S3, calculating the eigenvalue and the eigenvector of the discrimination matrix D, comparing whether the eigenvalue and the eigenvector of each motion chain are the same or not, if different, the motion chains are heterogeneous, and if the motion chains are the same, the motion chains are isomorphic.
Further, when comparing whether the feature codes of the motion chains are identical, firstly, according to M in the feature codes M of the motion chains n The value size is relative to M in each feature code M n Reordered and then compared.
Further, according to m n The value of M in the feature code M of each motion chain is from small to large or from large to small n And (5) rearrangement is carried out.
Further, in comparing whether the sum arrays of the discrimination matrices D of the respective kinematic chains are identical, the respective values in the sum arrays of the discrimination matrices D are rearranged from small to large or from large to small according to the respective values in the sum arrays of the discrimination matrices D, and then compared.
Further, when comparing whether the eigenvalues and eigenvectors of the motion chains are the same, the eigenvalues of the discrimination matrix D of each motion chain are rearranged in order from small to large or from large to small, and then compared.
Compared with the prior art, the invention has at least the following beneficial effects: the invention uses prime numbers to construct feature codes and corresponding adjacency matrix to clearly describe the adjacency relation difference between components in topological graph. The method realizes the one-to-one correspondence between the kinematic chain diagram and the components, is visual in expression, can judge whether isomorphism is realized without complex operation, and is very simple and efficient in judgment.
Drawings
FIG. 1 is a schematic view of an exercise training a1 according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of an exercise training a2 according to an embodiment of the present invention;
fig. 3 is a schematic view of the structure of the exercise amount a3 according to the embodiment of the present invention.
Detailed Description
For a better understanding of the present invention, the following examples are further illustrative of the present invention, but the contents of the present invention are not limited to the following examples only.
The invention provides a kinematic chain isomorphism judging method based on a prime asymmetric adjacency matrix, which comprises the following steps:
s1: three ten-bar kinematic chains with one degree of freedom in the double twist of 1 as shown in fig. 1-3 were analyzed. Firstly, the components in each motion chain are numbered, the number of the components is not limited, and the components only need to be numbered in sequence. And then, according to the number of kinematic pairs of each component in the three kinematic chains, corresponding prime numbers are given to the kinematic pairs, so that the feature codes Ma1, ma2 and Ma3 of the three kinematic chains are obtained, and the feature codes corresponding to the kinematic chain a1, the kinematic chain a2 and the kinematic chain a3 are as follows:
Ma1=[5,3,2,2,2,2,2,2,2,3]
Ma2=[5,3,2,2,2,2,2,2,2,3]
Ma3=[5,3,2,2,2,2,2,3,2,2];
since the components in the three kinematic chains are arbitrarily numbered, rearranging the values in the feature codes does not affect the feature codes. In order to facilitate comparing whether the feature codes of the three motion chains are the same, the values in the feature codes of the three motion chains can be rearranged from small to large or from large to small, and in this embodiment, the feature codes are rearranged from large to small to obtain the following feature codes corresponding to the three motion chains:
Ma1=[5,3,3,2,2,2,2,2,2,2]
Ma2=[5,3,3,2,2,2,2,2,2,2]
Ma3=[5,3,3,2,2,2,2,2,2,2];
by comparison, it can be found that the feature codes of the three motion chains are the same, and whether isomorphism exists cannot be determined, so that the next determination is needed.
S2: the 3 kinematic chains containing the compound strands are directly converted into an asymmetric matrix, corresponding prime numbers are given to the kinematic pairs of each component in each kinematic chain according to the number of the kinematic pairs, and corresponding element values in the matrix are determined according to whether every two components are adjacent to each other, so that an n multiplied by n asymmetric adjacent matrix of each kinematic chain is obtained:
wherein i, j represents the number of members in the kinematic chain, n represents the number of members in the kinematic chain, and the element a of the diagonal of the matrix i,j (i=j) is 0; other elements a of the matrix i,j (i+.j, i= … n, j= … n) is the prime number corresponding to the number of j kinematic pairs constructed when the component i and the construction j are contiguous; when building block i and build j are not contiguous then a i,j The value of (i+.j, i= … n, j= … n) is zero;
in this embodiment, the asymmetric adjacency matrix of exercise a1 is:
asymmetric adjacency matrix A a1 Element a of (2) ij (i.noteq.j) represents the structure in the motion chain a1Connection between member i and member j, e.g. element a 12 The value of (2) is 3, which means that the first member is adjacent to the second member and the element value of the second member is 3 (i.e., the number of kinematic pairs of the second member). Element a 21 A value of 5 indicates that the member two is adjacent to the member one and the element value of the member one is 5 (i.e., the number of kinematic pairs of the member 1).
Similarly, the asymmetric adjacency matrix of the kinematic chain a2 is:
the asymmetric adjacency matrix of the kinematic chain a3 is:
s3: constructing a discrimination matrix D for the matrix A constructed by each moving chain in the step S2, wherein D=A, calculating the sum array of the discrimination matrices D of each moving chain, comparing whether the sum arrays of the discrimination matrices D of each moving chain are identical or not, if not, isomerism is realized, and if so, the next step is carried out;
in this embodiment, the discrimination matrices of a1, a2 and a3 and the sum array of the discrimination matrices are respectively as follows:
in order to facilitate comparison of whether the sum arrays of the three motion chains are identical, the numerical values in the sum arrays of the discrimination matrices D of the 3 motion chains are rearranged from large to small to obtain a sum array H as follows. Of course, the values in the sum array of the 3 kinematic chains can also be rearranged from small to large.
H a1 =[136 128 120 108 83 79 72 44 39 36]
H a2 =[136 128 120 108 83 79 72 44 39 36]
H a3 =[136 124 112 108 79 79 72 55 44 36]
Thus, it can be seen from the above information that H a1 And H is a2 Identical to H a3 Different. Thus, both kinematic chains a1 and a2 are different from kinematic chain a 3. Due to H a1 And H is a2 It is therefore impossible to determine whether the moving chains a1 and a2 are isomorphic, and next the next step is to determine whether the moving chains a1 and a2 are isomorphic.
S4: and (3) calculating the discrimination matrix D of the motion chain which is the same as the array in the step (S3), calculating the eigenvalue and the eigenvector of the discrimination matrix D, comparing whether the eigenvalue and the eigenvector of the discrimination matrix D of each motion chain are the same or not, if different, the discrimination matrix D is heterogeneous, and if the eigenvalue and the eigenvector are the same, the discrimination matrix D is isomorphic.
In the present embodiment, the discrimination matrix D of the motion chains a1 and a2 is calculated a1 And D a2 Can be calculated by matlab to obtain the following characteristic value Ta1 and characteristic vector Sa1 of the motion chain a1 and the following characteristic value Ta2 and characteristic vector Sa2 of the motion chain a2:
as can be seen from the table, the eigenvalues and eigenvectors of the discrimination matrices D of the kinematic chain a1 and the kinematic chain a2 are the same, so the kinematic chain a1 and the kinematic chain a2 are isomorphic. So far, the determination of whether the kinematic chain is isomorphic is ended.
While the invention has been described with respect to the preferred embodiments, it will be understood that the invention is not limited thereto, but is capable of modification and variation without departing from the spirit of the invention, as will be apparent to those skilled in the art.

Claims (5)

1. A kinematic chain isomorphism judging method based on a prime number asymmetric adjacency matrix is characterized by comprising the following steps:
s1: numbering each component in each moving chain, assigning corresponding prime numbers to each component in each moving chain according to the number of kinematic pairs of each component, and constructing feature codes M, M= [ M ] of each moving chain 1 ,m 2 ,m 3 …m n ]Wherein m is n Comparing whether the feature codes of all the motion chains are identical or not for the prime numbers corresponding to the number of the motion pairs of the component with the number of n, if not, isomerism is achieved, and if so, performing the next step;
s2: constructing an asymmetric adjacency matrix A for the kinematic chains with the same feature codes in the step S1, assigning corresponding prime numbers to the kinematic pairs of each component in each kinematic chain according to the number of the kinematic pairs, and determining corresponding element values in the matrix according to whether the two components are adjacent to each other or not to obtain an n multiplied by n asymmetric adjacency matrix of each kinematic chain:
wherein i, j represents the number of members in the kinematic chain, n represents the number of members in the kinematic chain, and the element a of the diagonal of the matrix i,j (i=j) is 0; other elements a of the matrix i,j (i+.j, i= … n, j= … n) is the prime number corresponding to the number of j kinematic pairs constructed when the component i and the construction j are contiguous; when building block i and build j are not contiguous then a i,j The value of (i+.j, i= … n, j= … n) is zero;
s3: constructing a corresponding discrimination matrix D for each motion chain asymmetric adjacent matrix A constructed in the step S2, wherein D=A' and calculating the sum array of the discrimination matrices D of each motion chain, and then comparing whether the sum array of the discrimination matrices D of each motion chain is identical or not, if not, the discrimination matrices D of each motion chain are heterogeneous, and if so, the next step is carried out;
s4: and (3) calculating a discrimination matrix D of the motion chains with the same discrimination matrix D as the array in the step S3, calculating the eigenvalue and the eigenvector of the discrimination matrix D, comparing whether the eigenvalue and the eigenvector of each motion chain are the same or not, if different, the motion chains are heterogeneous, and if the motion chains are the same, the motion chains are isomorphic.
2. The method for determining isomorphism of moving chains based on prime number asymmetric adjacency matrix as recited in claim 1, wherein M in the feature code M of each moving chain is first determined according to whether the feature codes of each moving chain are identical or not n The value size is relative to M in each feature code M n Reordered in a certain order and compared.
3. The method for determining isomorphism of a moving chain based on a prime number asymmetric adjacency matrix as claimed in claim 2, characterized in that according to m n The value of M in the feature code M of each motion chain is from small to large or from large to small n And (5) rearrangement is carried out.
4. The method for determining the isomorphism of a moving chain based on a prime number asymmetric adjacency matrix according to claim 1, wherein, when comparing whether the sum array of the discrimination matrix D of each moving chain is the same, the values in the sum array of the discrimination matrix D are rearranged from small to large or from large to small according to the values in the sum array of the discrimination matrix D, and then compared.
5. The method for determining the isomorphism of a moving chain based on a prime number asymmetric adjacency matrix according to claim 1, wherein when comparing whether the eigenvalues and eigenvectors of the moving chains are identical, the eigenvalues of the determination matrix D of each moving chain are rearranged in order from small to large or from large to small, and then compared.
CN201911094538.2A 2019-11-11 2019-11-11 Method for determining isomorphism of kinematic chain based on prime number asymmetric adjacency matrix Active CN111104563B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911094538.2A CN111104563B (en) 2019-11-11 2019-11-11 Method for determining isomorphism of kinematic chain based on prime number asymmetric adjacency matrix

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911094538.2A CN111104563B (en) 2019-11-11 2019-11-11 Method for determining isomorphism of kinematic chain based on prime number asymmetric adjacency matrix

Publications (2)

Publication Number Publication Date
CN111104563A CN111104563A (en) 2020-05-05
CN111104563B true CN111104563B (en) 2024-04-16

Family

ID=70420410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911094538.2A Active CN111104563B (en) 2019-11-11 2019-11-11 Method for determining isomorphism of kinematic chain based on prime number asymmetric adjacency matrix

Country Status (1)

Country Link
CN (1) CN111104563B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002236446A (en) * 2001-02-09 2002-08-23 Nippon Steel Corp Method and device for information processing regarding prime number
CN105447277A (en) * 2015-12-28 2016-03-30 泉州装备制造研究所 Isomorph identification method for complex-hinge-containing kinematic chains based on topological characteristic loop codes
WO2017008043A1 (en) * 2015-07-08 2017-01-12 Brown University Homomorphic encryption
CN107391861A (en) * 2017-07-28 2017-11-24 珞石(山东)智能科技有限公司 Industrial robot loading kinetics parameter identification method independent of body kinetic parameter
CN107944097A (en) * 2017-11-07 2018-04-20 武汉科技大学 Plane enclosed cuts with scissors the numbering matrix description method of kinematic chain again
CN108170642A (en) * 2017-12-25 2018-06-15 武汉科技大学 Kinematic chain isomorphic products method based on number matrix
CN108628809A (en) * 2018-05-16 2018-10-09 武汉科技大学 Using the isomorphic products method of the kinematic chain of joint-joint adjacency matrix description

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3944566B2 (en) * 2001-07-30 2007-07-11 国立大学法人 東京大学 A method of dynamically determining joint acceleration in a link mechanism.
US9262380B2 (en) * 2013-05-22 2016-02-16 International Business Machines Corporation Calculating node centralities in large networks and graphs

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002236446A (en) * 2001-02-09 2002-08-23 Nippon Steel Corp Method and device for information processing regarding prime number
WO2017008043A1 (en) * 2015-07-08 2017-01-12 Brown University Homomorphic encryption
CN105447277A (en) * 2015-12-28 2016-03-30 泉州装备制造研究所 Isomorph identification method for complex-hinge-containing kinematic chains based on topological characteristic loop codes
CN107391861A (en) * 2017-07-28 2017-11-24 珞石(山东)智能科技有限公司 Industrial robot loading kinetics parameter identification method independent of body kinetic parameter
CN107944097A (en) * 2017-11-07 2018-04-20 武汉科技大学 Plane enclosed cuts with scissors the numbering matrix description method of kinematic chain again
CN108170642A (en) * 2017-12-25 2018-06-15 武汉科技大学 Kinematic chain isomorphic products method based on number matrix
CN108628809A (en) * 2018-05-16 2018-10-09 武汉科技大学 Using the isomorphic products method of the kinematic chain of joint-joint adjacency matrix description

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A kinematic chain isomorphism identification algorithm using optimized circuit simulation method;Huiliang Shang;2016 35th Chinese Control Conference (CCC);20160829;全文 *
单春成 ; 曹毅 ; 刘凯 ; .基于改进的杆组邻接矩阵的平面运动链类型自动综合.机械传动.2016,(09),全文. *
平面连杆变胞机构构型综合推演方法及应用;孙伟;中国博士学位论文全文数据库工程科技Ⅱ辑;20190815(第8期);C029-5 *
罗贤海 ; .机构运动链邻接矩阵的素数表示与同构判别.机械工程学报.2013,(05),第28-35页. *

Also Published As

Publication number Publication date
CN111104563A (en) 2020-05-05

Similar Documents

Publication Publication Date Title
CN108170642B (en) Kinematic chain isomorphism judgment method based on numbering matrix
Zhou et al. Spatio-temporal tensor completion for imputing missing internet traffic data
Kim et al. A symbiotic evolutionary algorithm for the integration of process planning and job shop scheduling
CN110677284B (en) Heterogeneous network link prediction method based on meta path
US20090028465A1 (en) Image upscaling technique
KR20070011304A (en) Radio communication system, reception device, demodulation method used for them, and program thereof
CN111104563B (en) Method for determining isomorphism of kinematic chain based on prime number asymmetric adjacency matrix
CN111008359B (en) Hamming number kinematic chain isomorphism judging method based on prime number asymmetric adjacency matrix
US20150323512A1 (en) Computer-Implemented Systems and Methods for Generating Generalized Fractional Designs
Filipiak et al. Optimality and efficiency of circular neighbor balanced designs for correlated observations
Chye et al. Selection criteria of the compression algorithm in information-measuring system
Dudin et al. A Retrial BMAP/PH/N queueing system with Markov modulated retrials
CN1714535A (en) Method and device for designing a data network
CN110866293B (en) Isomorphism judging method, system, device and medium based on prime number layering matrix
Lin Designs of variable resolution
CN110533159A (en) A method of deep learning data set quality and model interpretation are improved based on Jie's scientific guidance
CN109768890A (en) The oriented Weighted Complex Networks networking method of symbolism based on STL decomposition method
Durand et al. Interaction terms in non-linear PLS via additive spline transformation
Weishaupl et al. Parallelization of cellular neural networks for image processing on cluster architectures
CN113933263B (en) Near infrared spectrum model building method based on index fitting
CN107958027A (en) A kind of Sensor Network data capture method ensured with QoS
Rajih et al. Alternating least squares identification of under-determined mixtures based on the characteristic function
CN106374936A (en) Power grid real-time control method based on compressed sensing technology
Camacho Páez et al. Group-wise Partial Least Square Regression
Ron et al. A fast multigrid algorithm for energy minimization under planar density constraints

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
GR01 Patent grant
GR01 Patent grant