1.Introduction

An aircraft hydraulic power supply system provides high-pressure fluid for the actuation system,braking system,landing gear system,and other sub-function systems.As the power source of an aircraft hydraulic system,an aviation piston pump’s performance in fluences flight safety directly.Therefore,an aircraft prognostics and health management(PHM)system appears to keep high reliability and long life of an aerial piston pump.In PHM technologies,an accurate estimation of the remaining useful life(RUL)is the most difficult issue because it is related to the failure physics and stress spectrum imposed on a hydraulic pump.Since the structure of an aviation piston pump is very complicated,its failure generation and development are comprehensively affected by inner frictional pairs with uncertain characteristics.Statistically,an aviation piston pump shows variant degradation paths under divers iform operating conditions,which unavoidably brings about a great deal of uncertainties and difficulties in the analytic solution of the RUL.Although the life of a certain type of aviation piston pump can be obtained through tens of thousands of hours of experiments under a fixed spectrum,it is difficult to give the exact RUL under an arbitrary condition.A prognostic estimation method of RUL is imminently needed,which will highly benefit the reduction of costs by providing the possibility to de fine predictive maintenance strategies and prolonging useful life.

Over past decades,a lot of research has been conducted in estimating the RUL of machinery.The methods can be generally divided into two categories1:data-driven methods and model-based methods.Typical data-driven methods based on machine learning are artificial neural networks(ANNs)and the hidden semi-Markov model(HSMM).Zangenehmadar and Moselhi2used an ANN to assess the RUL of pipelines successfully in which more than 80000 groups of data were used for training.Dong et al.3–5applied the HSMM function for machine health prognosis and veri fied the method by using data from a real hydraulic pump health monitoring application case study.In fact,the data used for training were far more than the sample size of a certain type of aviation piston pump.Several thousand of hours were taken to obtain only one set of lift-cycle data of an aviation piston pump.6He et al.7presented a health monitoring and prognostic method using the PSOSVM to predict the RUL for an axial piston pump.The small sample problem is what makes an aviation piston pump distinguish from traditional machinery.Model-based methods like physics-of-failure(POF)and filter-based methods are also widely utilized for prediction of component life.Liu et al.8developed a failure physics model for the creep fatigue of a piston,and the degradation mechanism was analyzed for prognosis.Lamoureux et al.9defined a health indicator to describe the degradation of an aircraft engine fuel pumping unit by using a linear regression method.To build a model by POF,the mechanism needs to be known firstly.Thereby,this kind of method can hardly be used in a complicated system like an aviation piston pump whose failure mechanism is still under research.Filter-based functions take advantages in combining system models with experimental data.For linear systems with Gaussian noise,Kalman filter(KF)is a commonly used prognostic technique,and its effectiveness has been proven in many works.10–13Extended KF(EKF)and unscented KF(UKF)methods are modified KFs to cope with non-linear systems while limitation isshown in somesy stems with high nonlinearity.

According to the fact that the sample size of an aviation piston pump is very small,a data-driven method would not be a good choice for life prediction.Among model-based methods,particle filter(PF)has shown great advantages as an efficient prognostics tool in handling the uncertainty and noise affecting measurements.14A dual-particle- filter method was used to estimate the state of charge for power Li-ion batteries.15To address the particle impoverishment problem,a modified particle filter,named intelligent particle filter(IPF),was proposed by Yin and Zhu.16Miao et al.17introduced an improved PF algorithm–unscented particle filter(UPF)into battery RUL prediction,and the analytical results showed that UPF could predict the actual RUL with an error less than 5%.Zio et al.18–21improved the method a lot by applying PF functions in different degrading systems,and the framework to estimate the RUL of nonlinear components provides ideas for the prognosis of pump systems.However,to adopt the method,a degrading model should be built.Some parts of a piston pump has been modeled like a friction mechanism model of oil between the valve plate and the cylinder block in axial piston pumps22and a wear mechanism model of friction pairs23while a physical model that can be used for prediction has not been proposed.24,25For a system that is difficult to be modeled,an empirical model or a model built according to historical data is used to describe the degrading process.Fagogenis et al.26proposed an auto-regressive(AR)model with an RUSBoost classifier,and a CMAPSS dataset provided by the NASA AMES research center was used to verify the performance of computing the RUL of turbofan engines.A gray prognostic model based on the Markov process was used for a gas turbine compressor’s state estimation.27These kinds of models compromise the merits of data-driven models and have less demand in data quantity.

