CN107220497A - A kind of Circularity error evaluation method based on packet learning aid algorithm - Google Patents

A kind of Circularity error evaluation method based on packet learning aid algorithm Download PDF

Info

Publication number
CN107220497A
CN107220497A CN201710387364.3A CN201710387364A CN107220497A CN 107220497 A CN107220497 A CN 107220497A CN 201710387364 A CN201710387364 A CN 201710387364A CN 107220497 A CN107220497 A CN 107220497A
Authority
CN
China
Prior art keywords
student
mrow
algorithm
group
achievement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710387364.3A
Other languages
Chinese (zh)
Inventor
杨洋
李明
韦庆玥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201710387364.3A priority Critical patent/CN107220497A/en
Publication of CN107220497A publication Critical patent/CN107220497A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Abstract

Involved in the present invention is a kind of Circularity error evaluation method based on packet learning aid algorithm, first, according to lowest area principal, sets up the solution mathematical modeling of Circularity error evaluation;Secondly, by obtaining the measurement data of tested circularity key element, and set up mathematical modeling is combined, sets up the object function of problem;Finally, object function is solved using packet learning aid algorithm, key step includes the step such as mutual study between the Parameter Initialization procedure of algorithm, the teaching phase of student performance, student, and it is high for standard learning aid algorithmic procedure solving precision, the problems such as convergence rate is slower, while using population grouping strategy, strengthening the global optimizing ability of algorithm, so as to sufficiently improve computational accuracy, finally solved.Parameter is less needed for this method, algorithm stability is good, surveyed data can be sufficiently applied to, it is more novel in algorithm application, test function is more accurate relative to other algorithms most in use, faster, solution procedure complies fully with the Minimum Area principle in international standard to iteration speed, therefore computational solution precision is higher.

Description

