vpasolve() gives wrong solution

4 Ansichten (letzte 30 Tage)
John Arbolino
John Arbolino am 13 Sep. 2020
Kommentiert: Walter Roberson am 17 Sep. 2020
I am trying to find the time at which my equation equals a value, but whenever I use vpasolve I am given a completely wrong answer of a negative time value. When I also add t>0 it gives no answer even though the function does equal that value at some positive t value.
syms s t
y(s) = 36/(s*(s^2+2*s+36));
response(t) = ilaplace(y(s))
s = vpasolve(response == 0.1)
st = vpasolve(response == 0.1, t>0)
response(s)
response(t) =
1 - exp(-t)*(cos(5.9161*t) + 0.1690*sin(5.9161*t))
s =
-226.9853
st =
Empty sym: 0-by-1
ans =
-1.0328e+62
You can see that when I try and use solve() without a specific time domain, I receieve a value S that when I plug back into response() is wildly off the value I am looking for. I know for a fact that the function response(t) does cross the value 0.1 in the positve time domain, I have graphed it.

Akzeptierte Antwort

Walter Roberson
Walter Roberson am 13 Sep. 2020
s = vpasolve(response == 0.1, 0.07)
  2 Kommentare
John Arbolino
John Arbolino am 17 Sep. 2020
This works thanks! Where is the default starting value, because I assumed it would be zero.
Walter Roberson
Walter Roberson am 17 Sep. 2020
I think the default start is 0. However, 0 is not a good starting point for this particular equation.

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Kategorien

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

Tags

Produkte


Version

R2019a

Community Treasure Hunt

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

Start Hunting!

Translated by