To address the problem mentioned above,a novel particl efilter based prognostic method for an aviation piston pump’s RUL prediction,named adaptive-order particle filter(AOPF),is proposed in this work.Though the wear mechanism of the pump is complex,the return oil flow has been proven to be a logical characteristic of the pump internal wear status.23The degradation of the pump is re flected by the increase of the flow while the flow shows a non-smooth characteristic that violent fluctuation occurs randomly.The main task is to handle the uncertainties which are classi fied in three categories.28The first one is the uncertainty of future degradation progress of an aviation piston pump which may be caused by unknown load spectrum and random environment factors and will result in different degrading paths,because of which model should be modi fied timely in order to maintain accuracy.Secondly,in modeling,an incomplete data set and some sufficient but not necessary assumptions and simpli fications may prevent a prognostic model from precision.A model with biased parameters will in fluence the performance in prediction.Thirdly,data collected by sensors and acquisition systems are often accompanied by measurement noise.29–31

The order number influences the accuracy of the prediction.Although the AOPF only takes the model which fits best,the relationship between the order number and the performance of that is expected to be known so that an optimization could be applied for a reduction of time cost.By fixing the order number of the AOPF,the results of a 7-order model and a 20-order model are shown in Fig.10,and the performance is listed in Table 1.The 20-order model shows a better performance than that of the 7-order model.The errors of the 7-order model and the 20-order model are 68 hours and 55 hours,respectively.As it is analyzed above,a higher-order model may fit the actual data better.

2.Adaptive-order particle filter based prognostics

2.1.Degrading characteristics of aviation piston pumps

There are four main friction pairs in a typical aviation piston pump:pair of the cylinder block and the valve plate,pair of the swash plate and the slipper,pair of the slipper and the piston,and pair of the piston and the cylinder block.Under normal circumstances,oil film exists between the friction pairs.It is shown that wear occurs and becomes serious with an increase of the serving time.22The wear process is accompanied by emerging abrasives.The accumulative amount of wear particles can indicate the degree of wear to some extent.However,the sizes of the abrasives are usually too small to be detected.There is no sensor installed to detect wear particles as well.Consequently,abrasives’increase cannot be utilized as a degrading index of the pump.The leakage caused by wear can also lead to a decrease of the outlet flow rate.In an actual hydraulic system of an aircraft,the outlet flow is controlled to be at a stable value,which means that the decline of the outlet flow cannot be monitored.In fact,the leakage oil is usually discharged out of the pump through the return oil outlet.By collecting the data of the return oil flow,the leakage of the pump can therefore be reflected.

(4)表达式可以出现于等号的左右两边，从而构成一句完整的句子。[注][美]斯鲁格：《弗雷格》，江怡译，北京：中国社会科学出版社，1989年，第267页。

The mechanism of wear is very complicated that it cannot be described only by establishing the model of each friction pair.As is mentioned above,abrasives occur randomly during the wear.The abrasives produced by one friction pair will transmit to another friction pair.A mutual reinforcement of wear then happens which leads to the uncertainties of the degrading path.It is obvious that the trend of the return oil flow is increased for wear can only become more serious.What should be taken care of is that short-term declines often appear in the increasing trend,which cannot be explained by measurement noises.The phenomenon has been described as a process with independent incensement by Wang et al.23On this basis,a pump can be modeled with a system noise whose distribution obeys the process with independent in censement.

2.2.Basic algorithm of particle filter based prognostics

Particle filter has become an attractive state estimation method for its capability of accounting for the randomness of the process and the noise affecting measurements.The main idea of PF is that particles are used to describe a probability distribution instead of an analytical expression.To de fine a system,consider a state sequence x0:k={x0，x1，...，xk}and its corresponding measurement sequence y0:k={y0，y1，...，yk},where k denotes a time metric.The state space model is commonly defined as