A kind of Circularity error evaluation method based on packet learning aid algorithm
Technical field
It is more particularly to a kind of based on packet learning aid algorithm the invention belongs to the digital measuring field of machine components Circularity error evaluation method.
Background technology
Deviation from circular from is one of important Form and position error of revolving parts, affects the assembly precision of part and uses the longevity Life.Therefore, how accurately to obtain the control information of part circularity, it is ensured that the designing quality of engineering goods, be important at this stage One of research topic.
The national Specification error evaluation method of four kinds of circularity, including minimum circumscribed circle method, maximum inscribed circle method, most Small square law and minimum area method, wherein least square method and minimum area method is most widely used.For machine components circle The evaluation algorithm of error is spent, Wang Dongxia etc. is in its document《Application of the differential evolution algorithm in Circularity error evaluation》In employ Differential evolution algorithm is evaluated to deviation from circular from, and Yue Wuling etc. is in its document《Deviation from circular from based on imitative delta algorithm is fast Fast accurate evaluation》In employ imitative delta algorithm deviation from circular from evaluated, Chen Sa et al. is in its document《Lost based on improving The Circularity error evaluation of propagation algorithm》Middle use genetic algorithm is evaluated to flatness error.
In summary, continuing to develop with computer technology and artificial intelligence technology, Circularity error evaluation method is also got over Come more perfect, precision also more and more higher, and wherein topmost method is exactly that minimum area method is entered with reference to all kinds of optimized algorithms Row constantly research.Learning aid optimized algorithm is a kind of new global intelligent optimization algorithm proposed in 2012, its algorithm parameter Seldom, algorithm robustness very well, can further improve arithmetic accuracy by introducing corresponding strategies, be calculated because learning aid optimizes The calculating performance of method well, has been widely used in all kinds of engineering fields at this stage.
The present invention is on the basis of basic learning aid algorithm, to add grouping strategy to algorithm initialization process, enter one Step improves the search capability of algorithm, it is to avoid algorithm is absorbed in local optimum, to improve the evaluating precision of part flatness.
The content of the invention
It is an object of the invention to provide a kind of Circularity error evaluation method based on packet learning aid algorithm, to improve circularity The evaluating precision and computational stability of error.
A kind of Circularity error evaluation method based on packet learning aid algorithm, mainly including herein below:
Step one:The measured circle key element of tested part is determined, the measurement data points of part are obtained by three coordinate measuring machine For Pi(xi, yi) (i=1,2 ..., n);
Step 2:According to the lowest area principal of deviation from circular from, shown in the range formula of point to the center of circle, such as formula (1);Build The mathematical modeling of vertical deviation from circular from, shown in such as formula (2).Wherein RijFor the distance of measuring point to the center of circle, (xi, yi) sat for measuring point Mark, (aj, bj) it is concentric circles central coordinate of circle value, f (a`, b`) is deviation from circular from;
F (a ', b ')=min (max (Rij)-min(Rij)) (2)
Step 3:Read measurement data Pi(xi, yi) (i=1,2 ..., n), bring into formula (2), to learning aid algorithm Parameter initialized, it is main to include initialization student X=(x1, x2..., xi), i=1, wherein 2 ..., N, xiRepresent class Student i, N in level are student's quantity.For some student xi, there is x againi=(xi1, xi2..., xij), j=1,2 ..., D, xijThe subject j, D for representing student's i are total subject, and optimal solution is xt, i.e. fitness value is optimal for class's teacher's achievement, religion Iterations W, ramping constraint iterations w, into step 4 with learning optimized algorithm;
Step 4:It is to strengthen the information interchange ability between student using the purpose of population packet and strategy of shuffling, by building Vertical new group and the teacher of each group, so as to increase the ability of searching optimum of algorithm, it is to avoid algorithm is precocious.Population point is defined first Group policy:Determine group's number m in class, every group of number of student n, so that it is determined that produce initial student's number for F (F=m × N) fitness of each student is calculated, student is ranked up according to fitness by quality, and record the adaptation of whole population The optimal student x of angle valuet.Student is grouped by group number m, criteria for classifying is that the student that will be ranked first is put into the 1st Group, the 2nd student is put into the 2nd group, and m-th of student is put into m-th of group, and m+1 student then is put into the 1st again Group, by that analogy, the 2m student are put into m-th of group, wherein, q-th of group YqExpression formula is:Yq=x (q+m (p- 1)), p=1 ..., n;Q=1 ... m).Global optimum student x is recorded simultaneouslytWith Straight A x in groupT, m
Step 5:The fitness function value of each student is calculated, and according to achievement of the formula (3), (4) and (5) to student It is updated, in formula:X (i, j) ' is the achievement after student i subject j renewals;Before x (i, j) updates for student i subject j Achievement;Rand is the random number between [0,1];xt(j) it is the subject j achievements of teacher;T is the teaching factor;X (j) is student section Mesh j achievement;M (j) is j-th of section's purpose average achievement in class;Round is the bracket function rounded up.Teaching phase After the completion of, by contrasting the quality of fitness value, the renewal of solution is completed, into step 5;
X (i, j) '=x (i, j)+rand × (xt(j)-tM(j)) (3)
T=round [1+rand (0,1)] (4)
Step 6:Student performance is carried out again " to learn ", concrete mode is by randomly selecting two students, contrasting it Fitness value obtains size, the carry out study renewal small to fitness value for making fitness value big, and more new formula is such as shown in (6), formula In:x′aAchievement after learning for student a;xaAchievement before learning for student a;xbAchievement before learning for student b.When learning After the completion of journey, further with the achievement x of teachertContrasted, the renewal of solution is completed, into step 7;
x′a=xb+rand×|xa-xb| (6)
Step 7:Whether evaluation algorithm is when iteration completion, and after W iteration is met, algorithm is terminated, now globally optimal solution xtFitness function value be required flatness error value, if algorithm fails to complete W iteration, this iteration is obtained Globally optimal solution brings next iteration process into as class teacher.
Compared with prior art, the present invention has the advantage that:
By setting up the minimum area method parametric equation of Flatness error evaluation, asking for flatness is more intuitively reflected Mathematical modeling is solved, is a kind of general mathematical modeling, without the coordinate transformation process in computational geometry and measuring point preprocessing process Etc. cumbersome modeling process, surveyed data can be sufficiently applied to, and be can apply among substantial amounts of measuring point data, In algorithm application, the algorithm mutually introduces search by hill climbing strategy twice, entered on the basis of standard learning aid optimized algorithm One step improves the precision of algorithm, and faster, solution procedure complies fully with the Minimum Area principle in international standard to convergence rate, because This computational solution precision is higher.
Brief description of the drawings
Fig. 1 is the algorithm flow of the present invention;
Fig. 2 is the deviation from circular from iterative curve map of the present invention.
Embodiment
In order to become apparent from the design and advantage of the specific expression present invention, following processes will be with reference to accompanying drawing to whole The protocol procedures of evaluation are described in detail.
The invention provides a kind of Flatness error evaluation method based on secondary learning aid algorithm of climbing the mountain, it is main include with Lower content:
Step one:The measured circle key element of tested part is determined, the measurement data points of part are obtained by three coordinate measuring machine For Pi(xi, yi) (i=1,2 ..., n);
Step 2:According to the lowest area principal of deviation from circular from, shown in the range formula of point to the center of circle, such as formula (1);Build The mathematical modeling of vertical deviation from circular from, shown in such as formula (2).Wherein RijFor the distance of measuring point to the center of circle, (xi, yi) sat for measuring point Mark, (aj, bj) it is concentric circles central coordinate of circle value, f (a`, b`) is deviation from circular from;
F (a ', b ')=min (max (Rij)-min(Rij)) (2)
Step 3:Read measurement data Pi(xi, yi) (i=1,2 ..., n), bring into formula (2), to learning aid algorithm Parameter initialized, it is main to include initialization student X=(x1, x2..., xi), i=1, wherein 2 ..., N, xiRepresent class Student i, N in level are student's quantity.For some student xi, there is x againi=(xi1, xi2..., xij), j=1,2 ..., D, xijThe subject j, D for representing student's i are total subject, and optimal solution is xt, i.e. fitness value is optimal for class's teacher's achievement, religion Iterations W, ramping constraint iterations w, into step 4 with learning optimized algorithm;
Step 4:It is to strengthen the information interchange ability between student using the purpose of population packet and strategy of shuffling, by building Vertical new group and the teacher of each group, so as to increase the ability of searching optimum of algorithm, it is to avoid algorithm is precocious.Population point is defined first Group policy:Determine group's number m in class, every group of number of student n, so that it is determined that produce initial student's number for F (F=m × N) fitness of each student is calculated, student is ranked up according to fitness by quality, and record the adaptation of whole population The optimal student x of angle valuet.Student is grouped by group number m, criteria for classifying is that the student that will be ranked first is put into the 1st Group, the 2nd student is put into the 2nd group, and m-th of student is put into m-th of group, and m+1 student then is put into the 1st again Group, by that analogy, the 2m student are put into m-th of group, wherein, q-th of group YqExpression formula is:Yq=x (q+m (p- 1)), p=1 ..., n;Q=1 ... m).Global optimum student x is recorded simultaneouslytWith Straight A x in groupT, m
Step 5:The fitness function value of each student is calculated, and according to achievement of the formula (3), (4) and (5) to student It is updated, in formula:X (i, j) ' is the achievement after student i subject j renewals;Before x (i, j) updates for student i subject j Achievement;Rand is the random number between [0,1];xt(j) it is the subject j achievements of teacher;T is the teaching factor;X (j) is student section Mesh j achievement;M (j) is j-th of section's purpose average achievement in class;Round is the bracket function rounded up.Teaching phase After the completion of, by contrasting the quality of fitness value, the renewal of solution is completed, into step 5;
X (i, j) '=x (i, j)+rand × (xt(j)-tM(j)) (3)
T=round [1+rand (0,1)] (4)
Step 6:Student performance is carried out again " to learn ", concrete mode is by randomly selecting two students, contrasting it Fitness value obtains size, the carry out study renewal small to fitness value for making fitness value big, and more new formula is such as shown in (6), formula In:x′aAchievement after learning for student a;xaAchievement before learning for student a;xbAchievement before learning for student b.When learning After the completion of journey, further with the achievement x of teachertContrasted, the renewal of solution is completed, into step 7;
x′a=xb+rand×|xa-xb| (6)
Step 7:Whether evaluation algorithm is when iteration completion, and after W iteration is met, algorithm is terminated, now globally optimal solution xtFitness function value be required roundness error, if algorithm fails to complete W iteration, by this iteration obtain it is complete Office's optimal solution brings next iteration process into as class teacher.
It is described above, it is only that the present invention teaches good embodiment, but protection scope of the present invention is not limited thereto, Any one skilled in the art the invention discloses technical scope in, the variations and alternatives that can be readily occurred in, It should all be included within the scope of the present invention, therefore, protection scope of the present invention should be with scope of the claims It is defined.

