CN110750848A - Method for estimating remaining life of software-hardware degradation system by considering software operation - Google Patents

Method for estimating remaining life of software-hardware degradation system by considering software operation Download PDF

Info

Publication number
CN110750848A
CN110750848A CN201810815534.8A CN201810815534A CN110750848A CN 110750848 A CN110750848 A CN 110750848A CN 201810815534 A CN201810815534 A CN 201810815534A CN 110750848 A CN110750848 A CN 110750848A
Authority
CN
China
Prior art keywords
degradation
software
hardware
performance
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201810815534.8A
Other languages
Chinese (zh)
Inventor
张德平
韩佳佳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201810815534.8A priority Critical patent/CN110750848A/en
Publication of CN110750848A publication Critical patent/CN110750848A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Debugging And Monitoring (AREA)

Abstract

The invention discloses a method for estimating the residual life of a software-hardware degradation system by considering software operation, which comprises the following steps: constructing and screening the performance degradation indexes of the software-hardware system; constructing a system degradation model; providing a residual life estimation algorithm; analyzing a performance degradation model based on the hidden Markov model; and optimizing the system segmental degradation model by adjusting the state transition parameters of the software, further respectively constructing the system degradation model by the performance indexes, and extrapolating the residual life estimation of the whole software-hardware system according to the system structure of the system. The invention has the advantages that: aiming at the traditional problem of estimating the residual life of a software and hardware system by independently considering the software reliability or the hardware reliability, a discrete hidden Markov process is adopted to describe the transfer between software modules, the relation between a hidden state and actual degradation is established, different degradation models are established according to the number of inflection points in a system performance degradation index sample under different software running conditions, and the residual life is estimated according to a system structure.

Description

