中易网

怎样用matlab实现主成分分析?

答案:1  悬赏:0  
解决时间 2021-11-08 14:12
怎样用matlab实现主成分分析?
最佳答案
function [lambda,T,fai]=MSA2(A)
%求标准化后的协差矩阵,再求特征根和特征向量
%标准化处理
[p,n]=size(A);
for j=1:n
mju(j)=mean(A(:,j));
sigma(j)=sqrt(cov(A(:,j)));
end
for i=1:p
for j=1:n
Y(i,j)=(A(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
% disp('方差贡献率:');
% disp(fai);
u=T(:,n);
B=[];
h=length(A(:,1));
for k=1:n
m1=mean(A(:,k));
t=(A(:,k)-m1).^2;
m2=sqrt(sum(t))/(h-1);
B=[B,(A(:,k)-m1)./m2];
end
y=B*u;
x1=1:1:length(y);
plot(x1,y);
xlabel('时间/小时')
ylabel('综合指标')
title('综合指标-时间曲线')
我要举报
如以上问答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
华硕24寸的宽屏液晶显示器加底座有多高?
丽江古城到玉龙雪山要多长时间
如何瘦肚肚?产后
瑞士提契诺州的花岗岩石在哪里
聚苯防火保温板价格高吗?
说说手摇式碎纸机好不好?
日立用电梯三 二楼高的多少钱
求电影《成事在人》英文影评,有些原创,不能
本溪如家宾馆小市店地址在什么地方,想过去办
电信100m光纤宽带用什么路由器好
头发稍短如何扎头发简单好看图解
直埋保温管由几个部分组成?
雅歌丽舍美发地址在什么地方,我要处理点事
怎样报名沁阳一中
蒙牛 纯牛奶利乐枕是低脂吗
推荐资讯
空压机有哪些故障呢?
深圳哪裡好玩啊?
今天流行长款羽绒服还是短款的
求心理罪前传1-4部全集txt和周瞳探案1-4部全
诚信手机店地址在哪,我要去那里办事
南京装修网 好点的,麻烦推荐下
打算怀孕前,梦见摘茄子玉米香瓜子把香瓜吃了
玻璃体混沌
荣威350怎么评价?
激光能彻底祛斑吗,有没有副作用
南水北调东线沟通哪4大水系
房屋建筑施工总承包包括不包括钢结构
手机登qq时,显示手机磁盘不足,清理后重新登
刺客的套装怎么选啊?