Ваш ответ правильный и общий.Однако обратите внимание, что при условии, что ваши внутренние массивы имеют одинаковый размер (а не просто одинаковую размерность), существует также следующий более быстрый способ:
julia> matrix = [rand(1,2,3) for _ in 1:4]; # some test data
julia> @btime a = cat($matrix..., dims=4); # your solution
11.519 μs (80 allocations: 3.83 KiB)
julia> @btime b = reshape(collect(Iterators.flatten($matrix)), (1,2,3,4)); # much faster solution
611.960 ns (55 allocations: 2.27 KiB)
julia> a == b
true