角度
的相位角。
库::`工程师`
语法
函数调用
-
[参数:theta]=角度(<参数:z>>)— 返回间隔中的相位角[−π, π]对于复数数组的每个元素[参数:z]. 角度在[参数:theta]这样，z=abs(z)。*exp(i*theta).
例子：
复数的模量和相位
Details
让我们创建一个复数 z 并计算其模量 r 和相 西塔.
import EngeeDSP.Functions: angle
import EngeeDSP.Functions: abs
z = 2 * exp(im * 0.5)
r = abs(z)
theta = angle(z)
println("z=",z)
println("r=",r)
println("theta=",theta)z=1.7551651237807455 + 0.958851077208406im
r=2.0
theta=0.5FFT阶段
Details
让我们创建一个由两个频率的正弦曲线组成的信号 15 Hz和 40 赫兹。 第一正弦曲线具有相位 −π/4，而第二个 — π/2. 用频率对信号进行采样 100 赫兹一秒。
import EngeeDSP.Functions: angle, fft, fftshift
fs = 100
t = 0:1/fs:1-1/fs
x = cos.(2*pi*15*t .- pi/4) .- sin.(2*pi*40*t);算信号的傅立叶变换。 让我们绘制傅立叶变换模量对频率的依赖性。
y = fft(x)
z = fftshift(y)
ly = length(y)
f = (-ly/2:ly/2-1)/ly*fs;
plot(f, abs(z), xlabel="Frequency (Hz)", ylabel="|y|")
通过去除低幅度变换的值来计算傅立叶变换的相位。 让我们绘制相位与频率的相关性。
tol = 1e-6
z[abs(z) .< tol] .= 0
theta = angle(z)
plot(f, theta/pi, xlabel="Frequency (Hz)", ylabel="Phase / pi")