exm0651_2.m

书名应为《掌握和精通MATLAB》附有完整的源程序

%exm0651_2.m

clear
A=sym([1 1/2 1/2 -1;1 1 -1 1;1 -1/4 -1 1;-8 -1 1 1]);
b=sym([0;10;0;1]);

syms k
A2=A(1:3,:);
X2=A2\b(1:3,1)					%求一个特解：最少非零元素的最小二乘解
XX2=X2+k*null(A2)					%构成通解
A2*XX2								%验算