Laplace transform not getting Value

1 view (last 30 days)
I tried to solve lapalce but not getting the value (Xs). Find the code below, Is there any issues in the code?
syms x(t) Xs
eqn = diff(x,t,2)+2*diff(x,t)+26*x(t) == 10*cos(t)*(heaviside(t-pi));
eqnLT = laplace(eqn)
eqnLT = 
eqnLT = subs(eqnLT,laplace(x(t)), Xs);
eqnLT = subs(eqnLT, {x(0), diff(x(t), t, 0)}, {1/2,1});
Xs = solve(eqnLT, Xs)
Xs = 

Accepted Answer

Star Strider
Star Strider on 16 Jan 2022
I do not see any specific problems.
syms x(t) Xs s
eqn = diff(x,t,2)+2*diff(x,t)+26*x(t) == 10*cos(t)*(heaviside(t-pi));
eqnLT = laplace(eqn)
eqnLT = 
eqnLT = subs(eqnLT,laplace(x(t)), Xs);
eqnLT = subs(eqnLT, {x(0), diff(x(t), t, 0)}, {1/2,1})
eqnLT = 
Xs = isolate(eqnLT, Xs)
Xs = 
X(s) = rhs(Xs)
X(s) = 
Perhaps the isolate function will do what you want.
.
  4 Comments

Sign in to comment.

More Answers (1)

Walter Roberson
Walter Roberson on 16 Jan 2022
syms x(t) Xs
eqn = diff(x,t,2)+2*diff(x,t)+26*x(t) == 10*cos(t)*(heaviside(t-pi));
eqnLT = laplace(eqn)
eqnLT = 
eqnLT = subs(eqnLT,laplace(x(t)), Xs);
dx = diff(x,t);
eqnLT = subs(eqnLT, {x(0), dx(0)}, {1/2,1});
Xs = solve(eqnLT, Xs)
Xs = 

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by