high-frequency projection operator
Phf
Returns the NxN matrix of rank N-M that projects splines coefficients xi into the high-frequency subspace. This is used to filter the signal.