e0504.m

来自「《MATLAB 7.X程序设计》-王建卫-源代码及习题答案」· M 代码 · 共 14 行

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]);
%创建方程系数的符号常数矩阵                                                    
X=A\b         % X=inv(A)*b等价于X=A\b。
%如果求解前3个方程所构成的“欠定”方程组呢？
A=sym([1 1/2 1/2 -1;1 1 -1 1;1 -1/4 -1 1;-8 -1 1 1]);A2=A(1:3,:);b=sym([0;10;0;1]);
X=A2\b(1:3,1) ;         %求特解
syms k;
X2=X+k*null(A2)         %构成通解
A2*X2;