Circuit Theory/Convolution Integral/Examples/example49/VL

From testwiki
Jump to navigation Jump to search
series LRC circuit ... find voltage across the resistor

Given that the source voltage is (2t-3t2), find voltage across the resistor.

This is the VL solution.

Outline:

Transfer Function

H(s)=VLVS=s4+s+10.25s
simplify(s/(4 + s + 1/(0.25*s)))
H(s)=s2s2+4s+4

Homogeneous Solution

solve(s^2 + 4.0*s + 4.0,s)

There are two equal roots at s = -2, so the solution has the form:

VLh(t)=Ae2t+Bte2t+C1

Particular Solution

After a long time attached to a unit step function source, the inductor has shorted and the capacitor has opened. All the drop is across the capacitor.

VLp=0

This also means that C1 has to be zero.

Initial Conditions

So far the full equation is:

VL(t)=Ae2t+Bte2t

Initial voltage is all across the inductor.

VL(0)=1=A
A=1

At this point will have to do integral .. to get to the current. There is no other way to use the known initial conditions: current (initially zero), and VC (initially zero). Will have to introduce integration constant and then evaluate that. More chance of mistakes, more complex, so start over with something else.

i(t)=1L0tVL(t)dt=0t(e2xB(x)e2x)dx
f := (exp(-2*x) - B*x*exp(-2*x));
S :=int(f,x=0..t)
i(t)=B*(e2t(2t+1)414)e2t2+12+C1
i(0)=0=C1

Ok so C1 is zero. Now need to find B. Find B by doing another integral to get VC:

VC(t)=1C0ti(t)dt=4*0t(B*(e2t(2t+1)414)e2t2+12)dx
f := (4*(B*(exp(-2*x)*(2*x+1)/4 -1/4) - exp(-2*x)/2 + 1/2));
S :=int(f,x=0..t)
VC(t)=B+2t+e2tBtBe2tBte2t1+C1
VC(0)=0=B+1B1+C1
C1=0

Still doesn't help us find B. Guess B = 2 since (2t-Bt) has to equal zero if VC is going to converge on 1. Then see if VC(∞) = 1:

VC(t)=2+e2t2e2t2te2t1

At t = ∞ what is the 2te-2t term's value?

limit(B*t*exp(-t),t = infinity)

Mupad says 0.

VC()=21=1

Yes! B = 2 works ... looks like the only thing that works .... So:

i(t)=2*(e2t(2t+1)414)e2t2+12=te2t
simplify(2*(exp(-2*t)(2*t+1)/4-1/4) - exp(-2*t)/2 + 1/2)

This means that VR is:

VR(t)=4*i(t)=4te2t

Impulse Solution

Taking the derivative of the above get:

VRδ(t)=4e2t8te2t

Convolution Integral

VR(t)=0t(4e2(tx)8(tx)e2(tx))(2x3x2)dx
f := (4*exp(-2*(t-x)) - 8*(t-x)exp(-2*(t-x)))*(2*x-3*x^2);
S :=int(f,x=0..t)
VR(t)=88e2t10te2t6t

There will not be any constant since again, V_R(t) = 0 after a long time ... and the capacitor opens. Template:BookCat