How to plot a impulse unit function using Plots.jl and SymPy.jl?

General Usage
,
1

Hello. Good evening.

I’m doing all exercises of the Sadiku’s book in julia programming.

The code is in github

Im stucked to plot the impulse unit function after get it differentiating the unit step function.
The plot I get is only a line in zero, as plots.jl cant understand that the interval is all zero for values != t0. May someone help?

code:

let
	@syms t t₀

	u = sympy.Piecewise(
		(0, Lt(t, 0)),
		(1, Gt(t, 0))
	)
	@info u
	uₜ₀ = subs(u, t=> t - t₀)
	vₜ = 10(subs(uₜ₀, t₀=>2) - subs(uₜ₀, t₀=>5))
	@info vₜ
	p1 = plot(vₜ, xlims=(0, 6), lw=3, legend=:none, size=(400, 300))
	xlabel!(L"t")
	ylabel!(L"v(t)")
	title!(L"Pulso de Tensão $v(t)$")

	δ = diff(vₜ)
	# Ainda não sei como plotar a função de impulso unitário.
	@info δ
    plot( δ)
end

Results issues

image
image600×400 6.98 KB

Correct

image
image595×321 3.58 KB

2

Hi,

First of all, mathematically speaking a step function simply is not differentiable at the location of the step (according to the usual notion of differentiability), so in particular the derivative is not a dirac impulse. At every other point the derivative does exist and is zero.

I’m not very familiar with SymPy (and had to rely on PythonCall to get it working), but if I understand correctly, you’re explicitly not defining the function u: t \mapsto u(t) at t = 0. Consequently, δ will not be defined in 2 or 5. This is not related to the fact that the shifted u is non-differentiable here, just that you did not include it in the domain.

julia> δ
Python: 10*Piecewise((0, (t > 2) | (t < 2))) - 10*Piecewise((0, (t > 5) | (t < 5)))

julia> δ.subs(t, 2)
Python: nan

julia> δ.subs(t, 5)
Python: nan

If you would define u at t = 0 via Ge(t, 0), I get

julia> δ
Python: 0

so we still have issues with the non-differentiability.

In any case, if you want to recreate the plot, you could just manually add arrows

using Plots, LaTeXStrings

plot([0., 6.], [0., 0.], color=:black,
     xlabel=L"t", ylabel=L"$\frac{dv}{dt}$", legend=:none)
plot!([2., 2.], [0., 10.], arrow=true, color=:black, )
plot!([5., 5.], [0., -10.], arrow=true, color=:black)

arrows
arrows600×400 3.2 KB

(I’m not sure how to get the axes to intersect at the origin in Plots.jl.)

3

Your derivative returns nan at both 2 and 5. This function gives it a value and then plots (the arrow works with GR at least)

## δ(2) is nan; so this adjusts
f(x) = isnan(δ(x)) ? float(vₜ(x+.001) - vₜ(x-.001)) : float(δ(x))
plot(0:8, f; linetype=:sticks, arrow=true, framestyle=:origin)

(Strangely, I thought defining u so it is defined at 0 might make this work, but that gives back a 0 derivative.)