QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 12162|回复: 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 微信微信
    Kingwax        

    1

    主题

    3

    听众

    98

    积分

    升级  97.89%

  • TA的每日心情

    2014-9-21 10:10
  • 签到天数: 20 天

    [LV.4]偶尔看看III

    回复

    使用道具 举报

    zhangbq        

    0

    主题

    13

    听众

    208

    积分

    升级  54%

  • TA的每日心情
    奋斗
    2015-7-15 09:40
  • 签到天数: 76 天

    [LV.6]常住居民II

    自我介绍
    一名知识薄弱的大学生,12年入学
    回复

    使用道具 举报

    0

    主题

    9

    听众

    127

    积分

    升级  13.5%

  • TA的每日心情
    奋斗
    2016-4-13 14:16
  • 签到天数: 74 天

    [LV.6]常住居民II

    自我介绍
    我来自中北大学

    社区QQ达人

    群组2015年美赛冲刺

    群组物联网工程师培训

    回复

    使用道具 举报

    Viming        

    1

    主题

    9

    听众

    285

    积分

    升级  92.5%

  • TA的每日心情

    2015-2-8 15:24
  • 签到天数: 104 天

    [LV.6]常住居民II

    自我介绍
    MI

    社区QQ达人 新人进步奖

    群组华南理工大学

    回复

    使用道具 举报

    swjtuwl        

    0

    主题

    13

    听众

    32

    积分

    升级  28.42%

  • TA的每日心情
    难过
    2015-4-26 11:36
  • 签到天数: 23 天

    [LV.4]偶尔看看III

    自我介绍
    参加数学建模国赛
    回复

    使用道具 举报

    swjtuwl        

    0

    主题

    13

    听众

    32

    积分

    升级  28.42%

  • TA的每日心情
    难过
    2015-4-26 11:36
  • 签到天数: 23 天

    [LV.4]偶尔看看III

    自我介绍
    参加数学建模国赛
    回复

    使用道具 举报

    1

    主题

    11

    听众

    586

    积分

    升级  95.33%

  • TA的每日心情
    开心
    2014-9-26 14:03
  • 签到天数: 58 天

    [LV.5]常住居民I

    群组Matlab讨论组

    回复

    使用道具 举报

    陈锦欣        

    2

    主题

    14

    听众

    440

    积分

    升级  46.67%

  • TA的每日心情
    开心
    2016-3-11 14:32
  • 签到天数: 175 天

    [LV.7]常住居民III

    自我介绍
    努力奋斗

    社区QQ达人

    回复

    使用道具 举报

    0

    主题

    11

    听众

    175

    积分

    升级  37.5%

  • TA的每日心情
    慵懒
    2014-11-12 23:58
  • 签到天数: 28 天

    [LV.4]偶尔看看III

    自我介绍
    路过

    群组国赛讨论

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

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2025-8-17 01:02 , Processed in 0.900566 second(s), 104 queries .

    回顶部