Circuit Theory/Convolution Integral/Examples/Example43

From testwiki
Jump to navigation Jump to search
2R 1L 1C circuit to solve using the convolution integral .. for wikibook circuit analysis

Given that is = 1 + cos(t), find io using the convolution integral.

Outline:

Transfer Function

H(s)=𝕀o𝕀s=𝕍T𝕀s*𝕀o𝕍T=111/s+11+1s+1*1s+1
simplify(1/((s+1)*(s + 1 + 1/(s+1))))
H(s)=1s2+2s+2

Homogeneous Solution

Setting the denominator to zero and find the values of s:

solve(s^2 + 2*s + 2)
s=1±i

This means the solution has the form:

ioh(t)=et(Acost+Bsint)+C1

Particular Solution

After a long time the cap opens and the inductor shorts. The current is split between the two resistors. Each will get 1/2 of the unit step function source which would be 1/2 amp:

iop(t)=12

Initial Conditions

io=12+et(Acost+Bsint)+C1

The particular solution still has to apply so at t= ∞:

io()=12=12+C1
C1=0

Initially the current has to be zero in this leg so:

io(0)=12+A=0
A=12

The initial voltage across the cap is zero, and across the leg is zero and across the inductor is zero. So:

VL(t)=L*dio(t)dt
f := 1/2 + exp(-t)*((-1/2)*cos(t) + B*sin(t));
g = diff(f,t)
VL(t)=et(B(cos(t)sin(t))+12(cos(t)+sin(t))
VL(0)=B+12=0
B=12

So now:

io=12(1et(cost+sint))

Impulse Response

Taking the derivative of the above

f := 1/2*(1-exp(-t)*(cos(t) + sin(t)));
g = diff(f,t)

get:

ioδ(t)=etsint

Convolution Integral

io(t)=0tet+xsin(tx)*(1+cos(x))dx
f := exp(x-t)*sin(t-x)*(1 + cos(x));
S =int(f,x=0..t);
io(t)=0.2cos(t)+0.4sin(t)0.7cos(t)et1.1sin(t)e1+0.5

Template:BookCat