solving diffrential equation systems

13 Ansichten (letzte 30 Tage)
mohsen safaei
mohsen safaei am 7 Nov. 2014
Beantwortet: VBBV am 21 Mär. 2022
Hello everybody I want to solve this series of diffrential equation and everytime matlab types an eror what should I do please somebody help me v1=1000; v2=1500; v3=100; k1=.1; k2=.2; k3=.4; syms f(t) g(t) y(t) S=dsolve(diff(f)==((1/v1)*(1000+100*g-1100*f-k1*v1*f)),diff(g)==((1/v2)*(1100*f-1100*g-k2*v2*g)),diff(y)==((1/v3)*(1000*g-1000y-k3*v3*y)) and f(0) and g(0) and y(0) are all zero

Antworten (2)

Torsten
Torsten am 7 Nov. 2014
What exactly is the error message ?
Best wishes
Torsten.
  2 Kommentare
mohsen safaei
mohsen safaei am 7 Nov. 2014
hello this is the eror | Error: Unexpected MATLAB expression.
Torsten
Torsten am 10 Nov. 2014
Maybe because there is a parenthesis missing behind the expression
S=dsolve(diff(f)==((1/v1)*(1000+100*g-1100*f-k1*v1*f)),diff(g)==((1/v2)*(1100*f-1100*g-k2*v2*g)),diff(y)==((1/v3)*(1000*g-1000y-k3*v3*y))
?
Best wishes
Torsten.

Melden Sie sich an, um zu kommentieren.


VBBV
VBBV am 21 Mär. 2022
v1=1000;
v2=1500;
v3=100;
k1=.1; k2=.2; k3=.4;
syms f(t) g(t) y(t)
eq1 = diff(f)==((1/v1)*(1000+100*g-1100*f-k1*v1*f))
eq1(t) = 
eq2 = diff(g)==((1/v2)*(1100*f-1100*g-k2*v2*g))
eq2(t) = 
eq3 = diff(y)==((1/v3)*(1000*g-1000*y-k3*v3*y)) % use the multiplication operator
eq3(t) = 
S=dsolve([eq1 eq2 eq3])
S = struct with fields:
g: exp((t*(82^(1/2) - 32))/30)*(82^(1/2)/300 + 14/15)*(C2 - (5500*exp(-(t*(82^(1/2) - 32))/30)*(140*82^(1/2) - 41))/(178391*(82^(1/2) - 32))) - exp(-(t*(82^(1/2) + 32))/30)*(82^(1/2)/300 - 14/15)*(C3 - (5500*exp((t*(82^(1/2) + 32))/30)*(140*82^(1/… f: exp((t*(82^(1/2) - 32))/30)*((23*82^(1/2))/550 - 173/1100)*(C2 - (5500*exp(-(t*(82^(1/2) - 32))/30)*(140*82^(1/2) - 41))/(178391*(82^(1/2) - 32))) - exp(-(t*(82^(1/2) + 32))/30)*((23*82^(1/2))/550 + 173/1100)*(C3 - (5500*exp((t*(82^(1/2) + 32))… y: exp((t*(82^(1/2) - 32))/30)*(C2 - (5500*exp(-(t*(82^(1/2) - 32))/30)*(140*82^(1/2) - 41))/(178391*(82^(1/2) - 32))) + exp(-(t*(82^(1/2) + 32))/30)*(C3 - (5500*exp((t*(82^(1/2) + 32))/30)*(140*82^(1/2) + 41))/(178391*(82^(1/2) + 32))) + exp(-(52…

Kategorien

Mehr zu Mathematics finden Sie in Help Center und File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by