我想使用标量函数的点分形式来计算整个输入矩阵的值。该函数定义如下:
function SpecificCharge(U, r, B)
( 2 * U ) / ( r^2 * B^2 )
end
,我想输入以下内容:
const Us = [ # voltages U
228.9, 243, 257, 271, 285, 300
]
const Rs = [ # radii r
9 7 5.8 5 4.3 3.9;
10.5 8.1 6.7 5.7 4.9 4.5;
10.1 8.3 7.1 6.5 6 4.8;
11.1 9.0 7.2 6.5 5.6 5.1;
10.5 8.3 7.8 6.5 5.7 5.4;
11.9 8.9 8.0 7.4 6.2 5.5
]
const Bs = [ # flux densities B
0.0007332339999999999, 0.00089328, 0.0010421599999999999, 0.00119104, 0.00133992, 0.0014873112
]
# Calling the dotted version of the function
specific_charges = SpecificCharge.(Us, Rs, Bs)
结果应为与Rs
相同尺寸的矩阵,每个半径由SpecificCharge
给出的相应特定电荷代替。每对U
-B
对应一个磁通密度,Rs
中的行对应于电压U
,列对应磁通密度B
。
问题是,我不确定茱莉亚广播的工作方式。如果我只是适当地对向量SpecificCharge.
和Us
进行转置,函数Bs
是否照原样工作?或者我需要使函数本身复杂化才能完成这项工作?
答案 0 :(得分:14)
执行功能
SpecificCharge
。如果我只是适当地对向量Us
和Bs
进行转置,或者我需要使函数本身复杂化才能完成这项工作,便可以按原样工作?
您尝试过吗?它按书面形式工作,但是,我猜您想换位Us
或Bs
以获得正确的“耦合”。
当您使用具有不同数组维度的广播(在示例中为矢量和矩阵)时,Julia将在单例维度中“扩展”数组。例如。向量将扩展为矩阵,该矩阵在每一列中均具有该向量的副本。这是一个说明这一点的例子
julia> A = ["A1", "A2"];
julia> B = ["B11" "B12" "B13";
"B21" "B22" "B23"];
julia> f(args...) = join(args, ", ");
julia> f.(A, B)
2×3 Array{String,2}:
"A1, B11" "A1, B12" "A1, B13"
"A2, B21" "A2, B22" "A2, B23"
如您所见,向量A
在单例维度上已扩展为
julia> A2 = [A A A]
2×3 Array{String,2}:
"A1" "A1" "A1"
"A2" "A2" "A2"
,使其具有与B
相同的大小,然后将f
应用于每对。要在另一个维度上扩展向量,您需要对其进行置换;
julia> C = ["C1", "C2", "C3"];
julia> f.(B, permutedims(C))
2×3 Array{String,2}:
"B11, C1" "B12, C2" "B13, C3"
"B21, C1" "B22, C2" "B23, C3"
现在C
基本上已经在第一维(单例维)中扩展了
julia> C2 = [permutedims(C); permutedims(C); permutedims(C)]
3×3 Array{String,2}:
"C1" "C2" "C3"
"C1" "C2" "C3"
"C1" "C2" "C3"