subspacea（F，G，A）查找A-中矩阵F和G的列所跨越的两个子空间之间的所有min（size（orth（F），2），size（orth（G），2））主角基于标量积x'* A * y，其中A为Hermitian且为正定。主角的COS在统计中称为规范相关。
[theta，U，V] = subspacea（F，G，A）还计算左和右
主（规范）向量-分别为U和V的列。

如果F和G是单位长度的向量，且A = I，则在精确算法中，角度为ACOS（F'* G）。如果未提供A作为第三个参数，则A = I，并且函数给出的最大角度与Andrew Knyazev的SUBSPACE.m相同。MATLAB的SUBSPACE.m函数仍然设计不当，无法准确计算某些角度。

该算法在AV Knyazev和ME Argentati的《基于A的标量积中子空间之间的主角：算法和摄动估计》中进行了描述。SIAM J.Scientific Computing，23（2002），否。2009年6月6日。http://dx.doi.org/10.1137/S1064827500377332