where f and h represent the state transition function and measurement function,respectively.xk∈Rnxis the state to be estimated,and yk∈Rnyis the observation,where nxand nyare the dimensions of the states and the observations.uk∈Rnxand vk∈Rnyare independent,identically distributed noises with known probability densities.

According to the Bayes’rule,the posterior probability density can be described as

where the requirement is that p(y)＞0,and p(x)is the probability of x.The state estimation problem is to solve the probability density function(PDF)contains all the information about the state xk,which is inferred from the observationsis the initial distribution of the states and is known previously.The PDFis obtained recursively from the prediction step and the update step.

In the prediction step,the system model is used to obtain the prior PDF of the state at time k by the Chapman–Kolmogorov equation as follows:

in whichis defined by the system model Eq.(1).has been obtained by the update step in last loop,andis known previously.

which means that the importance weight is the ratio between two distributions.It can be deduced as

whereis the normalizing constant,

and p(yk|xk)can be obtained by Eq.(2).

世界范围内资本账户开放的浪潮大致经历了80年代各国争相放开，90年代发展中国家陆续发生货币危机从而产生质疑，再到目前注重时序性和渐进性政策三个阶段。历史经验表明，资本项目开放的收益多体现在发达国家，而较少见诸发展中国家，而一些经济、金融欠发达的中小发展中国家，往往因为缺少资本管制的屏障保护而成为国际投机资本攻击的牺牲品。因此，如何充分享受资本项目开放的收益成为政策面临的重大挑战。前国际货币基金组织总裁Camdessus提出，亚洲金融危机的根源在于既定的宏观经济环境和制度背景下，资本项目开放措施与经济金融等改革的衔接出现了问题[1]。

For the system defined by Eqs.(1)and(2),f and h are usually nonlinear,which leads to a difficulty in obtaining an analytic solution for the posterior distributionA PF approximates it with a set of samples (particles)where N is the particle number.The initial particles are drawn from p(x0).To address the problem of sampling from the posterior distribution,importance sampling is employed.Samples can be easily drawn from an importance distribution q(x0:k|y0:k),and the importance weight for each particle can be calculated by

60例甲状腺肿瘤患者经腺叶切除术治疗,除2例全切患者出现暂时性的甲状旁腺功能低下,2例出现皮下出血的情况之外,其余患者均无喉返神经、喉上神经损伤和永久性甲状旁腺和甲状腺功能减退和低下的问题。经随访,无患者出现再发性肿瘤。

wheredenotes the weight of particle i at time k.Thereby,the posterior distribution p(xk|y0:k)can be approximated by particles as

where δ is the Dirac function.

3.深港股市的动态相关系数具有显著的时变特征，并且整体上呈现“N”型趋势。且“深港通”启动初期，深成指收益率和恒生指数收益率的动态相关系数呈现波动下降趋势，随着政策的深入推行逐渐呈现跃升趋势，说明“深港通”对深港两地联动效应的加强作用需要一定的过渡期才能实现。

2.3.Adaptive-order particle filter

2.3.1.Particle filter with a high-order model

Note that the state space model given by Eqs.(1)and(2)is described by a Markov model of first order.Actually,the system model may be given by

in which the system is described as a p order model.The PF framework still takes hold and is proven as follows.

In the prediction step,the prior PDF of the state at time k is

基于中国能源强度在地理空间上集聚分布的客观事实，本文借助探索性空间数据分析方法进一步考察变量的空间自相关性。对能源强度的全域自相关性进行检验，如表3所示，ln ei、ln gdp、ln fdi的全局Moran's I指数值均显著为正，即在全域范围内存在显著的空间正自相关性。因此，在建立模型时应充分考虑地区观测值之间的空间溢出效应。

wherecan be obtained by Eq.(9),and p(xk-p|y0:k-p)is obtained by the update step in loop k-p.p(x0)=p(x0|y0:0)is known previously.

In the update step,with a new observation yk,the posterior PDF of the state can be calculated by Eqs.(5)and(6).Thus,the Bayesian solution is capable to be used in a high-order model.The proof of the applicability of importance sampling is as follows.

The original definition of the importance weight of a particle is

where σ(x0:t)is the joint probability distribution as

尼康为这款产品开发和制造的快门单元，其快门叶片是由碳素纤维和复合材料制成。对组装好的相机快门进行循环测试证明，D700相机的快门在极端环境中经过约150000次的测试。自我诊断式快门可持续地监控并保持快门的精确性。

