今天继续,尝试加入一些范例
依然是简单的内容:
%%%%%%%%%%%%%%%%%%%%%说明 %%%%%%%%%%%%%%%%%%%
行空间的基:按行的角度来看待矩阵,更多介绍在代码说明里,简单的利用了昨天的代码。
%%%%%%%%%%%%%%%%%%%%%CODE %%%%%%%%%%%%%%%%%%
function B = rowbasis(A)
% rowbasis 行空间的基
%
% B = rowbasis(A) 返回A的行空间的一组基
% 以B的列表示
% A的行空间相当于A的转置的列空间
% rowbasis 找出A'中线性无关的
% 前r列,r代表A的秩
%
B = colbasis(A');
%%%%%%%%%%%%%%%%%%%%%范例%%%%%%%%%%%%%%%%%%
>> A=[1,1,1;1,3,2;2,3,4]
A =
1 1 1
1 3 2
2 3 4
>> rowbasis(A)
ans =
1 1 2
1 3 3
1 2 4
%%%%%%%%%%%%%%%%%%%%%说明 %%%%%%%%%%%%%%%%%%%
零空间的基:这个零空间的代码值得好好看看,体现了线性代数里关于列意义的基本思想。
直接看去你可能会不知道这些代码在干什么。因为一个简单的语句可能包含了很多操作。<