基于PCA的故障诊断方法(matlab)

该博客介绍了PCA(主成分分析)在故障诊断中的应用,通过T2和SPE统计量来识别系统故障。首先解释PCA的原理和数据预处理步骤,接着详细阐述如何计算协方差矩阵、特征值和特征向量,并选择合适的特征进行降维。然后,定义了T2和SPE统计量的限值计算方法,以及如何利用这些统计量判断系统是否发生故障。最后,通过Matlab代码展示了整个过程,并分析了贡献图以定位故障源。

1. PCA原理分析

PCA的原理主要是将原始数据进行降维。其具体工作原理参照:CodingLabs - PCA的数学原理

2. 数据预处理

 训练数据集(只有正样本)为X_{n*m}维数据,即有n个采样值,每个采样值有m个特征。

\begin{pmatrix} x_{11} &x_{12} & \cdots & x_{1m}\\ x_{21}& x_{22} &\cdots &x_{2m} \\ \vdots&\vdots &\ddots &\vdots \\ x_{n1}&x_{n2} &\cdots &x_{nm} \end{pmatrix}

2.1 数据归一化

将数据X针对每个特征归一化为均值为0,均方根为1的数据。

x_{i,j}^*=\frac{x_{i,j}-\bar{x_i}}{\sqrt{s_{i}}}

其中:

\bar {x_i}=\frac{1}{n}\sum_{j=1}^n{x_{ji}},i=1,2,...,n

s_{i}=\frac{1}{n-1}\sum_{j=1}^n{(x_{ji}-\bar {x_i})^2}, i=1,2,...,n

3. PCA降维

3.1 首先求取协方差矩阵

协方差矩阵的公式为:

R=\frac{1}{n-1}X^TX

计算出来的协方差矩阵为特征m*m维矩阵。

3.2 求取特征值和特征向量

求取协方差矩阵R的特征值和特征向量,并将特征值按照从大到小的顺序排列

评论 24
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值