In the update step,a new observation ykis available and can be used to update the prior PDF of the state via Bayes’rule as follows:

Substitute Eq.(12)into Eq.(13),and then

Take an optimal ofasand then Eq.(14)can be described as

For each particle,the importance can be obtained by

Then the importance sampling process is proven to be effective in a high-order model.The res ampling step changes the importance sampling weight to be equal by adjusting the number of the particles like

where j is the new particle index after resampling,and xij∈ {xi，i=1，2，...，N}.Substituting Eq.(17)into Eq.(16)yields

where u0is the obeyed uniform distribution.(2)For k=p，p+1，···

By a condensation filter,32take a sub-optimal of q asThe weighting function would be like

The posterior distribution can be approximated by Eq.(8).

2.3.2.Adaptive-order particle filter prediction framework

A general PF uses a fixed dynamic model to estimate system states.In fact,it is difficult to build a dynamic model for a piston pump analytically.The complicated system’s degrading mechanism is still unknown whose model is usually built empirically or built by data-driven models.Models built like these might not be guaranteed to accord with the actual degrading process,which means that updating the model is necessary under such a circumstance.

According to the return oil flow data shown in Fig.1,the curves of return oil flow are not monotonously increasing so that sometimes a decline appears and lasts for a few hours.The most importance in a PF is to ensure the accuracy of the prior distribution.However,a first-order data-driven model could not deal with such un smooth data.Take the hypothetical sequence shown in Fig.2 for an example.If the dynamic model is first-order,the estimated value is more likely to show at xk+1.If a second-order model is taken into consideration,seems to be a more reasonable value.A highorder model could contain more recent information.A shortterm trend by the latest states is more likely to reflect the possible states of next time steps.Note that when the parameters of all the high-order variables are equal to zero,the model is a first-order model.

Based on the theory mentioned above,we propose an adaptive-order particle filter method,which can be applied to a dynamic model with a current time state related to the states at not only last time step but also several time steps previously which can be described by Eqs.(9)and(2).The adaptive order means that the order number of the dynamic model is updated at each time step adaptively.The main difference between a general PF and a high-order particle filter is the prediction step during that to predict the state at time k,if the order of the dynamic model is Ok=p,then the states of the particlesshould be stored without aresampling step.{O0:k}is the order sequence of the updated dynamic model.Thus a recursive function instead of a onestep function is applied to sampling particles from the suboptimal of q,which is

An empirical or data-driven model is built according to the historical data firstly.To fit a new observation sequence,parameters of the model are usually modi fied accordingly.The model may already be accurate enough to describe the trend until time k,when a new data comes.Without changing the form of the original model,the adaptive-order framework re-determines the order number of the model by calculating the accumulated errors.If a higher-order model is calculated to be more fitted,the order number is then updated.By rebuilding the model in a higher-order form,the short-term trend is likely to be more accurate.The adaptive-order framework is shown in Fig.3.

Fig.1 Return oil flow.

Fig.2 Hypothetical sequence.

Fig.6 shows the short-term prediction results by the GM(1,1)and the proposed AOPF method.Both of the algorithms show good performance in the one-step prediction.Each point of the curves stands for a prediction obtained by the prior data.The gray method fits well while the AOPF shows a better result which is demonstrated in Fig.7.The red curve and the blue curve indicate the mean square errors(MSEs)of gray prediction and AOPF prediction,respectively.To get Fig.7,each point is calculated by accumulating all the MSE of the previous points.Obviously,the MSE of the AOPF prediction is lower than that of the gray prediction after 47 time steps.Raw data of the return oil flow illustrate a large fluctuation at the beginning during that from 20 hours the return oil flow continues to fail which impacts the performance of the method.However,after a few steps,the AOPF performs a good anti-interference characteristic.The blue curve keeps lower than the red curve,which means that the error of the AOPF prediction is smaller than that of the gray prediction and the noise is depressed after enough prior knowledge accumulated.

李栋，山西太谷县人，1969年出生，中国诗歌学会会员，山西省作协会员。1992年发起创办太阳谷文学会。2010年起触网写诗，组诗偶见于《人民文学》《诗刊》《星星》等处。

