CN111784193A - Product performance consistency inspection method based on normal distribution - Google Patents

Product performance consistency inspection method based on normal distribution Download PDF

Info

Publication number
CN111784193A
CN111784193A CN202010689731.7A CN202010689731A CN111784193A CN 111784193 A CN111784193 A CN 111784193A CN 202010689731 A CN202010689731 A CN 202010689731A CN 111784193 A CN111784193 A CN 111784193A
Authority
CN
China
Prior art keywords
normal distribution
groups
samples
probability density
density function
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.)
Granted
Application number
CN202010689731.7A
Other languages
Chinese (zh)
Other versions
CN111784193B (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.)
National University of Defense Technology
Original Assignee
National University of Defense 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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202010689731.7A priority Critical patent/CN111784193B/en
Publication of CN111784193A publication Critical patent/CN111784193A/en
Application granted granted Critical
Publication of CN111784193B publication Critical patent/CN111784193B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00Computing arrangements based on specific mathematical models
    • G06N7/01Probabilistic graphical models, e.g. probabilistic networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/04Manufacturing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Abstract

The invention provides a product performance consistency testing method based on normal distribution, which comprises the following steps: acquiring performance parameters of a plurality of products; selecting two groups of measurement data samples from the performance parameters, wherein the two groups of samples are subjected to normal distribution; respectively estimating normal distribution parameters of the two groups of measurement data samples according to the two groups of samples; further obtaining the intersection point of the normal distribution probability density function curves of the two groups of samples; then obtaining the area of the overlapping part of the probability density function of the two groups of samples and the area surrounded by the abscissa axis as consistency measurement; finally, the two groups of samples are compared with a given significance level or confidence level to determine whether the two groups of samples have consistency at the given confidence level. The consistency test method is different from a classical consistency test method, does not need to construct statistics, defines consistency measurement by starting from the concept and the mathematical meaning of a probability density function, and avoids the difficulty of constructing sufficient statistics.

Description