Claims (1)

1. a kind of Flatness error evaluation method based on secondary learning aid algorithm of climbing the mountain, it is characterised in that use following steps:
Step one:The measured circle key element of tested part is determined, the measurement data points for obtaining part by three coordinate measuring machine are Pi (xi, yi) (i=1,2 ..., n);
Step 2:According to the lowest area principal of deviation from circular from, shown in the range formula of point to the center of circle, such as formula (1);Set up circle The mathematical modeling of error is spent, shown in such as formula (2).Wherein RijFor the distance of measuring point to the center of circle, (xi, yi) it is measuring point coordinate, (aj, bj) it is concentric circles central coordinate of circle value, f (a`, b`) is deviation from circular from;
<mrow> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>a</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>b</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
F (a ', b ')=min (max (Rij)-min(Rij)) (2)
Step 3:Read measurement data Pi(xi, yi) (i=1,2 ..., n), bring into formula (2), to the ginseng of learning aid algorithm Number is initialized, main to include initialization student X=(x1, x2..., xi), i=1, wherein 2 ..., N, xiRepresent in class Student i, N is student's quantity.For some student xi, there is x againi=(xi1, xi2..., xij), j=1,2 ..., D, xijTable The subject j, D of dendrography life i are total subject, and optimal solution is xt, i.e. fitness value is optimal for class's teacher's achievement, learning aid The iterations W of optimized algorithm, ramping constraint iterations w, into step 4;
Step 4:It is to strengthen the information interchange ability between student using the purpose of population packet and strategy of shuffling, it is new by setting up Group and each group teacher, so as to increase the ability of searching optimum of algorithm, it is to avoid algorithm is precocious.Population packet plan is defined first Slightly:Group's number m in class, every group of number of student n are determined, is counted so that it is determined that producing initial student's number for F (F=m × n) The fitness of each student is calculated, student is ranked up according to fitness by quality, and records the fitness value of whole population Optimal student xt.Student is grouped by group number m, it is the 1st small that criteria for classifying is that the student that will be ranked first is put into Group, the 2nd student is put into the 2nd group, and m-th of student is put into m-th of group, is then again put into m+1 student the 1st small Group, by that analogy, the 2m student are put into m-th of group, wherein, q-th of group YqExpression formula is:Yq=x (q+m (p-1)), p =1 ..., n;Q=1 ... m).Global optimum student x is recorded simultaneouslytWith Straight A x in groupT, m
Step 5:The fitness function value of each student is calculated, and the achievement of student is carried out according to formula (3), (4) and (5) Update, in formula:X (i, j) ' is the achievement after student i subject j renewals;X (i, j) is the achievement before student i subject j renewals; Rand is the random number between [0,1];xt(j) it is the subject j achievements of teacher;T is the teaching factor;X (j) is student's subject j's Achievement;M (j) is j-th of section's purpose average achievement in class;Round is the bracket function rounded up.Teaching phase is completed Afterwards, by contrasting the quality of fitness value, the renewal of solution is completed, into step 5;
X (i, j) '=x (i, j)+rand × (xt(j)-tM(j)) (3)
T=round [1+rand (0,1)] (4)
<mrow> <mi>M</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
Step 6:Student performance is carried out again " to learn ", concrete mode is by randomly selecting two students, contrasting its adaptation Angle value obtains size, the carry out study renewal small to fitness value for making fitness value big, and more new formula is such as shown in (6), in formula: x′aAchievement after learning for student a;xaAchievement before learning for student a;xbAchievement before learning for student b.Work as learning process After the completion of, further with the achievement x of teachertContrasted, the renewal of solution is completed, into step 7;
x′a=xb+rand×|xa-xb| (6)
Step 7:Whether evaluation algorithm is when iteration completion, and after W iteration is met, algorithm is terminated, now globally optimal solution xt's Fitness function value is required flatness error value, if algorithm fails to complete W iteration, the overall situation that this iteration is obtained Optimal solution brings next iteration process into as class teacher.
CN201710387364.3A 2017-05-26 2017-05-26 A kind of Circularity error evaluation method based on packet learning aid algorithm Pending CN107220497A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710387364.3A CN107220497A (en) 2017-05-26 2017-05-26 A kind of Circularity error evaluation method based on packet learning aid algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710387364.3A CN107220497A (en) 2017-05-26 2017-05-26 A kind of Circularity error evaluation method based on packet learning aid algorithm

