Circuit Theory/2Source Excitement/Example45

From testwiki
Revision as of 14:31, 22 December 2022 by imported>Ameisenigel (broken file links: file was deleted)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
Problem L = 1H, C=1F, R=1Ω, find ir

Particular/Steady State solution

Inductor short, cap open, Vs = 5 μ(t),find ir

Homogeneous/Transient Solution

Loop equation:

Vs(t)=VL(t)+VCR(t)
Vs(t)=Lditdt+VRC
VCR=iR*R
iC=CdVCRdt
it=iR+iC=VCRR+CdVCRdt
Vs(t)=L(d(VCRR)dt+Cd2(VCR)dt2)+VCR

Differential equation that needs to be solved:

0=Ld(VCRR)dt+LCd2(VCR)dt2+VCR

Guess:

VCR=Aest

Substitute to check if possible:

0=Ld(AestR)dt+LCd2(Aest)dt2+Aest
0=L(sAestR)+LC(s2Aest)+Aest
0=L(sAR)+LC(s2A)+A
0=LRs+LCs2+1

So the answer is going to be second order, thus guess was wrong, but can guess more accurately now by computing roots of the above second order equation:

s2+2s+1
s1,2=1,1

Both roots are negative and equal, so the new guess is:

VCR=Aet+Btet

Checking again by plugging into s2 + 2s + 1 = 0:

(Aet+BtetBetBet)+2(AetBtet+Bet)+Aet+Btet=0

Yes it equals zero now! So can go on. Have to add a constant to the differential equation solution so Vcr is:

VCR=Aet+Btet+C1

Without Initial Conditions .. Finding the Constants

Have initial conditions: VCR(0+) = 0 since initially cap is a short and impedance times the derivative of the inductor current it(0+) = 5. Turning this into an equation:

VCR(0+)=0=A+C1

The final voltage across the parallel RC combination is going to be 5 volts (after a very long time) because the capacitor opens and the inductor shorts.

VCR()=C1=5A=5

This is the matlab code that computes the limit:

syms A B C1 t
f = A*exp(-t) + B*t*exp(-t) + C1;
limit(f,t,inf)

Only B is unknown now:

VCR=5et+Btet+5

The initial voltage across the inductor is going to be 5 volts. But this does not lead to the value of B. Another initial condition is that the initial current through the capacitor (even though it is initially a short) is zero because the inductor is initially an open. This leads to B:

ic=CdVCRdt=5et+BetBtet
ic(0+)=5+B=0B=5

Now VCR is:

VCR=5(1ettet)

Which means that ir is:

iR=VCRR=10(1ettet)

Without C_1 constant

Trying to do this problem without the C_1 constant ends in something like this:

VL(0+)=5=B(22)

Which has no solution. Or it can lead to 5=5 where the constant disappears from the equation without finding a number for it. Or it can lead to A or B equaling infinity. Any of these non-answers means a mistake was made somewhere. Template:BookCat