PRACTICAL MATLAB® FOR ENGINEERS PRACTICAL MATLAB

(sharon) #1

Alternating Current Analysis 303


title(‘angle(VC) vs. w’,’FontSize’,9)
xlabel(‘frequency in w (rad/sec)’), ylabel(‘angle (degrees)’)
figure(3)
IC = VC./XC; IL = VC./ZS;
subplot(2,2,1)
plot(w,abs(IC))
grid on
set(gca,’fontsize’,9)
title(‘abs[IC(w)] vs. w’,’FontSize’,9)
ylabel(‘Current (amps)’)
subplot(2,2,2)
plot(w,angle(IC)*deg)
grid on
set(gca,’FontSize’,9)
title(‘angle[IC(w)] vs. w’,’FontSize’,9)
ylabel(‘angle (degrees)’)
subplot(2,2,3)
plot(w, abs(IL))
grid on
set(gca,’FontSize’,9)
title(‘abs(IL) vs. w’,’FontSize’,9)
xlabel(‘frequency w (rad/sec)’), ylabel(‘Current (amps)’)
subplot(2,2,4)
plot(w,angle(IL)*deg)
grid on
set(gca,’FontSize’,9)
title(‘angle(IL) vs. w’,’FontSize’,9)
xlabel(‘frequency in w (rad/sec)’), ylabel(‘angle (degrees)’)
The script fi le RLC_parallel_ analysis is executed and the results are shown in Figures 3.81
through 3.83.

0 5 10 15 0 5 10 15

(^051015051015)
× 104
× 10 −^5 ×^10
(^4) × 104
abs[Z(w)] versus w angle[Z(w)] versus w
abs(YT) versus w angle(YT) versus w
Amplitude (Ohms)
Amplitude (Ohms
−^1 )
Angle (degrees)
− 10
− 5
0
5
10
Angle (degrees)
− 10
− 5
0
5
10
6
5.5
5
4.5
frequency w (rad/sec) × 104 frequency w (rad/sec) × 104
2.2
2.1
2
1.9
1.8
1.7
FIGURE 3.81
Plots of Example 3.19.

Free download pdf