CN112884355B - Proportional electromagnet electromagnetic force linear characteristic evaluation method based on multiple correlation coefficients - Google Patents

Proportional electromagnet electromagnetic force linear characteristic evaluation method based on multiple correlation coefficients Download PDF

Info

Publication number
CN112884355B
CN112884355B CN202110279685.8A CN202110279685A CN112884355B CN 112884355 B CN112884355 B CN 112884355B CN 202110279685 A CN202110279685 A CN 202110279685A CN 112884355 B CN112884355 B CN 112884355B
Authority
CN
China
Prior art keywords
current
working
electromagnetic force
operating
discrete
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
CN202110279685.8A
Other languages
Chinese (zh)
Other versions
CN112884355A (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.)
Changsha University of Science and Technology
Original Assignee
Changsha University of 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 Changsha University of Science and Technology filed Critical Changsha University of Science and Technology
Priority to CN202110279685.8A priority Critical patent/CN112884355B/en
Publication of CN112884355A publication Critical patent/CN112884355A/en
Application granted granted Critical
Publication of CN112884355B publication Critical patent/CN112884355B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

The invention discloses a method for evaluating the linear characteristic of electromagnetic force of a proportional electromagnet based on a complex correlation coefficient, and belongs to the field of performance evaluation and analysis of the proportional electromagnet. Firstly, dividing a full working condition plane formed by the working current of an electromagnet and the working stroke of an armature; then obtaining the electromagnetic force output by all discrete working condition points through simulation or experiment means; calculating the average value of the electromagnetic force of the discrete working condition points under the same working current to obtain a regression sample point; performing linear regression fitting on the sample points to obtain a regression equation; obtaining the overall complex correlation coefficient between the current and the electromagnetic force according to the overall average calculation or the weighted average calculation of the working area under different working conditionsR 2(ii) a According to the overall complex correlation coefficientR 2And (4) judging the linear characteristic of the electromagnetic force of the proportional electromagnet. The method is simple and reliable, can objectively and quantitatively evaluate the linear characteristic of the electromagnetic force of the proportional electromagnet, and can be used for measuring the product performance of the proportional electromagnet and comparing and analyzing the performance of different proportional electromagnets.

Description

