矩阵正交补的求法：
利用QR分解即可.
[m,n]=size(A);
[Q,R]=qr(A);
那么Q(:,n+1:m)就是你所要的.