Publications (1)

Publication Number Publication Date
CN107220497A true CN107220497A (en) 2017-09-29

Family

ID=59946725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710387364.3A Pending CN107220497A (en) 2017-05-26 2017-05-26 A kind of Circularity error evaluation method based on packet learning aid algorithm

Country Status (1)

Country Link
CN (1) CN107220497A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108319764A (en) * 2018-01-15 2018-07-24 湖北汽车工业学院 Evaluation method for spatial straightness errors method based on longicorn palpus searching algorithm
CN109323677A (en) * 2018-08-21 2019-02-12 上海隧道工程有限公司 Improve the Circularity error evaluation algorithm of cuckoo searching algorithm
CN109917647A (en) * 2019-03-06 2019-06-21 南京航空航天大学 One kind optimizing sliding-mode control based on the improved learning aid algorithm of instructional strategies and filled Spacecraft
CN110579201A (en) * 2019-07-25 2019-12-17 北京航空航天大学 Flatness evaluation method based on differential evolution algorithm

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108319764A (en) * 2018-01-15 2018-07-24 湖北汽车工业学院 Evaluation method for spatial straightness errors method based on longicorn palpus searching algorithm
CN109323677A (en) * 2018-08-21 2019-02-12 上海隧道工程有限公司 Improve the Circularity error evaluation algorithm of cuckoo searching algorithm
CN109323677B (en) * 2018-08-21 2020-08-11 上海隧道工程有限公司 Roundness error evaluation algorithm for improved cuckoo search algorithm
CN109917647A (en) * 2019-03-06 2019-06-21 南京航空航天大学 One kind optimizing sliding-mode control based on the improved learning aid algorithm of instructional strategies and filled Spacecraft
CN109917647B (en) * 2019-03-06 2020-12-11 南京航空航天大学 Teaching and learning algorithm improved based on teaching strategy and liquid-filled spacecraft optimization sliding mode control method
CN110579201A (en) * 2019-07-25 2019-12-17 北京航空航天大学 Flatness evaluation method based on differential evolution algorithm
CN110579201B (en) * 2019-07-25 2021-06-01 北京航空航天大学 Flatness evaluation method based on differential evolution algorithm

