Julia 之初体验（五）复数

复数

Julia用全局常数im代表复数i，代表-1的平方根。

julia> 1+2im
1 + 2im

julia> (1 + 2im)*(2 - 3im)
8 + 1im

julia> (1 + 2im)/(1 - 2im)
-0.6 + 0.8im

julia> (1 + 2im) + (1 - 2im)
2 + 0im

julia> (-3 + 2im) - (5 - 1im)
-8 + 3im

julia> (-1 + 2im)^2
-3 - 4im

julia> (-1 + 2im)^2.5
2.729624464784009 - 6.9606644595719im

julia> (-1 + 2im)^(1 + 1im)
-0.27910381075826657 + 0.08708053414102428im

julia> 3(2 - 5im)
6 - 15im

julia> 3(2 - 5im)^2
-63 - 60im

julia> 3(2 - 5im)^-1.0
0.20689655172413796 + 0.5172413793103449im

julia> 2(1 - 1im)
2 - 2im

julia> (2 + 3im) - 1
1 + 3im

julia> (1 + 2im) + 0.5
1.5 + 2.0im

julia> (2 + 3im) - 0.5im
2.0 + 2.5im

julia> 0.75(1 + 2im)
0.75 + 1.5im

julia> (2 + 3im) / 2
1.0 + 1.5im

julia> (1 - 3im) / (2 + 2im)
-0.5 - 1.0im

julia> 2im^2
-2 + 0im

julia> 1 + 3/4im
1.0 - 0.75im

注意：3/4im = 3/(4*im) =3im/(4im*im)= -(3/4*im),Julia把复数看成一个数字。

Julia还提供了用于处理复杂值的标准函数：取实部、取虚部系数、共轭复数、弧度等

julia> z = 1 + 2im
1 + 2im

julia> real(1 + 2im) # real part of z
1

julia> imag(1 + 2im) # imaginary part of z
2

julia> conj(1 + 2im) # complex conjugate of z
1 - 2im

julia> abs(1 + 2im) # absolute value of z
2.23606797749979

julia> abs2(1 + 2im) # squared absolute value
5

julia> angle(1 + 2im) # phase angle in radians
1.1071487177940904

通常，复数的绝对值（abs）是其到零的距离。 abs2给出绝对值的平方，并且对于复数特别有用，因为它避免了平方根。 angle返回以弧度表示的相位角（也称为自变量或arg函数）。 还为复数定义了其他基本功能，如：

julia> sqrt(1im)
0.7071067811865476 + 0.7071067811865475im

julia> sqrt(1 + 2im)
1.272019649514069 + 0.7861513777574233im

julia> cos(1 + 2im)
2.0327230070196656 - 3.0518977991518im

julia> exp(1 + 2im)
-1.1312043837568135 + 2.4717266720048188im

julia> sinh(1 + 2im)
-0.4890562590412937 + 1.4031192506220405im

注意，数学函数通常在应用于实数时返回实数值，而在应用于复数时返回复数值。 例如，即使-1 == -1 + 0im，sqrt在应用于-1和-1 + 0im时的行为也会有所不同：sqrt(-1)会报错，因为实数范围内，-1的平方根是不存在的。只能写sqrt(-1 + 0im)，这样返回复数值就是im了。

julia> sqrt(-1)
ERROR: DomainError with -1.0:
sqrt will only return a complex result if called with a complex argument. Try sqrt(Complex(x)).
Stacktrace:
[...]

julia> sqrt(-1 + 0im)
0.0 + 1.0im

julia> a = 1; b = 2; a + b*im
1 + 2im

Julia推荐使用复数函数，避免了乘法和加法运算。

julia> a = 1; b = 2; complex(a, b)
1 + 2im

Inf和NaN也允许用在复数中。

julia> 1 + Inf*im
1.0 + Inf*im

julia> 1 + NaN*im
1.0 + NaN*im