QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 12165|回复: 45
打印 上一主题 下一主题

主成分分析(PCA)Matlab源码

[复制链接]
字体大小: 正常 放大

2620

主题

162

听众

1万

积分

升级  0%

  • TA的每日心情
    开心
    2015-3-12 15:35
  • 签到天数: 207 天

    [LV.7]常住居民III

    社区QQ达人 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    群组第六届国赛赛前冲刺培

    群组国赛讨论

    群组2014美赛讨论

    群组2014研究生数学建模竞

    群组数学中国试看培训视频

    跳转到指定楼层
    1#
    发表于 2014-8-30 11:55 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    主成分分析(PCA)Matlab源码
    principal component analysis(PCA)主成分分析,Matlab实例代码

    function main()
    %*************主成份分析************
    %
    %see also http://www.matlabsky.com
    %
    %读入文件数据
    X=load('data.txt');
    %==========方法1:求标准化后的协差矩阵,再求特征根和特征向量=================
    %标准化处理
    [p,n]=size(X);
    for j=1:n
        mju(j)=mean(X(:,j));
        sigma(j)=sqrt(cov(X(:,j)));
    end
    for i=1:p
        for j=1:n
            Y(i,j)=(X(i,j)-mju(j))/sigma(j);
        end
    end
    sigmaY=cov(Y);
    %求X标准化的协差矩阵的特征根和特征向量
    [T,lambda]=eig(sigmaY);
    disp('特征根(由小到大):');
    disp(lambda);
    disp('特征向量:');
    disp(T);
    %方差贡献率;累计方差贡献率
    Xsum=sum(sum(lambda,2),1);
    for i=1:n
        fai(i)=lambda(i,i)/Xsum;
    end
    for i=1:n
        psai(i)= sum(sum(lambda(1:i,1:i),2),1)/Xsum;
    end
    disp('方差贡献率:');
    disp(fai);
    disp('累计方差贡献率:');
    disp(psai);
    %综合评价....略
    复制代码

    游客,如果您要查看本帖隐藏内容请回复

    点评

    dddddolores  1234567890-1234567890  发表于 2015-2-5 21:09
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏1 支持支持0 反对反对0 微信微信

    1

    主题

    2

    听众

    1024

    积分

    升级  2.4%

  • TA的每日心情
    擦汗
    2023-5-2 09:40
  • 签到天数: 298 天

    [LV.8]以坛为家I

    国际赛参赛者

    网络挑战赛参赛者

    网络挑战赛参赛者

    邮箱绑定达人

    回复

    使用道具 举报

    3

    主题

    5

    听众

    32

    积分

    升级  28.42%

  • TA的每日心情
    开心
    2020-5-31 08:31
  • 签到天数: 18 天

    [LV.4]偶尔看看III

    自我介绍
    热爱数学,更爱数模
    回复

    使用道具 举报

    3

    主题

    5

    听众

    32

    积分

    升级  28.42%

  • TA的每日心情
    开心
    2020-5-31 08:31
  • 签到天数: 18 天

    [LV.4]偶尔看看III

    自我介绍
    热爱数学,更爱数模
    回复

    使用道具 举报

    0

    主题

    2

    听众

    44

    积分

    升级  41.05%

    该用户从未签到

    回复

    使用道具 举报

    csmsh123        

    0

    主题

    2

    听众

    12

    积分

    升级  7.37%

  • TA的每日心情
    开心
    2019-1-4 19:08
  • 签到天数: 1 天

    [LV.1]初来乍到

    回复

    使用道具 举报

    0

    主题

    3

    听众

    66

    积分

    升级  64.21%

  • TA的每日心情
    无聊
    2018-11-4 16:06
  • 签到天数: 1 天

    [LV.1]初来乍到

    国际赛参赛者

    回复

    使用道具 举报

    0

    主题

    3

    听众

    43

    积分

    升级  40%

  • TA的每日心情
    开心
    2019-5-13 18:07
  • 签到天数: 3 天

    [LV.2]偶尔看看I

    回复

    使用道具 举报

    0

    主题

    5

    听众

    53

    积分

    升级  50.53%

  • TA的每日心情
    开心
    2018-2-8 12:36
  • 签到天数: 18 天

    [LV.4]偶尔看看III

    自我介绍

    邮箱绑定达人

    回复

    使用道具 举报

    0

    主题

    13

    听众

    78

    积分

    升级  76.84%

  • TA的每日心情

    2017-9-16 10:08
  • 签到天数: 9 天

    [LV.3]偶尔看看II

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2025-8-17 09:50 , Processed in 0.928154 second(s), 108 queries .

    回顶部