Laplace Transform of Given Differential Equation
357 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hello, I have the differential equation with initial condtions: y'' + 2y' + y = 0, y(-1) = 0, y'(0) = 0.
I need to use MATLAB to find the need Laplace transforms and inverse Laplace transforms.
I'm not sure if what I have so far is correct, here is what I have...
syms s t Y;
f = 0;
F = laplace(f,t,s);
Y1 = s*Y - 0;
Y2 = s*Y1 - 0;
laplaceSol = solve(Y2 + 2*Y1 + Y - F, Y) %Laplace Transform
invlaplaceSol = ilaplace(laplaceSol,s,t) %Inverse Laplace Transform
I get the following as output.
laplaceSol = 0
invlaplaceSol = 0
I also have the following code in an m-file.
function myplot(f,interv)
% myplot(f,[a,b])
% plot f for interval [a,b]
% here f is a symbolic expression, or a string
%
% example:
% myplot('x^2',[-1,1])
% syms x; myplot(x^2,[-1,1])
f = sym(f);
tv = linspace(interv(1),interv(2),300);
T = findsym(f,1);
plot(tv,double(subs(f,T,tv)))
Thank you,
4 commentaires
MD.AL-AMIN
le 2 Mar 2024
Déplacé(e) : Sam Chak
le 3 Mar 2024
y''(t) + 4y'(t) + 8y(t) = x' (t)+x(t) with x(t) = e ^ (- 4t) * u(t) y(0) = 0 and y'(0) = 0 matlab code
MD.AL-AMIN
le 2 Mar 2024
Déplacé(e) : Sam Chak
le 3 Mar 2024
I don't know how to solve by using MATLAB?
Réponses (2)
Sulaymon Eshkabilov
le 25 Avr 2022
Laplace transform does not work at t ~0 initial conditions and thus, here dsolve() might be a better option, e.g.:
syms y(t)
Dy=diff(y,t);
D2y = diff(y,t,2);
Eqn = D2y == -2*Dy-y;
ICs = [y(-1)==0, Dy(0)==0];
S = dsolve(Eqn, ICs)
7 commentaires
Walter Roberson
le 25 Avr 2022
Though, considering that dsolve() using Dy(0) == 0 gives the same solution as dsolve() with y(-1)==0, then you could potentially ignore the y(-1) and go ahead with laplace.
Sulaymon Eshkabilov
le 25 Avr 2022
Note that if your system has "zero" ICs and not excitation force; therefore, your system solution (response) will be zero. If you set one of your ICs, non-zero varlue and then you'll see something, e.g.:
syms s Y t
y0=0;
dy0=-1; %
Y1 = s*Y - y0;
Y2 = s*Y1- dy0;
Sol = solve(Y2 + 2*Y1 + Y, Y)
Sol = ilaplace(Sol,s,t)
fplot(Sol, [-1, 1])
% Verify: alternative solution with dsolve() gives the same result
syms y(t)
Dy=diff(y,t);
D2y = diff(y,t,2);
Eqn = D2y == -2*Dy-y;
ICs = [y(0)==0, Dy(0)==-1];
S = dsolve(Eqn, ICs);
fplot(S, [-1, 1])
1 commentaire
Walter Roberson
le 25 Avr 2022
However... the posters have been clear that they have an initial condition at y(-1) not an initial condition at y(0) . Which is a problem for laplace transforms.
Voir également
Catégories
En savoir plus sur Calculus dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!