[发明专利]一种载荷谱分布估计的拟合方法

说明书

一种载荷谱分布估计的拟合方法，包括：对载荷谱分布进行聚类分析；利用EM算法进行混合高斯参数估计；混合高斯拟合。本发明的一种载荷谱分布估计的拟合方法，简单实用、操作方便、易于编程的优点，除此之外可以通过调整高斯分项的个数可以对任意载荷循环分布进行拟合，大大提高的该方法的灵活性与拟合精度，为农机产品可靠性设计分析与验证以及制造材料与工艺提供理论支撑。

技术领域

本发明涉及一种载荷数据的处理。特别是涉及一种用于联合收割机田间作业载荷数据处理的一种载荷谱分布估计的拟合方法。

背景技术

我国耕地地形及地貌复杂多样、作物种植结构及区域差异较大、作物种类及品种数量众多，亟需田间作业可靠性高、工况环境适应性强、作业质量满足要求的农机装备。但由于农机装备在设计研发、生产制造及使用维护过程中缺乏完善的田间作业基础数据支撑以及可靠有效的田间试验验证技术和方法，产品质量、田间作业可靠性和环境适应性始终难以保证。

目前，对于联合收割机关键零部件的可靠性分析常用的方法是根据疲劳寿命曲线研究疲劳寿命与材料强度之间的关系。然而，农机装备在田间工况下所受的载荷是随机的，具有显著的复杂非线性特征。基于此，需要研究农田复杂开放工况环境下作业载荷特点及载荷分布，以达到精确再现整机及关键零部件在服役寿命内所承受的载荷特性。描述随机载荷信号常用的方法是建立一种概率密度模型对载荷参数进行拟合。

发明内容

本发明所要解决的技术问题是，提供一种既简单实用、操作方便，又能实现对任意载荷参数分布进行拟合的载荷谱分布估计的拟合方法。

本发明所采用的技术方案是：一种载荷谱分布估计的拟合方法，包括如下步骤：

1)对载荷谱分布进行聚类分析；

2)利用EM算法进行混合高斯参数估计；

3)混合高斯拟合，所述混合高斯拟合的数学表达式为：

其中

式中：α_k为高斯分项p_k的权值；d为随机变量的维数，对于二维数据空间S＝(S_a,S_m)，d＝2；μ_k为第k个高斯分项的均值向量；∑_k为第k个高斯分项的协方差矩阵；S_a表示载荷循环的幅值；S_m表示载荷循环的均值。

步骤1)包括：

(1)对原始载荷-时间信号进行双参数雨流统计计数处理，生成二维数据空间S＝(S_a,S_m)，式中S_a表示载荷循环的幅值；S_m表示载荷循环的均值；

(2)对二维数据空间S＝(S_a,S_m)，选取K个点作为K个聚类的初始点，计算二维数据空间所有点到K个聚类初始点的距离，将每个数据点归到距离它最近的那个聚类初始点所形成的类别中；

(3)用公式计算出K个类别新的聚类初始点，重新计算二维数据空间所有点到K个新的聚类初始点的距离，n_w为第w个类别数据点的个数，μ^w为第w个类别的均值向量，s_w表示第w个类别中所有数据点；1≤w≤K；

(4)利用τ_nw||s-μ^w||²求出所有数据点分别到每一个类别的聚类初始点的距离，对得到的所有距离用公式进行求和，式中，J为求和结果，N为二维数据空间数据点个数，τ_nw表示数据点n被归类到第w类别的时候为1，否则为0；