Proportional electromagnet electromagnetic force linear characteristic evaluation method based on multiple correlation coefficients
Technical Field
The invention belongs to the field of performance evaluation and analysis of proportional electromagnets, and particularly relates to a method for evaluating the linear characteristics of the electromagnetic force of a proportional electromagnet based on a complex correlation coefficient.
Background
The proportional electromagnet as the electro-mechanical converter of electro-hydraulic proportional control element is an automatic control element with wide application, can make the liquid pressure and flow change continuously and proportionally with the control signal, and has the advantages of low cost, simple structure, good universality and strong anti-pollution capability. In order to realize the proportional control characteristic of the proportional electromagnet, the electromagnetic force of the proportional electromagnet is required to have good linear characteristic, that is, the control current and the output electromagnetic force have good linear relation. However, currently, the evaluation of the linear characteristic of the electromagnetic force of the proportional electromagnet is mostly based on subjective qualitative judgment of a current-force curve diagram (chengda. research of novel disc-type proportional electromagnet [ D ]. university of zhejiang, 2009 ]), and an objective and quantitative evaluation index is lacked, so that the product performance of the proportional electromagnet is not easy to be measured, and the performance of different proportional electromagnets is not easy to compare and analyze.
Disclosure of Invention
In order to solve the problems of the existing evaluation method, the invention provides a simple, reliable, objective and quantitative evaluation method for the electromagnetic force linear characteristic of the proportional electromagnet.
The purpose of the invention is realized as follows:
the method for evaluating the electromagnetic force linear characteristic of the proportional electromagnet based on the complex correlation coefficient comprises the following steps:
step 1, dividing a full-working-condition plane;
step 2, obtaining the electromagnetic force of the discrete working condition points;
step 3, calculating the overall complex correlation coefficient R between the current and the electromagnetic force2；
Step 4, according to the integral complex correlation coefficient R between the current and the electromagnetic force2The linear characteristic of the electromagnetic force of the proportional electromagnet is judged according to the numerical value of the proportional electromagnet.
As a further explanation of the above evaluation method:
further, the dividing method of the full working condition plane in the step 1 comprises the following steps:
1.1, determining the working current of the proportional electromagnet and the working stroke range of the armature, wherein the working current range is marked as [ ia,id]The range of the armature stroke is denoted as [ x ]a,xd]；
1.2, equally dividing and dispersing the full working condition plane formed by the working current and the working stroke, and further obtaining corresponding discrete working condition points (i) in the full working condition planen,xm) Wherein inExpressed as the operating current range ia,id]The nth operating current value, x, corresponding to the discrete number of equal divisionsmIs a working stroke range [ xa,xd]The corresponding mth working stroke position is equally divided and dispersed.
Further, the method for obtaining the electromagnetic force of the discrete operating point in the step 2 is a simulation or experiment means, and specifically comprises the following steps:
2.1, adjusting the working stroke to a certain discrete working condition point value and then keeping the working stroke unchanged;
2.2, sequentially adjusting the working current to different discrete working condition point values, and respectively measuring corresponding electromagnetic force;
and 2.3, adjusting the working stroke to another discrete working condition point value, keeping the working stroke unchanged, and repeating the process to obtain the electromagnetic force of all the discrete working condition points.
Further, the overall complex correlation coefficient R between the current and the electromagnetic force in step 32The calculation method comprises the following steps:
3.1, respectively calculating the average value F (i) of the electromagnetic force of the discrete working points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure GDA0003372891680000021
In the formula F (i)n,xm) Representing discrete operating points (i)n,xm) Corresponding electromagnetic force, f represents the number of equally divided working strokes;
3.2 sample point (i) by least squares or partial least squaresn,F(in)a) Linear regression is carried out to obtain a linear regression equation of current-electromagnetic force
Figure GDA0003372891680000022
Wherein
Figure GDA0003372891680000023
Is inCorresponding regression average electromagnetic force predicted values, wherein k and b are regression coefficients;
3.3 calculating the difference between the same currentsMean value of the electromagnetic force in the operating stroke
Figure GDA0003372891680000024
The expression is
Figure GDA0003372891680000025
F(in)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current, and e is the part for equally dividing the working current;
3.4 calculating the sum of the squares of the total deviations SST, which is expressed as
Figure GDA0003372891680000026
Wherein F (i)t)aAverage value of electromagnetic force of discrete operating points having different operating strokes for equal operating currents itRepresenting the operating current range [ i ]a,id]The t-th working current value;
3.5, calculating the regression square sum SSR, wherein the expression is
Figure GDA0003372891680000027
Wherein
Figure GDA0003372891680000028
The regression average electromagnetic force predicted value corresponding to the discrete operating point with the same working current and different working strokes;
3.6, calculating to obtain an integral complex correlation coefficient R2The expression is
Figure GDA0003372891680000029
Further, the overall complex correlation coefficient R between the current and the electromagnetic force in step 32The calculation method comprises the following steps:
3.1, respectively calculating the average value F (i) of the electromagnetic force of the discrete working points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure GDA0003372891680000031
In the formula F (i)n,xm) Representing discrete operating points (i)n,xm) Corresponding electromagnetic force inRepresenting the operating current range [ i ]a,id]Inner nth operating current value, xmRepresenting the working stroke range [ xa,xd]The mth working stroke position in the inner part, and f represents the number of equally divided working strokes;
3.2, defining the working area of primary and secondary current, and the interval [ ia,ib)、(ic,id]For the current minor operating region, interval [ ib,ic]Is a current main working area, wherein ia、ib、ic、id∈[ia,id]And i isa<ib<ic<id；
3.3 sample points (i) of each current working area by using least square method or partial least square methodn,F(in)a) Performing linear regression to obtain three groups of current-electromagnetic force linear regression equations:
Figure GDA0003372891680000032
Figure GDA0003372891680000033
Figure GDA0003372891680000034
wherein
Figure GDA0003372891680000035
Respectively i in different current working regionsα、iβ、iγCorresponding regression mean electromagnetic force prediction value, k1、b1、k2、b2、k3、b3Is a regression coefficient, where iαIs a current working region [ ia,ib) Of the alpha operating current value, iβIs a current working region [ ib,ic]Of (d) the beta-th operating current value, iγIs a current working region (i)c,id]The γ -th operating current value of (1);
3.4, respectively calculating the average value of the average electromagnetic force in the same current working area under the same current working condition
Figure GDA0003372891680000036
Figure GDA0003372891680000037
Wherein
Figure GDA0003372891680000038
iα∈[ia,ib) And g represents that the operating point falls in the current operating region [ i ]a,ib) The number of (2);
Figure GDA0003372891680000039
iβ∈[ib,ic]h represents that the operating point falls in the current operating region [ i ]b,ic) The number of (2);
Figure GDA00033728916800000310
iγ∈(ic,id]and j represents that the operating point falls in the current operating region (i)c,id]The number of (2);
3.5, respectivelyCalculating the sum of squared deviations SST of each current working areaab、SSTbc、SSTcdWherein
Figure GDA00033728916800000311
iα∈[ia,ib)；
Figure GDA00033728916800000312
iβ∈[ib,ic]；
Figure GDA00033728916800000313
iγ∈(ic,id]；
Wherein F (i)α)a、F(iβ)a、F(iγ)aThe average value of the electromagnetic force of discrete working condition points with equal working current and different working strokes in each current working area;
3.6, calculating the regression square sum SST of each current working area respectivelyab、SSTbc、SSTcdWherein
Figure GDA0003372891680000041
iα∈[ia,ib)；
Figure GDA0003372891680000042
iβ∈[ib,ic]；
Figure GDA0003372891680000043
iγ∈(ic,id]；
Wherein
Figure GDA0003372891680000044
Respectively is itA regression average electromagnetic force predicted value corresponding to any current working condition point in the current working area;
3.7 respectively calculating the complex correlation coefficient R of each current working area2 ab、R2 bc、R2 cdWherein
Figure GDA0003372891680000045
Figure GDA0003372891680000046
Figure GDA0003372891680000047
3.8 setting weighting coefficients of each current working area, ia,ib)、[ib,ic]、(ic,id]The weighting coefficients of the corresponding complex correlation coefficients are respectively K1、K2、K3；
3.9, carrying out weighted calculation on the electromagnetic force complex correlation coefficient of each working current region according to the current region to obtain the electromagnetic force integral complex correlation coefficient R2The expression is R2＝K1·R2 ab+K2·R2 bc+K3·R2 cd。
Further, the overall complex correlation coefficient R between the current and the electromagnetic force in step 32The calculation method comprises the following steps:
3.1, define the working area of the primary and secondary travel, the interval [ xa,xb)、(xc,xd]For the secondary working area of the journey, interval [ xb,xc]Is the main working area of the stroke, wherein xa、xb、xc、xd∈[xa,xd]And x isa<xb<xc<xd；
3.2, settingWeighting factor of each trip working area, trip working area xa,xb)、[xb,xc]、(xc,xd]The weighting coefficients of the corresponding complex correlation coefficients are S1、S2、S3；
3.3, respectively calculating the average value F (i) of the electromagnetic force of the discrete working condition points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure GDA0003372891680000048
In the formula xδ∈[xa,xb)，xε∈[xb,xc]，xζ∈(xc,xd]，F(in,xδ) Is an operating current inTime travel working area [ x ]a,xb) Any working stroke position xδElectromagnetic force of F (i)n,xε) Is an operating current inTime travel working area [ x ]b,xc]Any working stroke position xεElectromagnetic force of F (i)n,xζ) Is an operating current inTime travel working area (x)c,xd]Any working stroke position xζU, v, w respectively indicate that the working point falls in the stroke working area [ x ]a,xb)、[xb,xc]、(xc,xd]The number of (2);
3.4 sample point (i) by least squares or partial least squaresn,F(in)a) Linear regression is carried out to obtain a linear regression equation of current-electromagnetic force
Figure GDA0003372891680000051
Wherein
Figure GDA0003372891680000052
Is inCorresponding regression average powerPredicting the magnetic force, wherein k and b are regression coefficients;
3.5, calculating the average value of the electromagnetic force under different working strokes of the same current
Figure GDA0003372891680000053
The expression is
Figure GDA0003372891680000054
Wherein F (i)n)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current, and e is the part for equally dividing the working current;
3.6 calculating the sum of squared deviations SST, which is expressed as
Figure GDA0003372891680000055
Wherein F (i)t)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current;
3.7 calculating regression Square sum SSR with the expression of
Figure GDA0003372891680000056
Wherein
Figure GDA0003372891680000057
The regression average electromagnetic force predicted value corresponding to the discrete operating point with the same working current and different working strokes;
3.8, calculating to obtain an integral complex correlation coefficient R2The expression is
Figure GDA0003372891680000058
Further, the overall complex correlation coefficient R between the current and the electromagnetic force in step 32Calculation methodThe method comprises the following steps:
3.1, define the working area of the primary and secondary travel, the interval [ xa,xb)、(xc,xd]For the secondary working area of the journey, interval [ xb,xc]Is the main working area of the stroke, wherein xa、xb、xc、xd∈[xa,xd]And x isa<xb<xc<xd；
3.2 setting weighting coefficients of each stroke working area, the stroke working area [ xa,xb)、[xb,xc]、(xc,xd]The weighting coefficients of the corresponding complex correlation coefficients are S1、S2、S3；
3.3, respectively calculating the average value F (i) of the electromagnetic force of the discrete working condition points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) The expression is
Figure GDA0003372891680000059
In the formula xδ∈[xa,xb)，xε∈[xb,xc]，xζ∈(xc,xd]U, v, w respectively indicate that the operating point falls in the stroke operating region [ x ]a,xb)、[xb,xc]、(xc,xd]The number of (2);
3.4, defining the working area of primary and secondary current, and the interval [ ia,ib)、(ic,id]For the current minor operating region, interval [ ib,ic]Is a current main working area, wherein ia、ib、ic、id∈[ia,id]And i isa<ib<ic<id；
3.5 sample points (i) of each current working area by using least square method or partial least square methodn,F(in)a) Performing linear regression to obtain three groups of current-electromagnetic force linear regression equations:
Figure GDA0003372891680000061
Figure GDA0003372891680000062
Figure GDA0003372891680000063
wherein
Figure GDA0003372891680000064
Respectively i in different current working regionsα、iβ、iγCorresponding regression mean electromagnetic force prediction value, k1、b1、k2、b2、k3、b3Is a regression coefficient;
3.6, respectively calculating the average value of the average electromagnetic force in the same current working area under the same current working condition
Figure GDA0003372891680000065
Figure GDA0003372891680000066
Wherein
Figure GDA0003372891680000067
iα∈[ia,ib) And g represents that the operating point falls in the current operating region [ i ]a,ib) The number of (2);
Figure GDA0003372891680000068
iβ∈[ib,ic]h represents that the operating point falls in the current operating region [ i ]b,ic) The number of (2);
Figure GDA0003372891680000069
iγ∈(ic,id]and j represents that the operating point falls in the current operating region (i)c,id]The number of (2);
3.7, respectively calculating the sum of squared deviations SST of each current working areaab、SSTbc、SSTcdWherein
Figure GDA00033728916800000610
iα∈[ia,ib)；
Figure GDA00033728916800000611
iβ∈[ib,ic]；
Figure GDA00033728916800000612
iγ∈(ic,id]；
Wherein F (i)α)a、F(iβ)a、F(iγ)aThe average value of the electromagnetic force of discrete working condition points with equal working current and different working strokes in each current working area;
3.8, respectively calculating the regression square sum SST of each current working areaab、SSTbc、SSTcdWherein
Figure GDA00033728916800000613
iα∈[ia,ib)；
Figure GDA00033728916800000614
iβ∈[ib,ic]；
Figure GDA00033728916800000615
iγ∈(ic,id]；
Wherein
Figure GDA00033728916800000616
Respectively obtaining regression average electromagnetic force predicted values corresponding to any current working condition point in the current working area;
3.9 respectively calculating the complex correlation coefficient R of each current working area2 ab、R2 bc、R2 cdWherein
Figure GDA0003372891680000071
Figure GDA0003372891680000072
Figure GDA0003372891680000073
3.10 setting weighting coefficients of each current working area, current working area [ ia,ib)、[ib,ic]、(ic,id]The weighting coefficients of the corresponding complex correlation coefficients are respectively K1、K2、K3；
3.11, carrying out weighted calculation on the electromagnetic force complex correlation coefficient of each working current region according to the current region to obtain the electromagnetic force integral complex correlation coefficient R2The expression is R2＝K1·R2 ab+K2·R2 bc+K3·R2 cd。
Further, the method for determining the linear characteristic of the electromagnetic force of the proportional electromagnet in step 4 includes: integral complex correlation coefficient R between current and electromagnetic force2The closer to 1, the better the linear characteristic of the electromagnetic force of the proportional electromagnet is; integral complex correlation coefficient R between current and electromagnetic force2The closer to 0, the worse the linear characteristic of the proportional electromagnet electromagnetic force.
The invention has the advantages that: the method for evaluating the linear characteristic of the electromagnetic force of the proportional electromagnet quantitatively judges the linear characteristic of the electromagnetic force of the proportional electromagnet by adopting the overall complex correlation coefficient between the current and the electromagnetic force, considers the influence of the proportional electromagnet under different working currents and working stroke working conditions, combines a method of regional weighting calculation, and can comprehensively, objectively and quantitatively evaluate the linear characteristic of the electromagnetic force of the proportional electromagnet. The method can be widely used for measuring the product performance of the proportional electromagnet and comparing and analyzing the performance of electromagnets with different proportions.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a flow chart of calculating a complex correlation coefficient from a ensemble averaged electromagnetic force;
FIG. 3 is a schematic view of the full operating face of a proportional electromagnet;
FIG. 4 is a flow chart of evaluation based on current partition weighting;
FIG. 5 is a schematic diagram of current-partition-based electromagnetic force complex correlation coefficient weighting;
FIG. 6 is a flow chart of evaluation based on trip partition weighting;
FIG. 7 is a schematic diagram of electromagnetic force weighting based on trip zones;
fig. 8 is a flow chart of evaluation based on current and trip partition quadratic weighting.
Detailed Description
In order to more specifically describe the present invention, the following detailed description is provided for the technical solution of the present invention with reference to the accompanying drawings and the specific embodiments.
As shown in fig. 1, an embodiment of the present invention provides a method for evaluating linear characteristics of electromagnetic force of a proportional electromagnet based on a complex correlation coefficient, and the specific implementation manner is as follows.
The first embodiment is as follows:
as shown in fig. 2, the present embodiment specifically includes the following steps:
step one, dividing a full-working-condition plane
1.1 determining the working current of the proportional electromagnet and the working stroke range of the armature, the working current range being denoted as [ ia,id]The range of the armature stroke is denoted as [ x ]a,xd]As shown in fig. 3;
1.2 equally dividing and dispersing the whole working condition plane formed by the working current and the working stroke, and further obtaining corresponding discrete working condition points (i) in the whole working condition planen,xm) I.e. discrete points in the working area shown in FIG. 3, where inExpressed as the operating current range ia,id]The nth operating current value, x, corresponding to the discrete number of equal divisionsmIs a working stroke range [ xa,xd]The corresponding mth working stroke position is equally divided and dispersed.
Step two, obtaining the electromagnetic force of the discrete operating point
The method for obtaining the electromagnetic force of the discrete working condition point by means of simulation or experiment specifically comprises the following steps:
2.1 adjusting the working stroke to a certain discrete working condition point value and then keeping the working stroke unchanged;
2.2 adjusting the working current to different discrete working condition point values in sequence, and respectively measuring corresponding electromagnetic force;
and 2.3, adjusting the working stroke to another discrete operating point value, keeping the working stroke unchanged, and repeating the process to obtain the electromagnetic force of all the discrete operating points.
Step three, calculating the integral complex correlation coefficient R between the current and the electromagnetic force2
3.1 calculating the average value F (i) of the electromagnetic force of the discrete operating points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure GDA0003372891680000081
In the formula F (i)n,xm) Representing discrete operating conditionsPoint (i)n,xm) Corresponding electromagnetic force inRepresenting the operating current range [ i ]a,id]Inner nth operating current value, xmRepresenting the working stroke range [ xa,xd]The mth working stroke position in the inner part, and f represents the number of equally divided working strokes;
3.2 applying least squares or partial least squares to the sample points (i)n,F(in)a) Linear regression is carried out to obtain a linear regression equation of current-electromagnetic force
Figure GDA0003372891680000082
Wherein
Figure GDA0003372891680000083
Is inCorresponding regression average electromagnetic force predicted values, wherein k and b are regression coefficients;
3.3 calculating the average value of the electromagnetic force under different working strokes of the same current
Figure GDA0003372891680000084
The expression is
Figure GDA0003372891680000085
Wherein F (i)n)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current, and e is the part for equally dividing the working current;
3.4 calculate the sum of the squares of the total deviations SST, expressed as
Figure GDA0003372891680000091
Wherein F (i)t)aAverage value of electromagnetic force of discrete operating points having different operating strokes for equal operating currents itRepresenting the operating current range [ i ]a,id]The t-th working current value;
3.5 calculating the regression Square sum SSR, the expression is
Figure GDA0003372891680000092
Wherein
Figure GDA0003372891680000093
The regression average electromagnetic force predicted value corresponding to the discrete operating point with the same working current and different working strokes;
3.6 calculating to obtain an integral complex correlation coefficient R2The expression is
Figure GDA0003372891680000094
Step four, according to the integral complex correlation coefficient R between the current and the electromagnetic force2The linear characteristic of the electromagnetic force of the proportional electromagnet is judged
Integral complex correlation coefficient R between current and electromagnetic force2The closer to 1, the better the linear characteristic of the electromagnetic force of the proportional electromagnet is; integral complex correlation coefficient R between current and electromagnetic force2The closer to 0, the worse the linear characteristic of the proportional electromagnet electromagnetic force.
The second embodiment is as follows:
referring to fig. 4, the difference between the present embodiment and the first embodiment is that the overall complex correlation coefficient R between the current and the electromagnetic force is calculated in the third step2The method specifically comprises the following steps: 3.1 calculating the average value F (i) of the electromagnetic force of the discrete operating points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure GDA0003372891680000095
In the formula F (i)n,xm) Representing discrete operating points (i)n,xm) Corresponding electromagnetic force, f representing the working strokeThe number of divided portions;
3.2 defining the working area of primary and secondary current, the interval [ ia,ib)、(ic,id]For the current minor operating region, interval [ ib,ic]Is a current main working area, wherein ia、ib、ic、id∈[ia,id]And i isa<ib<ic<id；
3.3 sample points (i) for each current operating region by using least squares or partial least squaresn,F(in)a) Performing linear regression to obtain three groups of current-electromagnetic force linear regression equations:
Figure GDA0003372891680000096
Figure GDA0003372891680000097
Figure GDA0003372891680000098
wherein
Figure GDA0003372891680000099
Respectively i in different current working regionsα、iβ、iγCorresponding regression mean electromagnetic force prediction value, k1、b1、k2、b2、k3、b3Is a regression coefficient, where iαIs a current working region [ ia,ib) Of the alpha operating current value, iβIs a current working region [ ib,ic]Of (d) the beta-th operating current value, iγIs a current working region (i)c,id]The γ -th operating current value of (1);
3.4 calculating the average value of the average electromagnetic force in the same current working area under the same current working condition
Figure GDA0003372891680000101
Figure GDA0003372891680000102
Wherein
Figure GDA0003372891680000103
iα∈[ia,ib) And g represents that the operating point falls in the current operating region [ i ]a,ib) The number of (2);
Figure GDA0003372891680000104
iβ∈[ib,ic]h represents that the operating point falls in the current operating region [ i ]b,ic) The number of (2);
Figure GDA0003372891680000105
iγ∈(ic,id]and j represents that the operating point falls in the current operating region (i)c,id]The number of (2);
3.5 calculating the sum of squared deviations SST of each current working area respectivelyab、SSTbc、SSTcdWherein
Figure GDA0003372891680000106
iα∈[ia,ib)；
Figure GDA0003372891680000107
iβ∈[ib,ic]；
Figure GDA0003372891680000108
iγ∈(ic,id]；
Wherein F (i)α)a、F(iβ)a、F(iγ)aThe average value of the electromagnetic force of discrete working condition points with equal working current and different working strokes in each current working area;
3.6 calculating the regression Square sum SST of each current working areaab、SSTbc、SSTcdWherein
Figure GDA0003372891680000109
iα∈[ia,ib)；
Figure GDA00033728916800001010
iβ∈[ib,ic]；
Figure GDA00033728916800001011
iγ∈(ic,id]；
Wherein
Figure GDA00033728916800001012
Respectively obtaining regression average electromagnetic force predicted values corresponding to any current working condition point in the current working area;
3.7 respectively calculating the complex correlation coefficient R of each current working area2 ab、R2 bc、R2 cdWherein
Figure GDA00033728916800001013
Figure GDA00033728916800001014
Figure GDA00033728916800001015
3.8 setting weighting coefficients for each current working region, current working region [ ia,ib)、[ib,ic]、(ic,id]Weighting factor K of the corresponding complex correlation coefficient1、K2、K3As shown in fig. 5;
3.9 weighting the complex correlation coefficient of the electromagnetic force in each working current region according to the current region to obtain the overall complex correlation coefficient R between the current and the electromagnetic force2The expression is R2＝K1·R2 ab+K2·R2 bc+K3·R2 cd。
The third concrete implementation mode:
referring to fig. 6, the difference between the present embodiment and the first embodiment is that the overall complex correlation coefficient R between the current and the electromagnetic force is calculated in the third step2The method specifically comprises the following steps:
3.1 defining the working area of the primary and secondary strokes, the interval [ xa,xb)、(xc,xd]For the secondary working area of the journey, interval [ xb,xc]Is the main working area of the stroke, wherein xa、xb、xc、xd∈[xa,xd]And x isa<xb<xc<xd；
3.2 setting weighting coefficients of each travel working area, travel working area [ xa,xb)、[xb,xc]、(xc,xd]Weighting coefficient S of corresponding complex correlation coefficient1、S2、S3As shown in fig. 7;
3.3 calculating the average value F (i) of the electromagnetic force of the discrete working points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure GDA0003372891680000111
In the formula xδ∈[xa,xb)，xε∈[xb,xc]，xζ∈(xc,xd]，F(in,xδ) Is an operating current inTime travel working area [ x ]a,xb) Any working stroke position xδElectromagnetic force of F (i)n,xε) Is an operating current inTime travel working area [ x ]b,xc]Any working stroke position xεElectromagnetic force of F (i)n,xζ) Is an operating current inTime travel working area (x)c,xd]Any working stroke position xζU, v, w respectively indicate that the operating point falls in the stroke operating region [ x ]a,xb)、[xb,xc]、(xc,xd]The number of (2);
3.4 sample points (i) are aligned using least squares or partial least squaresn,F(in)a) Linear regression is carried out to obtain a linear regression equation of current-electromagnetic force
Figure GDA0003372891680000112
Wherein
Figure GDA0003372891680000113
Is inCorresponding regression average electromagnetic force predicted values, wherein k and b are regression coefficients;
3.5 calculating the average value of the electromagnetic force under different working strokes of the same current
Figure GDA0003372891680000114
The expression is
Figure GDA0003372891680000115
Wherein F (i)n)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current, and e is the part for equally dividing the working current;
3.6 calculate the sum of the squares of the total deviations SST, expressed as
Figure GDA0003372891680000116
Wherein F (i)t)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current;
3.7 calculating the regression Square sum SSR, the expression is
Figure GDA0003372891680000117
Wherein
Figure GDA0003372891680000121
The regression average electromagnetic force predicted value corresponding to the discrete operating point with the same working current and different working strokes;
3.8 calculating to obtain an integral complex correlation coefficient R2The expression is
Figure GDA0003372891680000122
The fourth concrete implementation mode:
referring to fig. 8, the difference between the present embodiment and the first embodiment is that the overall complex correlation coefficient R between the current and the electromagnetic force is calculated in the third step2The method specifically comprises the following steps:
3.1 defining the working area of the primary and secondary strokes, the interval [ xa,xb)、(xc,xd]For the secondary working area of the journey, interval [ xb,xc]Is the main working area of the stroke, wherein xa、xb、xc、xd∈[xa,xd]And x isa<xb<xc<xd；
3.2 setting the weighting coefficient of each stroke working area, and falling into the stroke secondary working area [ x ] for the working pointa,xb)、(xc,xd]According to a weighting coefficient S1、S3Performing weighted calculation to the working point falling in the travel main working area [ x ]b,xc]By a weighting factor S2Performing a weighting calculation, as shown in fig. 7;
3.3 calculating the average value F (i) of the electromagnetic force of the discrete working points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure GDA0003372891680000123
In the formula xδ∈[xa,xb)，xε∈[xb,xc]，xζ∈(xc,xd]U, v, w respectively indicate that the operating point falls in the stroke operating region [ x ]a,xb)、[xb,xc]、(xc,xd]The number of (2);
3.4 define the working area of primary and secondary current, interval [ ia,ib)、(ic,id]For the current minor operating region, interval [ ib,ic]Is a current main working area, wherein ia、ib、ic、id∈[ia,id]And i isa<ib<ic<id；
3.5 sample points (i) of each current working area by using least square method or partial least square methodn,F(in)a) Performing linear regression to obtain three groups of current-electromagnetic force linear regression equations:
Figure GDA0003372891680000124
Figure GDA0003372891680000125
Figure GDA0003372891680000126
wherein
Figure GDA0003372891680000127
Respectively i in different current working regionsα、iβ、iγCorresponding regression mean electromagnetic force prediction value, k1、b1、k2、b2、k3、b3Is a regression coefficient;
3.6 calculating average value of average electromagnetic force in same current working region
Figure GDA0003372891680000128
Figure GDA0003372891680000129
Wherein
Figure GDA0003372891680000131
iα∈[ia,ib) And g represents that the operating point falls in the current operating region [ i ]a,ib) The number of (2);
Figure GDA0003372891680000132
iβ∈[ib,ic]h represents that the operating point falls in the current operating region [ i ]b,ic) The number of (2);
Figure GDA0003372891680000133
iγ∈(ic,id]and j represents that the operating point falls in the current operating region (i)c,id]The number of (2);
3.7 calculating the sum of squared deviations SST of each current working area respectivelyab、SSTbc、SSTcdWherein
Figure GDA0003372891680000134
iα∈[ia,ib)；
Figure GDA0003372891680000135
iβ∈[ib,ic]；
Figure GDA0003372891680000136
iγ∈(ic,id]；
Wherein F (i)α)a、F(iβ)a、F(iγ)aThe average value of the electromagnetic force of discrete working condition points with equal working current and different working strokes in each current working area;
3.8 calculating the regression Square sum SST of each current working areaab、SSTbc、SSTcdWherein
Figure GDA0003372891680000137
iα∈[ia,ib)；
Figure GDA0003372891680000138
iβ∈[ib,ic]；
Figure GDA0003372891680000139
iγ∈(ic,id]；
Wherein
Figure GDA00033728916800001310
Respectively obtaining regression average electromagnetic force predicted values corresponding to any current working condition point in the current working area;
3.9 respectively calculating the complex correlation coefficient R of each current working area2 ab、R2 bc、R2 cdWherein
Figure GDA00033728916800001311
Figure GDA00033728916800001312
Figure GDA00033728916800001313
3.10 setting weighting coefficients of each current working area, current working area [ ia,ib)、[ib,ic]、(ic,id]The weighting coefficients of the corresponding complex correlation coefficients are respectively K1、K2、K3As shown in fig. 5;
3.11 weighting the complex correlation coefficient of the electromagnetic force in each working current region according to the current region to obtain the overall complex correlation coefficient R between the current and the electromagnetic force2The expression is R2＝K1·R2 ab+K2·R2 bc+K3·R2 cd。
The present invention is capable of other embodiments and its several details are capable of modifications in various obvious respects, all without departing from the spirit and scope of the present invention.

