如何在Julia中实现类似Matlab的动态数组填充？

一维数组的动态填充 #

你目前用append!的方式是Julia中一维数组动态填充的常规写法，另外也可以用push!更简洁：

x = Float64[]
for a in 1:10
    push!(x, a)
end

不过Julia不像Matlab那样支持直接通过x[a] = a自动扩展空数组（会因索引超出数组长度报错），push!/append!是内置的标准动态填充方式。

多维数组的动态扩展方案 #

Julia默认的多维数组是固定大小的，直接通过索引赋值会报错，以下是几种可行的替代方案：

1. 手动实现动态扩展逻辑

可以写一个辅助函数，在赋值前检查数组大小，自动扩展并填充默认值（比如0）：

function set_dynamic!(A, val, idx...)
    current_size = size(A)
    new_size = Tuple(max(c, i) for (c,i) in zip(current_size, idx))
    if new_size != current_size
        # 创建新数组并复制原有元素，其余位置填充0
        new_A = zeros(eltype(A), new_size)
        new_A[[1:s for s in current_size]...] = A
        A = new_A
    end
    A[idx...] = val
    return A
end

# 使用示例
x = Matrix{Int}(undef, 0, 0)
x = set_dynamic!(x, 5, 2, 3)  # 扩展为2×3数组，x[2,3] = 5
x = set_dynamic!(x, 10, 1, 1) # 赋值x[1,1] = 10

2. 用字典模拟稀疏式赋值

如果大部分位置是空值，可以用字典存储索引与值的映射，最后再转换为数组：

# 用元组作为键存储多维索引
x_dict = Dict{Tuple{Int,Int}, Float64}()
x_dict[(2,3)] = 5.0
x_dict[(1,1)] = 10.0

# 转换为完整数组
max_i = maximum(k[1] for k in keys(x_dict))
max_j = maximum(k[2] for k in keys(x_dict))
x = zeros(Float64, max_i, max_j)
for ((i,j), val) in x_dict
    x[i,j] = val
end

3. 使用第三方扩展数组包

ExtensibleArrays.jl专门提供了支持自动扩展的数组类型，行为和Matlab非常接近：

# 先安装包：using Pkg; Pkg.add("ExtensibleArrays")
using ExtensibleArrays

# 初始化空的可扩展数组，未赋值位置默认是undef
x = ExtArray{Float64}(undef, 0, 0)
x[2,3] = 5.0  # 自动扩展为2×3数组
x[1,1] = 10.0

# 也可以初始化时填充默认值（比如0）
x = ExtArray{Float64}(zeros, 0, 0)
x[3,4] = 7.0

注意事项 #

和Matlab一样，这种动态扩展的方式效率较低，因为每次扩展都需要重新分配内存并复制数据。如果提前知道数组的最终大小，优先直接初始化固定大小的数组（比如x = zeros(Float64, 10, 10)）再赋值，能大幅提升性能。

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