Laplace transform not getting Value
显示 更早的评论
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 = subs(eqnLT,laplace(x(t)), Xs);
eqnLT = subs(eqnLT, {x(0), diff(x(t), t, 0)}, {1/2,1});
Xs = solve(eqnLT, Xs)
采纳的回答
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 = subs(eqnLT,laplace(x(t)), Xs);
eqnLT = subs(eqnLT, {x(0), diff(x(t), t, 0)}, {1/2,1})
Xs = isolate(eqnLT, Xs)
X(s) = rhs(Xs)
.
4 个评论
BEENA BASHEER
2022-1-16
The equation needs to be set up correctly, and here it is not.
There is a typographical error in:
eqnLT = subs(eqnLT, {x(0), diff(x(t), t, 0)}, {1/2,1})
since it should be:
eqnLT = subs(eqnLT, {x(0), subs(diff(x(t), t), t, 0)}, {1/2,1});
and with that change:
syms x(t) Xs s X(s)
eqn = diff(x,t,2)+2*diff(x,t)+26*x(t) == 10*cos(t)*(heaviside(t-pi));
eqnLT = laplace(eqn);
eqnLT = subs(eqnLT,laplace(x(t)), X(s));
eqnLT = subs(eqnLT, {x(0), subs(diff(x(t), t), t, 0)}, {1/2,1});
eqnLT = isolate(eqnLT, X(s))
It does not eimplify further. (I checked.)
.
BEENA BASHEER
2022-1-16
Star Strider
2022-1-16
As always, my pleasure!
.
更多回答（1 个）
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 = subs(eqnLT,laplace(x(t)), Xs);
dx = diff(x,t);
eqnLT = subs(eqnLT, {x(0), dx(0)}, {1/2,1});
Xs = solve(eqnLT, Xs)
1 个评论
Sarthak
2024-5-27
%write a matlab code to find the solution of following differential
%equation using Laplace transform method
%d^2y/dx^2-3dy/dx+2y=sinx, y(0)=2, y'(0)=2
类别
在 帮助中心 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
