如何用Matlab将矩阵零空间的线性无关向量扩展为零空间基

核心思路 #

已知矩阵A和其零空间中k个线性无关向量x₁,…,xₖ，要扩展得到完整零空间基，关键是找到与已知向量线性无关且属于A零空间的向量，可通过正交补+零空间求解的组合方式实现。

步骤与代码实现 #

1. 先验证已知向量的合法性

首先必须确保输入的向量确实满足要求：

• 属于A的零空间：A*xᵢ ≈ 0（考虑浮点误差，用范数判断）
• 线性无关：向量组的秩等于k

% 示例输入（可替换为你的矩阵和向量）
A = [1 2 3 4; 5 6 7 8]; % 2×4矩阵，零空间维数为4-2=2
x1 = [-1; 1; 0; 0];     % 已知的零空间线性无关向量
X = x1;                  % 将已知向量按列拼接成矩阵
k = size(X, 2);          % 已知向量的数量

% 验证：已知向量在A的零空间
if norm(A*X, 'fro') > 1e-10
    error('输入的向量不全在A的零空间中');
end
% 验证：已知向量线性无关
if rank(X) ~= k
    error('输入的向量线性相关');
end

2. 扩展得到完整零空间基

通过以下步骤生成补充的基向量：

1. 求已知向量组X的正交补空间，得到基矩阵Z（Z的列与X的所有列正交）
2. 在Z张成的空间内，求解A的零空间，得到的向量就是与X线性无关的零空间向量

% 步骤1：求X的正交补空间基
Z = null(X');

% 步骤2：求A在Z空间上的零空间，得到补充向量的系数矩阵Y
Y = null(A*Z);

% 生成补充的零空间向量
X_extend = Z*Y;

% 合并得到完整的零空间基
full_null_basis = [X, X_extend];

3. 验证结果正确性

可以通过两个指标确认结果：

• full_null_basis的列数等于size(A,2) - rank(A)（零空间的维数）
• A*full_null_basis的Frobenius范数远小于1（浮点误差范围内的零矩阵）

% 验证基的维数是否正确
expected_dim = size(A,2) - rank(A);
if size(full_null_basis,2) ~= expected_dim
    error('扩展后的基维数不正确');
end

% 验证基向量都在零空间
disp('A乘以完整基的Frobenius范数（应接近0）：');
disp(norm(A*full_null_basis, 'fro'));

注意事项 #

• 如果已知向量已经是完整的零空间基，X_extend会是空矩阵，无需额外处理
• 浮点运算中，判断零矩阵时建议使用1e-10这类小阈值，避免因精度问题误判

内容的提问来源于stack exchange，提问作者district9