(1)Initialization step. The initial modelis built by the historical data,in which p0is the order number of the model.The number of particles N is set.The particles and their weights are initialized as

转变临时应对的想法，完善事件舆情的预警体系，使得事件到来时有相应的应对措施引导舆论。事件管理如今应该常态化，与其在事件发生时手忙脚乱最后事件也不能够妥善解决，不如早作打算，做好事件预警工作。随着自媒体的产生，在如今自媒体也能够对突发性公共卫生事件舆情进行预警和检测。

(a) Calculate the prior estimation by the current modelwhich is described by the particles:is the updated state model.

五是垦地合作取得新进展。我们分别与湛江、茂名、阳江、揭阳市政府签订战略合作框架协议，在统一优化整合区域资源、加速生产要素合理流动、加快发展现代农业、推进农垦土地开发利用、建设产业园区、加强农场基础建设、推进城镇化进程、发展公共事业等方面开展深度合作，促进农垦与地方经济社会融合发展。垦区各产业集团结合实际采取订单收购、“公司+农户”等方式，带动农户发展种养业，促进周边农村经济发展和农民增收致富。

Fig.3 Adaptive-order framework.

Fig.4 Flow chart of the AOPF.

(c) Normalize the weight of each particle:

(d) Resampling the particles by the weight.

(e) Update the order number of the model and get a new current model

(3)Output the predicted states sequence:whereis the estimated value of

The parameters of the initial model will not be as important as the form of the model for that the parameters from the historical data will be replaced by the parameters for the new data sequence by the model updating procedure.The model form f0:kis not modified and will be used in the prediction step.When there is no observation arriving at time k+1,the model is fixed,and the weights of the particles will not change.To get the RUL,the current dynamic model is the only dependence,which means that the recursive state in the future xk+tis obtained by

where t is the step of time that xk+t≥Threshold(return oil flow).The distribution of a piston pump’s RUL TRULis defined as

3.Example verification

In this section,data collected by a full-life experiment for a certain type of piston pump are used to demonstrate the forecasting performance of the proposed AOPF algorithm.The return flow data were measured by a turbine flow sensor installed on the return oil pipe of a pump.The work operation of the pump was full output flow by which an accelerated wear process was adopted so that the cost could be reduced.The flow data were recorded every hour by the same measurement system.

Based on the historical data,a gray model GM(1,1)is commonly used as the data-driven model for a prediction of the piston pump.At least four points are needed to initial the model which means that the minimum order number of a GM(1,1)is 4.By accumulating the raw sequence X(0)={x0，x1，...，xk},the accumulated sequencek=0，1，2，...is used to fitby

where

and

The data shown in Fig.5 are the experimental results of the return oil flow of a certain pump.When the return oil flow reaches 2.8 L/min,the wear condition of the pump is considered to be the threshold of the failure of wear according to the design parameters.The trend of the return oil flow is increasing except some parts of reduction,which means that the leakage of the pump is increasing and gradually leads to an incapability for the pump to provide high-pressure hydraulic power for the actuation system.

Fig.5 Experimental results of the return oil flow.

In the proposed model,the adaptive-order model is embedded in the PF measurement equation.The flow chart of the proposed updated-order particle filter is depicted in Fig.4,and its main steps are as follows:

The adaptive-order number of each step is displayed in Fig.8.The order numbers show a polarization that most of the order numbers are very small and large order numbers usually appear and last for some time.Comparing Figs.6–8,a strong connection appears between the large order number and the decline of the raw data.When the trend is rising,the adaptive order number is stabilized at a low level.To adapt the unusual trend,the AOPF matches with a higher-order model adaptively.

The characteristics of the adaptive-order numbers shown in Fig.8 provide an optimization strategy.In a rising stage,the best order number occurs between 0 and 40.Then if the trend of the sequence is determined,the computation burden can be largely declined.The strategy shows a greater significance with an increase of the time.In a later period,the resolving time is impacted mainly by the time on MSE calculation.The strategy could keep the computing time rise steadily.