Claims (5)

1. A method for evaluating the linear characteristic of electromagnetic force of a proportional electromagnet based on multiple correlation coefficients is characterized by comprising the following steps:
step 1, dividing a full-working-condition plane;
step 2, obtaining the electromagnetic force of the discrete working condition points;
step 3, calculating the overall complex correlation coefficient R between the current and the electromagnetic force2；
Step 4, according to the integral complex correlation coefficient R between the current and the electromagnetic force2The linear characteristic of the electromagnetic force of the proportional electromagnet is judged according to the numerical value of the proportional electromagnet;
the dividing method of the full-working-condition plane in the step 1 comprises the following steps:
1.1, determining the working current of the proportional electromagnet and the working stroke range of the armature, wherein the working current range is marked as [ ia,id]The range of the armature stroke is denoted as [ x ]a,xd]；
1.2, equally dividing and dispersing the full working condition plane formed by the working current and the working stroke, and further obtaining corresponding discrete working condition points (i) in the full working condition planen,xm) Wherein inExpressed as the operating current range ia,id]The nth operating current value, x, corresponding to the discrete number of equal divisionsmIs a working stroke range [ xa,xd]The m-th working stroke position corresponding to the equal division and dispersion;
the method for obtaining the electromagnetic force of the discrete working condition points in the step 2 is a simulation or experiment means, and specifically comprises the following steps:
2.1, adjusting the working stroke to a certain discrete working condition point value and then keeping the working stroke unchanged;
2.2, sequentially adjusting the working current to different discrete working condition point values, and respectively measuring corresponding electromagnetic force;
2.3, adjusting the working stroke to another discrete working point value, keeping the working stroke unchanged, and repeating the process to obtain the electromagnetic force of all the discrete working points;
the overall complex correlation coefficient R between the current and the electromagnetic force in the step 32The calculating method specifically comprises the following steps:
3.1, define the working area of the primary and secondary travel, the interval [ xa,xb)、(xc,xd]For the secondary working area of the journey, interval [ xb,xc]Is the main working area of the stroke, wherein xa、xb、xc、xd∈[xa,xd]And x isa<xb<xc<xd；
3.2 setting weighting coefficients of each stroke working area, the stroke working area [ xa,xb)、[xb,xc]、(xc,xd]The weighting coefficients of the corresponding complex correlation coefficients are S1、S2、S3；
3.3, respectively calculating the average value F (i) of the electromagnetic force of the discrete working condition points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure FDA0003372891670000011
xδ∈[xa,xb)，xε∈[xb,xc]，xζ∈(xc,xd]，F(in,xδ) Is an operating current inTime travel working area [ x ]a,xb) Any working stroke position xδElectromagnetic force of F (i)n,xε) Is an operating current inTime travel working area [ x ]b,xc]Any working stroke position xεElectromagnetic force of F (i)n,xζ) Is an operating current inTime travel working area (x)c,xd]Any working stroke position xζU, v, w respectively indicate that the operating point falls in the stroke operating region [ x ]a,xb)、[xb,xc]、(xc,xd]The number of (2);
3.4, defining the working area of primary and secondary current, and the interval [ ia,ib)、(ic,id]For the current minor operating region, interval [ ib,ic]Is a current main working area, wherein ia、ib、ic、id∈[ia,id]And i isa<ib<ic<id；
3.5 using least squares or partial minimaTwo multiplication for each current operating region sample point (i)n,F(in)a) Performing linear regression to obtain three groups of current-electromagnetic force linear regression equations:
Figure FDA0003372891670000021
Figure FDA0003372891670000022
Figure FDA0003372891670000023
3.6, respectively calculating the average value of the average electromagnetic force in the same current working area under the same current working condition
Figure FDA0003372891670000024
Figure FDA0003372891670000025
Wherein
Figure FDA0003372891670000026
g represents that the operating point falls in the current operating region [ i ]a,ib) The number of (2);
Figure FDA0003372891670000027
h represents that the operating point falls in the current operating region [ i ]b,ic) The number of (2);
Figure FDA0003372891670000028
j represents that the operating point falls in the current operating region (i)c,id]The number of (2);
3.7 minutesRespectively calculating the sum of squared deviations SST of each current working areaab、SSTbc、SSTcdWherein
Figure FDA0003372891670000029
Figure FDA00033728916700000210
Figure FDA00033728916700000211
Wherein F (i)α)a、F(iβ)a、F(iγ)aThe average value of the electromagnetic force of discrete working condition points with equal working current and different working strokes in each current working area;
3.8, respectively calculating the regression square sum SST of each current working areaab、SSTbc、SSTcdWherein
Figure FDA00033728916700000212
Figure FDA00033728916700000213
Figure FDA00033728916700000214
Wherein
Figure FDA00033728916700000215
Respectively obtaining regression average electromagnetic force predicted values corresponding to any current working condition point in the current working area;
3.9 respectively calculating the complex correlation coefficient R of each current working area2 ab、R2 bc、R2 cdWherein
Figure FDA0003372891670000031
Figure FDA0003372891670000032
Figure FDA0003372891670000033
3.10 setting weighting coefficients of each current working area, current working area [ ia,ib)、[ib,ic]、(ic,id]The weighting coefficients of the corresponding complex correlation coefficients are respectively K1、K2、K3；
3.11, carrying out weighting calculation on the complex correlation coefficient under each working current according to the current region to obtain an integral complex correlation coefficient R2The expression is R2＝K1·R2 ab+K2·R2 bc+K3·R2 cd。
2. The method for evaluating the electromagnetic force linearity characteristics of proportional electromagnets based on complex correlation coefficient as claimed in claim 1, wherein the overall complex correlation coefficient R between the current and the electromagnetic force in step 3 is2The calculating method specifically comprises the following steps:
3.1, respectively calculating the average value F (i) of the electromagnetic force of the discrete working points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure FDA0003372891670000034
In the formula F (i)n,xm) Representing discrete operating points (i)n,xm) Corresponding electromagnetic force, f represents the number of equally divided working strokes;
3.2 sample point (i) by least squares or partial least squaresn,F(in)a) Linear regression is carried out to obtain a linear regression equation of current-electromagnetic force
Figure FDA0003372891670000035
Wherein
Figure FDA0003372891670000036
Is inCorresponding regression average electromagnetic force predicted values, wherein k and b are regression coefficients;
3.3, calculating the average value of the electromagnetic force under different working strokes of the same current
Figure FDA0003372891670000037
The expression is
Figure FDA0003372891670000038
Figure FDA0003372891670000039
Wherein F (i)n)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current, and e is the part for equally dividing the working current;
3.4 calculating the sum of the squares of the total deviations SST, which is expressed as
Figure FDA00033728916700000310
Wherein F (i)t)aAverage value of electromagnetic force of discrete operating points having different operating strokes for equal operating currents itRepresenting the operating current range [ i ]a,id]The t-th working current value;
3.5 calculating regression SquareAnd SSR, expressed as
Figure FDA00033728916700000311
Wherein
Figure FDA00033728916700000312
The regression average electromagnetic force predicted value corresponding to the discrete operating point with the same working current and different working strokes;
3.6, calculating to obtain an integral complex correlation coefficient R2The expression is
Figure FDA00033728916700000313
3. The method for evaluating the electromagnetic force linearity characteristics of proportional electromagnets based on complex correlation coefficient as claimed in claim 1, wherein the overall complex correlation coefficient R between the current and the electromagnetic force in step 3 is2The calculating method specifically comprises the following steps:
3.1, respectively calculating the average value F (i) of the electromagnetic force of the discrete working points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure FDA0003372891670000041
In the formula F (i)n,xm) Representing discrete operating points (i)n,xm) Corresponding electromagnetic force, xm∈[xa,xd]F represents the number of equally divided working strokes;
3.2, defining the working area of primary and secondary current, and the interval [ ia,ib)、(ic,id]For the current minor operating region, interval [ ib,ic]Is a current main working area, wherein ia、ib、ic、id∈[ia,id]And i isa<ib<ic<id；
3.3, miningUsing least square method or partial least square method to sample point (i) of each current working arean,F(in)a) Performing linear regression to obtain three groups of current-electromagnetic force linear regression equations:
Figure FDA0003372891670000042
Figure FDA0003372891670000043
Figure FDA0003372891670000044
wherein
Figure FDA0003372891670000045
Respectively i in different current working regionsα、iβ、iγCorresponding regression average electromagnetic force predicted value, wherein iαIs a current working region [ ia,ib) Of the alpha operating current value, iβIs a current working region [ ib,ic]Of (d) the beta-th operating current value, iγIs a current working region (i)c,id]Of (1) the gamma-th operating current value, k1、b1、k2、b2、k3、b3Is a regression coefficient;
3.4, respectively calculating the average value of the average electromagnetic force in the same current working area under the same current working condition
Figure FDA0003372891670000046
Figure FDA0003372891670000047
Wherein
Figure FDA0003372891670000048
g represents that the operating point falls in the current operating region [ i ]a,ib) The number of (2);
Figure FDA0003372891670000049
h represents that the operating point falls in the current operating region [ i ]b,ic) The number of (2);
Figure FDA00033728916700000410
j represents that the operating point falls in the current operating region (i)c,id]The number of (2);
3.5, respectively calculating the sum of squared deviations SST of each current working areaab、SSTbc、SSTcdWherein
Figure FDA00033728916700000411
Figure FDA00033728916700000412
Figure FDA00033728916700000413
Wherein F (i)α)a、F(iβ)a、F(iγ)aThe average value of the electromagnetic force of discrete working condition points with equal working current and different working strokes in each current working area;
3.6, calculating the regression square sum SST of each current working area respectivelyab、SSTbc、SSTcdWherein
Figure FDA0003372891670000051
Figure FDA0003372891670000052
Figure FDA0003372891670000053
Wherein
Figure FDA0003372891670000054
Respectively obtaining regression average electromagnetic force predicted values corresponding to any current working condition point in the current working area;
3.7 respectively calculating the complex correlation coefficient R of each current working area2 ab、R2 bc、R2 cdWherein
Figure FDA0003372891670000055
Figure FDA0003372891670000056
Figure FDA0003372891670000057
3.8 setting weighting coefficients of each current working area, ia,ib)、[ib,ic]、(ic,id]The weighting coefficients of the corresponding complex correlation coefficients are respectively K1、K2、K3；
3.9, carrying out weighted calculation on the complex correlation coefficient of each working current region according to the current region to obtain the integral complex correlation coefficient R2The expression is R2＝K1·R2 ab+K2·R2 bc+K3·R2 cd。
4. The method for evaluating the electromagnetic force linearity characteristics of proportional electromagnets based on complex correlation coefficient as claimed in claim 1, wherein the overall complex correlation coefficient R between the current and the electromagnetic force in step 3 is2The calculating method specifically comprises the following steps:
3.1, define the working area of the primary and secondary travel, the interval [ xa,xb)、(xc,xd]For the secondary working area of the journey, interval [ xb,xc]Is the main working area of the stroke, wherein xa、xb、xc、xd∈[xa,xd]And x isa<xb<xc<xd；
3.2 setting weighting coefficients of each stroke working area, the stroke working area [ xa,xb)、[xb,xc]、(xc,xd]The weighting coefficients of the corresponding complex correlation coefficients are S1、S2、S3；
3.3, respectively calculating the average value F (i) of the electromagnetic force of the discrete working condition points with the same working current and different working strokesn)aObtaining a series of current and electromagnetic force linear regression sample points (i)n,F(in)a) Wherein
Figure FDA0003372891670000058
xδ∈[xa,xb)，xε∈[xb,xc]，xζ∈(xc,xd]，F(in,xδ) Is an operating current inTime travel working area [ x ]a,xb) Any working stroke position xδElectromagnetic force of F (i)n,xε) Is an operating current inTime travel working area [ x ]b,xc]Any working stroke position xεElectromagnetic force of F (i)nX ζ) is the operating current inTime travel working area (x)c,xd]Any working stroke position xζU, v, w respectively indicate that the operating point falls in the stroke operating region [ x ]a,xb)、[xb,xc]、(xc,xd]The number of (2);
3.4 sample point (i) by least squares or partial least squaresn,F(in)a) Linear regression is carried out to obtain a linear regression equation of current-electromagnetic force
Figure FDA0003372891670000061
Wherein
Figure FDA0003372891670000062
Is inCorresponding regression average electromagnetic force predicted values, wherein k and b are regression coefficients;
3.5, calculating the average value of the electromagnetic force under different working strokes of the same current
Figure FDA0003372891670000063
The expression is
Figure FDA0003372891670000064
Figure FDA0003372891670000065
F(in)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current, and e is the part for equally dividing the working current;
3.6 calculating the sum of squared deviations SST, which is expressed as
Figure FDA0003372891670000066
Wherein F (i)t)aThe average value of the electromagnetic force of discrete working condition points with different working strokes for equal working current;
3.7 calculating regression Square sum SSR with the expression of
Figure FDA0003372891670000067
Wherein
Figure FDA0003372891670000068
The regression average electromagnetic force predicted value corresponding to the discrete operating point with the same working current and different working strokes;
3.8, calculating to obtain an integral complex correlation coefficient R2The expression is
Figure FDA0003372891670000069
5. The method for evaluating the electromagnetic force linearity characteristics of a proportional electromagnet based on complex correlation coefficient as claimed in claim 1, wherein the step 4 is performed according to the overall complex correlation coefficient R between the current and the electromagnetic force2The method for judging the linear characteristic of the electromagnetic force of the proportional electromagnet comprises the following steps: integral complex correlation coefficient R between current and electromagnetic force2The closer to 1, the better the linear characteristic of the electromagnetic force of the proportional electromagnet is; integral complex correlation coefficient R between current and electromagnetic force2The closer to 0, the worse the linear characteristic of the proportional electromagnet electromagnetic force.
CN202110279685.8A 2021-03-16 2021-03-16 Proportional electromagnet electromagnetic force linear characteristic evaluation method based on multiple correlation coefficients Active CN112884355B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110279685.8A CN112884355B (en) 2021-03-16 2021-03-16 Proportional electromagnet electromagnetic force linear characteristic evaluation method based on multiple correlation coefficients

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110279685.8A CN112884355B (en) 2021-03-16 2021-03-16 Proportional electromagnet electromagnetic force linear characteristic evaluation method based on multiple correlation coefficients