Method for estimating remaining life of software-hardware degradation system by considering software operation
Technical Field
The invention relates to an analysis method for efficiently estimating system-level residual life, and particularly provides a simple and easy-to-use mathematical analysis method aiming at the problems of complicated simulation experiment process and huge manpower and material consumption of a large software-hardware integrated complex system.
Background
With the continuous improvement of the informatization level, the functions of various control systems are gradually realized by software, and the characteristics of high density of the software and continuous improvement of the complexity are presented. Integration of hardware and software has become an inherent form of large systems of all kinds. The system reliability of the system depends on software and hardware respectively, and further depends on the mutual influence in the process of combining the software and the hardware. Therefore, the reliability analysis of software and hardware interaction for software and hardware systems is more and more widely regarded. The traditional residual life estimation theory needs a large amount of life data, and with the technical progress, the life of modern large-scale systems is longer and longer, so that the failure is rarely generated in a limited development period, enough life data is difficult to obtain, and the existing reliability assessment is difficult to implement. In fact, failure mechanism analysis is the basis for reliability assessment and life prediction. For most systems, the failure mechanism is related to the performance degradation process, which therefore contains a large amount of reliability information. The traditional reliability theory only focuses on the information such as the service life data and the failure times at the failure moment, neglects the relation between the system failure and the performance degradation process, and has the defects in the reliability and residual life evaluation of a high-reliability long-service-life system. In fact, whether the system is failed or not, a cluster of performance degradation data which changes along with time can be measured, the performance degradation data provides a rich and key information source for product life prediction and reliability evaluation, and life characteristics and reliability of the system can be extrapolated, so that a new feasible way is provided for reliability evaluation.
The demand of modern military industry and the development of scientific technology determine that a mechanical system is continuously developed towards high speed, high precision and high reliability; on the other hand, many military system devices are bulky in structure and complex in function. Even if the device is excellent, its performance and state gradually deteriorate with continuous operation and cause occurrence of a failure. Failure of the system not only reduces production efficiency and system reliability, but also may cause serious safety problems and even cause personal injury and death. If the faults occurring in the system operation process can be accurately predicted and corresponding maintenance strategies can be formulated, the method is very meaningful for improving the system performance and avoiding the occurrence of major accidents and casualties. In addition, along with the aggravation of market competition, the design and manufacturing technology of the system is continuously improved, the performance of hardware materials is continuously improved, the reliability of the system is higher and higher, and the service life is longer and longer. In order to ensure the availability of system functions, the reliability index of the system needs to be determined, predicted and evaluated, otherwise, the functional design loses significance. Therefore, how to quickly and accurately estimate the reliability of a high-reliability long-life system within a limited research and development period and test expenses is a problem to be solved in the current reliability engineering field. The existing comprehensive test method for the reliability of the software and hardware hybrid system simply avoids the difference of the failure mechanisms of the software and hardware and does not solve the influence caused by the difference of the failure mechanisms of the software and hardware.
Reliability assessment modeling performance indicators of degraded systems is more advantageous and has been extensively studied by many scholars. The model building method of the degradation quantity generally adopts a random process method, namely a certain random process function is selected to describe the degradation process, and then the relevant parameters of the random process characteristic are solved and determined according to the degradation data. Through distribution fitting inspection, if the distribution rule of the degradation amount is not obvious, the relative increment of the degradation data meets a certain distribution rule, or a random accumulated change rule exists in the performance index, a random process with an independent increment characteristic is used, but modeling on the degradation data is a more appropriate choice.
However, problems associated with large system evaluations still exist, particularly the dynamic modeling of the operational process over time and its complexity. Based on the method, a new viewpoint is provided from a modeling object, the problem of estimation of the residual service life of a software-hardware system influenced by the operation of software modules is considered, the software operation is considered as one kind of impact or influence of hardware performance degradation, the transition between different software modules is described by adopting a discrete hidden Markov process, the relation between the actual degradation and the hidden state is constructed, and the degradation process of the hardware performance is analyzed on the basis.
Disclosure of Invention
The invention aims to invent a method for more accurately estimating the residual life of a system in limited budget and time, namely a degradation quantity segmented modeling method based on a hidden Markov model. And the traditional method for estimating the residual life based on failure data is improved, and the used degradation data samples are easier to obtain and have more key information.
The specific technical scheme of the invention comprises the following steps:
the method comprises the following steps: and (3) defining a software-hardware degradation system, constructing a degradation performance index, and extracting and screening the index.
Step two: and constructing a system degradation model by using a hidden Markov model according to the degradation performance index.
Step three: and designing a system-level residual life estimation algorithm according to the system architecture.
Step four: and (4) experimental design, wherein feature extraction and screening are carried out on the sample data set, the model is constructed, the residual life is predicted, a performance evaluation index is given, and the precision of the estimated value of the residual life is obtained.
The technical scheme provided by the invention has the following beneficial effects:
the hidden Markov model-based degradation quantity segmented modeling method provided by the invention aims at the problems of complicated simulation experiment process and huge manpower and material consumption of a large software-hardware integrated complex system, provides a simple and easy-to-use mathematical analysis method, combines the hidden Markov model to perform segmented modeling of degradation indexes, efficiently improves the accuracy of residual life estimation, improves the traditional failure data-based residual life estimation method, and uses degradation data samples which are easier to obtain and have more key information.
Drawings
FIG. 1 is a general flow diagram of the process of the present invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings and related algorithms.
The invention is based on the degradation data of the system to analyze, and obtains the performance degradation index as a data sample by a signal decomposition and dimension reduction method. And selecting a hidden Markov model to describe the influence of the operation of the software module on the hardware performance, predicting a hidden software state sequence, and solving a path of software operation state transition. Therefore, the degradation data of the hardware performance is modeled in a segmented mode, and the estimation of the residual service life at different stages is closer to a true value. The specific implementation steps are as follows:
1. index extraction and dimensionality reduction
The estimation of the remaining life of the software-hardware system first needs to determine key characteristic indexes representing the operation state of the system. The characteristic indexes not only comprise functional performance parameters of the system, but also comprise indirect characteristic parameters extracted from sensor monitoring signals of vibration, noise, temperature and the like. These features are greatly affected by software shock. During the operation of the software-hardware system, various detection signals of the system are information carriers reflecting whether the system is normally operated or not. Because the sensor monitoring information can not be directly used for degradation modeling, in order to improve the degradation evaluation precision, a signal processing method is required to be utilized to deeply research a life characteristic extraction method, so that the adopted method and the extracted characteristics can fully reflect the running state of the system. At present, the characteristic domains for representing the performance degradation indexes of the vibration signals are 3 types of time domains, frequency domains and time-frequency domains. The time domain features include an average value, a root mean square, an RMS, a kurtosis index, and the like. The frequency domain characteristics may represent the dispersion or concentration of the spectrum, the variation of the main band, etc. The frequency domain and time-frequency domain features are processed by methods such as wavelet transformation, empirical mode decomposition and the like. The specific process is shown as algorithm 1 and algorithm 2. In addition, after feature extraction is performed on the system monitoring signal data, some feature vectors may be very correlated, resulting in overlapping information. Therefore, in order to reduce data redundancy, the experiment adopts a Principal Component Analysis (PCA) method to reduce the dimension. The PCA protocol is shown in algorithm 3.
Algorithm 1: extraction of features by wavelet decomposition
In wavelet decomposition, a Hilbert space is decomposed into orthogonal sums of a plurality of subspaces according to different scale factors j, and further decomposed into a closed-packet subspace, and a wavelet packet base sequence is constructed to meet the following conditions:
wherein h and g are high-pass filter and low-pass filter coefficients respectively; k is 1, 2, …, j is a time position parameter and n is a frequency parameter.
After decomposition, reconstructing the decomposition coefficients of each frequency band to the 0 th layer by layer, and extracting signals in each frequency band, wherein if the data length of the original signal x (k) is N, the signal energy of each frequency band after decomposition is represented as:
Figure BSA0000167636500000032
wherein the content of the first and second substances,
Figure BSA0000167636500000033
representing the magnitude of a discrete signal located in a subspace in a wavelet decomposition at a resolution j. The feature vector corresponding to the energy sequence can be determined therefrom.
And 2, algorithm: EMD decomposition method for extracting features
The EMD decomposition involves the concept of eigenmode functions, which must satisfy the following two conditions:
1) the number of local extreme points and zero crossings must be equal or differ by at most one over the entire data range.
2) At any time, the average value of the upper envelope formed by the local maximum point and the lower envelope formed by the local minimum point is zero, that is, the upper envelope and the lower envelope are locally symmetrical with respect to the time axis.
EMD decomposition process:
1) all local maximum points of the signal data x (t) are found, all the local maximum points are connected by a cubic spline curve to form an upper envelope line, and all the local minimum points are connected to form a lower envelope line.
2) Mean value m of upper and lower envelope lines1Subtracting the mean value from the original signal data to obtain a new data sequence h1
x(t)-m1=h1
If h is1Two conditions for the eigenmode function are satisfied, then h1Is the first IMF component c decomposed1. Otherwise h1Instead, h is required to be1Repeating step 1) as the original signal, and continuing the screening until a first IMF component c is generated1
3) And (3) calculating:
r1=x(t)-c1
will r is1Repeating the steps 1) to 3) as the original signal, and so on to obtain n components of the original signal and a residual component r which can not be decomposed any morenNamely:
Figure BSA0000167636500000041
algorithm 3: principal component analysis method for reducing dimension
(1) Raw data indices were normalized:
(2) calculating a correlation coefficient matrix R ═ (R)ij):
Figure BSA0000167636500000043
rij(i, j ═ 1, 2.. times, p) is the correlation coefficient of the original variable, rij=rji
(3) Calculating a characteristic value and a characteristic vector:
solving a characteristic equation | λ I-R | ═ 0, solving characteristic values by a common Jacobian method, and arranging the characteristic values in the order of magnitude: lambda [ alpha ]1≥λ2≥...≥λp≥0。
Respectively determining corresponding characteristic values lambdaiCharacteristic vector e ofi(i ═ 1, 2,. p.) such that | | | ei||=1;
(4) Calculating the principal component contribution rate and the accumulated contribution rate:
contribution toRate of change(i=1,2,...,p)
Cumulative contribution rate
Figure BSA0000167636500000051
(i=1,2,...,p)。
2. Hidden Markov model
The method comprises the steps of constructing a system degradation and failure domain model by using a hidden Markov model, estimating parameters of the hidden Markov model by using a Baum-Welch algorithm, predicting a hidden state sequence by using a Viterbi algorithm, and solving a path of software running state transition. And (3) assuming that K inflection points exist in the hidden sequence, determining an optimal value of K by adopting a Bayesian Information Criterion (BIC), segmenting the degraded data according to the value of K, and respectively constructing a hardware performance degradation model for each segment.
The hidden Markov model is specifically processed as described in algorithm 4.
And algorithm 4: based on hidden Markov model
Inputting: observing the sequence of states, the parameter K, and the termination condition
Initialization: initial state probability vector pi, state transition probability matrix A
And (3) outputting: software running state transition path
# set initial hidden State
z<-runif(121,min=1,max=3)
# Bayesian Information Criterion (BIC) determination of k
k<-3
Probability of transfer
p<-matrix(c(.4,.3,.3,.4,.5,.1,0,0,1),3,3)
Initial probability initial value of # hidden state
pi=c(1,0,0,0)
# known hidden State sequence, time series of varying points T
Figure BSA0000167636500000052
Figure BSA0000167636500000061
3. Model of degenerate domain
Y (n) is at time tnThe predicted sensor index value d (n), Y (n) -Y (n-1) is set to a period [ (n-1) h, nh]The performance degradation model of the software-hardware system can be described as:
Figure BSA0000167636500000062
the performance degradation models given by the above equations may further be combined, assuming that under a given software operating state Z (n-1) ═ i, Z (n) ═ j, the system performance degradation increment d (n) has a probability density function G (d | i, j, Θ)ij) I.e. by
D(n)~G(d|i,j,Θij),i,j∈I
Wherein, thetaijIs a parameter vector determined by Z (n-1) ═ i, Z (n) ═ j. Thus, of all possible I, j e I, the probability density set G ═ G (d | I, j, Θ) of the system performance degradation increments d (n) in each stateij),i,j∈I}。
Under the condition of determining the distribution of the degradation increment in different software running states, the future performance degradation amount of the system in different software running states can be easily analyzed, and the future performance degradation amount can be specifically determined by the following formula:
wherein n iscIs the current number of cycles, Y (n'c) To be at the moment of timeThe amount of degradation in the performance index observed.
At a given significance level of α of 0.05, T-tests were chosen for goodness of fit tests, all with P values greater than 0.05, so the original hypothesis was accepted.According to the set hardware performance failure threshold D, the performance degradation data obeys normal distribution
Figure BSA0000167636500000064
The probability distribution function of (a) is:
Figure BSA0000167636500000065
and estimating parameters of the performance degradation model by using a nonlinear least square method to obtain normal distribution model parameters of each section of data. Given different time nodes, the remaining life of each piece of hardware can be estimated, assuming a failure threshold.
The evaluation index is measured by Root Mean Square Error (RMSE), mean error (AE) and Mean Square Percentage Error (MSPE) between the estimated residual life and the real residual life, and the smaller the value, the better the prediction result.
Root Mean Square Error (RMSE)
Figure BSA0000167636500000066
Mean error (AE)
Figure BSA0000167636500000067
Mean Square Percent Error (MSPE)
Wherein, yiDenotes true value, y'iIndicating the predicted value.
And algorithm 5: construction of degenerate domain model
Inputting: degradation index data, parameter K
And (3) outputting: k degradation model functions
% plotted initial data histogram
histfit(data);
% normal distribution test
normplot(data);
% parameter estimation
[mu,sigma,muci,sigmaci]=normfit(data);
% T test
[h,sig,ci]=ttest(data,mu);
4. Remaining life estimation algorithm
Assuming that the software-hardware system failure is defined as the time when the system performance index value enters the failure area for the first time under the operation state of different software modules, and the failure time is the time when the system performance index value enters the failure area for the first time. According to the hidden Markov process of software operation, the time t is generated by simulationNNh software execution state sequence { z (n): n iscN is less than or equal to N. Each software execution state sequence generated by simulation can be regarded as a determined software running path, and the residual life of the system can be used as the residual life estimation of the whole software-hardware system by solving the average residual life of all paths by utilizing a law of large numbers. The specific estimation process is as follows:
(1) the maximum possible life (cycle) of the software-hardware system is set to N, and the number of simulation generated path samples is set to M (for example, M is 1000). Where the maximum lifetime N may be determined from historical data observed by the actual software-hardware system, such as 557 maximum system lifetime in the observed system performance data set for a certain equipment system, and thus N may be set to 600.
(2) According to the observed system performance index data set, index data with obvious performance index change in different software running states are selected, and failure domain models of different performance indexes in different software running states are constructed.
(3) Generating M software execution transfer paths by simulation, and generating corresponding system performance index degradation amount in each software execution transfer, thereby forming corresponding degradation path data pair { (i, v)i),i=nc,nc+1,...,N}。
(4) Degenerate path data pair { (i, v) generated using simulationi),i=nc,nc+ 1.. ang.N }, at each occurrenceA point in time tnCalculate the expected degradation value of performance E [ Y (n) for each state]Then, the remaining life of z (n) k, k e I in different software execution states is determined. Here, the remaining life is determined by the time when the system performance index first enters the failure region:
Figure BSA0000167636500000081
from the M paths, an estimate of the remaining life of the system in M +1 states can thus be obtained.
(5) In order to reduce the data pollution of the collected performance degradation indexes, the simulation is repeatedly carried out by adopting b different simulation starting points, namely, the simulation is carried out by adopting
Figure BSA0000167636500000082
Taking the time point as a new current time point, repeatedly carrying out the simulation of (3) and (4), and estimating
Figure BSA0000167636500000083
The estimates are then weighted averaged
Figure BSA0000167636500000084
Wherein the weight ω isiAnd assigning according to the reliability of each sample data.
(6) For a sufficiently large integer M, the remaining life of the different subsystems can be estimated:
(7) the remaining life of the different subsystems is estimated separately.
(8) And estimating the system-level residual life of the software-hardware system by using the system structure.
5. Residual life estimation accuracy improvement
The method comprises the steps of extracting and screening indexes by using a signal decomposition algorithm and a principal component analysis method to obtain performance degradation indexes, describing the influence of the operation of a software module on the hardware performance by using a hidden Markov model, predicting a hidden software state sequence and solving the path of software operation state transition. Therefore, the degradation data of the hardware performance is subjected to segmented modeling, so that the estimation of the residual service life at different stages is closer to a true value, and the prediction precision is improved. And the traditional method for estimating the residual life based on failure data is improved, and the used degradation data samples are easier to obtain and have more key information.