However,many functions can provide a high precision in the short-term prediction.What challenges is the long-term prediction.In order to verify the performance of a long-term prediction,the raw data are divided into two sets.One is used for model adaption and the other for proving.Taking the first 600 points as the training data,a gray model sequence is generated.Meanwhile,a gray model initialized by the first 4 points are updated adaptively by the AOPF.Long-term prediction results by GM(1,1)and the proposed AOPF method are shown in Fig.9.The actual failure time of the pump is 1186 hours.By gray forecasting with GM(1,1),the time that reaches the threshold is 1103 hours,which means that the error is 73 hours.Whereas the life predicted by the AOPF is 1144 hours and the error is 42 hours,which means that the error is reduced by 42.5%.Note that the precision of the AOPF is affected by the number of particles.More particles lead to a heavier computation burden.In the proposed stage,the number of the particles is 300.In the long-term prediction,the measurement error and the system error are accumulated in the datadriven model.The AOPF can promote the long-term prediction accuracy through filtering the noise of each measurement point.

全球移动通信协会（GSMA）首席执行官洪曜庄看到，在新技术的浪潮中，传统媒体转型势在必行。未来的数字世界，移动设备将会是数字变革中的核心，传统媒体要做好快速迁移的准备。没有一成不变的格局，也没有一劳永逸的创新。在这个“新”与“旧”快速更迭的时代，谁会是下一个领跑者，是一个待解的谜题。拥抱变化、随时出发，是人们给出的题解，当然，也有人在用想象让人们走在时间之前，感触未来。

Fig.6 Short-term prediction comparison.

Fig.7 MSEs of return oil flow comparison between grey prediction and AOPF prediction.

Fig.8 Adaptive-order number of the AOPF.

Fig.9 Long-term prediction comparison.

The rest of this paper is organized as follows:Section 2 describes the particle filtering framework for estimating the RUL and the proposed AOPF technique is described in detail;in Section 3,the application of the method is verified by an experimental test,and a comparison with traditional methods is discussed;in Section 4,some conclusions and remarks are drawn.

Fig.10 Different order models’prediction results comparison.

Table 1 Performances of 7-order and 20-order models.

Type 7-order model 20-order model MSE of return oil flow(L/min) 0.0091 0.0060 RMSE of return oil flow(L/min) 0.0953 0.0775

Fig.11 Different order number models’MSEs of the return oil flow.

For an intuitive understanding of the relationship between the order number and the performance,Fig.11 shows the MSEs of different order number models at time 600 h.The order number is from 7 to 60.With the order number increasing,the trend of the MSE is decreasing.The lowest point of the curve occurs near 30 after which the curve shows a fluctuation around 5.5×10-3L/min.It is not that the higher the order number is,the lower the error is.By taking an appropriate order,the error can be reduced by more than 50%which will increase the precision of the RUL prediction accordingly.

4.Conclusions

An AOPF method is proposed in this paper to improve the prediction accuracy of an aviation piston pump’s RUL.In the AOPF,the accuracy of long-term prediction is promoted via changing the order number of a model adaptively.With this strategy,the model is modified with a new observation arriving at every time step,so that the information of the new point and the empirical knowledge can be well fused.Compared to a general PF,the state prior distribution in the AOPF is calculated recursively.In the AOPF,the order number is determined by the minimum mean square error between estimated states and raw data.It shows that the MSE between future data and predicted states is reduced by the proposed method.The experimental results indicate that the AOPF reduces the error by 42.5%in contrast with the traditional grey forecasting method.

两种针刺工艺路线主要区别在于进入预刺之前，PPS含基布滤料会作为夹心层铺在净气面（底层）上，然后迎尘面 （面层）纤网会铺在PPS基布上，形成三层结构，PPS无基布滤料在进入预刺之前不会有放卷PPS基布在两层纤网中间，针刺环节只会发生纤维与纤维之间的抱合。由于无基布滤料没有基布作为加强层，在针刺工艺参数设置上，会适当降低针刺密度，增加针刺深度，增加针刺机工作辊速度，增加梳理机喂入量，以降低对纤维的损伤，同时保障纤维之间的抱合力。

However,there are still some improvements that could be made on the AOPF.Firstly,the experimental results show that the overall trend of the MSE is decreasing.When the order number is near 30,the MSE reaches a minimum value and after that,the MSE begins to stabilize.To get a minimum MSE,the MSEs of all order number models should be computed,which leads to a heavy computation burden.A method of optimization may be applied to improve the adaptive strategy.

