peak limiting current mode control
determine d
(%i1) | e1: jc = j0 + m1 * d; |
determine j1
(%i2) | e2: j1 = jc + m2 * (1 - d); |
solve for j1 and d
(%i3) | s1: linsolve([e1, e2], [j1, d]); |
(%i4) | ej1: s1[1]; |
find the steady state j0, labeled J0
(%i5) | e3: ev(ej1, j1 = j0); |
(%i6) | s2: linsolve(e3, j0); |
(%i7) | J0: rhs(s2[1]); |
introduce perturbations
(%i8) | e4: ev(ej1, j1 = J0 + dj1, j0 = J0 + dj0); |
solve for dj1
(%i9) | e4: ratsimp(e4); |
final result
(%i10) | solve(e4, dj1); |