Matlab是否有内置函数生成指定稀疏度的对角占优随机矩阵？

嗨！作为刚接触Matlab的新手，你的需求其实挺具体的——目前Matlab并没有直接一步到位的内置命令能完全满足你说的所有要求，但我们可以用几个基础内置函数组合起来轻松实现，代码逻辑也很容易理解，我一步步给你拆解：

一、生成稠密的目标矩阵 #

如果你的矩阵阶数n不大，用稠密矩阵处理更直观：

1. 生成初始n阶[0,1]随机矩阵
   用Matlab内置的rand()函数就能生成元素均匀分布在[0,1]的随机矩阵：

   n = 5;          % 替换成你需要的矩阵阶数
   x = 30;         % 非对角元素置0的比例（百分比，比如30代表30%）
   A = rand(n);    % 初始随机矩阵
   

2. 将指定比例的非对角元素置0
   我们需要先定位所有非对角元素的位置，再随机挑选x%的位置设为0：

   % 获取所有非对角元素的线性索引
   non_diag_idx = find(~eye(n));
   % 随机选择x%的非对角元素索引
   zero_count = round(x / 100 * length(non_diag_idx));
   zero_idx = randsample(non_diag_idx, zero_count);
   % 将选中的位置置0
   A(zero_idx) = 0;
   

3. 设置对角元素为对应列的和，实现对角占优
   计算每一列的总和，然后把对角元素替换成这个总和——这样对角元素必然大于等于该列其他元素的和（因为所有元素都是非负的），自然满足对角占优：

   col_sums = sum(A);                  % 计算每列的总和
   A(logical(eye(n))) = col_sums;      % 替换对角元素
   

二、生成稀疏的目标矩阵 #

如果n很大（比如几百上千阶），稀疏矩阵能大幅节省内存、提升计算效率，用Matlab的稀疏矩阵工具链实现：

1. 生成初始稀疏随机矩阵
   用sprand()生成非对角的稀疏随机矩阵，再单独添加初始对角元素：

   n = 100;        % 大阶数示例
   x = 30;
   % 生成非对角的稀疏随机矩阵（保留(100-x)%的非对角元素）
   non_diag_sparse = sprand(n, n, (100 - x)/100);
   non_diag_sparse(logical(eye(n))) = 0;  % 清空对角位置，避免重复
   % 生成对角元素的稀疏矩阵
   diag_sparse = sparse(1:n, 1:n, rand(n,1));
   % 合并得到初始稀疏矩阵
   A_sparse = non_diag_sparse + diag_sparse;
   

2. 更新对角元素为列和
   逻辑和稠密矩阵一致，只是针对稀疏矩阵操作：

   col_sums = sum(A_sparse);
   A_sparse(logical(eye(n))) = col_sums;
   

补充说明 #

• 上述代码里的n和x都可以替换成你实际需要的数值；
• 因为最终对角元素是对应列的总和（包含原对角元素），所以不管原矩阵的非零元素分布如何，都能保证严格的对角占优（所有元素非负，对角元等于列内所有元素之和）。

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