Circuit Theory/TF Examples/Example14

From Wikibooks, open books for an open world
Jump to: navigation, search
Finding io using transfer functions

Find a partical solution to one of the re-occurring examples in this book. Given VS, find io. This has to be split into two transfer functions:

\frac{\mathbb{I}_o}{\mathbb{V}_s} = \frac{\mathbb{I}_o}{\mathbb{I}_t} \frac{\mathbb{I}_t}{\mathbb{V}_s}

Which is of the form:

\frac{\mathbb{I}_o}{\mathbb{V}_s} = \text{current divider} * \text{admittance}

The current divider is going to be:

\frac{\mathbb{I}_o}{\mathbb{I}_t} = \frac{\frac{1}{R_3}}{\frac{1}{R_3} + \frac{1}{sL_2 + \frac{1}{sC_2}}} = \frac{sL_2 + \frac{1}{sC_2}}{R_3 + sL_2 + \frac{1}{sC_2}}

Substituting into the transfer function:

\frac{\mathbb{I}_o}{\mathbb{V}_s} = \frac{sL_2 + \frac{1}{sC_2}}{R_3 + sL_2 + \frac{1}{sC_2}} * \frac{1}{R_1 + \frac{1}{sC_1} + sL_1 + \frac{1}{\frac{1}{R_3} + \frac{1}{sL_2 + \frac{1}{sC_2}}}}

Substituting and solving using this mupad code:

R1 :=1000;C1 :=.0001;C2 :=.0002;R3 := 2000;L1 := .001;L2 := .001;
f := (s*L2 + 1/(s*C2))/(R3 + s*L2 + 1/(s*C2))*1/(R1 + 1/(s*C1) + s*L1 + 1/(1/R3 + 1/(s*L2 + 1/(s*C2)))):
factor(f);
\frac{\mathbb{I}_o}{\mathbb{V}_s} = \frac{1000*(s^2 + 5000000.0)*s}{s^4 + 5000000.0*s^3 + 2.000015*10^{12}*s^2 + 3.5*10^{13}*s + 5.0*10^{13}}

The differential equation for this would be:

1000({d^2 V_s \over dt^2} + 5000000.0){d V_s \over dt} = {d^4 I_o \over dt^4} + 5000000.0{d^3 I_o \over dt^3} + 2.000015*10^{12}{d^2 I_o \over dt^2} + 3.5*10^{13}{d I_o \over dt} + 5.0*10^{13}I_o

Set the source Vs to zero and then find Io by adding together Second order techniques (Real roots, are underdamped. Equal roots are critically damped and imaginary roots are overdamped.

Solving for the zeros (setting denominator to zero):

solve(s^4 + 5000000.0*s^3 + 2.000015*10^12*s^2 + 3.5*10^13*s + 5.0*10^13 = 0, s)
\text{zeros} = -4561551.036, -438431.4637, -15.93127461, -1.569297305

So in this case the solution would be:

i_o = A e^{-4561551.036t} + B e^{-438431.4637t} + C e^{-15.93127461t} + D e^{-1.569297305} + E

The next step would be to find four initial conditions of i0 (and final condition). Most texts don't attempt to solve this complex of a problem. The steady state particular solution was easily found long ago with phasors. Transfer functions help to find a general solution to any type of voltage source. And it helps find the homogeneous solution that describes a ramp up when the circuit is first turned on. What circuits do engineers design that is worth this effort? Filters. The fact that today most of this is done digitally just means that all this needs to be repeated in a digital context. (Simulation becomes reality!). So using math tools to solve these problems in prepration for understanding more advanced tools is reason for continuing this thread of content.

So on to filters (Bode Plots .. Fourier Analysis) where the constants disappear and we stay in the complex domain. We will need the poles for to analyze as a filter (setting numerator to zero):

solve(s*(1000.0*s^2 + 5000000000.0)=0,s)
\text{poles} = 0, -2236.067977i, 2236.067977i