数模·降维

0xMayL

已于 2024-07-27 12:25:51 修改

阅读量446

点赞数 6

分类专栏：数学建模 # 降维文章标签：数学建模

于 2024-07-25 14:28:17 首次发布

本文链接：https://blog.csdn.net/2301_80132162/article/details/140686901

版权

数学建模同时被 2 个专栏收录

11 篇文章 0 订阅

订阅专栏

降维

1 篇文章 0 订阅

订阅专栏

降维

用较少的新变量代替较多的旧变量，并且最大限度保留原来变量的特征，叫做降维

主成分分析

1.标准化处理

也可以不进行标准化处理，到时候使用协方差矩阵即可，但是如果数据进行标准化后，一定要使用相关系数矩阵

在这里插入图片描

2.计算协方差矩阵或者相关系数矩阵

在这里插入图片描述

3.特征值

这里的特征值大小蕴含着方差大小

在这里插入图片描述

4.计算贡献率

将特征值进行排序，特征值越大的贡献率越大，越应该成为主成分

在这里插入图片描述

5.主成分确定

计算累加贡献率时应该对特征值列表进行排序(注意排序后下标变了)
对累计贡献率取一个阈值85%

在这里插入图片描述

6++ 分析结果

PCA+改进代码

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt


def pca():
    data = pd.read_excel('girl3.xlsx')  # 有可能依赖openxyl
    print('请确保数据进行标准化！')
    covariance_matrix = data.cov(ddof=1)  ## 参数ddof=1，表示结果除以N-1
    # 打印协方差矩阵
    print(F'协方差矩阵如下\n{covariance_matrix}')

    # 数据标准化后使用相关系数矩阵，数据未标准化可以使用协方差矩阵
    # 计算数据集中各列之间的相关系数
    correlation_matrix = data.corr()
    print(f'相关系数矩阵如下\n{correlation_matrix}')

    featValue, featVec = np.linalg.eig(correlation_matrix)  # 求解协方差矩阵的特征值和特征向量
    # print(f'特征值如下{featValue}长度为{len(featValue)}')
    dict={}
    for i in range(len(featValue)):
        dict[featValue[i]]=i
    print(f'特征值如下{featValue}长度为{len(featValue)}')
    featValue = sorted(featValue)[::-1]  # 对特征值进行排序，特征值大的贡献多
    print(f'特征向量矩阵如下\n{featVec}')

    x = np.linspace(0, 1, len(featValue))
    plt.plot(x, featValue)
    plt.show()

    # 主成分贡献和累加贡献
    featValue_contribute = featValue / np.sum(featValue)
    cumsum_contribute = np.cumsum(featValue_contribute)
    print(f'累加贡献如下{cumsum_contribute}')

    # 选出主成分
    k =[] # 选择85%作为阈值
    for i in range(len(cumsum_contribute)):
        if cumsum_contribute[i]<0.85:
            k.append(dict[featValue[i]])
    k=sorted(k)
    print(f'选择的主成分如下{k}')

    # 选出主成分对应的特征向量矩阵
    selectVec = np.matrix(featVec.T[k]).T
    finalData = np.dot(data, selectVec)
    # finalData = data @ selectVec
    print(F'特征向量矩阵如下\n{finalData}')
if __name__ == '__main__':
    pca()