using FFTW
using TyPlot

定义时间范围和采样点数

t = 0:0.01:5 # 从0到5，步长为0.01
n = length(t)

定义单位阶跃信号

u = [x >= 0 ? 1 : 0 for x in t]

计算FFT

U = fft(u)

计算频率轴

frequencies = fftfreq(n, 1/0.01)

使用 TyPlot 绘制频谱

TyPlot.figure()
TyPlot.plot(frequencies, abs.(U))

只返数据.png

您好，图形已经完成绘制，因为abs.(U)的返回值均在0附近，因此绘制的图形为一条竖直的直线。

abs.(U)
501-element Vector{Float64}:
 500.9999999999999
   1.7319371752068542e-14
   3.3465790649875366e-15
   1.8226809783403305e-14
   1.213306252814783e-14
   2.7107982737991706e-14
   1.648024610214602e-14
   1.439755064796734e-14
   3.4216651175203884e-15
   6.512801026707881e-15
   1.7928982842403996e-14
   1.4637206250852028e-14
   4.7508602430004135e-15
   ⋮
   1.5798211851355722e-14
   6.182793562880679e-15
   1.504244928654353e-14
   1.7627699830891913e-14
   5.047936739424075e-15
   1.4096511187343402e-14
   1.2132869306563273e-14
   1.602160649646797e-14
   2.4939647858786302e-14
   2.1687736366240954e-14
   2.430707493483532e-14
   5.4216521744726745e-15

定义时间范围和采样点数

定义单位阶跃信号

计算FFT

计算频率轴

使用 TyPlot 绘制频谱

全部回答 1