Publications (2)

Publication Number Publication Date
CN112884355A CN112884355A (en) 2021-06-01
CN112884355B true CN112884355B (en) 2021-12-31

Family

ID=76042416

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110279685.8A Active CN112884355B (en) 2021-03-16 2021-03-16 Proportional electromagnet electromagnetic force linear characteristic evaluation method based on multiple correlation coefficients

Country Status (1)

Country Link
CN (1) CN112884355B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012105469A (en) * 2010-11-11 2012-05-31 Panasonic Corp Method for estimating magnetic pole position of motor
CN110929437A (en) * 2019-10-28 2020-03-27 温州大学 Moving-iron type proportional electromagnet constant force prediction method based on response surface

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010113250A1 (en) * 2009-03-31 2010-10-07 トヨタ自動車株式会社 Joint quality examining device and joint quality examining method
CN106095646A (en) * 2016-06-27 2016-11-09 江苏迪纳数字科技股份有限公司 Hadoop performance cluster computational methods based on multiple linear regression model
CN108805343A (en) * 2018-05-29 2018-11-13 祝恩元 A kind of Scientech Service Development horizontal forecast method based on multiple linear regression
CN109992857A (en) * 2019-03-19 2019-07-09 浙江大学 Automobile tail light shape design evaluation and prediction technique

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012105469A (en) * 2010-11-11 2012-05-31 Panasonic Corp Method for estimating magnetic pole position of motor
CN110929437A (en) * 2019-10-28 2020-03-27 温州大学 Moving-iron type proportional electromagnet constant force prediction method based on response surface

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
高速电磁阀电磁力全工况关键参数相关性分析;刘鹏等;《农业工程学报》;20150331;第89-96页 *
高速电磁阀电磁力近似模型的构建与分析;刘鹏等;《农业工程学报》;20150830;第96-102页 *