Product performance consistency inspection method based on normal distribution
Technical Field
The invention relates to the technical field of product performance consistency inspection methods, in particular to a normal distribution-based product performance consistency inspection method.
Background
The normal distribution is a distribution form most commonly used in the engineering field, and occurs in many fields, such as thermal noise current in semiconductor devices, weight of bagged cement, filling amount of liquid rocket propellant, temperature in constant-temperature storage containers in industrial production, and the like, and measurement data thereof are subject to the normal distribution. When the physical quantities are measured in multiple batches, two or more groups of different data are obtained, and whether the two or more groups of data are subjected to the same normal distribution needs to be judged so as to analyze the consistency of the measured data and judge the stability of the production process and technical means, such as the stability of temperature control in a constant-temperature storage container, the stability of a bagged cement production line and the like.
In engineering practice, two different sets of measurement data which are subject to normal distribution are obtained, such as the magnitude of thermal noise current in two batches of semiconductor devices, the weight of two batches of bagged cement produced in different time periods, the temperature in a constant-temperature storage container under two power-on starts and the like, theoretical analysis and actual measurement data show that the two sets of measurement data are both subject to normal distribution, but whether the two sets of data are subject to the same normal distribution, namely whether normal distribution parameters are the same or not, further analysis is needed, and an analysis conclusion can be used for judging the stability of the production process. For two sets of samples that obey a normal distribution (referred to as measured data in the aforementioned practice), a consistency test is performed on the two sets of samples to determine whether the two sets of samples obey the same distribution at a given confidence level or significance level. The normal distribution consistency test method adopts a classical mean-variance test method, firstly, the F distribution is used for carrying out consistency test on the variances of two groups of data, if the variance of two groups of data passes the test, the t distribution is continuously used for carrying out the consistency test on the means of two groups of data, if the variance of two groups of data passes the test, the two groups of data are considered to be obeyed by the same normal distribution, and if any parameter of the variance and the mean does not pass the consistency test, the two groups of normal distribution data are considered to be not obeyed by the same normal distribution. This method requires the construction of test statistics that are subject to F-distribution and t-distribution, respectively, to be implemented.
Disclosure of Invention
The invention aims to provide a method for testing the consistency of product performance based on normal distribution, which solves the problems.
The invention provides a product performance consistency testing method based on normal distribution, which comprises the following steps:
acquiring performance parameters of a plurality of products; the performance parameters are extracted from different batches of products or extracted from different moments of the products;
selecting two groups of measurement data samples from the performance parameters, wherein the two groups of measurement data samples are subjected to normal distribution;
respectively estimating the mean value and the mean square error of normal distribution of two groups of measurement data samples by utilizing a maximum likelihood estimation method according to the normal distribution probability density function;
obtaining the intersection point of the normal distribution probability density function curves of the two groups of measurement data samples according to the mean value and the mean square error of the normal distribution of the two groups of measurement data samples;
according to the intersection point of the probability density function curves of the two groups of measurement data samples, the area of the overlapping part of the normal distribution probability density function curves of the two groups of measurement data samples and the region surrounded by the abscissa axis is obtained;
and judging whether the two groups of measurement data samples have consistency or not according to the area of the overlapping part of the normal distribution probability density function curves of the two groups of measurement data samples and the area surrounded by the abscissa axis and a preset significance level or confidence coefficient.
Further, the estimating the mean and mean square error of the normal distribution of the two groups of measurement data samples by using a maximum likelihood estimation method according to the normal distribution probability density function includes:
the normal distribution probability density function is:
Figure BDA0002588894900000021
wherein mu is the mean value of normal distribution, and sigma is the mean square error of normal distribution;
according to the maximum likelihood estimation method, the estimated values of two groups of data normal distribution parameters are respectively
Figure BDA0002588894900000022
Figure BDA0002588894900000023
Wherein n is the number of samples in the sample set X, and m is the number of samples in the sample set Y;
then the normal distribution probability density function obeyed by the two groups of data is respectively:
Figure BDA0002588894900000024
Figure BDA0002588894900000025
wherein the content of the first and second substances,
Figure BDA00025888949000000310
is the mean of a normal distribution under the first set of samples,
Figure BDA00025888949000000311
is the mean square error of a normal distribution under the first set of samples,
Figure BDA00025888949000000312
is the mean of a normal distribution for the second set of samples,
Figure BDA00025888949000000313
is the mean square error of the normal distribution for the second set of samples.
Further, the obtaining an intersection point of the normal distribution probability density function curves of the two groups of measurement data samples according to the normal distribution probability density functions of the two groups of measurement data samples includes:
making the two probability density functions equal, i.e. fx(r)=fy(r) then
Figure BDA0002588894900000031
Wherein the first set of measurement data samples is X ═ { X1,x2,…,xnY ═ Y in the second set of measurement data samples1,y2,…,ymN is the number of the first group of samples, and m is the number of the second group of samples; f. ofx(r) probability density function corresponding to the first set of samples, fy(r) represents a probability density function corresponding to the second set of samples, r representing a variable;
Figure BDA0002588894900000032
namely, it is
Figure BDA0002588894900000033
An equation can be obtained
Figure BDA0002588894900000034
When in use
Figure BDA0002588894900000035
When there is a crossing point
Figure BDA0002588894900000036
When in use
Figure BDA0002588894900000037
When, there are two intersections, respectively
Figure BDA0002588894900000038
Wherein
Figure BDA0002588894900000039
Further, according to the intersection point of the probability density function curves of the two groups of measurement data samples,obtaining the area of the overlapping part of the normal distribution probability density function curves of the two groups of measurement data samples and the region surrounded by the abscissa axis, wherein the area comprises the following steps: when in use
Figure BDA0002588894900000041
When there is a crossing point
Figure BDA0002588894900000042
If it is not
Figure BDA0002588894900000043
The calculation formula of the area of the overlapped part is
Figure BDA0002588894900000044
If it is not
Figure BDA0002588894900000045
The calculation formula of the area of the overlapped part is
Figure BDA0002588894900000046
If it is not
Figure BDA0002588894900000047
The two probability density function curves are completely overlapped, and the area of the overlapped part
cr=1
Wherein c isrIndicating the size of the consistency.
Further, the obtaining of the area of the overlapping portion of the normal distribution probability density function curves of the two groups of measurement data samples and the region surrounded by the abscissa axis according to the intersection point of the probability density function curves of the two groups of measurement data samples further includes:
when in use
Figure BDA0002588894900000048
When, there are two intersections, respectively
Figure BDA0002588894900000049
Order to
Figure BDA00025888949000000410
The calculation of the area of the overlapped part can be divided into three parts
When in use
Figure BDA00025888949000000411
When, let f1(r)=min(fx(r),fy(r)) to obtain
Figure BDA00025888949000000412
When in use
Figure BDA00025888949000000413
When, let f2(r)=min(fx(r),fy(r)) to obtain
Figure BDA00025888949000000414
When in use
Figure BDA00025888949000000415
When, let f3(r)=min(fx(r),fy(r)) to obtain
Figure BDA00025888949000000416
The area of the overlapped part of the two normally distributed probability density function curves and the area surrounded by the abscissa axis is
cr=c3+c2+c1
Further, the determining whether the two groups of measurement data samples have consistency according to the area of the overlapping portion of the normal distribution probability density function curves of the two groups of measurement data samples and the area surrounded by the abscissa axis and a preset significance level or confidence degree includes:
if the preset significance level is alpha, the judgment rule of whether the two groups of samples are consistent is
If c isr1- α, the two groups of samples were considered consistent at a significance level of α;
if c isr< 1- α, the two groups of samples were considered to be inconsistent at a significance level of α.
The invention has the technical effects that:
(1) the invention provides a product performance consistency test method based on normal distribution, which utilizes the area of the overlapping part of two groups of sample experience probability density functions and the area surrounded by the abscissa axis as the consistency measurement of two samples and provides a corresponding calculation method and a judgment criterion. Since the probability density function and the cumulative distribution function are two basic concepts in statistics, the probability density function has an integral of 1 throughout the domain of the random variable. Under the condition of a given random variable interval, the area enclosed by the probability density function and the interval is the cumulative distribution function under the interval, namely the probability. According to the concept and the meaning of the probability density function, if the overlapping area of the two normal distribution probability density function curves and the area surrounded by the abscissa axis is larger, the approximation degree of the two probability density functions is larger, and the consistency of the two groups of corresponding samples is better. The innovation point of the method is that the area of the overlapping part of the two groups of normal distribution sample empirical probability density function curves and the area surrounded by the abscissa axis is defined as the measure of the consistency of the two groups of normal distribution samples, a new method for judging the consistency of the two groups of normal distribution samples is provided, the calculation process of the consistency is provided, the basic problem in the consistency test of the two groups of normal distribution samples is solved, and a feasible method is provided for the consistency test of the normal distribution.
(2) The product performance consistency testing method based on normal distribution provided by the invention has the advantages of clear mathematical concept, clear calculation steps, reasonability and feasibility. Moreover, the product performance consistency test method based on normal distribution provided by the invention does not need to construct sample-based statistics and even know the distribution form of the statistics, but starts with the concept and mathematical meaning of the probability density function, defines consistency measurement, and provides a feasible calculation method for solving the consistency test problem of two groups of normal distribution samples.
The above and other aspects of the present invention will become apparent from the following description, which refers in particular to various embodiments of the method for checking the consistency of product properties based on normal distribution according to the present invention.
Drawings
FIG. 1 is a schematic flow diagram of a preferred embodiment method of the present invention;
fig. 2 is a schematic area diagram of an overlapping portion of two normal distribution probability density function curves and a region surrounded by an abscissa axis, and an intersection point of the two probability density function curves is given, and a diagonally shaded portion in the diagram is the overlapping portion.
Detailed Description
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the invention and, together with the description, serve to explain the invention and not to limit the invention.
Referring to fig. 1, the method for checking the consistency of product performance based on normal distribution includes the following steps:
step one, acquiring performance parameters of a plurality of products; the performance parameters are extracted from different batches of products or extracted from different moments of the products; selecting two groups of measurement data samples from the performance parameters, wherein the two groups of measurement data samples are subjected to normal distribution;
secondly, respectively estimating the mean value and the mean square error of normal distribution of two groups of measurement data samples by using a maximum likelihood estimation method according to the probability density function of the normal distribution;
step three, obtaining the intersection point of the normal distribution probability density function curves of the two groups of measurement data samples according to the mean value and the mean square error of the normal distribution of the two groups of measurement data samples;
obtaining the area (shown in figure 2) of the overlapping part of the normal distribution probability density function curves of the two groups of measurement data samples and the region surrounded by the abscissa axis according to the intersection point of the probability density function curves of the two groups of measurement data samples, and taking the area as the measure of the consistency of the two groups of samples;
and fifthly, judging whether the two groups of measurement data samples have consistency or not according to the area of the overlapping part of the normal distribution probability density function curves of the two groups of measurement data samples and the area surrounded by the abscissa axis and a preset significance level or confidence coefficient.
Since the probability density function and the cumulative distribution function are two basic concepts in statistics, the probability density function has an integral of 1 throughout the domain of the random variable. Under the condition of a given random variable interval, the area enclosed by the probability density function and the interval is the cumulative distribution function under the interval, namely the probability. According to the concept and the meaning of the probability density function, if the overlapping area of the two normal distribution probability density function curves and the area surrounded by the abscissa axis is larger, the approximation degree of the two probability density functions is larger, and the consistency of the two groups of corresponding samples is better. The innovation point of the method is that the area of the overlapping part of the two groups of normal distribution sample empirical probability density function curves and the area surrounded by the abscissa axis is defined as the measure of the consistency of the two groups of normal distribution samples, a new method for judging the consistency of the two groups of normal distribution samples is provided, the calculation process of the consistency is provided, the basic problem in the consistency test of the two groups of normal distribution samples is solved, and a feasible method is provided for the consistency test of the normal distribution.
Specifically, the second step further comprises:
the normal distribution probability density function is:
Figure BDA0002588894900000071
wherein mu is the mean value of normal distribution, and sigma is the mean square error of normal distribution;
according to the maximum likelihood estimation method, the estimated values of two groups of data normal distribution parameters are respectively
Figure BDA0002588894900000072
Figure BDA0002588894900000073
Wherein n is the number of samples in the sample set X, and m is the number of samples in the sample set Y;
then the normal distribution probability density function obeyed by the two groups of data is respectively:
Figure BDA0002588894900000074
Figure BDA0002588894900000075
wherein the content of the first and second substances,
Figure BDA0002588894900000076
is the mean of a normal distribution under the first set of samples
Figure BDA0002588894900000077
Figure BDA0002588894900000078
Is the mean square error of a normal distribution under the first set of samples,
Figure BDA0002588894900000079
is the mean of a normal distribution under the first set of samples
Figure BDA00025888949000000710
Figure BDA00025888949000000711
Is the mean square error of the normal distribution under the first set of samples.
Specifically, the third step further comprises:
making the two probability density functions equal, i.e. fx(r)=fy(r) then
Figure BDA00025888949000000712
Wherein the first set of measurement data samples is X ═ { X1,x2,…,xnY ═ Y in the second set of measurement data samples1,y2,…,ymN is the number of the first group of samples, and m is the number of the second group of samples; f. ofx(r) probability density function corresponding to the first set of samples, fy(r) represents a probability density function corresponding to the second set of samples, r representing a variable;
Figure BDA00025888949000000713
namely, it is
Figure BDA00025888949000000714
An equation can be obtained
Figure BDA0002588894900000081
When in use
Figure BDA0002588894900000082
When there is a crossing point
Figure BDA0002588894900000083
When in use
Figure BDA0002588894900000084
When, there are two intersections, respectively
Figure BDA0002588894900000085
Wherein
Figure BDA0002588894900000086
Specifically, the fourth step further comprises:
when in use
Figure BDA0002588894900000087
When there is a crossing point
Figure BDA0002588894900000088
If it is not
Figure BDA0002588894900000089
The calculation formula of the area of the overlapped part is
Figure BDA00025888949000000810
If it is not
Figure BDA00025888949000000811
The calculation formula of the area of the overlapped part is
Figure BDA00025888949000000812
If it is not
Figure BDA00025888949000000813
The two probability density function curves are completely overlapped, and the area of the overlapped part
cr=1
Wherein c isrIndicating the size of the consistency.
When in use
Figure BDA00025888949000000814
When, there are two intersections, respectively
Figure BDA00025888949000000815
Order to
Figure BDA00025888949000000816
The calculation of the area of the overlapped part can be divided into three parts
When in use
Figure BDA00025888949000000817
When, let f1(r)=min(fx(r),fy(r)) to obtain
Figure BDA00025888949000000818
When in use
Figure BDA00025888949000000819
When, let f2(r)=min(fx(r),fy(r)) to obtain
Figure BDA0002588894900000091
When in use
Figure BDA0002588894900000092
When, let f3(r)=min(fx(r),fy(r)) to obtain
Figure BDA0002588894900000093
The area of the overlapped part of the two normally distributed probability density function curves and the area surrounded by the abscissa axis is
cr=c3+c2+c1
Specifically, the fifth step further comprises:
if the preset significance level is alpha, the judgment rule of whether the two groups of samples are consistent is
If c isr1- α, the two groups of samples were considered consistent at a significance level of α;
if c isr< 1- α, the two groups of samples were considered to be inconsistent at a significance level of α.
In order to better explain the technical scheme provided by the invention, the following description is made in conjunction with specific examples.
(1) Assuming that two different batches of semiconductor devices are produced in a factory, the thermal noise current of the two batches of devices is measured, and two sets of measurement data samples are obtained, wherein X is ═ { X ═ respectively1,x2,…,xn},Y={y1,y2,…,ymAnd the two sets of measurement data are subjected to consistency test to judge whether the thermal noise performance of the semiconductor products in two batches is the same or not. Where the data set X has a total of 18 samples, i.e. n 18. The data set Y has 15 samples, i.e., m 15.
X={4.61 -4.63 0.19 -0.22 1.46 1.44 -2.09 0.41 0.01 2.38 3.78 3.83 -2.09 0.73 -3.14 -2.84 0.48 5.10},
Y={-2.08 2.49 0.10 5.47 -3.36 1.13 3.21 5.40 7.18 1.34 -4.97 -1.97 -3.25 10.40 -1.46}。
(2) According to the normal distribution probability density function, respectively estimating normal distribution parameter parameters under two groups of data by utilizing a maximum likelihood estimation method
Figure BDA00025888949000000910
And
Figure BDA00025888949000000911
. For the first set of data X, the distribution parameter estimate is
Figure BDA0002588894900000094
For the second set of data Y, the distribution parameter estimate is
Figure BDA0002588894900000095
Wherein n is the number of samples in the sample set X, and m is the number of samples in the sample set Y.
Then the normal distribution probability density function obeyed by the two groups of data is respectively:
Figure BDA0002588894900000101
Figure BDA0002588894900000102
(3) after the normal distribution probability density functions of the two groups of data are obtained, the intersection point of the two probability density function curves is calculated. Let fx(r)=fy(r) due to
Figure BDA0002588894900000103
The equation can be obtained
Figure BDA0002588894900000104
Can obtain
Figure BDA0002588894900000105
Figure BDA0002588894900000106
Then
Figure BDA0002588894900000107
(4) And calculating the area of the overlapped part of the two normal distribution probability density function curves and the region surrounded by the abscissa axis according to the intersection point characteristics of the two probability density function curves, and taking the area as the measure of the consistency of the two groups of samples. Order to
Figure BDA0002588894900000108
Due to the fact that
Figure BDA0002588894900000109
The calculation of the area of the overlapped part can be divided into three parts
1) When in use
Figure BDA00025888949000001010
When, let f1(r)=min(fx(r),fy(r))=fx(r) obtaining
Figure BDA00025888949000001011
2) When in use
Figure BDA00025888949000001012
When, let f2(r)=min(fx(r),fy(r))=fy(r) obtaining
Figure BDA00025888949000001013
3) When in use
Figure BDA00025888949000001014
When, let f3(r)=min(fx(r),fy(r))=fx(r) obtaining
Figure BDA00025888949000001015
The area of the overlapped part of the two normally distributed probability density function curves and the area surrounded by the abscissa axis is
cr=c3+c2+c1=0.7688
(5) Judging whether the two groups of data samples have consistency or not according to a given significance level or confidence;
assuming a significance level of α of 0.2 and a confidence level of 1- α of 0.8, c is assignedr< 1- α ═ 0.8, i.e., the two sets of samples X and Y were considered to have no agreement at significance level 0.2.
The product performance consistency testing method based on normal distribution provided by the invention has clear mathematical concept, clear calculation steps, and is reasonable and feasible. Moreover, the product performance consistency test method based on normal distribution provided by the invention does not need to construct sample-based statistics and even know the distribution form of the statistics, but starts with the concept and mathematical meaning of the probability density function, defines consistency measurement, and provides a feasible calculation method for solving the consistency test problem of two groups of normal distribution samples.
It will be clear to a person skilled in the art that the scope of the present invention is not limited to the examples discussed in the foregoing, but that several amendments and modifications thereof are possible without deviating from the scope of the present invention as defined in the attached claims. While the invention has been illustrated and described in detail in the drawings and the description, such illustration and description are to be considered illustrative or exemplary and not restrictive. The invention is not limited to the disclosed embodiments.
Variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims. In the claims, the term "comprising" does not exclude other steps or elements, and the indefinite article "a" or "an" does not exclude a plurality. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measures cannot be used to advantage. Any reference signs in the claims shall not be construed as limiting the scope of the invention.

