Sinusoidal Steady State Response of Linear Circuits The circuit shown on Figure 1 is driven by a sinusoidal voltage source vs(t) of the form vtso()= vcos(ωt) (1.1) i(t) R +-vR(t) + vs(t) vc(t) C - Figure 1. Series RC circuit driven by a sinusoidal forcing function Our goal is to determine the voltages vc(t) and the current i(t) which will completely characterize the “Steady State” response of the circuit. The equation that describes the behavior of this circuit is obtained by applying KVL around the mesh. vtRc()+ v()t= vs(t) (1.2) Using the current voltage relationship of the resistor and the capacitor, Equation (1.2) becomes dv ()t RCvc +=()tvcos(ωt) (1.3) dt co Note that the coefficient RC has the unit of time. (Ohm)(Farad) → seconds Before proceeding with the solution of this differential equation let’s explore its physical significance. This linear circuit is driven (forced) by an independent sinusoidal source, vs(t). We may view its response (its effect on the circuit) as the superposition of the response with the source set equal to zero (source-free or natural response (vch(t)) ) and the forced response (vcp(t)). vtcc()= vh()t+ vcp(t) (1.4) 6.071/22.071 Spring 2006, Chaniotakis and Cory 1 Schematically the superposition is shown on Figure 2(a) and (b). R R + + vch(t) C vs(t) vcp(t) C - - (a) (b) Figure 2 In mathematical language we call these two responses the homogeneous solution (vch(t)) ) and the particular solution (vcp(t)) of the equation that characterizes the system (Equation (1.3) in our case) The homogeneous solution corresponds to the differential equation dv ()t RC ch + v ()t = 0 (1.5) dt ch And the particular solution to the equation dv ()t RCvcp +=()tvcos(ωt) (1.6) dt cp o The homogeneous solution (or the natural response of the system) has the form ⎡ −t ⎤ vtch ()= Bexp (1.7) ⎣⎢ RC ⎦⎥ The particular solution (or the forced response of the system) is the cosine function of amplitude A, frequency ω, and phase φ. vtcp ()= Acos(ωt+φ) (1.8) And so the general form of the system response is ⎡ −t ⎤ vtc ()=+Acos(ωφt )+Bexp (1.9) ⎣⎢ RC ⎦⎥ 6.071/22.071 Spring 2006, Chaniotakis and Cory 2 At this time let’s recall the problem statement which says that we are interested in obtaining the Steady State response of the system. This is equivalent to saying that the source vs(t) was connected to the system long time ago and all transient phenomena are gone. In order to be more precise, if the time t RC then the exponential term in Equation (1.9) would go to zero. In this case the observable and thus the important response of the system is the Steady State response which is given by vtc ()= Acos(ωt+φ) (1.10) We may now proceed to determine the details of the solution by calculating the amplitude A and the phase φ. To do this we substitute the form of the solution (Equation (1.10) ) into the differential Equation (1.3). Before we make that substitution lets use the trigonometric identities to expand Equation (1.10). vt()=+Acos(ωt φ) c (1.11) =−Atcosφ cos(ωφ) Asin sin(ωt) And the corresponding derivative dv c =−Aω cosφωsin( tA) − ωsinφcos(ωt) (1.12) dt Substituting Equations (1.11) and (1.12) into Equation (1.3) we have −−Atω cosφωsin()Aωsinφcos()ωt+ 1 v (1.13) [Atcosφ cos()ωφ−=Asinsin()ωt]0 cos(ωt) RCRC Collecting the coefficients of cos(ωt) and sin(ωt) we have the following equations for the unknowns A and φ. A −Aωφcos −sinφ=0 (1.14) RC A v −+Aωφsin cosφ=0 (1.15) RCRC These equations are independent and thus they may be solved simultaneously for A and φ. 6.071/22.071 Spring 2006, Chaniotakis and Cory 3 From Equation (1.14) we obtain the phase φ =arctan (−RC ω) (1.16) and A becomes v A = o (1.17) cosφ −ωφRC sin Therefore, the Steady State response of the system is v vt()=o cos(ωt+φ ) (1.18) c cosφω− RC sinφ whereφω=arctan ()−RC Figure 3 shows a plot of the phase and the amplitude ratio Av/ o as a function of the dimensionless parameter ωRC . Observe the frequency selectivity of this system. It passes “low” frequencies while it attenuates “higher” frequencies. 0.0 -10.0 -20.0 -30.0 -40.0 -50.0 -60.0 -70.0 -80.0 -90.0 0 1 2 3 4 5 6 7 8 9 10 ωRC 1.0 0.8 0.6 0.4 0.2 0.0 0 1 2 3 4 5 6 7 8 9 10 ωRC Figure 3 6.071/22.071 Spring 2006, Chaniotakis and Cory 4 The voltage across the resistor, vR (t) , may also be determined by calculating the current i(t) and multiplying it by R. dvco()t v Cω ⎛π ⎞ it()==C cos⎜ωφt++⎟ (1.19) dt cosφω− RC sinφ ⎝⎠2 vRoω C ⎛π ⎞ vtR ()==i()tR cos⎜ωφt++⎟ (1.20) cosφω− RC sinφ ⎝⎠2 Figure 4 shows the plot of amplitude ratio vR/ vo as a function of the dimensionless parameterωRC . Here we see the complementary behavior to that shown on Figure 3. Observe again the frequency selectivity of this system. If the output is thus taken across the resistor it passes “high” frequencies while it attenuates “lower” frequencies. 1.0 0.8 0.6 0.4 0.2 0.0 0 1 2 3 4 5 6 7 8 9 10 ωRC Figure 4. 6.071/22.071 Spring 2006, Chaniotakis and Cory 5 Now let’s look at a few frequencies of interest 1. For ω = 0 (dc signal) the phase φ = 0 and the amplitude A= vo . The voltage across the capacitor is constant. No current flows in the circuit The capacitor behaves as an open circuit. The circuit equivalent when ω = 0 is shown on Figure 5. R R ≡ vo C vo C Figure 5 2. For ωRC =1 vv The phase φ =−450 and the amplitude A ==oo 1/ 2 +1/ 2 2 R vo vtcos(ω ) C cos(ωπt − / 4) o 2 Here current is flowing and thus some power is dissipated in R. Also the capacitor stores some energy. 6.071/22.071 Spring 2006, Chaniotakis and Cory 6 3. For ωRC 1 vv The phase φ =−900 and the amplitude A = oo=→0 0(−−ωωRC 1)RC As ω increases more power is dissipated in the resistor R. When ω →∞ the capacitor acts as a short circuit and all power is dissipated in R. The circuit equivalent when ω →∞ is shown on Figure 6. R R C C vto cos(ω ) ⎯⎯ω→∞⎯→ vto cos(ω ) Figure 6 6.071/22.071 Spring 2006, Chaniotakis and Cory 7 Now let’s consider the RL circuit i(t) R + vs(t) vL(t) L - We would like to characterize this circuit by obtaining the current i(t) and the voltage vL(t). Apply KVL around the mesh di()t vt()=+i()tR L (1.21) s dt With vs(t) of the form vtso( ) = vcos(ωt) , Equation (1.21) becomes di()t R v +=it() 0 cos(ωt) (1.22) dt L L L Henry The coefficient is a time constant. → seconds R Ohm Here again we are interested in the behavior of the system for times that are long L compared to the time constant . In this scenario the only contribution to the solution is R again the one forced by the sinusoidal source voltage vtso( ) = vcos(ωt) . And the form of the forced response solution is it()=+Acos(ωt φ) (1.23) =−Acosφ cos()ωφtAsinsin()ωt di()t =−Aω cosφωsin( tA) − ωsinφcos(ωt) (1.24) dt Substituting back into the differential equation (1.22) we get 6.071/22.071 Spring 2006, Chaniotakis and Cory 8 −−Atω cosφωsin()Aωsinφcos()ωt+ R v (1.25) [Atcosφ cos(ωφ) −=Asin sin(ωt)]0 cos(ωt) LL We now separate the cosine and sine functional forms of the solution since they are independent contributions. The coefficients of the sin(ωt) terms are: RA −Aωφcos −sinφ=0 (1.26) L And the coefficients of the cos(ωt) terms are: RA v −+Aωφsin cosφ=0 (1.27) L L Equation (1.26) gives us the expression for the phase φ. −1 ⎛ωL ⎞ φ =tan ⎜−⎟ (1.28) ⎝⎠R And the amplitude A becomes vR/ A = o (1.29) ωL cosφ − sinφ R And the solution for i(t) becomes vR/ it()=o cosωt+φ (1.30) L () cosφω− sinφ R −1 ⎛⎞ωL Where the phase φ =−tan ⎜⎟ ⎝⎠R Figure 7 (a) and (b) show the plots of the phase and the amplitude as a function of the ωL parameter . R 6.071/22.071 Spring 2006, Chaniotakis and Cory 9 0.0 -10.0 -20.0 -30.0 -40.0 -50.0 -60.0 -70.0 -80.0 -90.0 0 1 2 3 4 5 6 7 8 9 10 ωL/R (a) 1.0 0.8 0.6 0.4 0.2 0.0 0 1 2 3 4 5 6 7 8 9 10 ωL/R (b) Figure 7 6.071/22.071 Spring 2006, Chaniotakis and Cory 10 Let’s look at a few important frequencies as before 1. For ω = 0 (dc signal) the phase φ = 0 and the amplitude Av= o / R. The voltage across the inductor is zero and the current flowing in the circuit is vRo / .