Claims (5)

1. The method for improving the residual life estimation precision of the software and hardware degradation system under the software running condition is characterized by comprising the following steps of:
(1) and (3) defining a software-hardware degradation system, constructing a degradation performance index, and extracting and screening the index.
(2) And constructing a system hardware performance degradation model by using the hidden Markov model according to the degradation performance index.
(3) And designing a system-level residual life estimation algorithm according to the system architecture.
(4) And (4) experimental design, wherein feature extraction and screening are carried out on the sample data set, the model is constructed, the residual life is predicted, a performance evaluation index is given, and the precision of the estimated value of the residual life is obtained.
2. The method of claim 1, wherein the step (1) of establishing and screening the degradation performance index comprises the following steps:
by using the signal processing method, the adopted method and the extracted characteristics can fully reflect the running state of the system. At present, the characteristic domains for representing the performance degradation indexes of the vibration signals are 3 types of time domains, frequency domains and time-frequency domains. The time domain features include an average value, a root mean square, an RMS, a kurtosis index, and the like. The frequency domain and time-frequency domain features are processed by methods such as wavelet transformation, empirical mode decomposition and the like. In addition, after feature extraction is performed on the system monitoring signal data, some feature vectors may be very correlated, resulting in overlapping information. Therefore, in order to reduce data redundancy, the experiment adopts a Principal Component Analysis (PCA) method to reduce the dimension.
Wavelet decomposition process:
in wavelet decomposition, a Hilbert space is decomposed into orthogonal sums of a plurality of subspaces according to different scale factors j, and further decomposed into a closed-packet subspace, and a wavelet packet base sequence is constructed to meet the following conditions:
Figure FSA0000167636490000011
wherein h and g are high-pass filter and low-pass filter coefficients respectively; k is 1, 2, …, j is a time position parameter and n is a frequency parameter.
After decomposition, reconstructing the decomposition coefficients of each frequency band to the 0 th layer by layer, and extracting signals in each frequency band, wherein if the data length of the original signal x (k) is N, the signal energy of each frequency band after decomposition is represented as:
wherein the content of the first and second substances,representing the magnitude of a discrete signal located in a subspace in a wavelet decomposition at a resolution j. The feature vector corresponding to the energy sequence can be determined therefrom.
The EMD decomposition involves the concept of eigenmode functions, which must satisfy the following two conditions:
1) the number of local extreme points and zero crossings must be equal or differ by at most one over the entire data range.
2) At any time, the average value of the upper envelope formed by the local maximum point and the lower envelope formed by the local minimum point is zero, that is, the upper envelope and the lower envelope are locally symmetrical with respect to the time axis.
EMD decomposition process:
1) all local maximum points of the signal data x (t) are found, all the local maximum points are connected by a cubic spline curve to form an upper envelope line, and all the local minimum points are connected to form a lower envelope line.
2) Mean value m of upper and lower envelope lines1Subtracting the mean value from the original signal data to obtain a new data sequence h1
x(t)-m1=h1
If h is1Two conditions for the eigenmode function are satisfied, then h1Is the first IMF component c decomposed1. Otherwise h1Instead, h is required to be1Repeating step 1) as the original signal, and continuing the screening until a first IMF component c is generated1
3) And (3) calculating:
r1=x(t)-c1
will r is1Repeating the steps 1) to 3) as the original signal, and so on to obtain n components of the original signal and a residual component r which can not be decomposed any morenNamely:
Figure FSA0000167636490000021
the PCA algorithm is summarized as follows:
raw data were normalized:
Figure FSA0000167636490000022
calculating a correlation coefficient matrix R ═ (R)ij):
Figure FSA0000167636490000023
Is the correlation coefficient of the original variable, rij=rji
Calculating a characteristic value and a characteristic vector:
solving a characteristic equation | λ I-R | ═ 0, solving characteristic values by a common Jacobian method, and arranging the characteristic values in the order of magnitude:
λ1≥λ2≥...≥λp≥0。
respectively finding corresponding characteristic valuesλiCharacteristic vector e ofi(i ═ 1, 2,. p.) such that | | | ei||=1。
Calculating the principal component contribution rate and the accumulated contribution rate:
rate of contribution
Figure FSA0000167636490000031
Cumulative contribution rate
Figure FSA0000167636490000032
3. The method according to claim 1, wherein the hidden Markov model is used to construct the system degradation model in step (2), and the method is divided into two steps, namely, the hidden Markov model is estimated by using Baum-Welch algorithm, and the hidden state sequence is predicted by using Viterbi algorithm to obtain the path of the software running state transition. Assuming that K inflection points exist in the hidden sequence, determining an optimal value of K by adopting a Bayesian Information Criterion (BIC), segmenting the degraded data according to the value of K, and respectively constructing a hardware performance degradation model for each segment; and secondly, a distribution function of the degradation amount can be obtained by fitting nonparametric distribution of degradation index data when the hardware is degraded under different software running states. Assuming that, given a software operating state Z (n-1) ═ i, Z (n) ═ j, the system performance degradation increment d (n) has a probability density function G (d | i, j, Θ)ij) I.e. by
D(n)~G(d|i,j,Θij),i,j∈I
Wherein, thetaijIs a parameter vector determined by Z (n-1) ═ i, Z (n) ═ j. Thus, of all possible I, j e I, the probability density set G ═ G (d | I, j, Θ) of the system performance degradation increments d (n) in each stateij),i,j∈I}。
Under the condition of determining the distribution of the degradation increment in different software running states, the future performance degradation amount of the system in different software running states can be easily analyzed, and the future performance degradation amount can be specifically determined by the following formula:
Figure FSA0000167636490000033
wherein n iscFor the current number of cycles, Y (n)c) To be at the moment of time
Figure FSA0000167636490000034
The amount of degradation in the performance index observed.
4. The method of claim 1, wherein step (3) designs a system level remaining life estimation algorithm. According to the observed system performance index data set, index data with obvious performance index change in different software running states are selected, and degradation models of different performance indexes in different software running states are constructed. At each time point tnCalculate the expected degradation value of performance E [ Y (n) for each state]Then, the remaining life of z (n) k, k e I in different software execution states is determined. Here, the remaining life is determined by the time when the system performance index first enters the failure region:
Figure FSA0000167636490000035
for a sufficiently large integer M, the remaining life of the different subsystems can be estimated:
Figure FSA0000167636490000036
the remaining life of the different subsystems is estimated separately. And estimating the system-level residual life of the software-hardware system by using the system structure.
5. The method according to claim 1, wherein the step (4) predicts the remaining life based on the constructed model, and derives the accuracy of the estimated value of the remaining life based on the evaluation index.
Obtaining a state sequence of hidden software operation according to the state transition probability by applying a hidden Markov model, and obtaining a hidden software operation state sequence according to a transition nodeObtaining the degradation amount of a hardware system in different states, obtaining that each section of data belongs to a normal distribution model by Matlab fitting, selecting T test to carry out goodness-of-fit test under the condition that the given significance level α is 0.05, and obtaining that the P values of the test are all larger than 0.05, so that the original hypothesis is accepted
Figure FSA0000167636490000041
The probability distribution function of (a) is:
Figure FSA0000167636490000042
and estimating parameters of the performance degradation model by using a nonlinear least square method to obtain normal distribution model parameters of each section of data. Given different time nodes, the remaining life of each piece of hardware can be estimated, assuming a failure threshold.
The evaluation index is measured by Root Mean Square Error (RMSE), mean error (AE) and Mean Square Percentage Error (MSPE) between the predicted residual life and the real residual life, and the smaller the value, the better the prediction result.
Root Mean Square Error (RMSE)
Figure FSA0000167636490000043
Mean error (AE)
Figure FSA0000167636490000044
Mean Square Percent Error (MSPE)
Figure FSA0000167636490000045
Wherein, yiDenotes true value, y'iIndicating the predicted value.
CN201810815534.8A 2018-07-20 2018-07-20 Method for estimating remaining life of software-hardware degradation system by considering software operation Pending CN110750848A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810815534.8A CN110750848A (en) 2018-07-20 2018-07-20 Method for estimating remaining life of software-hardware degradation system by considering software operation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810815534.8A CN110750848A (en) 2018-07-20 2018-07-20 Method for estimating remaining life of software-hardware degradation system by considering software operation

