CN113420916A - Design method of multivariate quality characteristic dynamic monitoring strategy - Google Patents

Design method of multivariate quality characteristic dynamic monitoring strategy Download PDF

Info

Publication number
CN113420916A
CN113420916A CN202110638721.5A CN202110638721A CN113420916A CN 113420916 A CN113420916 A CN 113420916A CN 202110638721 A CN202110638721 A CN 202110638721A CN 113420916 A CN113420916 A CN 113420916A
Authority
CN
China
Prior art keywords
monitoring
state
production process
scheme
statistics
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
CN202110638721.5A
Other languages
Chinese (zh)
Other versions
CN113420916B (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.)
Beijing Institute of Electronic System Engineering
Original Assignee
Beijing Institute of Electronic System Engineering
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 Beijing Institute of Electronic System Engineering filed Critical Beijing Institute of Electronic System Engineering
Priority to CN202110638721.5A priority Critical patent/CN113420916B/en
Publication of CN113420916A publication Critical patent/CN113420916A/en
Application granted granted Critical
Publication of CN113420916B publication Critical patent/CN113420916B/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/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Business, Economics & Management (AREA)
  • Databases & Information Systems (AREA)
  • Economics (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Software Systems (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Business, Economics & Management (AREA)
  • Probability & Statistics with Applications (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Computing Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • General Factory Administration (AREA)

Abstract

One embodiment of the invention discloses a design method of a multivariate quality characteristic dynamic monitoring strategy, which comprises the following steps: designing a loose monitoring scheme and a strict monitoring scheme criterion: in the initial stage of the production process, a loose monitoring scheme is adopted to collect data, monitoring statistics is obtained, the state of the production process is judged, if the production process is in a controlled state, the loose monitoring scheme is still adopted when the monitoring statistics is obtained next time, if the production process is in an out-of-control state, an alarm signal is given, and if the production process is in an alert state, the strict monitoring scheme is adopted when the monitoring statistics is obtained next time. Based on the design criteria, an optimization model with the shortest average alarm time is established, and design parameters of the optimal loose monitoring scheme and the optimal strict monitoring scheme are solved.

Description

Design method of multivariate quality characteristic dynamic monitoring strategy
Technical Field
The present invention relates to the field of production monitoring. And more particularly, to a design method of a multivariate quality characteristic dynamic monitoring strategy.
Background
With the gradual development of the traditional production process towards the direction of integration and diversification, more than one key quality characteristic often exists in the production process, so that the multivariate quality characteristic monitoring method based on the multivariate statistical process inspection theory is widely applied to the field of production process monitoring, can timely discover the quality drift of the multivariate production process, and provides decision information support for production process quality managers.
The design parameters of the traditional multivariate quality characteristic monitoring method comprise sampling sample size, sampling interval and control limit, the design parameters are fixed, and the corresponding scheme can be calculated by giving specific false alarm rate requirements. It is characterized by simple design, but this method often requires a large number of samples. In practical engineering application, sometimes a sample to be detected is very expensive, the detection cost can be obviously increased by the requirement of a large sample, and meanwhile, in order to meet the requirement of practical production management, the detection time often has a specific upper limit requirement. When the limitation of resource conditions such as sample size, time, cost and the like cannot meet the requirement of a large sample of the traditional method, once the quality of a multivariate production process drifts, if the quality characteristic monitoring is carried out by continuously adopting the traditional design method with fixed parameters, the average alarm time is obviously increased, so that a manager cannot quickly find the quality defect of the production process. Under the condition that the quality defect of the production process cannot be found in time, the production process is in a quality drifting state in the last period, the defective rate of products is further improved, the quality loss of the production process is aggravated, and unnecessary operation cost is brought to the production process.
Disclosure of Invention
In view of this, a first embodiment of the present invention provides a method for designing a multivariate quality characteristic dynamic monitoring strategy, including:
designing a loose monitoring scheme and a strict monitoring scheme criterion:
in the initial stage of the production process, a loose monitoring scheme is adopted to collect data, a first monitoring statistic is obtained, the state of the production process is judged, if the production process is in a controlled state, the loose monitoring scheme is still adopted when the monitoring statistic is obtained next time, if the production process is in an out-of-control state, an alarm signal is given, and if the production process is in an alert state, a strict monitoring scheme is adopted when the monitoring statistic is obtained next time.
In one embodiment, the loosely monitored scenario parameters include: first sample size n1First sampling interval t1A first control limit CL1And a first warning limit WL1
And acquiring monitoring statistics according to the first sampling sample amount, if the monitoring statistics is smaller than a first warning limit, judging that the production process is in a controlled state, otherwise judging whether the monitoring statistics is larger than the first control limit, if so, judging that the production process is in an out-of-control state, giving an alarm signal, otherwise, judging that the production process is in a warning state, and acquiring the monitoring statistics according to the strict scheme when acquiring the monitoring statistics next time.
In one embodiment, the rigorous monitoring regimen data comprises: second sample size n2Second sampling interval t2A second control limit CL2And a second warning limit WL2
And acquiring monitoring statistics according to the second sampling sample size, if the monitoring statistics is smaller than a second warning limit, judging that the production process is in a controlled state, acquiring the monitoring statistics according to a loose monitoring scheme next time, otherwise, judging whether the monitoring statistics is larger than the second control limit, if so, judging that the production process is in an out-of-control state, and giving an alarm signal, otherwise, judging that the production process is in a warning state, and when the monitoring statistics is acquired next time, continuing to adopt the strict scheme.
In a specific embodiment, the optimal design parameters of the loose monitoring scheme and the strict monitoring scheme are selected, so that the alarm time of the monitoring method is shortest.
In a specific embodiment, the step of solving the optimal design parameters of the loose monitoring scheme and the strict monitoring scheme includes:
the type of distribution of the monitoring statistics is obtained,
a state transition probability matrix is established and,
establishing an average alarm time expression of the dynamic monitoring scheme,
the shortest average alarm time is taken as an optimization criterion to establish the optimumModeling to obtain optimal loose monitoring scheme and strict monitoring scheme
Figure BDA0003106258620000021
In one embodiment, obtaining the monitoring statistic expression comprises:
obtaining n v-element quality characteristic samples X from production line1,X2,...,Xn,XiI is more than or equal to 1 and less than or equal to n, constructing sample statistics
Figure BDA0003106258620000022
Establishing expressions for monitoring statistics
Figure BDA0003106258620000023
Wherein the sample X is samplediSubject to a target value of mu0And a multivariate normal distribution of the covariance matrix sigma,
calculating the size of the drift
When a mass characteristic shift of Δ occurs, calculation is made
Figure BDA0003106258620000031
When no mass drift occurs, the first monitoring statistic obeys chi-square distribution with the degree of freedom v; when a quality characteristic shift of Δ occurs, the first monitoring statistic obeys a degree of freedom of v and the position parameter is nd2Non-central chi-square distribution of (1) is recorded as
Figure BDA0003106258620000032
In a specific embodiment, the state transition probability matrix is:
Figure BDA0003106258620000033
wherein ,R*Indicating that the monitoring statistic is in a controlled state, T*Indicating that the monitored statistic is armed, pR*R*Indicating that the state remains at R between two samples*The probability of (c) is:
Figure BDA0003106258620000034
pR*T*indicates that the current state is R*The next state is T*The probability of (c) is:
Figure BDA0003106258620000035
pT*R*indicates that the current state is T*The next state is R*The probability of (c) is:
Figure BDA0003106258620000036
pT*T*indicates that the current state is T*The next state is T*The probability of (c) is:
Figure BDA0003106258620000037
Figure BDA0003106258620000038
a cumulative probability distribution function representing a non-central chi-squared distribution.
In one embodiment, the average alarm time expression of the dynamic monitoring scheme is:
ATS(n1,t1,WL1,CL1;n2,t2,WL2,CL2|Δ)
=[p0,1-p0]×[I2×2-QΔ]-1×[t1,t2]'
wherein ,
Figure BDA0003106258620000039
p0=P(Y<WL1|Y<CL1)=P(Y<WL2|Y<CL2)。
in one embodiment, an optimization model is established and an optimization scheme is solved
Figure BDA0003106258620000041
Comprises the following steps:
minATS(n1,t1,WL1,CL1,n2,t2,WL2,CL2|Δ)
s.t.n1p0+n2(1-p0)≤n*
t1p0+t2(1-p0)≤t*
n1≤n2
t2≤t1
p0=P(Y<WL1|Y<CL1)=P(Y<WL2|Y<CL2)
Figure BDA0003106258620000042
in one embodiment, the average alarm time of the optimal monitoring method is calculated as
Figure BDA0003106258620000043
The invention has the following beneficial effects:
the invention solves the problems that the traditional multivariate quality characteristic monitoring method needs large sample amount and low monitoring efficiency and can not meet the actual requirements of high-efficiency monitoring and quick feedback in the modern production process, can flexibly meet the requirements of actual monitoring constraint conditions and provide a dynamic optimization design scheme with the shortest average alarm time, can find the problem of quality drift in shorter time compared with the traditional multivariate quality characteristic monitoring method, further improves the monitoring efficiency, is convenient for engineering technicians to use, and has good application value.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 shows a flow diagram of a method for dynamic monitoring of quality characteristics according to an embodiment of the invention.
Fig. 2 shows a flowchart of a method for obtaining optimal loose monitoring solution data and strict monitoring solution data according to an embodiment of the invention.
Detailed Description
In order to make the technical solutions and advantages of the present invention clearer, the following will describe embodiments of the present invention in further detail with reference to the accompanying drawings.
As shown in fig. 1, a method for designing a multivariate quality characteristic dynamic monitoring strategy includes:
designing a loose monitoring scheme and a strict monitoring scheme criterion:
in the initial stage of the production process, a loose monitoring scheme is adopted to collect data, monitoring statistics is obtained, the state of the production process is judged, if the production process is in a controlled state, the loose monitoring scheme is still adopted when the monitoring statistics is obtained next time, if the production process is in an out-of-control state, an alarm signal is given, and if the production process is in an alert state, a strict monitoring scheme is adopted when the monitoring statistics is obtained next time.
In a specific embodimentIn an example, the loose monitoring scheme parameters include: first sample size n1First sampling interval t1A first control limit CL1And a first warning limit WL1
And acquiring a monitoring statistic according to the first sampling sample amount, if the monitoring statistic is smaller than a first warning limit, judging that the production process is in a controlled state, otherwise, judging whether the monitoring statistic is larger than the first control limit, if so, judging that the production process is in an out-of-control state, giving an alarm signal, otherwise, judging that the production process is in a warning state, and adopting the strict scheme when the monitoring statistic is acquired next time.
In one embodiment, the rigorous monitoring regimen data comprises: second sample size n2Second sampling interval t2A second control limit CL2And a second warning limit WL2
And acquiring monitoring statistics according to the second sampling sample amount, if the monitoring statistics is smaller than a second warning limit, judging that the production process is in a controlled state, acquiring the monitoring statistics according to a loose monitoring scheme when acquiring the monitoring statistics next time, otherwise judging whether the monitoring statistics is larger than the second control limit, if so, judging that the production process is out of control, and giving an alarm signal, otherwise, judging that the production process is in a warning state, and adopting the strict scheme when acquiring the monitoring statistics next time.
More preferably, the design parameters of the optimal loose monitoring scheme and the optimal strict monitoring scheme are selected, so that the alarm time of the monitoring method is shortest.
In another embodiment, as shown in fig. 2, solving the design parameters of the optimal relaxed monitoring scheme and the optimal strict monitoring scheme includes:
obtaining a distribution type of the monitoring statistics, comprising:
obtaining n sampling samples X from production line1,X2,...,Xn,XiI is more than or equal to 1 and less than or equal to n, constructing sample statistics
Figure BDA0003106258620000051
Establishing expressions for monitoring statistics
Figure BDA0003106258620000052
Wherein the sample X is samplediSubject to a target value of mu0And a multivariate normal distribution of the covariance matrix sigma,
calculating the size of the drift
When a mass characteristic shift of Δ occurs, calculation is made
Figure BDA0003106258620000061
When no quality drift occurs, monitoring statistics obeying chi-square distribution with the degree of freedom v; when the quality characteristic drift is delta, the monitoring statistic obeys the degree of freedom v and the position parameter nd2Non-central chi-square distribution of (1) is recorded as
Figure BDA0003106258620000062
The distribution type of the obtained monitoring statistics can be obtained according to the step of obtaining the distribution type of the monitoring statistics, and will not be described in detail herein.
A state transition probability matrix is established and,
the state transition probability matrix is:
Figure BDA0003106258620000063
wherein ,R*Indicating that the monitoring statistic is in a controlled state, T*Indicating that the monitored statistic is armed, pR*R*Indicating that the state remains at R between two samples*The probability of (c) is:
Figure BDA0003106258620000064
pR*T*indicates that the current state is R*The next state is T*The probability of (c) is:
Figure BDA0003106258620000065
pT*R*indicates that the current state is T*The next state is R*The probability of (c) is:
Figure BDA0003106258620000066
pT*T*indicates that the current state is T*The next state is T*The probability of (c) is:
Figure BDA0003106258620000067
Figure BDA0003106258620000068
a cumulative probability distribution function representing a non-central chi-squared distribution.
Establishing an average alarm time expression of the dynamic monitoring scheme,
ATS(n1,t1,WL1,CL1;n2,t2,WL2,CL2|Δ)
=[p0,1-p0]×[I2×2-QΔ]-1×[t1,t2]'
wherein ,
Figure BDA0003106258620000069
p0=P(Y<WL1|Y<CL1)=P(Y<WL2|Y<CL2)
establishing an optimization model by taking the shortest average alarm time as an optimization criterion to obtain design parameters of an optimal loose monitoring scheme and a strict monitoring scheme
Figure BDA00031062586200000610
The optimization model is as follows:
minATS(n1,t1,WL1,CL1,n2,t2,WL2,CL2|Δ)
s.t.n1p0+n2(1-p0)≤n*
t1p0+t2(1-p0)≤t*
n1≤n2
t2≤t1
p0=P(Y<WL1|Y<CL1)=P(Y<WL2|Y<CL2)
Figure BDA0003106258620000071
according to the design parameters of the optimal loose monitoring scheme and the strict monitoring scheme, the average alarm time of the optimal monitoring method is calculated as
Figure BDA0003106258620000072
In one embodiment, let the target value be μ0=[0.7,1.2,1.4]The mass characteristic number v is 3, and the covariance matrix is:
Figure BDA0003106258620000073
the mass characteristic shift vector Δ is [0.04,0.15,0.3 ].
Resource bar according to actual monitoringPiece limitation and statistical property requirement, given an upper limit n of the average sample size*20, average sampling interval upper limit t*2 and an average false alarm rate upper limit p*The value is 0.005, so that,
constructing sample statistics
Figure BDA0003106258620000074
Establishing expressions for monitoring statistics
Figure BDA0003106258620000075
Calculating the size of the drift amount:
when the mass characteristic shift Δ ═ 0.04,0.15,0.3 occurred, the magnitude of the shift amount was calculated
Figure BDA0003106258620000081
When no quality drift occurs in the multi-element production process, monitoring the statistic Y to obey chi-square distribution with the degree of freedom of 3; when the multivariate mass characteristic shift Delta is [0.04,0.15,0.3]]Then, the monitoring statistic Y obeys a non-central chi-square distribution with a degree of freedom of 3 and a position parameter of 0.0991n, and is recorded as
Figure BDA0003106258620000082
According to the magnitude relation of chi-square statistic and control limit and warning limit to establish the transition probability matrix of self state of dynamic multi-element production process and monitoring scheme judging state,
Figure BDA0003106258620000083
wherein, the calculation steps of each element in the matrix are as follows:
Figure BDA0003106258620000084
Figure BDA0003106258620000085
Figure BDA0003106258620000086
Figure BDA0003106258620000087
establishing an average alarm time expression of the dynamic monitoring scheme,
Figure BDA0003106258620000088
wherein ,
Figure BDA0003106258620000089
Figure BDA00031062586200000810
determining the constraint conditions of the design of the monitoring optimization scheme according to the specific conditions of the monitoring resources, wherein the constraint conditions comprise an upper limit n of the average sample size*20, average sampling interval upper limit t*2 and false alarm rate upper bound p*0.005, and then solve the optimum
Figure BDA00031062586200000811
Figure BDA0003106258620000091
Finally, the optimal dynamic monitoring scheme is obtained
(n1,t1,WL1,CL1;n2,t2,WL2,CL2)=(18,2.08,6.37,13.02,40,0.66,6.02,9.84).
Estimating the monitoring efficiency of the calculated optimal dynamic monitoring scheme, and calculating the average alarm time of the optimal monitoring scheme:
ATS(18,2.08,6.37,13.02,40,0.66,6.02,9.84|[0.04,0.15,0.3])=20.53.
therefore, the monitoring statistic is obtained according to the loose monitoring scheme, the state of the production process is judged, if the monitoring statistic is smaller than the first warning limit of 6.37, the production process is judged to be in a controlled state, if not, whether the monitoring statistic is larger than the first control limit of 13.02 is judged, if yes, the production process is judged to be in an out-of-control state, an alarm signal is given, otherwise, the production process is judged to be in a warning state, and when the monitoring statistic is obtained next time, the monitoring statistic is obtained according to the strict scheme.
And acquiring a second monitoring statistic according to the strict scheme, if the monitoring statistic is less than a second warning limit of 6.02, judging that the production process is in a controlled state, acquiring the monitoring statistic next time according to a loose monitoring scheme, otherwise, judging whether the monitoring statistic is greater than a second control limit of 9.84, if so, judging that the production process is in an out-of-control state, giving an alarm signal, otherwise, judging that the production process is in a warning state, and executing according to the strict scheme when acquiring the monitoring statistic next time.
The alarm time given by the scheme is shortest, limited time and sample resources can be utilized to the maximum extent, resource waste is avoided, when parameter drift occurs in the multivariate quality characteristics, the drift can be found quickly by the provided method, and the quality loss in the production process is reduced.
As can be seen from the following table, at n*=10,h*At 2, p 0.005, and different amounts of drift, the results of comparing ATS of the conventional method and the proposed method,
Figure BDA0003106258620000092
Figure BDA0003106258620000101
therefore, under the condition of different drift amounts, the average alarm time ATS of the optimal dynamic monitoring method provided by the invention is obviously smaller than that of the traditional method. This shows that when parameter drift occurs in the multivariate quality characteristics, the proposed method can quickly find the drift, reducing the quality loss of the production process.
It should be understood that the above-mentioned embodiments of the present invention are only examples for clearly illustrating the present invention, and are not intended to limit the embodiments of the present invention, and it will be obvious to those skilled in the art that other variations or modifications may be made on the basis of the above description, and all embodiments may not be exhaustive, and all obvious variations or modifications may be included within the scope of the present invention.

Claims (10)

1. A multivariate quality characteristic dynamic monitoring strategy design method is characterized by comprising the following steps:
designing a loose monitoring scheme and a strict monitoring scheme criterion:
in the initial stage of the production process, a loose monitoring scheme is adopted to collect data, a first monitoring statistic is obtained, the state of the production process is judged, if the production process is in a controlled state, the loose monitoring scheme is still adopted when the monitoring statistic is obtained next time, if the production process is in an out-of-control state, an alarm signal is given, and if the production process is in an alert state, a strict monitoring scheme is adopted when the monitoring statistic is obtained next time.
2. The method of claim 1, wherein the loosely monitored scheme parameters comprise: first sample size n1First sampling interval t1A first control limit CL1And a first warning limit WL1
And acquiring monitoring statistics according to the first sampling sample amount, if the monitoring statistics is smaller than a first warning limit, judging that the production process is in a controlled state, otherwise judging whether the monitoring statistics is larger than the first control limit, if so, judging that the production process is in an out-of-control state, giving an alarm signal, otherwise, judging that the production process is in a warning state, and acquiring the monitoring statistics according to the strict scheme when acquiring the monitoring statistics next time.
3. The method of claim 2, wherein the rigorous monitoring scenario data comprises: second sample size n2Second sampling interval t2A second control limit CL2And a second warning limit WL2
And acquiring monitoring statistics according to the second sampling sample size, if the monitoring statistics is smaller than a second warning limit, judging that the production process is in a controlled state, acquiring the monitoring statistics according to a loose monitoring scheme next time, otherwise, judging whether the monitoring statistics is larger than the second control limit, if so, judging that the production process is in an out-of-control state, and giving an alarm signal, otherwise, judging that the production process is in a warning state, and when the monitoring statistics is acquired next time, continuing to adopt the strict scheme.
4. The method of claim 3, wherein the optimal design parameters for the loose monitoring solution and the strict monitoring solution are selected to minimize the alarm time of the monitoring method.
5. The method of claim 4, wherein the step of solving for optimal loosely and tightly monitored solution design parameters comprises:
the type of distribution of the monitoring statistics is obtained,
a state transition probability matrix is established and,
establishing an average alarm time expression of the dynamic monitoring scheme,
the shortest average alarm time is taken as an optimization criterion to establish the optimumModeling to obtain optimal loose monitoring scheme and strict monitoring scheme
Figure FDA0003106258610000011
6. The method of claim 5, wherein obtaining a monitoring statistic expression comprises:
obtaining n v-element quality characteristic samples X from production line1,X2,...,Xn,XiI is more than or equal to 1 and less than or equal to n, constructing sample statistics
Figure FDA0003106258610000021
Establishing expressions for monitoring statistics
Figure FDA0003106258610000022
Wherein the sample X is samplediSubject to a target value of mu0And a multivariate normal distribution of the covariance matrix sigma,
calculating the size of the drift
When a mass characteristic shift of Δ occurs, calculation is made
Figure FDA0003106258610000023
When no mass drift occurs, the first monitoring statistic obeys chi-square distribution with the degree of freedom v; when a quality characteristic shift of Δ occurs, the first monitoring statistic obeys a degree of freedom of v and the position parameter is nd2The non-central chi-square distribution of (a), noted as Y:
Figure FDA0003106258610000024
7. the method of claim 5, wherein the state transition probability matrix is:
Figure FDA0003106258610000025
wherein ,R*Indicating that the monitoring statistic is in a controlled state, T*Indicating that the monitored statistic is in an alert state,
Figure FDA0003106258610000026
indicating that the state remains at R between two samples*The probability of (c) is:
Figure FDA0003106258610000027
Figure FDA0003106258610000028
indicates that the current state is R*The next state is T*The probability of (c) is:
Figure FDA0003106258610000029
Figure FDA00031062586100000210
indicates that the current state is T*The next state is R*The probability of (c) is:
Figure FDA00031062586100000211
Figure FDA00031062586100000212
indicates that the current state is T*The next state is T*The probability of (c) is:
Figure FDA00031062586100000213
Figure FDA00031062586100000214
a cumulative probability distribution function representing a non-central chi-squared distribution.
8. The method of claim 5, wherein the average alarm time of the dynamic monitoring scheme is expressed as:
ATS(n1,t1,WL1,CL1;n2,t2,WL2,CL2|Δ)
=[p0,1-p0]×[I2×2-QΔ]-1×[t1,t2]'
wherein ,
Figure FDA0003106258610000031
p0=P(Y<WL1|Y<CL1)=P(Y<WL2|Y<CL2)。
9. the method of claim 5, wherein the optimization model is established and the optimization solution is solved
Figure FDA0003106258610000032
Comprises the following steps:
Figure FDA0003106258610000033
10. the method of claim 8, wherein the average alarm time of the optimal monitoring method is calculated as follows based on the optimal loose monitoring scheme data and the strict monitoring scheme data
Figure FDA0003106258610000034
CN202110638721.5A 2021-06-08 2021-06-08 Method for designing dynamic monitoring strategy of multiple quality characteristics Active CN113420916B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110638721.5A CN113420916B (en) 2021-06-08 2021-06-08 Method for designing dynamic monitoring strategy of multiple quality characteristics

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110638721.5A CN113420916B (en) 2021-06-08 2021-06-08 Method for designing dynamic monitoring strategy of multiple quality characteristics

Publications (2)

Publication Number Publication Date
CN113420916A true CN113420916A (en) 2021-09-21
CN113420916B CN113420916B (en) 2023-09-29

Family

ID=77788003

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110638721.5A Active CN113420916B (en) 2021-06-08 2021-06-08 Method for designing dynamic monitoring strategy of multiple quality characteristics

Country Status (1)

Country Link
CN (1) CN113420916B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005019808A (en) * 2003-06-27 2005-01-20 Trecenti Technologies Inc Quality control method, quality control system, and manufacturing method of semiconductor device
CN110794809A (en) * 2019-10-28 2020-02-14 郑州航空工业管理学院 VSI EWMA control chart economic design method based on preventive maintenance and quality loss function
CN111352848A (en) * 2020-03-09 2020-06-30 南京航空航天大学 Method for measuring monitorability probability of property in runtime verification
CN112377333A (en) * 2020-12-02 2021-02-19 西安航天动力研究所 Real-time fault monitoring method for turbopump of high-thrust liquid rocket engine

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005019808A (en) * 2003-06-27 2005-01-20 Trecenti Technologies Inc Quality control method, quality control system, and manufacturing method of semiconductor device
CN110794809A (en) * 2019-10-28 2020-02-14 郑州航空工业管理学院 VSI EWMA control chart economic design method based on preventive maintenance and quality loss function
CN111352848A (en) * 2020-03-09 2020-06-30 南京航空航天大学 Method for measuring monitorability probability of property in runtime verification
CN112377333A (en) * 2020-12-02 2021-02-19 西安航天动力研究所 Real-time fault monitoring method for turbopump of high-thrust liquid rocket engine

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
沈钰健: "发动机缸体加工研究和关键特性SPC质量监控应用", 《中国优秀硕士学位论文全文数据库 工 程科技Ⅱ辑》, no. 03, pages 035 - 216 *

Also Published As

Publication number Publication date
CN113420916B (en) 2023-09-29

Similar Documents

Publication Publication Date Title
Bukkapatnam et al. Machine learning and AI for long-term fault prognosis in complex manufacturing systems
JPWO2006114871A1 (en) Data collection device and gateway device
Riaz Monitoring process variability using auxiliary information
KR20230039873A (en) Method for detecting abnormality
CN111309502A (en) Solid state disk service life prediction method
CN112415437A (en) High-voltage switch cabinet cable contact fault detection method based on multiple linear regression model
CN113420916A (en) Design method of multivariate quality characteristic dynamic monitoring strategy
CN113210824B (en) Servo welding gun driving abnormity detection method and equipment
CN117033044B (en) Visual access system and method
CN111125195B (en) Data anomaly detection method and device
CN116781583A (en) Service quality variation detection method and device based on load variation
KR101538548B1 (en) Controling method of very small fraction of nonconforming
CN107832265A (en) The cpu load Forecasting Methodology of desktop based on state aware
CN115080342A (en) Method and system for automatically processing test data of power board card
JP7403408B2 (en) Monitoring equipment, monitoring methods and programs
CN114331688A (en) Method and device for detecting batch operation state of bank counter system business
CN110297140B (en) Fault prediction method and device of power distribution system
CN112364069A (en) Thermocouple fault early warning method and system based on time sequence and storage medium
CN116223012B (en) Intelligent detection method and system for electronic element
CN115437888B (en) Monitoring and early warning method based on micro-service
WO2022196666A1 (en) Information processing device and information processing method
CN112836898B (en) Cloud manufacturing workshop operation hour estimation method
CN114050981B (en) Automatic inspection method and system for stability of elastic network card
CN117113257A (en) Dynamic SPC control method based on multi-value scene
CN117789327A (en) Screw installation quality analysis, early warning and optimization method and system in product assembly process

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