PRACTICAL MATLAB® FOR ENGINEERS PRACTICAL MATLAB

Direct Current and Transient Analysis 171

5. 
   

   VR(t) versus t, where VR(t) = Ri(t)

   
   


6. 
   

   VL(t) versus t, where VL(t) = L
   di(t)
   ____dt

   
   


7. 
   

   VC(t) versus t, where VC(t) = (^) ( __C^1 ) ∫
   i (t) dt

   
   


8. 
   

   [VL(t) + VR(t)] ver su s t

   
   


9. 
   

   [VC(t) + VR(t)] ver su s t

   
   


10. 
    

    PR(t) versus t, where PR(t) = i^2 (t)R

    
    


11. 
    

    PL(t) versus t, where PL(t) = iL(t)VL(t)

    
    


12. 
    

    PC(t) versus t, where PC(t) = iC(t)VC(t)

    
    

FIGURE 2.68

Network of Example 2.22.

C = 0.5 μF

i(t)

+

−

−

R = 1000 Ω L = 30 mH

12

−0.5

i(t)

t

0.5

0

MATLAB Solution

% Script file: analysis

% Analysis of RLC Series circuit

% with non DC source

echo off;

syms t;

H1 = 2/pi*sin(pi*t);

H3 = 2/(3*pi)*sin(3*pi*t);

H5 = 2/(5*pi)*sin(5*pi*t);

H7 = 2/(7*pi)*sin(7*pi*t);

H9 = 2/(9*pi)*sin(9*pi*t);

H11 = 2/(11*pi)*sin(11*pi*t);

H13 = 2/(13*pi)*sin(13*pi*t);

H15 = 2/(15*pi)*sin(15*pi*t);

iamps = H1+H3+H5+H7+H9+H11+H13+H15;

figure(1)

ezplot(iamps,[0,3])

xlabel(‘ time ‘);ylabel(‘i(t)in amps’);

title(‘i(t) vs t’);grid on;

isquare=iamps^2;

figure(2);%part(c.2)

ezplot(isquare,[0,3])

xlabel(‘ time ‘);ylabel(‘i(t)^2’);

title(‘i(t)^2 vs t’);grid on;