ea562.m

来自「一本关于工程线性代数及其matlab编程方面的好书。」· M 代码 · 共 9 行

% 《工程线性代数（MATLAB版）》第五章例题5.6.2程序ea562
%
A=[1,1,1,1,1;3,2,1,1,-3;0,1,2,2,6;5,4,3,3,-1];
b=[7;-2;23;12];                  % 系数矩阵赋值
U0=rref([A,b])                   % 将增广矩阵化为最简行阶梯形式
[m,n]=size(A);x0=zeros(n,1);    % 将特解初始化为n×1向量
x0(1:m,:)=U0(:,end)               % 求原方程的一个特解           
Z=null(A,'r')                    % 求齐次方程组的整数基础解系