MATLAB中如何检查矩阵是否含全1对角线或为单位下三角矩阵

首先明确单位下三角矩阵的两个核心判定条件：

• 矩阵所有对角线元素必须为1
• 主对角线上方的所有元素必须为0（即矩阵是下三角矩阵）

现有代码的问题 #

你当前的代码逻辑存在两处错误：

1. all(diag(L)~=1)：这个条件只会在所有对角线元素都不等于1时触发报错，但实际我们需要的是只要有一个对角线元素不是1就报错，逻辑完全反了。
2. istril(L) == 1：istril(L)返回true时说明矩阵是下三角矩阵，这是符合要求的，但你现在的代码反而在符合要求时抛出错误，逻辑颠倒。

修正后的代码 #

% 检查对角线是否全为1
if ~all(diag(L) == 1)
    error('矩阵必须是单位下三角矩阵：对角线元素不全为1');
end

% 检查是否为下三角矩阵（上三角全0）
if ~istril(L)
    error('矩阵必须是单位下三角矩阵：不是下三角矩阵');
end

测试示例矩阵 #

你给出的示例矩阵：

L = [ 1 1 1 ; 0 1 1 ; 0 0 1 ]

这个矩阵的上三角部分有非零元素（第一行第二、三列，第二行第三列都是1），所以istril(L)会返回false，修正后的代码会触发第二个错误提示，符合预期——因为它不是下三角矩阵，自然也不是单位下三角矩阵。

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