How to Plot Inverse Laplace with Heaviside?

Question

Anthony Aquino 2023-11-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2046982-how-to-plot-inverse-laplace-with-heaviside

回答： Star Strider 2023-11-14

Hello. I want to plot an inverse laplace transform using simple MATLAB functions however, I get an error message when I try to. Does this have to do with the heaviside functions? The equation is provided below:

clc
syms s t
Q = 10/(s*(s^2+100*s+1600))-10*exp(-3*s)/(s*(s^2+100*s+1600));
pretty(Q)
q = ilaplace(Q);
pretty(q)
plot(t,q)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2023-11-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2046982-how-to-plot-inverse-laplace-with-heaviside#answer_1352127

在 MATLAB Online 中打开

Use the fplot function —

clc

syms s t

Q = 10/(s*(s^2+100*s+1600))-10*exp(-3*s)/(s*(s^2+100*s+1600));

pretty(Q)

10 exp(-3 s) 10 --------------------- - --------------------- 2 2 s (s + 100 s + 1600) s (s + 100 s + 1600)

q = ilaplace(Q)

q =

pretty(q)

exp(-80 t) exp(-20 t) / exp(240 - 80 t) exp(60 - 20 t) 1 \ 1 ---------- - ---------- - heaviside(t - 3) | --------------- - -------------- + ---- | 10 + --- 480 120 \ 4800 1200 1600 / 160

figure

fplot(q, [0 5])

grid

xlabel('t')

ylabel('q(t)')

axis('padded')

.