如何在MATLAB中构造6×6秩为3的对称矩阵？

以下是几种简单可靠的实现方法：

方法1：利用外积构造（最直观） #

对称低秩矩阵可通过满秩矩阵与其转置的乘积生成：若$B$是$6 \times 3$的满秩矩阵，则$A = B B^T$必为6×6对称矩阵，且秩等于$B$的秩（即3）。

MATLAB代码示例：

% 生成6×3的满秩随机矩阵
B = randn(6, 3);
% 构造对称矩阵A
A = B * B';
% 验证秩
disp('矩阵A的秩：');
disp(rank(A));

若需非随机的确定矩阵，可手动指定$B$的元素，比如：

B = [1 0 0; 0 1 0; 0 0 1; 1 1 0; 1 0 1; 0 1 1];
A = B * B';

方法2：正交相似变换构造 #

先构造秩为3的对角矩阵，再通过正交矩阵做相似变换——正交变换不改变矩阵的秩和对称性，最终得到的矩阵满足要求。

MATLAB代码示例：

% 生成6×6正交矩阵Q
Q = orth(randn(6));
% 构造秩为3的对角矩阵，前3个对角元非零，后3个为0
D = diag([2, 3, 5, 0, 0, 0]);
% 生成对称矩阵A
A = Q * D * Q';
% 验证对称性和秩
disp('矩阵A是否对称：');
disp(all(all(A == A')));
disp('矩阵A的秩：');
disp(rank(A));

方法3：分块矩阵构造 #

用一个3×3的满秩对称矩阵$S$，构建分块对称矩阵：
$$A = \begin{bmatrix} S & S \ S & S \end{bmatrix}$$
此时$A$是6×6对称矩阵，且秩等于$S$的秩（3）。

MATLAB代码示例：

% 生成3×3满秩对称矩阵S
S = randn(3);
S = S + S'; % 确保对称
% 构造分块对称矩阵A
A = [S, S; S, S];
% 验证秩
disp('矩阵A的秩：');
disp(rank(A));

注意：MATLAB浮点数计算存在微小误差，rank(A)会自动忽略数值误差，返回正确秩。

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