А по моему тут все гораздо проще
p=Pa,
p- искомая проекция
а- элемент
P=L*[(L^t)*L]^(-1)*(L^t) >>> ^t трансп ^(-1) обратная
L - матрица коэффициентов