基于python实现POD（本征正交分解）

前言

POD（本征正交分解），也称为PCA（主成分分析，是一种常用的数据降维方法，其具体的数学原理可以参考张洋老师的文章【PCA的数学原理-张洋】，本文基于具体算例，利用python实现POD方法。

POD算法步骤

假设有m条n维数据

1）将原始数据按列组成n行m列矩阵X

2）将X的每一行进行零均值化，即减去这一行的均值

3）求出协方差矩阵

4）求出协方差矩阵C的特征值即对应的特征向量

5）将特征向量按对应特征值大小从上到下按行排列成矩阵，取前k行组成矩阵P

6）即为降维到k维后的数据矩阵

python代码实现

import numpy as np
import scipy.io as io

def POD():
    #加载数据
    data = io.loadmat('./Data/POD_data')
    POD_data = data['