e0504.m

来自「matlab的一些学习资料和matlab的部分经典源码」· 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;