Publications (1)

Publication Number Publication Date
CN110750848A true CN110750848A (en) 2020-02-04

Family

ID=69275229

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810815534.8A Pending CN110750848A (en) 2018-07-20 2018-07-20 Method for estimating remaining life of software-hardware degradation system by considering software operation

Country Status (1)

Country Link
CN (1) CN110750848A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111832226A (en) * 2020-07-13 2020-10-27 中国南方电网有限责任公司超高压输电公司柳州局 IGBT residual life estimation method based on auxiliary particle filtering
CN113283533A (en) * 2020-06-12 2021-08-20 北京航空航天大学 Borrowable sample screening method and system for performance decline time series data
CN113791960A (en) * 2021-08-20 2021-12-14 西安理工大学 Method for predicting cloud server software aging based on EGGM model

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5465321A (en) * 1993-04-07 1995-11-07 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Hidden markov models for fault detection in dynamic systems
EP2172887A2 (en) * 2008-09-30 2010-04-07 Rockwell Automation Technologies, Inc. System and method for dynamic multi-objective optimization of machine selection, integration and utilization
CN104504296A (en) * 2015-01-16 2015-04-08 湖南科技大学 Gaussian mixture hidden Markov model and regression analysis remaining life prediction method
CN107729678A (en) * 2017-11-02 2018-02-23 中国科学院数学与系统科学研究院 A kind of satellite solar cell modeling of residual life and analysis method in orbit

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5465321A (en) * 1993-04-07 1995-11-07 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Hidden markov models for fault detection in dynamic systems
EP2172887A2 (en) * 2008-09-30 2010-04-07 Rockwell Automation Technologies, Inc. System and method for dynamic multi-objective optimization of machine selection, integration and utilization
CN104504296A (en) * 2015-01-16 2015-04-08 湖南科技大学 Gaussian mixture hidden Markov model and regression analysis remaining life prediction method
CN107729678A (en) * 2017-11-02 2018-02-23 中国科学院数学与系统科学研究院 A kind of satellite solar cell modeling of residual life and analysis method in orbit

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JEAN-BAPTISTE DURAND 等: "Software reliability modelling and prediction with hidden Markov chains" *
JINHEE PARK 等: "An Embedded Software Reliability Model with Consideration of Hardware Related Software Failures" *
刘河生 等: "隐马尔可夫模型的原理与实现" *
周建民 等: "基于隐马尔科夫模型的滚动轴承性能退化评估" *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283533A (en) * 2020-06-12 2021-08-20 北京航空航天大学 Borrowable sample screening method and system for performance decline time series data
CN113283533B (en) * 2020-06-12 2022-08-12 北京航空航天大学 Borrowable sample screening method and system for performance decline time series data
CN111832226A (en) * 2020-07-13 2020-10-27 中国南方电网有限责任公司超高压输电公司柳州局 IGBT residual life estimation method based on auxiliary particle filtering
CN111832226B (en) * 2020-07-13 2024-04-26 中国南方电网有限责任公司超高压输电公司柳州局 IGBT residual life estimation method based on auxiliary particle filtering
CN113791960A (en) * 2021-08-20 2021-12-14 西安理工大学 Method for predicting cloud server software aging based on EGGM model

