我用matlab中eig函数做了主成分分析,但是得到的特征值和原来的矩阵参量之间是什么对应关系那?比如我求了3*3的矩阵得出了一个从对角元素小到大排列的对角阵,但是这些特征值分别和原来矩
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/27 18:25:38
![我用matlab中eig函数做了主成分分析,但是得到的特征值和原来的矩阵参量之间是什么对应关系那?比如我求了3*3的矩阵得出了一个从对角元素小到大排列的对角阵,但是这些特征值分别和原来矩](/uploads/image/z/12355325-53-5.jpg?t=%E6%88%91%E7%94%A8matlab%E4%B8%ADeig%E5%87%BD%E6%95%B0%E5%81%9A%E4%BA%86%E4%B8%BB%E6%88%90%E5%88%86%E5%88%86%E6%9E%90%2C%E4%BD%86%E6%98%AF%E5%BE%97%E5%88%B0%E7%9A%84%E7%89%B9%E5%BE%81%E5%80%BC%E5%92%8C%E5%8E%9F%E6%9D%A5%E7%9A%84%E7%9F%A9%E9%98%B5%E5%8F%82%E9%87%8F%E4%B9%8B%E9%97%B4%E6%98%AF%E4%BB%80%E4%B9%88%E5%AF%B9%E5%BA%94%E5%85%B3%E7%B3%BB%E9%82%A3%3F%E6%AF%94%E5%A6%82%E6%88%91%E6%B1%82%E4%BA%863%2A3%E7%9A%84%E7%9F%A9%E9%98%B5%E5%BE%97%E5%87%BA%E4%BA%86%E4%B8%80%E4%B8%AA%E4%BB%8E%E5%AF%B9%E8%A7%92%E5%85%83%E7%B4%A0%E5%B0%8F%E5%88%B0%E5%A4%A7%E6%8E%92%E5%88%97%E7%9A%84%E5%AF%B9%E8%A7%92%E9%98%B5%2C%E4%BD%86%E6%98%AF%E8%BF%99%E4%BA%9B%E7%89%B9%E5%BE%81%E5%80%BC%E5%88%86%E5%88%AB%E5%92%8C%E5%8E%9F%E6%9D%A5%E7%9F%A9)
我用matlab中eig函数做了主成分分析,但是得到的特征值和原来的矩阵参量之间是什么对应关系那?比如我求了3*3的矩阵得出了一个从对角元素小到大排列的对角阵,但是这些特征值分别和原来矩
我用matlab中eig函数做了主成分分析,但是得到的特征值和原来的矩阵参量之间是什么对应关系那?
比如我求了3*3的矩阵得出了一个从对角元素小到大排列的对角阵,但是这些特征值分别和原来矩阵的哪列相对应那?
A =
1.0000 5.0000 0.5000
0.2000 1.0000 0.1429
2.0000 7.0000 1.0000
>> MSA2(A)
ans =
-0.0000 0 0
0 0.0487 0
0 0 2.9513
0,0.0487,2.9513被排序了,他们原来分别对应哪个?求高手解答
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('综合指标-时间曲线')
我用matlab中eig函数做了主成分分析,但是得到的特征值和原来的矩阵参量之间是什么对应关系那?比如我求了3*3的矩阵得出了一个从对角元素小到大排列的对角阵,但是这些特征值分别和原来矩
不知道你的MSA2是怎么回事
用[V,D]=eig(A)就好了,得到V就是特征向量,D得到是特征值对应的对角阵