【发布时间】:2014-02-06 06:31:23
【问题描述】:
我想写一个函数来计算矩阵的上三角性质。所以让我们说 a_i_j 是第 i^th 行和 j^th 列中的数字。一个矩阵是 如果所有 i > j 的 a_i_j = 0,则为上三角形。
【问题讨论】:
我想写一个函数来计算矩阵的上三角性质。所以让我们说 a_i_j 是第 i^th 行和 j^th 列中的数字。一个矩阵是 如果所有 i > j 的 a_i_j = 0,则为上三角形。
【问题讨论】:
试试这个:
(defn is-upper-triangular [m]
(->> (map-indexed vector m)
(mapcat (fn [[r v]] (take r v)))
(every? zero?)))
上面的代码从第一行取 0 个元素,从第二行取 1 个元素,从第三行取 2 个元素,等等……并检查所有取入的元素是否为零。如果全部为零,则为上三角形。
此代码不检查给定矩阵是否为正方形。如有必要,您可以添加此检查。
【讨论】:
upper-triangular? 尚未包含在 core.matrix API 中,但如果您使用 vectorz-clj,您可以通过 Java 互操作获得该函数:
(def a (array :vectorz [[1 2] [0 4]]))
(.isUpperTriangular a)
=> true
(.isUpperTriangular (transpose a))
=> false
【讨论】:
(use 'clojure.core.matrix) 才能将 core.matrix API 函数导入当前命名空间