Create periodic signals for simulating system response with lsim

[lucasmsoares96]

I need to generate a periodic quadratic wave to use as input to the lsim function. In Matlab I can generate this wave with the command:

z = tf('z',1);
H = (0.05*(z - 1)) / (z^2 - 1.85*z + 0.9);
[u,t] = gensig('square',20,120,1);%           <<<<
[y,t] = lsim(H,5*u,t);

How do I do the equivalent in Julia? I didn’t find about it in the ControlSystems.jl library.

[baggepinnen]

You can pass a function u(x, t), like this

julia> using ControlSystemsBase

julia> z = tf('z',1);

julia> H = (0.05*(z - 1)) / (z^2 - 1.85*z + 0.9);

julia> u(x, t) = sign(sin(t))

julia> res = lsim(H, u, 10);

julia> using Plots

julia> plot(res)

[lucasmsoares96]

to get the same output as equivalent code in matlab

I needed to do the following:

@time begin
    using Plots
    using ControlSystems
end

step(x::Float64) = x > 0 ? 1 : 0

@enum Type sine square pulse

function gensig(type::Type, tau::Int64)
    type == sine && return (x, t) -> sin(t * (2π / tau))
    type == square && return (x, t) -> step(-sin(t * (2π / tau)))
    type == pulse && return (x, t) -> (t % tau) == 0 ? 1 : 0
end

begin
    z = tf('z', 1)
    H = (0.05(z - 1)) / (z^2 - 1.85z + 0.9)
    u = gensig(square, 20)
    res = lsim(H, u, 120)
end

plot(res.t, [res.u' 5res.y'];
    lines=:stem,
    xticks=0:10:120
)

[baggepinnen]

If you’re bothered by the long load time of ControlSystems, load ControlSystemsBase instead, this example does not require full ControlSystems

[lucasmsoares96]

Ok. Thank you @baggepinnen.