先将xP=0两边取转置，得P'x'=0,求出x'再转置回来求出x. 当矩阵方程P'x'=0,rank(P)=r<n时,即P是奇异的，这时方程有无穷多个解,用MATLAB可求出它的一个基本解系， 基本解系的线性组合就是它的通解（一般解） 而求基本解系用matlab 中的命令 x=null(P', r )即可.其中:r=rank(P).
题主的问题可以用for双循环语句和solve（或vpasolve）函数求出FI(i)、GI(i)的值。
实现方法： k=1; for i=1:4 for j=1:4 syms FI0、GI0 [FI0,GI0]=vpasolve(FI0*dGI(i,j)-GI0*dSI(i,j)==0,FI0*dSI(i,j)-GI0*dGI(i,j)==0) FI(k)=FI0,GI(k)=GI0 k=k+1; end end
