How find symbolic solution matlab

I am trying to solve this simple autonomous differential equation using the dsolve function but MATLAB cannot find a symbolic solution to the function. Does anyone have any idea how I could solve the problem? The code is the following:
clc, clear all
syms y(t)
eqn = diff(y,t,2) == -9.8*sin(y) - 0.5*y'
Dy = diff(y,t)
cond = [y(0)==0.8, Dy(0)==0.8]
ySol(t) = dsolve(eqn,cond)
Thanks in advance to anyone who can help me.

4 Kommentare

Stephan
Stephan am 23 Nov. 2020
Bearbeitet: Stephan am 23 Nov. 2020
Sure about:
sin(y)
or did you mean:
sin(t)
Also does:
0.5*y'
mean the first derivative or the transpose of y - like Matlab will interpret this?
it's:
sin(y)
Because t is a sym var. And:
0.5*y'
Is the first derivative, wich i was also trying to declareas:
Dy = diff(y,t)
But it throws me the same warning, that MATLAB cannot find a symbolic solution.
John D'Errico
John D'Errico am 23 Nov. 2020
It cannot find a symbolic solution because one apparently does not exist. And that is why Stephan was asking if you really intended sin(t), not sin(y), as there may well be a solution in that case.
Juan Saa
Juan Saa am 23 Nov. 2020
I understand, thanks for your help

Melden Sie sich an, um zu kommentieren.

Antworten (1)

Stephan
Stephan am 23 Nov. 2020

1 Stimme

It appears that there either is no symbolic solution or Matlab is not powerful enough to find one. For sin(t) it would work. However Wolfram Alpha also does not provide an analytical solution for this ode - here is the numeric approach:
syms y(t)
Dy = diff(y,t)
eqn = diff(y,t,2) == -9.8*sin(y) - 0.5*Dy
% make numeric function from symbolic
[V,S] = odeToVectorField(eqn)
fun = matlabFunction(V,'Vars',{'t','Y'})
% solve ode numeric
[t,y] = ode45(fun,[0 10],[0.8, 0.8]);
plot(t,y)

2 Kommentare

Juan Saa
Juan Saa am 23 Nov. 2020
Exactly, I was also trying it with Wolfram Alpha and can't find a solution. I think it has to do with the differential equation being autonomous.
Juan Saa
Juan Saa am 23 Nov. 2020
Thank you very much for your help! :)

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Symbolic Math Toolbox finden Sie in Hilfe-Center und File Exchange

Gefragt:

am 23 Nov. 2020

Kommentiert:

am 23 Nov. 2020

Community Treasure Hunt

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

Start Hunting!

Translated by