离散时间信号积分问题

请教在Julia 中对离散信号求积分一般是怎么操作的呢?譬如就一定采样速率得到的加速度信号求一次或二次积分得到速度或位移信号

建议取slack问一下(搭梯子快一点)
workspace的url是julialang.slack.com
这里的论坛人有点少,可能没人做过这个

谢谢nesteiner,这个问题应该很简单,就是不知道有没有现成的函数或相关的库。
以下是尝试直接将离散信号逐个按梯形面积累加的算法,也不知道有没错或该怎么优化。

using Makie

#模拟信号生成
f1 = 5 #Hz
ω1 = 2pi*f1 #rad/s
a(t) = cos(ω1*t)

#采样
fs = 1000 #采样速率
ts = 1/fs 
t = 0:ts:1 #采样时间
L = length(t) #采样长度
A = a.(t) #离散区间信号

#待求积分
V = zeros(L-1)

#积分 v(t1) = v(t0) + a*δt
i = 2
V[1]=0.5*(A[1]+A[2])*ts
while  i< L
    V[i]=V[i-1]+0.5*(A[i]+A[i+1])*ts
    println(V[i])
    global i = i +1
end
scene1=Scene()
plot!(t,A,color=:red)
plot!(t[1:L-1],V,color=:blue)

京ICP备17009874号-2