Also Published As

Publication number Publication date
CN112884355A (en) 2021-06-01

Similar Documents

Publication Publication Date Title
Brünig et al. Damage and failure at negative stress triaxialities: Experiments, modeling and numerical simulations
DE10259609A1 (en) Method of determining fatigue sensitivity
Ali-Benyahia et al. Improvement of nondestructive assessment of on-site concrete strength: Influence of the selection process of cores location on the assessment quality for single and combined NDT techniques
Alayande et al. An overview and application of discriminant analysis in data analysis
CN108596369B (en) Space load prediction method based on multivariate model and blind number theory
CN103106813A (en) Traffic situation complexity assessment method based on index system
CN109632538A (en) Probabilistic Fatigue crack growth rate statistical analysis technique based on matched curve equivalency transform
CN112651119B (en) Multi-performance parameter acceleration degradation test evaluation method for space harmonic reducer
CN112884355B (en) Proportional electromagnet electromagnetic force linear characteristic evaluation method based on multiple correlation coefficients
US20230086774A1 (en) Method and system for predicting biological age on basis of various omics data analyses
Thanwane et al. The effect of measurement errors on the performance of the homogenously weighted moving average X¯ monitoring scheme with estimated parameters
Abbasi et al. Efficient CV control charts based on ranked set sampling
Cudney et al. Applying the Mahalanobis-Taguchi system to vehicle ride
CN113051733B (en) Method for evaluating electromagnetic force horizontal characteristics of proportional electromagnet
Bührer-Sekula et al. Simple dipstick assay for semi-quantitative detection of neopterin in sera
CN110489848B (en) Method for predicting corrosion fatigue crack propagation rates at different seawater flow velocities
CN107391962A (en) The method of gene or site to disease regulation relationship is analysed based on multigroup credit
CN109670278A (en) A kind of Probabilistic Fatigue crack growth rate statistical analysis technique based on Gaussian Profile
Petersen et al. A model for setting analytical quality specifications and design of control for measurements on the ordinal scale
CN102129515A (en) Method for evaluating fragility of water resource under climatic change
Zeng et al. Robust Phase I monitoring of profile data with application in low-E glass manufacturing processes
CN115169707A (en) Equipment energy consumption prediction method and device based on multiple linear regression
CN104933256A (en) Comprehensive parametric evaluation method for influence of ambient atmosphere condition on air quality
CN114938947A (en) Multi-level vaginal microecological assessment system construction method and device
Vineis et al. Technical variability in laboratory data

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