A 是需要分解的矩阵，p是最大特征值对应的特征向量：

% matlab代码如下<br>
% 给一个方阵A<br>
A=magic(5)<br>
% 对A特征分解，V为特征向量矩阵，D为对角矩阵<br>
[V,D]=eig(A)<br>
% 取D主对角线元素得向量d，d的每个元素都是A的特征值<br>
d=diag(D)<br>
% 降序排列d的元素，向量I记录了d中元素的顺序，即d(I(k))为A的第k大的特征值<br>
[~,I]=sort(d,\'descend\')<br>
% 取前三个最大的特征值<br>
k=3;<br>
% 前三个最大特征值对应的特征向量<br>
VI=V(:,I(1:k))<br>
% 前三个最大的特征值<br>
dI=d(I(1:k))<br>
% 前三个最大的特征值组成的对角矩阵<br>
DI=diag(dI)

matlab代码如下

您好，这段代码特别棒，可以请教一下如果用R写，应该怎样写吗？

您好，这段代码特别棒，可以请教一下用R应该怎么写吗？