orthcompletprojection.m

可以直接用

function P=OrthCompletProjection(Vectors)
%获取向量集Vectors的正交补投影算子P
[m,n]=size(Vectors);
SpectralPInverse=(inv(Vectors'*Vectors))*Vectors';%类数 X 波段数
E=eye(m);%
P=(E-(Vectors*SpectralPInverse));%光谱信号的正交空间的投影算子