Acknowledgements

This study was co-supported by the National Natural Science Foundation of China(Nos.51620105010,51575019),National Basic Research Program of China(No.2014CB046400),and Program 111 of China.

对于生产产品质量较高的农户来说，要不断推进农业合作社的改革创新，加强合作社的标准化管理生产模式和品牌建设。拓展农产品精深加工产业，不但可以提升产品质量和安全，产品的销路也能够得到保证。不同地区有不同的农业对策，根据不同地区的气候、地理条件、发展情况合理制定不同的生产销售模式。同时合作社管理者需要根据不同的经济发展条件及时变更销售策略，及时为农户找到适合的销路，不浪费管理资源，做到实事求是。根据各地的经济差异，规划农业产业化模式的普遍适用性也是有必要的，但同时农户又具有异质性，产品质量是销售的关键因素。要把产品质量当作农户异质性的标志，这是推进合作社工作发展的必要途径。

References

1.Sun J,Li H,Xu B.Prognostic for hydraulic pump based upon dctcomposite spectrum and the modi fied echo state network.Springerplus 2016;5(1):1293.

2.Zangenehmadar Z,Moselhi O.Assessment of remaining useful life of pipelines using different artificial neural networks models.J Perform Constru Facil 2016;30(5):04016032.

3.Dong M,Peng Y.Equipment PHM using non-stationary segmental hidden semi-markov model.Robot Comput Integr Manuf 2011;27(3):581–90.

4.Dong M,He D.Hidden semi-markov model-based methodology for multi-sensor equipment health diagnosis and prognosis.Eur J Oper Res 2007;178(3):858–78.

(b) When a new observation ykarrives,update the weight of each particle:

高职院校相对比较封闭，创新创业项目的开展应该在开放的环境中进行。创新团队外出交流，不仅可以学习到新的理念、技术，更可以认识技术同行，为自身项目完善、发展提供更多机会。外出参加的活动主题可以与创新团队的研究目标有关，也可以与学生专业相关，既可以是创新创业类活动，又可以是创新创业专项训练。既可以走出去，又可以引进来。创客空间可以利用自身场，举办讲座、沙龙等活动。

5.Dong M,He D.A segmental hidden semi-Markov model(HSMM)-based diagnostics and prognostics framework and methodology.Mech Syst Sig Proc 2007;21(5):2248–66.

6.Geng Y,Wang S,Zhang C.Life estimation based on unbalanced data for hydraulic pump.IEEE/CSAA international conference on aircraft utility systems;2016 Oct 10–12;Beijing,China.Piscataway(NJ):IEEE Press;2016.p.796–801.

7.He Z,Wang S,Wang K,Li K.Prognostic analysis based on hybrid prediction method foraxialpiston pump.IEEE international conference on industrial informatics;2012 Jul 25–27;Beijing,China.Piscataway(NJ):IEEE Press;2012.p.688–92.

8.Liu Y,Ji G,Li J,Zhao D.Physics-of-failure based approach for predicting reliability of mechanical products.International conference on reliability,maintainability and safety;2014 Aug 6–8;Guangzhou,China;2014.p.878–81.

9.Lamoureux B,MasséJR,Mechbal N.An approach to the health monitoring of a pumping unit in an aircraft engine fuel system.Proceedings of first European conference of the prognostics and health management society;2012 May 23–25;Beijing,China;2012.p.1–7.

10.Chen J,Ma C,Song D,Xu B.Failure prognosis of multiple uncertainty system based on kalman filter and its application to aircraft fuel system.Adv Mech Eng 2016;8(10).

11.Qian Y,Yan R,Hu S.Bearing degradation evaluation using recurrence quantification analysis and kalman filter.IEEE Trans Instrum Meas 2014;63(11):2599–610.

12.Mba D,Reuben LCK.Diagnostics and prognostics using switching kalman filters.Struct Heal Monit 2014;13(3):296–306.

13.Liu Y,Dong K,Wang H,Liu J,He Y,Pan L.Adaptive gaussian sum squared-root cubature kalman filter with split-merge scheme for state estimation.Chin J Aeronaut 2014;27(5):1242–50.

