Ask Your Question
1

Solving symbolically equation system

asked 2011-07-23 14:24:08 -0500

Ondra gravatar image

I have an equation system:

U(t) = R * I(t) + L * I'(t) + uC(t)

I(t) = C * uC'(t)

I want to know value of I'(t) and uC'(t), which is (-I(t) * R + U - uC(t))/L and I(t)/C respectively.

In sage I represent it this way:

R = 6; C = 10^(-4); L = 0.1
t = var('t')
U = function('U', t).function(t); I = function('I', t).function(t); uC = function('uC', t).function(t)

equations = [
U == R * I + L * I.diff(t) + uC,
I == C * uC.diff(t)
]

but

solve(equations, I.diff(t), uC.diff(t))

does't seem to be working. (TypeError: 'sage.symbolic.expression.Expression' object is not iterable)

why? how can I do this?

edit retag flag offensive close merge delete

1 answer

Sort by ยป oldest newest most voted
1

answered 2011-07-25 17:55:36 -0500

benjaminfjones gravatar image

As far as I know, solve only works with variables. So as a preliminary step you can substitute some new variables for I.diff(t) and uC.diff(t) and then solve. The other problem occuring is that to solve a system of equations for multiple variables, those variables should be specified as a list. Try this: (first substitute new variables DI and DuC for the two derivatives and then call solve)

sage: E2 = [ e.subs({I.diff(t): DI, uC.diff(t): DuC}) for e in equations ]
sage: E2
[t |--> U(t) == 0.100000000000000*DI + 6*I(t) + uC(t), t |--> I(t) == 1/10000*DuC]
sage: solve(E2, [DI, DuC])
[[DI == 10*U(t) - 60*I(t) - 10*uC(t), DuC == 10000*I(t)]]
edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 2011-07-23 14:24:08 -0500

Seen: 897 times

Last updated: Jul 27 '11