Similar Documents

Publication Publication Date Title
CN107941537B (en) A kind of mechanical equipment health state evaluation method
CN107292446B (en) Hybrid wind speed prediction method based on component relevance wavelet decomposition
CN102208028A (en) Fault predicting and diagnosing method suitable for dynamic complex system
CN110750848A (en) Method for estimating remaining life of software-hardware degradation system by considering software operation
CN110175541B (en) Method for extracting sea level change nonlinear trend
CN112116198A (en) Data-driven process industrial state perception network key node screening method
CN114548199A (en) Multi-sensor data fusion method based on deep migration network
CN115098962A (en) Method for predicting residual life of mechanical equipment in degradation state based on hidden half Markov model
CN114861349A (en) Rolling bearing RUL prediction method based on model migration and wiener process
CN111882120A (en) Power load prediction method based on VMD-MQPSO-BPn network
CN115051929A (en) Network fault prediction method and device based on self-supervision target perception neural network
CN114563671A (en) High-voltage cable partial discharge diagnosis method based on CNN-LSTM-Attention neural network
CN112836967B (en) New energy automobile battery safety risk assessment system
Jiang et al. Multi-objective optimal placement of sensors based on quantitative evaluation of fault diagnosability
Zhang et al. Remaining useful life prediction of lithium-ion batteries based on TCN-DCN fusion model combined with IRRS filtering
CN111858343A (en) Countermeasure sample generation method based on attack capability
CN113740671A (en) Fault arc identification method based on VMD and ELM
WO2023231374A1 (en) Semi-supervised fault detection and analysis method and apparatus for mechanical device, terminal, and medium
Niu et al. A hybrid bearing prognostic method with fault diagnosis and model fusion
CN116521495A (en) System performance bottleneck detection method based on reinforcement learning
CN111222203A (en) Method for establishing and predicting service life model of bearing
CN115935814A (en) Transformer fault prediction method based on ARIMA-SVM model
CN114841063A (en) Aero-engine residual life prediction method based on deep learning
Zhang et al. A novel hybrid model based on EMD-BPNN for forecasting US and UK stock indices
CN113821401A (en) WT-GA-GRU model-based cloud server fault diagnosis method

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200204