14.Zio E,Peloni G.Particle filtering prognostic estimation of the remaining useful life of nonlinear components.Reliab Eng Syst Saf 2011;96(3):403–9.

15.Liu X,Chen Z,Zhang C,Wu J.A novel temperature-compensated model for power li-ion batteries with dual-particle- filter state of charge estimation.Appl Energy 2014;123(3):263–72.

16.Yin S,Zhu X.Intelligent particle filter and its application to fault detection of nonlinear system.IEEE Trans Industr Electron 2015;62(6):3852–61.

17.Miao Q,Xie L,Cui H,Liang W,Pecht M.Remaining useful life prediction of lithium-ion battery with unscented particle filter technique.Microe lectron Reliab 2013;53(6):805–10.

18.Rigamonti M,Baraldi P,Zio E,Astigarraga D.Particle filter based prognostics for an electrolytic capacitor working in variable operating conditions.IEEE Trans Power Electron 2015;31(2):1567–75.

19.Hu Y,Baraldi P,Maio FD,Zio E.A particle filtering and kernel smoothing-based approach for new design component prognostics.Reliab Eng Syst Saf 2015;134:19–31.

20.Baraldi P,Mangili F,Zio E.Investigation of uncertainty treatment capability of model-based and data-driven prognostic methods using simulated data.Reliab Eng Syst Saf 2013;112(2):94–108.

21.Rigamonti M,Baraldi P,Zio E,Astigarraga D,Galarza A.A particle filtering-based approach for the prediction of the remaining useful life of an aluminum electrolytic capacitor.2nd European conference of the prognostics and health management society(PHME);2014 Jul 8–10;Nantes,France;2014.p.8–10.

22.Han L,Wang S,Zhang C.A partial lubrication model between valve plate and cylinder block in axial piston pumps.Proc Inst Mech Eng Part C J Mech Eng Sci 2015;229(17):3201–17.

23.Wang X,Lin S,Wang S,He Z,Zhang C.Remaining useful life prediction based on the wiener process for an aviation axial piston pump.Chin J Aeronaut 2016;29(3):779–88.

24.Guan C,Jiao Z,He S.Theoretical study of flow ripple for an aviation axial-piston pump with damping holes in the valve plate.Chin J Aeronaut 2014;27(1):169–81.

25.Guan C,Jiao Z.Modeling and optimal design of 3 degrees of freedom helmholtz resonator in hydraulic system.Chin J Aeronaut 2012;25(5):776–83.

26.Fagogenis G,Flynn D,Lane D.Novel rul prediction of assets based on the integration of auto-regressive models and an rusboost classifier.2014 IEEE conference on prognostics and health management(PHM);2014 Jun 22–25;Cheney,Washigton,USA.Piscataway(NJ):IEEE Press;2014.p.1–6.

27.Zhou D,Yu Z,Zhang H,Weng S.A novel grey prognostic model based on Markov process and grey incidence analysis for energy conversion equipment degradation.Energy 2016;109:420–9.

28.Zio E,Maio FD.A data-driven fuzzy approach for predicting the remaining useful life in dynamic failure scenarios of a nuclear system.Reliab Eng Syst Saf 2010;95(1):49–57.

29.Wang SP,Ma QX,Lu CQ.State recognition based on wavelet packet norm entropy for aircraft hydraulic pump.International conference on fluid power and mechatronics;2015 Aug 5–7;Harbin,China;2015.p.1318–23.

30.Lee J,Wu F,Zhao W,Ghaffari M,Liao L,Siegel D.Prognostics and health management design for rotary machinery systems—reviews,methodology and applications.Mech Syst Sig Proc 2014;42(1–2):314–34.

31.Mare JC,Vinson G,Prado T,Combacau M.Modelling and simulating the pump of an aerospace electro-hydrostatic module for fault detection and identification purposes.ASME/BATH 2014symposium on fluid power and motion control;2014 Sep 10–12;Bath,United Kingdom;2014.p.V001T01A2–VT01A2.

32.Isard M,Blake A.Condensation—conditional density propagation for visual tracking.Int J Comput Vision 1998;29(1):5–28.