Claims (6)

1. The method for checking the consistency of product performance based on normal distribution is characterized by comprising the following steps:
acquiring performance parameters of a plurality of products; the performance parameters are extracted from different batches of products or extracted from different moments of the products;
selecting two groups of measurement data samples from the performance parameters, wherein the two groups of measurement data samples are subjected to normal distribution;
respectively estimating the mean value and the mean square error of normal distribution of two groups of measurement data samples by utilizing a maximum likelihood estimation method according to the normal distribution probability density function;
obtaining the intersection point of the normal distribution probability density function curves of the two groups of measurement data samples according to the mean value and the mean square error of the normal distribution of the two groups of measurement data samples;
according to the intersection point of the probability density function curves of the two groups of measurement data samples, the area of the overlapping part of the normal distribution probability density function curves of the two groups of measurement data samples and the region surrounded by the abscissa axis is obtained;
and judging whether the two groups of measurement data samples have consistency or not according to the area of the overlapping part of the normal distribution probability density function curves of the two groups of measurement data samples and the area surrounded by the abscissa axis and a preset significance level or confidence coefficient.
2. The method for checking the product performance consistency based on normal distribution according to claim 1, wherein the step of respectively estimating the mean value and the mean square error of the normal distribution of two groups of measurement data samples by using a maximum likelihood estimation method according to the probability density function of the normal distribution comprises the following steps:
the normal distribution probability density function is:
Figure FDA0002588894890000011
wherein mu is the mean value of normal distribution, and sigma is the mean square error of normal distribution;
according to the maximum likelihood estimation method, the estimated values of two groups of data normal distribution parameters are respectively
Figure FDA0002588894890000012
Figure FDA0002588894890000013
Wherein n is the number of samples in the sample set X, and m is the number of samples in the sample set Y;
then the normal distribution probability density function obeyed by the two groups of data is respectively:
Figure FDA0002588894890000021
Figure FDA0002588894890000022
wherein the content of the first and second substances,
Figure FDA0002588894890000023
is the mean of a normal distribution under the first set of samples,
Figure FDA0002588894890000024
is the mean square error of a normal distribution under the first set of samples,
Figure FDA0002588894890000025
is the mean of a normal distribution for the second set of samples,
Figure FDA0002588894890000026
is the mean square error of the normal distribution for the second set of samples.
3. The method for checking the product performance consistency based on normal distribution according to claim 1, wherein obtaining the intersection point of the normal distribution probability density function curves of the two groups of measurement data samples according to the normal distribution probability density functions of the two groups of measurement data samples comprises:
making the two probability density functions equal, i.e. fx(r)=fy(r) then
Figure FDA0002588894890000027
Wherein the first set of measurement data samples is X ═ { X1,x2,…,xnY ═ Y in the second set of measurement data samples1,y2,…,ymN is the number of the first group of samples, and m is the number of the second group of samples; f. ofx(r) probability density function corresponding to the first set of samples, fy(r) represents a probability density function corresponding to the second set of samples, r representing a variable;
Figure FDA0002588894890000028
namely, it is
Figure FDA0002588894890000029
An equation can be obtained
Figure FDA00025888948900000210
When in use
Figure FDA00025888948900000211
When there is a crossing point
Figure FDA00025888948900000212
When in use
Figure FDA00025888948900000213
When, there are two intersections, respectively
Figure FDA00025888948900000214
Wherein
Figure FDA0002588894890000031
4. The method for checking the product performance consistency based on normal distribution according to claim 3, wherein the step of obtaining the area of the overlapping part of the normal distribution probability density function curves of the two groups of measurement data samples and the region surrounded by the abscissa axis according to the intersection point of the probability density function curves of the two groups of measurement data samples comprises the following steps:
when in use
Figure FDA0002588894890000032
When there is a crossing point
Figure FDA0002588894890000033
If it is not
Figure FDA0002588894890000034
The calculation formula of the area of the overlapped part is
Figure FDA0002588894890000035
If it is not
Figure FDA0002588894890000036
The calculation formula of the area of the overlapped part is
Figure FDA0002588894890000037
If it is not
Figure FDA0002588894890000038
The two probability density function curves are completely overlapped, and the area of the overlapped part
cr=1
Wherein c isrIndicating the size of the consistency.
5. The method for checking product performance consistency according to claim 3, wherein the area of the overlapping portion of the normal distribution probability density function curves of the two groups of measurement data samples and the region surrounded by the abscissa axis is obtained according to the intersection point of the probability density function curves of the two groups of measurement data samples, and the method further comprises:
when in use
Figure FDA0002588894890000039
When, there are two intersections, respectively
Figure FDA00025888948900000310
Order to
Figure FDA00025888948900000311
The calculation of the area of the overlapped part can be divided into three parts
When in use
Figure FDA00025888948900000312
When, let f1(r)=min(fx(r),fy(r)) to obtain
Figure FDA00025888948900000313
When in use
Figure FDA00025888948900000314
When, let f2(r)=min(fx(r),fy(r)) to obtain
Figure FDA00025888948900000315
When in use
Figure FDA0002588894890000041
When, let f3(r)=min(fx(r),fy(r)) to obtain
Figure FDA0002588894890000042
The area of the overlapped part of the two normally distributed probability density function curves and the area surrounded by the abscissa axis is
cr=c3+c2+c1
6. The method for checking the product performance consistency based on normal distribution according to claim 1, wherein the step of judging whether the two groups of measurement data samples have consistency according to the area of the overlapping part of the normal distribution probability density function curve of the two groups of measurement data samples and the area surrounded by the abscissa axis and a preset significance level or confidence level comprises the following steps:
if the preset significance level is alpha, the judgment rule of whether the two groups of samples are consistent is
If c isr1- α, the two groups of samples were considered consistent at a significance level of α;
if c isr< 1- α, the two groups of samples were considered to be inconsistent at a significance level of α.
CN202010689731.7A 2020-07-17 2020-07-17 Product performance consistency inspection method based on normal distribution Active CN111784193B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010689731.7A CN111784193B (en) 2020-07-17 2020-07-17 Product performance consistency inspection method based on normal distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010689731.7A CN111784193B (en) 2020-07-17 2020-07-17 Product performance consistency inspection method based on normal distribution

Publications (2)

Publication Number Publication Date
CN111784193A true CN111784193A (en) 2020-10-16
CN111784193B CN111784193B (en) 2024-03-26

Family

ID=72764570

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010689731.7A Active CN111784193B (en) 2020-07-17 2020-07-17 Product performance consistency inspection method based on normal distribution

Country Status (1)

Country Link
CN (1) CN111784193B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113743707A (en) * 2021-03-30 2021-12-03 中国人民解放军国防科技大学 Product credibility calculation method based on uniform distribution
CN113762981A (en) * 2021-03-30 2021-12-07 中国人民解放军国防科技大学 Product credibility calculation method based on exponential distribution

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001117906A (en) * 1999-10-14 2001-04-27 Building Research Inst Ministry Of Construction Method for estimating plane position accuracy of spatial information
JP2013178637A (en) * 2012-02-28 2013-09-09 Nippon Telegr & Teleph Corp <Ntt> Data verification apparatus, data verification method and program
CN106202002A (en) * 2016-06-28 2016-12-07 河海大学 A kind of for detecting the method whether series of hydrological parameter makes a variation
CN106777928A (en) * 2016-12-01 2017-05-31 中国人民解放军国防科学技术大学 Towards the Bayes prior distribation building methods of normal distribution data sample
WO2019015226A1 (en) * 2017-07-19 2019-01-24 厦门理工学院 Method for rapidly identifying wind speed distribution pattern
CN109858147A (en) * 2019-01-30 2019-06-07 西南石油大学 A kind of borehole well instability quantifying risk evaluation method based on Reliability Theory

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001117906A (en) * 1999-10-14 2001-04-27 Building Research Inst Ministry Of Construction Method for estimating plane position accuracy of spatial information
JP2013178637A (en) * 2012-02-28 2013-09-09 Nippon Telegr & Teleph Corp <Ntt> Data verification apparatus, data verification method and program
CN106202002A (en) * 2016-06-28 2016-12-07 河海大学 A kind of for detecting the method whether series of hydrological parameter makes a variation
CN106777928A (en) * 2016-12-01 2017-05-31 中国人民解放军国防科学技术大学 Towards the Bayes prior distribation building methods of normal distribution data sample
WO2019015226A1 (en) * 2017-07-19 2019-01-24 厦门理工学院 Method for rapidly identifying wind speed distribution pattern
CN109858147A (en) * 2019-01-30 2019-06-07 西南石油大学 A kind of borehole well instability quantifying risk evaluation method based on Reliability Theory

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
夏新涛;陈向峰;常振;: "模糊等价关系下滚动轴承振动性能变异分析", 航空动力学报, no. 11, 14 November 2018 (2018-11-14) *
杨华波;蔡洪;张士峰;: "基于混合验前分布的制导精度评定方法", 航空学报, no. 05, 25 May 2009 (2009-05-25) *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113743707A (en) * 2021-03-30 2021-12-03 中国人民解放军国防科技大学 Product credibility calculation method based on uniform distribution
CN113762981A (en) * 2021-03-30 2021-12-07 中国人民解放军国防科技大学 Product credibility calculation method based on exponential distribution
CN113762981B (en) * 2021-03-30 2023-11-28 中国人民解放军国防科技大学 Product credibility calculation method based on index distribution
CN113743707B (en) * 2021-03-30 2024-03-01 中国人民解放军国防科技大学 Product credibility calculation method based on uniform distribution

Also Published As

Publication number Publication date
CN111784193B (en) 2024-03-26

Similar Documents

Publication Publication Date Title
CN111784193B (en) Product performance consistency inspection method based on normal distribution
Mandelkern et al. Crystallization kinetics in high polymers. I. Bulk polymers
Al-Othman et al. Uncertainty modelling in power system state estimation
Morrison Obtaining uncertainty measures on slope and intercept of a least squares fit with Excel’s LINEST
Bertin et al. Energy dissipation rates, eddy diffusivity, and the Prandtl number: An in situ experimental approach and its consequences on radar estimate of turbulent parameters
CN106570286A (en) Bridge-ballastless track structure extreme temperature prediction method and system
CN106383501A (en) Carbon emission process monitoring and control method
CN107634516A (en) A kind of distribution method for estimating state based on Grey Markov Chain
CN109858112A (en) Numerical inversion analysis method based on structural stress monitoring result
CN115158076A (en) Metering error evaluation method, device and computer readable storage medium
CN111090907B (en) Flight test transition judgment method
US20090292493A1 (en) Probabilistic uncertainty estimator
Brown Evaluation of parametrization schemes for the convective boundary layer using large-eddy simulation results
CN113360983A (en) Slope reliability analysis and risk assessment method
Epstein et al. Stochastic analysis of meteorological fields
Keith et al. Interaction via diffusion of polymer crystals in growing spherulites
CN109145258B (en) Weibull distribution parameter confidence interval estimation method based on nonlinear fitting
CN113743707B (en) Product credibility calculation method based on uniform distribution
CN111811827B (en) Product performance consistency inspection method based on Rayleigh distribution
Hanaish et al. Daily rainfall disaggregation using HYETOS model for Peninsular Malaysia
EP4009023A1 (en) Method for simulation of a determination of a temperature of a medium
Monteiro et al. Asteroseismic inference for solar-type stars
CN113468795A (en) Universal building method for different parameter material partition pouring/filling dam deformation mixed model
Hermida et al. Analytical-numerical procedure to determine if a set of experimental curves can be superimposed to form a master curve
CN111210877A (en) Method and device for deducing physical property parameters

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