先将xP=0两边取转置，得P'x'=0,求出x'再转置回来求出x.当矩阵方程P'x'=0,rank(P)=r而求基本解系用matlab 中的命令 x=null(P', r )即可.其中:r=rank(P).
err = mean(sum((y-z).^2));
这样就行