Similar Documents

Publication Publication Date Title
CN107220497A (en) A kind of Circularity error evaluation method based on packet learning aid algorithm
CN106971087A (en) A kind of Flatness error evaluation method based on secondary learning aid algorithm of climbing the mountain
CN105630988A (en) Method and system for rapidly detecting space data changes and updating data
CN104166731A (en) Discovering system for social network overlapped community and method thereof
CN105550199A (en) Point position clustering method and point position clustering apparatus based on multi-source map
CN104881689A (en) Method and system for multi-label active learning classification
CN105225149A (en) A kind of reference scoring defining method and device
CN112884226A (en) Multi-agent algorithm-based territorial spatial pattern simulation planning method and system
CN103455612B (en) Based on two-stage policy non-overlapped with overlapping network community detection method
CN107492120A (en) Point cloud registration method
CN108376183A (en) A kind of city CA model construction method based on principle of maximum entropy
CN110097529A (en) A kind of farmland Grading unit division methods and system based on semantic rules
CN108171793A (en) A kind of method for detecting lamination area triangle gridding
CN106960102A (en) A kind of space linearity assessment method based on secondary glowworm swarm algorithm of climbing the mountain
CN112954623A (en) Resident occupancy rate estimation method based on mobile phone signaling big data
CN115688607A (en) Band-shaped cross-sea quasi-geoid refinement method based on multi-channel spectrum combination
CN104424276A (en) Method and device for self-updating fingerprint database based on manifold learning
CN107330233A (en) A kind of power transmission tower design wind speed analysis method and device
CN104462797A (en) Increment integration algorithm used for procedure parameter online testing
CN107436969A (en) A kind of three-dimensional multi-target orientation method based on genetic algorithm
CN104036024B (en) It is a kind of based on GACUC and the spatial clustering method of Delaunay triangulation network
CN107491841A (en) Nonlinear optimization method and storage medium
Hassan et al. Differentiated multiple aggregations in multidimensional databases
Wang et al. Fast hyperbolic mapping based on the hierarchical community structure in complex networks
Bonham et al. post-MORDM: Mapping policies to synthesize optimization and robustness results for decision-maker compromise

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
DD01 Delivery of document by public notice
DD01 Delivery of document by public notice

Addressee: Yang Yang

Document name: Notification of before Expiration of Request of Examination as to Substance

WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170929