274 Practical MATLAB® Applications for Engineers
MATLAB Solution
% Script file: curr _ volt
W = 300;
R = 19;
L = 60e-3;
Ang = W*L/R;
t = [0:.0001:.06];
V = (R^2+(W*L)^2)^.5*5*sin(W.*t+atan(Ang));
subplot(2,1,1)
I = 5*sin(W.*t);
plot (t,V), title(‘ v(t) vs. t’)
grid on
ylabel(‘Amplitude in volts’), title(‘v(t) vs.t’)
subplot(2,1,2)
plot(t,I)
grid on
title(‘i(t) vs. t’), ylabel(‘ Amplitude in amps’)
xlabel(‘time in sec.’)
The script fi le curr_volt is executed and the resulting plots are shown in Figure 3.55.
FIGURE 3.55
Plots of Example 3.6.
200
100
0
0
0
0.01 0.02 0.03 0.04 0.05 0.06
0 0.01 0.02 0.03
time in sec.
0.04 0.05 0.06
− 100
− 200
− 5
v(t) versus t
i(t) versus t
Amplitude in volts
Amplitude in amps
5