Circuit Theory/TF Examples/Example14

Finding i_o using transfer functions

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

$\frac{{𝕀}_o}{{𝕍}_s}=\frac{{𝕀}_o}{{𝕀}_t}\frac{{𝕀}_t}{{𝕍}_s}$

Which is of the form:

$\frac{{𝕀}_o}{{𝕍}_s}=\text{current divider}*\text{admittance}$

The current divider is going to be:

$\frac{{𝕀}_o}{{𝕀}_t}=\frac{\frac{1}{R_3}}{\frac{1}{R_3}+\frac{1}{s{L}_2+\frac{1}{s{C}_2}}}=\frac{s{L}_2+\frac{1}{s{C}_2}}{R_3+s{L}_2+\frac{1}{s{C}_2}}$

Substituting into the transfer function:

$\frac{{𝕀}_o}{{𝕍}_s}=\frac{s{L}_2+\frac{1}{s{C}_2}}{R_3+s{L}_2+\frac{1}{s{C}_2}}*\frac{1}{R_1+\frac{1}{s{C}_1}+s{L}_1+\frac{1}{\frac{1}{R_3}+\frac{1}{s{L}_2+\frac{1}{s{C}_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{{𝕀}_o}{{𝕍}_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(\frac{d^2{V}_s}{d{t}^2}+5000000.0)\frac{d{V}_s}{dt}=\frac{d^4{I}_o}{d{t}^4}+5000000.0\frac{d^3{I}_o}{d{t}^3}+2.000015*10^{12}\frac{d^2{I}_o}{d{t}^2}+3.5*10^{13}\frac{d{I}_o}{dt}+5.0*10^{13}{I}_o$

Set the source V_s to zero and then find I_o 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 i₀ (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 preparation 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$