Julia自定义Rank 2矩阵构造器报错：无匹配iterate(::Missing)方法

错误的核心原因是你自定义的RankOneMatrix和RankTwoMatrix虽然继承了AbstractMatrix{T}，但没有实现该接口的必需方法。虽然Julia官方文档确实说明无需手动定义iterate，但size和getindex是所有抽象矩阵类型必须实现的基础方法——REPL在尝试显示矩阵实例时，会依赖这两个方法获取矩阵的维度和元素，缺少它们会导致内部逻辑抛出Missing相关的迭代错误。

必须实现的接口方法 #

对于继承AbstractMatrix{T}的自定义类型，必须实现以下两个方法：

• size(M::YourMatrixType)：返回矩阵的行列维度元组
• getindex(M::YourMatrixType, i::Int, j::Int)：返回矩阵第i行第j列的元素

修正后的完整代码 #

# 定义RankOneMatrix并实现必需接口
struct RankOneMatrix{T} <: AbstractMatrix{T}
    v::AbstractVector{T}
    w::AbstractVector{T}
end

RankOneMatrix(v) = RankOneMatrix(v, v)

# 实现size方法：RankOneMatrix是v*w'，维度为(length(v), length(w))
Base.size(M::RankOneMatrix) = (length(M.v), length(M.w))

# 实现getindex方法：元素为v[i] * w[j]
Base.getindex(M::RankOneMatrix, i::Int, j::Int) = M.v[i] * M.w[j]

# 定义RankTwoMatrix并实现必需接口
struct RankTwoMatrix{T} <: AbstractMatrix{T}
    A::RankOneMatrix{T}
    B::RankOneMatrix{T}
end

RankTwoMatrix(v1, v2) = RankTwoMatrix(RankOneMatrix(v1), RankOneMatrix(v2))

# 实现size方法：需保证A和B维度一致，这里直接复用A的尺寸
Base.size(M::RankTwoMatrix) = size(M.A)

# 实现getindex方法：元素为A[i,j] + B[i,j]
Base.getindex(M::RankTwoMatrix, i::Int, j::Int) = M.A[i,j] + M.B[i,j]

# 测试初始化与显示
x = 1.0:5.0
y = 0.5:0.5:3.0
rank_2 = RankTwoMatrix(RankOneMatrix(collect(x)), RankOneMatrix(collect(y)))

# 验证可以正常显示和计算
println(rank_2)
println(rank_2[2,3])

补充说明 #

• 如果你需要支持矩阵的修改操作，还需要实现setindex!方法；如果不需要可变矩阵，可以忽略。
• 构造RankTwoMatrix时，建议添加维度检查，确保两个RankOneMatrix的行列数一致，避免后续计算出错：

  function RankTwoMatrix(A::RankOneMatrix{T}, B::RankOneMatrix{T}) where T
      @assert size(A) == size(B) "两个RankOneMatrix的维度必须一致"
      new{T}(A, B)
  end
  

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